Surface – particle

Particles can be sampled for unbiased estimates of surface in a number weighted, volume weighted, or surface-weighted manner.

Optical Rotator


thick section

isotropicTo estimate the surface of particles or spaces with the optical rotator (Tandrup, etal., 1997), isotropic, thick sections must be used. For details of the mechanics of the probe, please see optical rotator to estimate particle volume, but instead use this formula:

s = a ∑ lj g(lj)                (Tandrup, etal., 1997, equation 6)

s = estimate of surface

a = density of probing lines

lj = number of intersections of probing lines with the surface

OR surface formulas

d2 = the shortest distance from the grid line on the cross section to a unique point in the particle (Tandrup, etal., 1997, Fig. 2)

t = half of the thickness of the optical slice (Tandrup, etal., 1997, Fig. 1)





To estimate the surface area of particles or spaces with the surfactor (Jensen and Gundersen, 1987), isotropic sections must be used. Thick sections must be used if sampling is number weighted. If sampling is volume or surface weighted, thin sections may be used. As for all of the local particle volume and surface probes, we start by marking a unique point inside the particle (‘o’ in figure below).


Then a ray, l, is drawn from the unique point to the boundary. In the figure above two rays have been drawn one-hundred and eighty degrees from each other for efficiency. Three or four rays could have been used for more precision. The line segment from the unique point, o, to the boundary on the right has that distance, ‘l’, labeled. The thinner black line shows the tangent where ‘l’ hits the boundary, and the red line is perpendicular to the tangent. The angle, alpha, is between the normal and ‘l’. The second line segment extends from the unique point, ‘o’, and goes to the boundary on the left. Notice that the angle between the normal and the line from ‘o’ to the left boundary is greater than alpha on the right. This means there will be more surface area at the left intersection than the right intersection. If, unlike the figure, only one ray is used, the formula for the surface estimate is:

S = 4 π l2 c(β)                       (Jensen and Gundersen, 1987, equation 7)

S = surface area estimate

l = distance from unique point to boundary

β = the angle between the intersection of the surface and a plane that contains ‘l’, and ‘l’ (please see Jensen and Gundersen, 1987, Fig. 1 to understand how α is related to β)


surfactor equations

If more than one ray is used (e.g., in the figure above, two are used) the formula for the surface estimate is:

S = 2 π ∑l2 c(β)                    (Jensen and Gundersen, 1987, equation 8)

∑l2= sum from 1 to 2r of the distances from the unique point to the boundary squared

r = number of rays


Vertical Spatial Grid


verticalthick section


To get an estimate of the surface area (and volume) of a group of particles, the vertical spatial grid (Cruz-Orive and Howard, 1995) can be used. Unlike the surfactor and isotropic rotator (see above) that are local probes, the vertical spatial grid is a global probe. The local probes have an advantage over the vertical spatial grid in that for global probes there can be difficulty focusing on the tops and bottoms of the particles (Tandrup, etal., 1997, pg. 108, last paragraph).

The idea is to use two arrays of cycloid lines that are perpendicular to each other and have the minor axis of the cycloid lined up with the vertical axis (Cruz-Orive and Howard, 1995, fig. 3). Start with systematic random sampling so that you have a series of thin sampling planes. It will be easy to mark the intersections (I) of the cycloid lines that are in the XY plane with the line that is created by the cross section of the surface of the particle. Those intersections are called Ixy (Cruz-Orive and Howard, 1995, fig. 2b and 2c, where the cycloid intersects the boundary of the particle).

intersections focal plane one

In the figure above, the vertical axis direction had already been indicated and the cycloids are laid down with their minor axes aligned with the vertical direction. The red triangles indicate Ixy, the place where the surface of the particles crosses the lines. But we can’t see the cycloids in the XZ plane as lines because they are coming in and out the plane of section and would appear as points.

points focal plane one

So points along the XY cycloids (now dotted) are shown that represent the XZ cycloids (dashes along the dotted lines). If one of these goes through the particle it is marked (see red triangle above) indicating an intersection that will be used to calculate the number of Ixz. The points also can serve to estimate the volume of all of the particles.

Now the next sampling plane is brought into focus. In the images below the cross sections from the first sampling plane (green cross sections above) are shown in yellow, and the cross sections that are now in focus are in green.

intersections focal plane two

Again, it is easy to mark the intersections of the cycloids in this XY plane with the surface (red triangles above) to get Ixy. But to get Ixz we have to look at points, and to further complicate things, the lines are cycloids, not straight. This is why the XY grid of cycloids has been vertically shifted upward compared to where it was at the last focal plane (compare how the cycloids go through the yellow cross sections above to the green cross sections in the first image above). They are shifted (Cruz-Orive and Howard, 1995, pg. 149, second column) to follow the cycloids that are in the XZ plane as we focus up and down. Now we mark the points that will help calculate Ixz as well as be used for a volume estimate.

points focal plane two

In the case above, none of the dashes on the white dotted cycloid-lines fall inside the cross section, so nothing is marked. The program will keep track of how many of these points were the same points that were counted in the previous sampling plane. This is  Pi,i+1 (Cruz-Orive and Howard, 1995, Fig. 2), that in turn is used to calculate the number of intersections in the XZ plane:

Ixz,i = Pi + Pi+1 – 2Pi,i+1                (Cruz-Orive and Howard, 1995, fig. 2, figure legend)

In other words, the number of intersections of the cycloids going in and out of the plane with the surface of the particles is the number of points that hit a cross section on one focal plane plus the number of points that hit the cross section on the next focal plane minus twice the number of points that hit the cross section in both focal planes. Here is how to calculate the surface area estimate:

vsg equation                          (Cruz-Orive and Howard, 1995, equation 4)


S = the estimate of surface area for all the particles

a = test area per test point

dz = distance between sampling focal planes

l = length of cycloid per test point

Ixy = the sum of all the intersections counted over all of the focal planes

Ixz,i = Pi + Pi+1 – 2Pi,i+1


We are also reminded that we can calculate a volume estimate using the Cavalieri/point-counting estimator:

V = a (dz) ∑Pi               (Cruz-Orive and Howard, 1995, equation 9)

V = volume estimate

a = test area per test point

dz = distance between sampling focal planes

∑Pi = all of the points for all of the sampling-focal-planes added together



Cruz-Orive, L.M. and C.V. Howard (1994) Estimation of Individual Feature Surface Area with the Vertical Spatial Grid. J. of Microscopy, 178, 146 -151.

Jensen, E.B. and H.J. Gundersen (1987) Stereological Estimation of Surface Area of Arbitrary Particles. ACTA Stereol., 6/III, 25 – 30.

Tandrup, T., Gundersen, H.J.G. and E.B. Vedel Jensen (1997) The Optical Rotator. J. of Microscopy, 186, 108 -120.