## Abstract

We utilize Fourier-holographic light scattering angular spectroscopy to record the spatially resolved complex angular scattering spectra of samples over wide fields of view in a single or few image captures. Without resolving individual scatterers, we are able to generate spatially-resolved particle size maps for samples composed of spherical scatterers, by comparing generated spectra with Mie-theory predictions. We present a theoretical discussion of the fundamental principles of our technique and, in addition to the sphere samples, apply it experimentally to a biological sample which comprises red blood cells. Our method could possibly represent an efficient alternative to the time-consuming and laborious conventional procedure in light microscopy of image tiling and inspection, for the characterization of microscopic morphology over wide fields of view.

© 2006 Optical Society of America

## 1. Introduction

In biological systems, it is often necessary to determine structural information over a large range of size scales. For example, the size and shape of individual cells and the distribution of these cells throughout a sample both represent important features necessary for tissue characterization. Histological analysis of tissue sections is the most informative way to assess cellular composition and events in many situations [1]. All such examinations are laborious tasks, particularly when large tissue areas must be surveyed to obtain statistically significant cell populations, yet information on the microscopic scale must be discerned in order to provide accurate characterization. Recently, our group proposed a new approach to perform high-throughput structural characterization of a wide range of biological samples [2, 3]. Our technique is based on the intrinsic dependence of the angular distribution of elastic light scattering upon sample microstructure. By employing digital Fourier holography, we seek to map, with a single camera exposure, scatterer sizes (and potentially refractive-index contrast) over exceptionally large fields of view. Ultimately, this technique could enable single-image-capture ‘histological’ assessment of the microscopic properties of millimeter-scale samples without histological staining and without microscopic examinations such as cell counting - indeed, without even resolving the scatterers. Our publications thus far have entailed the first proposal and demonstration of this technique [2], and a preliminary experimental and theoretical analysis [3]. In this paper, we present a thorough theoretical treatment of the technique, and for the first time demonstrate experimentally the ability of our approach to generate fully spatially-resolved particle sizemaps for non-biological samples, and spatially-resolved angular spectra of a biological sample. We finally discuss outstanding issues and limitations in its implementation.

The technology of digital holographic imaging has developed considerably in recent years [4] due to the accessibility of high-power personal computers and advances in digital recording hardware. Its primary advantage over alternative techniques is the ability to form a digital representation of the full complex field distribution of an image, thus enabling the application of image processing and reconstruction algorithms to extract sample features, without recourse to complex optical setups. The renewed impetus in the field has been particularly striking with respect to the application of digital holographic microscopy, and a considerable body of literature has developed. Published journal papers have included treatments of quantitative phase contrast microscopy [5, 6, 7, 8], refractometry [9, 10], three-dimensional visualization [11], and synthetic-aperture optical imaging to achieve superresolution [12], or high resolution over large fields of view while utilizing low-numerical-aperture optics [13].

In parallel with developments in digital holographic microscopy, methods of characterizing a sample using the angular distribution [14] or spectrum [15] of elastically scattered light have been developed. These techniques have been termed ‘light scattering spectroscopy’ (LSS) and have been developed to quantify cellular or subcellular morphology (structure). LSS has shown promise in the detection of nuclear pleomorphism (the increase in size and shape variation of cell nuclei associated with dysplasia) in epithelial tumors *in vivo* [15, 16, 17, 18] and has been applied to bacteria size determination [19, 20]. A closely related microscopy technique has combined LSS with optical Fourier filtering to map morphology over a microscope’s field of view [21]. More recently, the importance of angular distributions has begun to be appreciated [22, 23], including the combined measurement of angular and spectral scattering distributions [24]. All such techniques are not intrinsically depth sectioning, which is a common important requirement for thick or turbid media. Both polarization gating [15, 16, 25, 26, 27] and coherence gating [28, 29] have been used to determine the depth-resolved spectral or angular distributions of scattered light.

In quantitative microscopy, the morphological parameters of samples are generally measured directly and accuracy depends largely on spatial resolution. The requirement for high spatial resolution typically restricts fields of view to substantially less than the region of interest. In LSS, microscopic morphology and spatial resolution are decoupled; measurements averaged over small but macroscopic areas of the sample are sensitive to microscopic morphology. Accurate measurements of the scatterer sizes can still be obtained, via the convenient use of lowresolution optics. The angular scattering spectra published elsewhere [14, 19, 20, 22, 23, 24, 28] were obtained point by point over a sample area. Our technique enables such spectra to be collected over millimeter-scale samples with a single image capture.

The remainder of the paper is organized as follows. In Section 2, we describe in detail the theory of the technique. In Section 3, we present its experimental realization and an analysis of the experimental setup. In Section 4, we present experimental results. We provide some discussion and draw some conclusions in Section 5.

## 2. Methodology

To determine the angular scattering spectra for each local area of a sample, it is sufficient to record the complex amplitude distribution of the scattered wave. Our approach is based on recording a digital hologram of the scattered light in a plane optically conjugate to the Fourier plane of the object. By performing filtering in either the Fourier plane or the reconstructed object plane, we are able to generate spatially-dependent angular scattering profiles for the sample.

