- Open Access
Quantitative phase imaging with scanning holographic microscopy: an experimental assesment
© Indebetouw et al; licensee BioMed Central Ltd. 2006
- Received: 19 September 2006
- Accepted: 28 November 2006
- Published: 28 November 2006
This paper demonstrates experimentally how quantitative phase information can be obtained in scanning holographic microscopy. Scanning holography can operate in both coherent and incoherent modes, simultaneously if desired, with different detector geometries. A spatially integrating detector provides an incoherent hologram of the object's intensity distribution (absorption and/or fluorescence, for example), while a point detector in a conjugate plane of the pupil provides a coherent hologram of the object's complex amplitude, from which a quantitative measure of its phase distribution can be extracted. The possibility of capturing simultaneously holograms of three-dimensional specimens, leading to three-dimensional reconstructions with absorption contrast, reflectance contrast, fluorescence contrast, as was previously demonstrated, and quantitative phase contrast, as shown here for the first time, opens up new avenues for multimodal imaging in biological studies.
- Spherical Wave
- Quantitative Phase
- Digital Holography
- Atomic Force Microscope Data
- Digital Holographic Microscopy
Microscopy is an essential tool in biological research, micromechanical testing, the integrated circuit industry, etc. The demand for higher resolution and contrast, shorter acquisition time, and multimodal imaging, among other desirable properties, has resulted in the recent invention, demonstration, and often rapid commercialization of a number of new technologies. In biological studies, two modalities appear to be of primary importance. They are fluorescence imaging for the specific identification of biomolecules in a labeled sample, and phase imaging for the determination of internal structures in unstained specimens. The conventional phase imaging methods (i.e. Zernike phase contrast, and Nomarski differential interference contrast) usually provide only the visualization of the phase of biological structures in a qualitative way, although it is possible to extract quantitative phase information with the differential interference contrast method . Recently, quantitative phase imaging, as provided by digital holographic microscopy  gave a new dimension to phase imaging by allowing the quantitative measurement of, for example, biomasses and fluid concentrations in cells. Quantitative phase imaging is also essential in measuring non-destructively the surface topography of biological samples, as well as of micromechanical systems, and integrated electrical circuits.
To our knowledge, no single instrument or imaging method can capture both fluorescence and quantitative phase information of 3D specimens simultaneously in holographic form. In this paper, we demonstrate experimentally that scanning holographic microscopy [3, 4], which was developed to obtain holograms of incoherent objects [5, 6], and has recently been shown to provide high resolution images of 3D fluorescent biological specimens , is also capable of providing quantitative phase information. The possibility of phase imaging by scanning holographic microscopy had been suggested earlier on theoretical grounds , but has not yet been demonstrated experimentally. This possibility opens up new potentials for multimodal imaging. For example, it would be possible to obtain, simultaneously if desired, absorption images carrying information of structural arrangements, fluorescence images of labeled specimens revealing functional activities, and quantitative phase images from which internal structures, biomasses, density, etc. could be measured.
The paper is organized as follows. Section two is a brief theoretical review of the two possible operation modes of scanning holography: the incoherent mode giving holograms of the object's intensity distribution (absorption, reflection, and fluorescence), and the coherent mode giving holograms of the object's complex amplitude distribution, and providing a quantitative measure of its phase distribution. Section three describes the experimental set up, and section 4 presents experimental results of phase images of unstained biological specimens, as well as of the quantitative profile of manufactured phase objects. The latter result is compared quantitatively with atomic force microscope measurements. Section five is a brief summary.
are the generalized defocused pupils . F-1 stands for inverse Fourier transform, and Ω is the frequency shift of one of the pupils. The transverse spatial frequency vector = P /λf 0 is proportional to the real space coordinate P in the pupil plane . λ is the wavelength of the illumination, and f 0 is the focal length of the objective. The two pupils used to obtain an in-line Fresnel hologram are, respectively, a spherical wave filling the pupil of the objective, and a point at the center of that pupil:
circ(x) is a disc function of unit radius. ρ MAX = sinα/λ is the cutoff frequency of the objective, where sinα = NA is its numerical aperture. The Fresnel pattern projected on the object is the interference of the Fourier transforms of the two pupils, namely a spherical wave and a plane wave in the paraxial approximation. The Fresnel number of this pattern is determined by the free parameter z 0, which is the distance from the objective's focal plane to the point where the spherical wave comes to a focus.
There are two possible modes of operation, depending on the detector geometry . With a spatially integrating detector, the resulting data is a convolution of the object's intensity distribution with the desired complex PSF (namely, a spherical wave with a radius of curvature z 0 + z, and a radius a = z 0sinα, in this case). This detection mode leads to a hologram from which the three-dimensional distribution of scattering intensity, absorption, and fluorescence intensity can be reconstructed. With a point detector at the center of a conjugate pupil plane, the resulting data is a convolution of the object's complex amplitude distribution with the same complex PSF. The reconstruction of this hologram gives the three-dimensional distribution of the specimen's complex amplitude transmittance. In particular, assuming that multiple scattering can be ignored, the phase of the reconstruction is a quantitative measure of the integrated optical path length through the specimen. These two modes of operation are similar to the usual coherent/incoherent imaging modes of a conventional system, which are obtained by using, respectively, a point source, or a large spatially incoherent source. In scanning holography, the detector size plays a similar role to that of the source size in conventional imaging.
For simplicity, let's assume an object with an amplitude transmittance T( , z). Note that for an incoherent or fluorescent object, the phase of the transmitted field is a random variable, and only the intensity I( , z) = |T( , z)|2 is measurable. For a quasi transparent object, the phase is equal to the integrated optical thickness of the object: Φ( ) = (2π/λ)∫dzn( , z), where n( , z) is the 3D distribution of refractive index. The amplitude distribution after the object is written as
A( , t) = ∫dzS( , z)T[ - S (t), z], (3)
where S (t) is the instantaneous position of the 2D raster scan.
The incoherent imaging mode is obtained with a spatially integrating detector, leading to a temporal signal proportional to the integrated intensity, ∫d2 r|A( , t)|2, which is stored in the computer. The data corresponding to each hologram line is cut from the signal, and band pass filtered to extract the term oscillating at Ω. The lines are then rearranged in a 2D format. The resulting hologram amplitude is found to be
H I ( ) = ∫dzI( , z) ⊕ [p 1( , z)p*2( , z)], (4)
where ⊕ symbolizes a convolution integral, I( , z) = |T( , z)|2, and p 1,2( , z) are the inverse Fourier transforms of the defocused pupil distributions (eq.1b). The superscript * stands for complex conjugate. The hologram is thus the convolution of the object intensity with a spherical wave, i.e. an in-line single-sideband Gabor hologram. In Fourier space, the hologram can be written as
where ⊕ symbolizes a correlation integral. With the pupils of eq.2, we find, in the paraxial approximation,
which is the Fourier transform of the Fresnel hologram of the object's intensity distribution.
The coherent imaging mode is obtained by using a pinhole or a point detector at the center of a conjugate plane of the pupil of the objective. This leads to the hologram amplitude
H C ( ) = ∫dz[T( , z) ⊕ p 1( , z)][T( , z) ⊕ p 2( , z)]*. (7)
Where again, ⊕ symbolizes a convolution product, and the superscript * stands for complex conjugate. In Fourier space,
With the pupils of eq.2, we find
Aside from an inconsequential complex constant (the first term under the integral), eq.9 is the Fourier transform of the Fresnel hologram of the object's complex amplitude distribution. Thus, the reconstruction of the hologram recorded in the coherent mode carries a quantitative measure of the object's phase distribution.
The experimental system sketched in fig. 1 uses a 20X objective with NA = 0.42 (Mitutoyo Plan Apo). The detectors are standard photomultiplier tubes (Hamamatsu). In the coherent mode, the signal is measured through a 10 μm pinhole placed in a conjugate plane of the objective's pupil. The specimen is mounted on a 2D piezo scanning stage (Physik Instrumente). The area scanned is 150 × 150 μm2 divided in 2000 lines with 75 nm line spacing. Each line is sampled at 10,000 samples per line. The data acquisition system is a Gage Scope with a sampling rate set at 105 samples per sec. The Fresnel pattern projected on the object has a diameter a ~ 50 μm, and a Fresnel number ~16. At the wavelength of 532 nm, this corresponds to a radius of curvature of the spherical wave z 0 ~ 75 μm in the focal plane of the objective, and an effective numerical aperture NA = sinα = a/z 0 ~ 0.34. The Rayleigh transverse resolution limit is thus λ/2NA ~ 0.8 μm, a number which has been verified in previous experiments .
The high resolution sampling is needed for two reasons: to insure an adequate representation of the hologram phase, and to satisfy the Shannon/Nyquist sampling theorem in the signal demodulation. The sampling theorem requires a minimum of two samples per resolution element. Thus, a single intensity image with size 150 × 150 μm2, and 0.8 μm resolution is adequately represented by an array of 375 × 375 samples. To represent a complex hologram of the same size, however, a finer sampling is needed to capture the hologram phase. The width of the outermost Fresnel zone of the pattern scanning the specimen is equal to the resolution limit (0.8 μm), and represents a phase excursion of π. We found experimentally that more than five samples per resolution element were needed to capture this phase excursion with sufficient resolution. In the present experiment, we chose a sampling interval of 75 nm, corresponding to ~10 samples per resolution element, and a phase sampling interval ~π/10. Note that in standard digital holography, a high resolution pixellated detector (e. g. CCD array) is needed to capture the hologram phase via interferences in the spatial domain. In this case, the detector's spatial resolution is often a factor limiting the resolution. In contrast, the scanning holographic method captures the hologram phase via the temporal modulation of the scanning pattern integrated on a non-imaging detector. The resolution of digital holography is limited by the spatial sampling of the detector, while that of scanning holography is limited by the temporal sampling rate.
From the Shannon/Nyquist sampling theorem, the signal modulation frequency must be smaller than 1/4 the sampling rate, and the band pass of the filter demodulating the signal must be smaller than the modulation frequency. In our set up, the sampling rate is 105 Hz, the modulation frequency is 25 kHz, obtained with a ramp-driven electro-optic phase modulator (Linos), and the filter band pass is 20 kHz. Note that this demodulation process is done digitally in the computer, since the scanning effects a mapping from time to space, temporal frequencies are directly related to sample numbers. In the spatial domain, each line of the hologram is chosen to have 10,000 samples, corresponding to the sampling rate of 105 Hz. The Fourier transform of each line, which also has 10,000 samples, has two modulation bands centered at sample 2500, and 7500 (corresponding to the two frequency sidebands at ± 25 kHz). After adjusting the initial phase of each line to that of a reference signal (see fig. 1), one of the modulation band is extracted to form a hologram line of 2000 samples (corresponding to the band width of 20 kHz). The final hologram is then assembled in a 2D array of 2000 × 2000 samples. The hologram is reconstructed by standard Fresnel back-propagation using a Matlab code (see ref. 6, 7 for more details). In the present experiment, the acquisition time of one hologram is limited by the speed of the scanning stage to about 4 min. With a Pentium D-820 dual processor, it takes less than 10s to construct the 2D hologram, and less than 50s to reconstruct a stack of 10 axial images.
Quantitative phase imaging of unstained biological specimens and relief optical surfaces using scanning holographic microscopy has been demonstrated experimentally for the first time. Quantitative phase measurement of biological specimens was first demonstrated with digital holographic microscopy , which is a coherent imaging process. The advantage of the scanning holographic method used in this work is that it is possible to obtain holograms in both the coherent and the incoherent imaging modes simultaneously. This opens up new possibilities for multimodal imaging, making it possible, for example, to acquire absorption images, fluorescence images, and quantitative phase images of three-dimensional specimens simultaneously using the same scanning holographic microscope.
This research was financially supported by The National Institute of Health, Office of Extramural Research, grant R21 RR018440. We thank Victoria Soghomonian, Daniel Davis, and Jean Heremans for providing the AFM data.
- Arnison MR, Larkin KG, Sheppard CJR, Smith NI, Cogswell CJ: Linear phase imaging using differential interference contrast microscopy. J Miroscopy 2004, 214: 7–12. 10.1111/j.0022-2720.2004.01293.xMathSciNetView ArticleGoogle Scholar
- Cuche E, Bevilacqua F, Depeursinge C: Digital holography for quantitative phase-contrast imaging. Opt Lett 1999, 24: 291–293.View ArticleGoogle Scholar
- Poon T-C, Korpel A: Optical transfer function of an acoustio-optic heterodyning image processor. Op Lett 1979, 4: 317–319.View ArticleGoogle Scholar
- Poon T-C: Scanning holography and two-dimensional image processing by acousto-optic two-pupil synthesis. J Opt Soc Am 1985, A2: 521–527.View ArticleGoogle Scholar
- Indebetouw G, Klysubun P, Kim T, Poon T-C: Imaging properties of scanning holographic microscopy. J Opt Soc Am 2000, A17: 380–390.View ArticleGoogle Scholar
- Indebetouw G, El Maghnouji A, Foster R: Scanning holographic microscopy with transverse resolution exceeding the Rayleigh limit, and extended depth of focus. J Opt Soc Am 2005, A22: 892–898. 10.1364/JOSAA.22.000892View ArticleGoogle Scholar
- Indebetouw G, Zhong W: Scanning holographic microscopy of three-dimensional fluorescent specimens. J Opt Soc Am 2006, A23: 1699–1707. 10.1364/JOSAA.23.001699View ArticleGoogle Scholar
- Lohmann AW, Rhode WT: Two-pupil synthesis of optical transfer functions. Appl Opt 1978, 17: 1141–1151.View ArticleGoogle Scholar
- Rosen J, Indebetouw G, Brooker G: Homodyne holographic microscopy. Opt Express 2006, 14: 4280–4285. 10.1364/OE.14.004280View ArticleGoogle Scholar
- Goodman JW: Introduction to Fourier Optics. 2nd edition. New York: McGraw Hill; 1996.Google Scholar
- McCutchen CW: Generalized aperture and the three-dimensional diffraction image. J Op Soc Am 1964, 54: 240–244.View ArticleGoogle Scholar
- Papouli A: Systems and Transforms with Applications to Optics. McGraw Hill; 1968:131.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.