Stereology: 3D Structure from 2D Sections

Why a flat slice lies about size

Everything you have measured so far in this rung came off a flat, polished micrograph — a two-dimensional section cut through a three-dimensional solid. Guide 2 already dropped a warning: the linear-intercept method gives a robust number, but a single section shows each grain only as the shape where the cutting plane happened to slice it, not as the grain itself. Picture slicing a salami, or a bag of oranges set in resin: your knife almost never passes through the fattest part of any one piece. Most cuts catch a sphere off-centre, so the little disc on the polished face is smaller — usually much smaller — than the sphere's true diameter.

How big is the lie? Take a sphere of diameter D and cut it with random parallel planes. The section is a circle whose diameter depends on where the plane fell, and averaged over all heights it comes out only about 0.79 D across — and often far smaller, because a cut near the pole gives a tiny circle. Draw a random line through the sphere instead, as the intercept method does, and its mean chord is exactly 4V/S = 2D/3, about 0.67 D (this is Cauchy's theorem for a convex body). Either way the section systematically shrinks what you see. So you cannot simply read a mean grain size off the apparent circles: the raw 2D picture is biased small, and its apparent size distribution is not the true 3D one.

The crown jewel: volume fraction, no shape assumed

The first and cleanest result is an identity discovered by the geologist Delesse in 1847. For a phase dispersed through a solid, its volume fraction equals its area fraction on a random section, which equals the fraction of random test lines lying in it, which equals the fraction of random test points that fall on it: Vv = Aa = Ll = Pp. Read that chain slowly, because it is remarkable — a volume you can never see is delivered exactly by counting dots on a flat picture, with no model of how the constituent is shaped. The easiest version to use is the last one, point counting: overlay a grid and count hits.

Work a case. Sprinkle a 100-point grid over a two-phase steel and 20 points land on the dark second phase, so its phase fraction is 20 percent by volume — full stop, no shape correction needed. The catch is statistics, not bias: with only 20 hits the relative uncertainty is roughly 1/sqrt(20) ~ 22 percent, so you actually have 20 plus-or-minus 4 percent. Precision improves only as 1/sqrt(count), so an accurate fraction needs hundreds or thousands of points spread over many fields of view. The identity is exact; the honest work is counting enough.

QUANTIFYING A 2D SECTION with a test grid   (matrix . , second phase #)

   o . . o . . o . . o . . o      o = test point   (25 points shown)
   . . . ##### . . . . . . .
   o . . ##### o . . o . . o      POINT COUNT (Delesse, model-free):
   . . . ### . . . . ## . .          points on #  /  all points
   o===o===#===o===o===#===o         = 5 / 25 = 0.20  ->  V_v = 20 %
   . . . . . . . . . ## . .          (a volume, exact, no shape assumed)
   o . . o . . o . . o . . o
                                  LINE INTERCEPT (count boundary crossings):
   the === row is a test line         crossings N on a line of true
   of true length L; count the        length L  ->  N_L = N / L
   interfaces N it crosses            S_v = 2 x N_L ,  L-bar = 1 / N_L

   MODEL-FREE (needs isotropic-uniform-random sampling):
      V_v = A_a = L_l = P_p        S_v = 2 N_L        L-bar = 1 / N_L

Counting crossings: boundary area and grain size

The same counting spirit reaches surfaces you cannot see edge-on. Lay a test line of known length across the section and count how often it crosses a grain boundary; call that count per unit true length N_L. Then the grain-boundary area per unit volume is exactly Sv = 2 N_L — again model-free, given fair sampling. This is the hidden engine under guide 2's linear-intercept method: the mean intercept length L-bar = 1/N_L is just the reciprocal of the same count, so measuring grain size and measuring boundary area are two readings of one tally.

Numbers make it concrete. Suppose a test line whose true length (after dividing by the magnification) is 5 mm crosses 250 boundaries. Then N_L = 50 per mm, the mean intercept is L-bar = 1/50 mm = 20 micrometres, and the grain-boundary area per volume is Sv = 2 x 50 = 100 mm^2 per mm^3. The factor of 2 is pure geometry: a flat interface tilted at every angle to the test line is, on average, crossed at a rate that undercounts its true area by exactly one half, and the 2 puts it back. That 20 micrometre intercept is what an ASTM grain-size number encodes — a coarser structure gives longer intercepts and a smaller ASTM number.

1. Prepare a plane, random section — grind, polish, and etch — so grains or phases show clear contrast, and record the true magnification so lengths on the image convert to real lengths.
2. For a volume or phase fraction, overlay a regular point grid and count the fraction of points landing on the feature: Pp = Vv, no shape model, just enough points for good statistics.
3. For grain-boundary area or grain size, overlay test lines of known true length and count boundary crossings N; then N_L = N/length, Sv = 2 N_L, and mean intercept L-bar = 1/N_L.
4. Repeat over many random fields and — for anything not clearly equiaxed — several section orientations, so the sampling is isotropic, uniform, and random.
5. Report the method with the number: an intercept size, an ASTM number, and an equivalent-circle diameter are different conventions and will not agree to the last digit.

What is exact, and what needs a model

Draw the honest line clearly. The mean quantities above — volume fraction, boundary area per volume, mean intercept — are exact and shape-free. But the moment you want the full 3D size distribution from the 2D one, you must assume a shape (usually spheres) and solve an inverse problem, the Saltykov or Wicksell unfolding. That inversion is ill-conditioned: small counting errors blow up, and it can even hand back negative numbers of grains in a size bin. So trust stereological means; treat a reconstructed 3D distribution as a fragile estimate, not gospel.

One prize genuinely lies beyond a single plane: the number of grains or particles per unit volume, Nv. From one section you can count particles per unit area, but converting to Nv needs their mean size — which you were trying to find — so the classic route is circular and shape-bound. Modern design-based stereology escapes this with the disector: look at two parallel sections a known small distance apart and count only the particles that appear in one but not the other (their tops). That count, divided by the sampled volume, gives Nv with no assumption about shape or size at all — the same unbiased spirit as Delesse, one dimension up.

From geometry to orientation

Stereology masters the geometry of the polycrystal — how much of each phase, how large the grains, how much boundary area — but it is completely blind to one thing: which way each grain's crystal lattice points. Two grains can have identical size and shape on the polished face while their atomic planes aim in utterly different directions, and no amount of point- or intercept-counting can tell them apart. Yet that orientation is exactly what makes a polycrystal anisotropic, and it is set by how the material was rolled, drawn, or grown.

Capturing that orientation needs a different probe and a different bookkeeping. Instead of a light micrograph you point a diffraction beam at each spot — electron backscatter diffraction reads a crystal orientation at every pixel of a grain map — and instead of a size distribution you summarise the result as a preferred orientation, or texture. The next two guides build exactly this language: the pole figure and inverse pole figure that plot which way the crystals face, and the orientation distribution function that captures the full statistics of orientation across the aggregate.