Figure 1 is a schematic showing the orientations of the object (*ξ*,*η*) and recording (*x*,*y*) planes, and the coordinate systems used to describe them. The complex representations of the Fourier transform of the scattered wave and the plane reference wave in the recording plane may be denoted *U _{S}*(

*x*,

*y*)=

*U*

_{0S}(

*x*,

*y*)exp[-

*jφ*(

_{S}*x*,

*y*)] and

*U*(

_{R}*x*,

*y*)=

*U*

_{0R}exp[-

*jφ*

*(*

_{R}*x*,

*y*)], respectively, where

*U*

_{0S}(

*x*,

*y*) is a real (positive) spatially-varying amplitude, and

*U*

_{0R}is a constant amplitude (assumed real and positive, without loss of generality). The term

*φ*(

_{S}*x*,

*y*) is the phase of the Fourier transform of the sample wave,

*φ*(

_{R}*x*,

*y*)=

*k*sin

*θ*(

_{r}*x*cos

*ϕ*+

_{r}*y*sin

*ϕ*) is the phase of the reference wave, linear with respect to spatial position,

_{r}*k*=2

*π*/λ represents optical wavenumber, and λ is the wavelength of the source. The angle

*θ*is the incident angle of the reference wave (relative to the

_{r}*z*-axis, which is normal to the Fourier plane) and

*ϕ*is the azimuthal angle, i.e., the angle between the reference plane of incidence and the

_{r}*x*coordinate axis.

The recorded intensity distribution is:

When this intensity is inverse Fourier-transformed, we obtain

$$\phantom{\rule{.9em}{0ex}}\phantom{\rule{.9em}{0ex}}\phantom{\rule{.9em}{0ex}}\phantom{\rule{.9em}{0ex}}\phantom{\rule{.9em}{0ex}}\phantom{\rule{.9em}{0ex}}+{U}_{0R}{u}_{S}^{*}(-{v}_{x}+\frac{\mathrm{sin}{\theta}_{r}\mathrm{cos}{\varphi}_{r}}{\lambda},-{v}_{y}+\frac{\mathrm{sin}{\theta}_{r}\mathrm{sin}{\varphi}_{r}}{\lambda}),$$

where 𝓕^{-1} denotes the inverse Fourier-transform operator, defined by the equation:

and (*ν _{x}*,

*ν*) represent the coordinates in the transform space. In Eq. (2),

_{y}*u*is the inverse Fourier transform of

_{S}*U*, and

_{S}*δ*represents the Dirac delta function. The first term on the right-hand side of the equation, Γ

_{u}(

*ν*,

_{x}*ν*)=

_{y}*u*(

_{S}*ν*,

_{x}*ν*)⊗

_{y}*u**

_{S}(-

*ν*, -

_{x}*ν*), is a zero-order autocorrelation and the second term is a high-intensity zero-order spot, located at the origin. The third and fourth terms are first-order twin images of the sample field and its (spatially inverted) complex conjugate, each translated in opposite directions due to the influence of the exponential carrier factors. If the first-order images are sufficiently spatially separated from the zero-order images, the squared magnitude of the entire transformed distribution yields reconstructed twin images of the scattered power from the sample. Otherwise, the first two terms can be removed by recording the sample and reference intensities separately, and subtracting both from Eq. (1).

_{y}The reconstructed field *u _{S}*(

*ν*,

_{x}*ν*) is a scaled version of the scattered field distribution in the object plane. We denote the latter distribution by

_{y}*V*(

_{O}*ξ*,

*η*) (utilizing the input plane coordinate system), such that:

where *M* is a constant (with dimensions of squared length) dependent on the optical elements of the setup and proportional to the illumination wavelength.

An important issue to consider in reconstructing the sample field distribution is that of defocusing, in both the object and recording planes. The recording plane may be deliberately defocused (i.e., located in a plane which is *not* conjugate to the back focal plane of the Fouriertransforming objective lens) when imaging strongly diffracting sample structures. This ensures that any tightly-focused bright spots in the recording plane are spread out over many pixels, improving signal-to-noise ratio, and avoiding detector saturation. For the present purposes, recording plane defocusing is unnecessary, since the detected sample wave is generated by scattering processes. However, for each recorded hologram, there remains the possibility that the sample plane is defocused (not located at the front focal plane of the objective lens).Within the assumption of Fresnel diffraction theory, we can correct for this offset by multiplying the recording plane field *U _{S}* by a quadratic phase factor, as follows. If the sample is located at a distance

*g*in front of the front focal plane, then the field in the front focal plane

*V*is related to the field in the sample plane

_{O}*V*via the convolution [30, p. 67]:

_{S}where:

Thus, *U _{S}*, the detected wave in the absence of defocus correction, is related to the defocuscorrected wave

*U*,

_{S}_{dc}via the equation:

where *K _{g}*, the Fourier transform of

*k*, may be represented:

_{g}If the recorded intensity distribution *I*(*x*,*y*) is divided by the quadratic correction factor, we obtain (noting that the inverse of *K _{g}* is equal to its complex conjugate):

where the arguments of *K _{g}* are taken to be (-

*x*/

*M*,-

*y*/

*M*) as in Eq. (7). It is clear from Eq. (9) that, after the recorded hologram is multiplied by the quadratic defocus correction factor, only one of the twin reconstructed images will be focused; the other will be doubly defocused.

The spatial (angular) frequency of the fringe pattern formed due to the interference between the sample and reference waves in the recording plane is calculated by projecting the difference between the sample and reference wavevectors, **k**
_{R} and **k**
* _{S}*, respectively, onto the recording plane. For an axial sample wave, this difference, the fringe vector

**k**

*=*

_{F}**k**

*-*

_{R}**k**

*, is depicted along with*

_{S}**k**

_{R}and

**k**

*in Fig. 1. The wavevectors both have magnitude*

_{S}*k*and

**k**

*has magnitude*

_{F}*k*=2

_{F}*k*sin(

*θ*/2). The fringe vector is oriented at angle

_{r}*θ*/2 to the recording plane, thus, its projection onto the plane has magnitude

_{r}*k*cos(

_{F}*θ*/2). The spatial period of the fringe pattern is therefore:

_{r}It should be noted that in Fourier holography, this fringe period is not dependent on the scattering angle of the light fromthe sample; instead, it depends only on the location of the scatterer in the sample plane. This is an important advantage of Fourier holography over other techniques, endowing it with the ability to record a large range of scattering angles on a low-resolution CCD sensor.

For a sample composed of discrete scatterers, angular scattering spectroscopy measurements may be sufficient to ascertain the scatterers’ sizes and their relative refractive indices with respect to the background medium. In the Fourier plane, there is a one-to-one correspondence between spatial position and scattered light direction. To determine the angular scattering distributions for selected sample regions, we have adopted two different (but essentially equivalent) techniques. In the first, we perform Fourier filtering using a selection of masks in the recording plane, each of which specifies (and weights) a particular scattering solid angle range. After filtering, a reconstructed image of the entire sample is formed only from the light scattered within that range. For each area of the sample, the variation in the reconstructed power corresponding to different mask positions can be used to determine sample structural characteristics. This technique, demonstrated in Refs. [2, 3], is useful when it is not convenient to select (manually or automatically) *a priori* sample regions of interest. In the second technique, we select such sample regions from the reconstructed object plane field distribution (where all detected scattering angles are utilized to form the reconstruction). The field within each region is Fourier-transformed to obtain a map of the scattering angle distribution. The principal advantage of this technique is the fact that it allows direct access to these two-dimensional maps. The distribution of the scattered power in each may be used to determine sample properties within the selected region.

The fundamental resolution limit of our method, which holds for both approaches above, is governed by the Fourier uncertainty relationship between angular resolution in the recording plane and spatial resolution in the sample plane. Assuming a (one-dimensional) Gaussian filter mask profile in the recording plane of full-width-at-half-maximum (FHWM) Δ*w*, then according to Eq. (4), the spatial resolution in the reconstruction plane is limited by a Gaussian coherent point-spread function FWHM of Δ*d*=4(ln2)*M*/(*φ*Δ*w*). In the paraxial approximation, the deflection angle *θ _{d}* (from the optical axis) is related to the focal point deflection distance in the recording plane

*r*by the equation

_{d}*θ*=λ

_{d}*r*/

_{d}*M*. Thus, the scattering angle resolution (FWHM in air) Δ

*θ*is related to spatial resolution in the reconstruction plane Δ

*d*by the equation:

It is clear from Eq. (11) that the only system parameter upon which the trade-off depends is the optical wavelength. For the red He-Ne laser wavelength *λ*=632.8nm and a FWHM angular resolution of Δ*θ*=1.0°, this yields a sample spatial resolution δ*d*=32*µ*m. Despite the fact that this value is clearly too large to directly resolve sample microstructure, it does not represent an impediment in our approach.

## 3. Experimental setup and procedure

Figure 2 presents a schematic diagram of our experimental setup. A light beam from a coherent source S (He-Ne laser) is split into reference and sample waves at the beamsplitter B1. A sample is placed at the object plane of the Fourier-transforming objective lens L1 (focal length *f*
_{1}=15mm) and illuminated by a plane wave. Its optical Fourier spectrum is imaged onto the recording plane via the lenses L2 and L3 (of focal lengths *f*
_{2}=15cm and *f*
_{3}=26cm, respectively). For this setup, the scaling constant *M*=*λ*
*f*
_{1}
*f*
_{3}/*f*
_{2}. Recording is performed using a charge-coupled device (CCD) matrix sensor (12 bit, 1392×1040 pixels, pixel length Δ*r*=4.65*µ*m). The reference wave is expanded using the telescopic system T and is directed offaxis onto the CCD matrix at angle *θ _{r}* of approximately 2.3°.

The direction of the light scattered by the sample is described by a polar (scattering) angle *θ _{s}*, the deviation from the forward direction, and an azimuthal angle. Its angular deviation

*θ*from the axis depends on both of these angles, as well as the direction of the illumination wave (described by

_{d}*θ*,

_{i}*ϕ*) and the refractive index of the sample medium

_{i}*n*

_{med}. The inset to Fig. 2 illustrates this relationship in the plane of incidence. Restricting our consideration to this plane, we may let

*θ*take both positive and negative values with respect to the optical axis, with the direction shown corresponding to a positive value. Within this plane,

_{d}*θ*is related to

_{s}*θ*via the equation:

_{d}The range of angles captured by the optical system is limited by the numerical aperture (NA) of the objective lens L1, the diameters and focal lengths of lenses L2 and L3, and the pixel size and count of the CCD sensor. In our case, the size of the CCD sensor is (for the most part) the limiting component of the system. For our experimental setup, in which *θ _{i}*≅49°, the angular deviation

*θ*range is about 17° (in air), centered on the axis. If the sample background medium is water (

_{d}*n*

_{med}=1.33), the scattering angles

*θ*which can be detected range from 139° to 152°. A rectangular field stop is placed in a plane conjugate to the sample, in order to restrict the field of view to a 1mm×2mm area, thus enabling clear identification of the first-order twin images.

_{s}The precise relationship between recording plane location and scattering direction can be determined. For our setup, the *ξ*, *η*, *x*, and *y* axes, depicted in Fig. 1, are aligned with the horizontal and vertical axes of the CCD recording area, respectively, and the azimuthal angle *ϕ*
_{i}≅35°. For our system parameters, curves in the recording plane corresponding to constant scattering angle *θ _{s}* (but varying azimuthal angle) are displayed in the left panel of Fig. 3. The curves are well approximated by straight lines perpendicular to the illumination-wave plane of incidence (in the paraxial approximation), i.e., oriented at angle

*ϕ*to the

_{i}*y*-axis. (The angular error associated with this straight-line approximation is less than 0.2° over virtually the entire recording plane.) The distance between two such lines corresponding to a scattering angle difference of

*θ*

*s*_{,diff}is approximately

*x̃*

_{diff}=(

*Mn*

_{med}/

*λ*)

*θ*

_{s}_{,diff}.

Vignetting due to the optical components between the sample and recording planes can limit the sample field of view. We assume that the objective lens most severely limits the range of scattered rays that can be detected by the system. Then the region of the Fourier spectrum (in the recording plane) contributed to by a scatterer at point (*ξ _{s}*,

*η*) in the sample plane is a circle centered at the point (

_{s}*x*,

_{s}*y*)=(

_{s}*f*

_{3}/

*f*

_{2})(

*ξ*,

_{s}*η*) of diameter (

_{s}*f*

_{3}/

*f*

_{2})

*D*, where

_{L}*D*is the diameter of the objective lens. Consider the sample area for which the CCD recording area lies entirely within this region. If the CCD recording area was circular with diameter

_{L}*D*, then this area would also be a circle with diameter

_{R}*D*=

_{S}*D*-(

_{L}*f*

_{2}/

*f*

_{3})

*D*(if

_{R}*D*>(

_{L}*f*

_{2}/

*f*

_{3})

*D*, and zero otherwise). Of course, for our setup, the recording area is rectangular, with dimensions 4.8×6.5 mm. The corresponding sample areas for a range of different objective lens diameters are shown in the right panel of Fig. 3. For our case, the objective lens diameter

_{R}*D*≅ 6mm, so that its NA was 0.2.

_{L}Let *d _{s}* be the diameter of the sample (or an alternative representative length). Then to ensure that the twin reconstructed images are spatially separated from the zero-order terms, it is necessary that [30, p. 309]:

To satisfy sampling restrictions due to the finite pixel size in the recording plane, it is necessary that *H _{f}*>2(D

*r*/√2), where the factor √2 arises due to the effective pixel size in a diagonal direction, and the factor 2 arises from Nyquist’s theorem, i.e., by Eq. (10):

For a 1-mm-diameter sample, Eqs. (13) and (14) imply that, for our system, it is necessary that 3.3°≤*θ _{r}*<5.5°. Our choice of

*θ*=2.3° clearly violates the lower limit, so there will be some overlap between the zero-order terms and the twin images. As explained above, however, this can be overcome by digitally subtracting the recorded reference and sample waves from the hologram. Equation (13) demonstrates that, for a given CCD sensor size and wavelength, if

_{r}*d*is increased, then

_{s}*M*must undergo a corresponding increase to satisfy the inequality. That is, there is a trade-off between measured scattering angle range and sample size.

We selected samples based on their scattering profiles in the range of scattering angles *θ _{s}* detectable by our system, from about 139° to 152° as described above. The choice of microsphere suspensions was natural since their angular scattering pattern follows a distinctive modulation (ripple) pattern described by Mie theory, with an (angular) period which (for the most part) decreases with increasing sphere diameter. Our samples comprised polysterene spheres suspended in distilled water. The microspheres were diluted to a volume concentration of 0.1% and a droplet was deposited into a 10×20 mm well on a microscope slide. To demonstrate the application of our approach to biological samples, we utilized a smear of erythrocytes, or red blood cells (RBCs). The RBCs were diluted with a droplet of NaCl solution (9%) and evenly smeared over a microscope slide. A coverslip was placed over the sample, and sealed at the edges. The normal RBC shape is a discocyte, an axially-symmetric disc indented on the axis [31]. This shape produces a distinctive angular scattering ripple pattern [31] which should be clearly observable over our angle range.

For spherical scatterers, the angular scattering distribution depends on scattering angle but not significantly on azimuthal angle, and this fact was utilized when processing the images. When applying the first technique of Section 2, the recording plane masks were chosen to be strips perpendicular to the dotted line in Fig. 3(a). Each thereby corresponded to a single scattering angle (and the full range of azimuthal angles). (The filters had a Gaussian cross-section.) By applying a range of such masks to the recorded image, it was possible to record the total power scattered at each scattering angle.We were thus able to generate a one-dimensional curve of scattered power vs. scattering angle for each selected sample region [2]. When applying the second technique, we sought to develop scatterer identification methods that were sensitive to structure parallel to the dotted line of Fig. 3(a) (in the scattering-angle-distribution map for each selected region), the direction of varying scattering angle.

## 4. Results

The scattered power angular distributions predicted by Mie theory are shown as blue curves in Fig. 4, assuming a refractive index ratio between sphere and background medium of *m*=*n*
_{sph}/*n*
_{med}=1.59/1.33=1.20, for a range of particle sizes. The ripple patterns are readily apparent, as well as the dependence of their angular periods on sphere size. It is this parameter of the curves that we utilize to determine sphere size in our samples. To assign a value to the apparent angular period in each case, a high-pass filter was applied to the curve in order to fit it accurately to a sinusoidal curve. The angular period was originally estimated from the average fringe spacing over the angle range 120°-170°, and the (4th-order Butterworth) filter cutoff frequency was equal to 0.75 times the estimated fringe frequency. The filtered curves are shown in red. They were fitted (in a minimum least-squares sense) to sinusoidal curves, and the frequencies recorded. By this process, it may be shown that the ripple angular frequency is almost linearly dependent on the Mie size parameter *α*=*πdn*
_{med}/*λ*, where *d* is the sphere diameter. Such an approximation is valid (at least) over the refractive index ratio range *m*=1.1-1.25, yielding a maximum error in detected diameter of less than 1*µ*m (for sphere sizes ranging from 1 to 20*µ*m) and a mean error of about 0.2*µ*m. A more direct Mie inversion procedure may be used for particle size/refractive index values outside these ranges. Thus, particle size can be recovered from the measured ripple angular frequency, to a degree of accuracy which should be sufficient for many applications. The minimum particle size measurable by our system can be estimated by determining the sphere diameter for which one full ripple cycle is visible over the angular range used. As is clear from Fig. 4, this minimum size is about 2*µ*m, assuming *m* is within the given range. In general, the particle-size sensitivity of our approach is limited by angular range and refractive index to the same extent as alternative angular-scattering-spectroscopy techniques [14]. Goniometric measurements published in the literature have clearly shown multiple Mie-theory ripples associated with sphere diameters of less than 1 or 2*µ*m [20, 22, 32].

Figure 5 demonstrates the application of Method 1 of Section 2 to a sample of 11.4-*µ*m spheres in water. The apparent brightness of the scattering regions of the sample (right-hand side) is clearly dependent on the recording plane strip-mask position (shown on the left-hand side) and, thus, on scattering angle.

Figures 6 and 7 show reconstructed images of samples of suspensions of 5.4 and 11.4-*µ*m spheres in water, respectively. In each case, five regions of the sample were selected, and the angular scattering (power) spectra corresponding to each is displayed. Note that for each, the speckle size in the two-dimensional spectra (and thus the angular resolution) appears to be inversely related to the size of the selected region (outlined in red), an observation consistent with the trade-off represented by Eq. (11). A two-dimensional (inverse) Fourier transform operation was applied to each power spectrum, and the results also displayed, with a dotted yellow line indicating spatial frequencies in the direction of varying scattering angle. Bright spots along this line correspond to detected ripple frequencies in the angular scattering spectra, and by detecting the peak position (indicated with a magenta cross), we can identify each region with a detected sphere diameter by invoking the linear relationship described in the first paragraph of this section. (Diffraction effects due to the shape of the recording array and low-spatial-frequency noise were mitigated by subtracting a CCD-sensor-shaped rectangle from the spectra before inverse Fourier-transforming them, so that their mean value was zero, and applying a radial square-root-profile mask to the Fourier-transformed spatial frequency spectra before peak detection.) The same process was applied to every region of the sample, and the detected sphere size in each displayed using a false color scale. The color scale shown indicates sphere size in micrometers.

The mean diameter of the detected spheres (ignoring spurious outliers) was 4.8*µ*m and 9.9*µ*m (with standard deviations 0.9*µ*m and 0.5*µ*m), respectively. The systematic error of 10–15% is due possibly to errors in the system scaling constant *M* (due to the large tolerances of the optical components used), leading to a smaller range of scattering angles being imaged onto the CCD detector than predicted theoretically. The relatively large variation in the values measured in the former case is due to the difficulty in precisely detecting the period of low-frequency fringes using a detector encompassing a limited angular range. In both cases, the scatterers on the far left of the reconstruction take on a ‘streaked’ appearance. This is due to vignetting at the boundary of the sample region; the scatterers are reconstructed using a reduced range of spatial frequencies, so they exhibit a resolution loss (and apparent broadening) in one direction.

Figure 8 presents a reconstruction and a false-color map showing the detected size distributions of a sample containing a mixture of both 5.4 and 11.4-*µ*m spheres. The two particle sizes are clearly distinguishable by their distinct hues used in the latter representation. Six different regions, three of each particle size, are highlighted as before. The detected sphere sizes are consistent with those measured in the previous two figures.

The reconstructed image and scattering spectra from five selected regions of a hologram of the red-blood-cells sample are shown in Fig. 9. The orientation of the incident beam with respect to the recording plane was different from the previous measurements. (In this case, *ϕ _{i}*≅140°.) Thus, the angle of the lines of constant scattering angle has been likewise varied. Clear ripple structure is apparent in the recorded power spectra (and highly visible peaks along the yellow dotted line in their inverse Fourier transforms), corresponding to a uniform spatial frequency. Since red blood cells are not spherical, it is not strictly appropriate to utilize Mie theory to determine their sizes. Nonetheless, for the purpose of comparison with the previous results, the same Mie inversion procedure was applied to the sample, and a false-color map of particle size generated, as before. (The fact that the relative refractive index

*m*between the blood cells and their background was outside the range specified earlier in this section was ignored for the purpose of this simple analysis.) For this sample, despite the fact that a majority (65%) of detected particle sizes were within 1

*µ*m of themean value (ignoring spurious outliers), there was much more variation in the detected blood cell sizes than in the sphere sizes of the earlier experiments. This is probably accounted for by the natural variation in the size and orientation of the particles, as well as the general inaccuracies inherent in applying aMie-theory inversion procedure to a distribution of non-spherical particles. However, the mean detected size of about 6

*µ*m correlates well with typical red-blood-cell sizes reported in the literature [31]. It is evident in the reconstructed image (as it had been in Figs. 6–8) that the focused regions corresponding to individual red blood cells take the form of a pair of closely-spaced bright points. This is unsurprising, since the far-field interference pattern generated by a pair of coherent point sources is a fringe pattern with a similar appearance to the two-dimensional spectra generated in these figures. We should therefore expect the focused scatterer images to be of this form when highly visible ripple patterns are present in the angular scattering spectra.

## 5. Discussion and Conclusion

The results presented in this paper clearly demonstrate the measurement of spatially resolved angular scattering spectra, and their use in estimating microscopic sphere sizes over large fields of view. We have confirmed that the spatial resolution required to resolve the sample’s microstructure is not necessary to characterize it. The theoretical (diffraction-limited) spatial/angular resolution trade-off for the present setup (using a Gaussian-profile spatial filter mask) in the direction of varying scattering angle may be determined using Eq. (11). A representative pair of values is Δ*θ*=1.3° (in air), and Δ*d*=25*µ*m. Utilizing the full mask length to reconstruct the sample, the angular/spatial resolutions in the direction of varying azimuthal angle are given by 8.9° (rectangular length)/6*µ*m (main lobe length). In principle, an image can be formed from a single exposure, in common with conventional microscopy, however, since we do not require spatial resolutions high enough to measure microparticle sizes directly, our system optical requirements are very modest. We can utilize a low-magnification, low-numerical-aperture objective, which allows us to form images over long working distances and millimeter-scale fields of view. This ability to form millimeter-scale images that provide the angular distribution of the scattered light in each local area is unique to our method.

Digitally recording the complex Fourier transformation of the scattered sample wave provides great flexibility for applying spatial filtering and, potentially, optical pattern recognition methods to select specific information from a sample. Besides Mie theory, other models to describe scattered light can be utilized in our technique. For example, for samples such as connective tissue that comprise non-spherical scatterers, alternative methods will be needed to quantify scatterer parameters.We are encouraged, in this regard, by the demonstrated sensitivity of light scattering to even small changes in the structure of a sample [21, 33]. In many instances, detection of change (e.g., differentiating normal from abnormal tissue structures) in itself will be sufficient to provide valuable capability. In our experiments, we used only the angular frequency of the Mie ripples to detect particle size. This parameter alone has low sensitivity to refractive index ratios (between spheres and background) over the range *m*=1.1-1.25. If we were able to improve our processing procedure to consider other parameters such as peak position or relative magnitude, our approach could have even greater discriminatory powers. Curiously, even for samples composed of identical particles (as in Figs. 6, 7), the peaks in the one-dimensional scattering spectra were not aligned over the five images. This observation is partially accounted for by the variation in the size of the particles used (the manufacturerspecified standard deviation was 0.14 and 0.21*µ*m, respectively, for the 5.4 and 11.4*µ*m sphere sizes; diameter variations of about 0.3 and 0.2*µ*m, respectively, would be sufficient to account for complete contrast reversal of the Mie ripples). Also, any variation in the particle shape from sphericity could account for these anomalies.

The accuracy of the technique could be improved by collecting scattered light over a larger solid angle range, and we are currently working to achieve this (see, for example, Ref. [13]). The backscattering geometry that we use has the advantage of not detecting specularly reflected or undeflected incident light, but our scattering-angle range of sensitivity could in this way be extended three- or four-fold.

Spatially-resolved Fourier-holographic angular scattering spectroscopy has great potential for application to the study of the microarchitecture of biological tissue. The technique could provide maps of the size ranges of key scatterers, including cells, cell nuclei, and organelles, as well as the fibrous components of tissues such as collagen and elastin matrices.

## Acknowledgments

SAA is supported by the Raine Medical Research Foundation.We gratefully acknowledge contributions to the work presented here by Guido Kudielka,Michael Same, and Abbey Trewenack.

## References and links

**1. **G. Sluder and D. E. Wolf, eds., *Digital Microscopy: a Second Edition of Video Microscopy* (Elsevier Academic Press, Amsterdam, 2003).

**2. **S. A. Alexandrov, T. R. Hillman, and D. D. Sampson, “Spatially resolved Fourier holographic light scattering angular spectroscopy,” Opt. Lett. **30**, 3305–3307 (2005). [CrossRef]

**3. **S. A. Alexandrov, T. R. Hillman, T. Gutzler, M. B. Same, and D. D. Sampson, “Particle sizing with spatiallyresolved Fourier-holographic light scattering angular spectroscopy,” in *BiOS 2006:Multimodal Biomedical Imaging*,
F. S. AzarD. N. Metaxas
, eds., Proc. SPIE6081, 608104 (2006). [CrossRef]

**4. **L. Yaroslavsky, *Principles of Digital Holography* (Kluwer Academic Publishers, Boston, 2003).

**5. **E. Cuche, F. Bevilacqua, and C. Depeursinge, “Digital holography for quantitative phase-contrast imaging,” Opt. Lett. **24**, 291–293 (1999). [CrossRef]

**6. **E. Cuche, P. Marquet, and C. Depeursinge, “Simultaneous amplitude-contrast and quantitative phase-contrast microscopy by numerical reconstruction of Fresnel off-axis holograms,” Appl. Opt. **38**, 6994–7001 (1999). [CrossRef]

**7. **G. Popescu, L. P. Deflores, J. C. Vaughan, K. Badizadegan, H. Iwai, R. R. Dasari, and M. S. Feld, “Fourier phase microscopy for investigation of biological structures and dynamics,” Opt. Lett. **29**, 2503–2505 (2004). [CrossRef] [PubMed]

**8. **P. Marquet, B. Rappaz, P. Magistretti, E. Cuche, Y. Emery, T. Colomb, and C. Depeursinge, “Digital holographic microscopy: a non-invasive contrast imaging technique allowing quantitative visualisation of living cells with subwavelength axial accuracy,” Opt. Lett. **30**, 468–470 (2005). [CrossRef] [PubMed]

**9. **M. Gustafsson and M. Sebesta, “Refractometry of microscopic objects with digital holography,” Appl. Opt. **43**, 4796–4801 (2004). [CrossRef] [PubMed]

**10. **M. Sebesta and M Gustafsson, “Object characterization with refractometric digital Fourier holography,” Opt. Lett. **30**, 471–473 (2005). [CrossRef] [PubMed]

**11. **B. Javidi, I. Moon, S. Yeom, and E. Carapezza, “Three-dimensional imaging and recognition of microorganism using single-exposure on-line (SEOL) digital holography,” Opt. Express **13**, 4492–4506 (2005). [CrossRef] [PubMed]

**12. **V. Mico, Z. Zalevsky, and J. Garcia, “Superresolution optical system by common-path interferometry,” Opt. Express **14**, 5168–5177 (2006). [CrossRef] [PubMed]

**13. **S. A. Alexandrov, T. R. Hillman, T. Gutzler, and D. D. Sampson, “Synthetic aperture Fourier holographic optical microscopy,” Phys. Rev. Lett. **97**, 168102 (2006). [CrossRef] [PubMed]

**14. **J. R. Mourant, T. M. Johnson, S. Carpenter, A. Guerra, T. Aida, and J. P. Freyer, “Polarized angular dependent spectroscopy of epithelial cells and epithelial cell nuclei to determine the size scale of scattering structures,” J. Biomed. Opt. **7**, 378–387 (2002). [CrossRef] [PubMed]

**15. **R. S. Gurjar, V. Backman, L. T. Perelman, I. Georgakoudi, K. Badizadegan, I. Itzkan, R. R. Dasari, and M. S. Feld, “Imaging human epithelial properties with polarized light-scattering spectroscopy,” Nature Med. **7**, 1245–1248 (2001). [CrossRef] [PubMed]

**16. **V. Backman, M. B. Wallace, L. T. Perelman, J. T. Arendt, R. Gurjar, M. G. Muller, Q. Zhang, G. Zonios, E. Kline, T. McGillican, S. Shapshay, T. Valdez, K. Badizadegan, J. M. Crawford, M. Fitzmaurice, S. Kabani, H. S. Levin, M. Seiler, R. R. Dasari, I. Itzkan, J. Van Dam, and M. S. Feld, “Detection of preinvasive cancer cells,” Nature (London) **406**, 35–36 (2000). [CrossRef] [PubMed]

**17. **K. Sokolov, R. Drezek, K. Gossage, and R. Richards-Kortum, “Reflectance spectroscopy with polarized light: Is it sensitive to cellular and nuclear morphology,” Opt. Express **5**, 302–317 (1999). [CrossRef] [PubMed]

**18. **H. Fang, M. Ollero, E. Vitkin, L. M. Kimerer, P. B. Cipolloni, M. M. Zaman, S. D. Freedman, I. J. Bigio, I. Itzkan, E. B. Hanlon, and L. T. Perelman, “Noninvasive sizing of subcellular organelles with light scattering spectroscopy,” IEEE J. Sel. Top. Quantum Electron. **9**, 267–276 (2003). [CrossRef]

**19. **A. Katz, A. Alimova, M. Xu, E. Rudolph, M. K. Shah, H. E. Savage, R. B. Rosen, S. A. McCormick, and R. R. Alfano, “Bacteria size determination by elastic light scattering,” IEEE J. Sel. Top. Quantum Electron. **9**, 277–287 (2003). [CrossRef]

**20. **W. P. Van De Merwe, J. Czégé,M. E Milham, and B. V. Bronk, “Rapid optically based measurements of diameter and length for spherical or rod-shaped bacteria *in vivo*,” Appl. Opt. **43**, 5295–5302 (2004). [CrossRef] [PubMed]

**21. **N. N. Boustany, S. C. Kuo, and N. V. Thakor, “Optical scatter imaging: subcellular morphometry *in situ* with Fourier filtering,” Opt. Lett. **26**, 1063–1065 (2001). [CrossRef]

**22. **J. D. Wilson, C. E. Bigelow, D. J. Calkins, and T. H. Foster, “Light scattering from intact cells reports oxidativestress-induced mitochondrial swelling,” Biophys. J. **88**, 2929–2939 (2005). [CrossRef] [PubMed]

**23. **J. D. Wilson and T. H. Foster, “Mie theory interpretations of light scattering from intact cells,” Opt. Lett. **30**, 2442–2444 (2005). [CrossRef] [PubMed]

**24. **Y. L. Kim, Y. Liu, R. K Wali, H. K. Roy, M. J. Goldberg, A. K. Kromin, K. Chen, and V. Backman, “Simultaneous measurement of angular and spectral properties of light scattering for characterization of tissue microarchitecture and its alteration in early precancer,” IEEE J. Sel. Top. Quantum Electron. **9**, 243–256 (2003). [CrossRef]

**25. **Y. Liu, Y. L. Kim, and V. Backman, “Development of a bioengineered tissue model and its application in the investigation of the depth selectivity of polarization gating,” Appl. Opt. **44**, 2288–2299 (2005). [CrossRef] [PubMed]

**26. **M. Bartlett, G. Huang, L. Larcom, and H. Jiang, “Measurement of particle size distribution in mammalian cells *in vitro* by use of polarized light spectroscopy,” Appl. Opt. **43**, 1296–1307 (2004). [CrossRef] [PubMed]

**27. **Y. Liu, Y. L. Kim, X. Li, and V. Backman, “Investigation of depth selectivity of polarization gating for tissue characterization,” Opt. Express **13**, 601–611 (2005). [CrossRef] [PubMed]

**28. **A. Wax, C. Yang, R. R. Dasari, and M. S. Feld, “Measurement of angular distributions by use of low-coherence interferometry for light-scattering spectroscopy,” Opt. Lett. **26**, 322–324 (2001). [CrossRef]

**29. **R. N. Graf and A. Wax, “Nuclear morphology measurements using Fourier domain low coherence interferometry,” Opt. Express **13**, 4693–4698 (2005). [CrossRef]

**30. **J. Goodman, *Introduction to Fourier optics* (McGraw-Hill, New York, 1996, 2nd ed.).

**31. **A. Karlsson, J. He, J. Swartling, and S. Andersson-Engels, “Numerical simulations of light scattering by red blood cells,” IEEE Trans. Biomed. Eng. **52**, 13–18 (2005). [CrossRef] [PubMed]

**32. **J. C. Ramella-Roman, P. R. Bargo, S. A. Prahl, and S. L. Jacques, “Evaluation of spherical particle sizes with an asymmetric illumination microscope,” IEEE J. Sel. Top. Quantum Electron. **9**, 301–306 (2003). [CrossRef]

**33. **M. T. Valentine, A. K. Popp, D. A. Weitz, and P. D. Kaplan, “Microscope-based static light-scattering instrument,” Opt. Lett. **26**, 890–892 (2001). [CrossRef]