3 April 2014 Review of digital holographic microscopy for three-dimensional profiling and tracking
Author Affiliations +
Digital holographic microscopy (DHM) is a potent tool to perform three-dimensional imaging and tracking. We present a review of the state-of-the-art of DHM for three-dimensional profiling and tracking with emphasis on DHM techniques, reconstruction criteria for three-dimensional profiling and tracking, and their applications in various branches of science, including biomedical microscopy, particle imaging velocimetry, micrometrology, and holographic tomography, to name but a few. First, several representative DHM configurations are summarized and brief descriptions of DHM processes are given. Then we describe and compare the reconstruction criteria to obtain three-dimensional profiles and four-dimensional trajectories of objects. Details of the simulated and experimental evidences of DHM techniques and related reconstruction algorithms on particles, biological cells, fibers, etc., with different shapes, sizes, and conditions are also provided. The review concludes with a summary of techniques and applications of three-dimensional imaging and four-dimensional tracking by DHM.



Nowadays, three-dimensional (3-D) profiling and tracking microscale objects have been receiving much attention due to their wide applications. For example, microsphere or microbubble counting and locating in a 3-D space occur in the fields of microfluidics, suspension rheology, crystallization, etc. Many microorganisms swim in 3-D and helical paths. Thus, measuring their 3-D trajectories in time is essential to obtain detailed information on biophysical processes, such as motile behavior and dynamic performance.1213.14 Quantitative analyses of cancer cell locomotion and shape change in a 3-D environment reveal the biological characteristics for clinical need.1516. 3-D imaging of randomly oriented microfibers and their interactions with surrounding free-swimming cells opens up new perspectives in cases where quantitative measurement of characteristics (size, length, orientation, speed, displacement, etc.) is of great interest.21

Various approaches have been demonstrated for 3-D profiling and tracking of micro-sized objects, and optical techniques and numerical localization algorithms have been widely chosen as remarkable tools since they have the advantages of being full-field, label-free, noncontact, noninvasive, etc.2223. Digital holographic microscopy (DHM) is an emerging technology of a new paradigm in general imaging and biomedical applications.44,45 Various techniques of DHM, including off-axis DHM, in-line DHM, quantitative phase microscopy by digital holography (DH-QPM), etc., have been proven to be potent tools to profile and track micro-sized objects in 3-D volumes. In conventional microscopy techniques, only two-dimensional focused images on a fixed plane can be recorded, while information not in the focal plane would be permanently missed. However, DHM provides a tool to refocus an object and is able to record a hologram containing all the real-time, 3-D structures of an object in the absence of mechanical focus adjustment. And the information is available in digital form for postprocessing. DHM has been shown to be the key to 3-D particle image velocimetry (PIV)4647.48.49.50 since the conventional PIV techniques have the inherent limitations of thin depth of field, instrument complexity, and impossibility of real-time imaging.5152.53 However, DHM is able to overcome these and has the simplicity of the microscope and sample preparation, maximum information and resolution, reduction in time and data amount, etc. Aberrations and background distortions of the optical field can be minimized by available DHM techniques. The DHM for 3-D profiling and tracking generally is a two-step process: first a hologram is recorded digitally, and then the hologram is numerically reconstructed in different depths to yield an image of the object by various numerical diffraction methods, such as angular spectrum method,44,45 Kirchhoff-Helmholtz transform,54 Fresnel transform,55,56 Huygens convolution,44,45 etc. This reconstruction results in complex field of object, and one can extract its amplitude or phase profile to represent the object. Investigations of the third dimension based on the reconstructed data are then performed to determine focal planes and complete the 3-D pattern recognition. For a moving object, the reconstruction of successive holographic data yields a complete four-dimensional space-time record of the dynamic process. DHM is demonstrated to have the capacity of monitoring the 3-D distribution and motion pattern of particles, living cells, and fibers with different shapes (spherical, needle shaped, and randomly oriented), size (few to hundred micrometer), and conditions (static, suspended, and flow-through) in real time.

Description of several representative DHM techniques is introduced in Sec. 2, and the holographic information processing methods for 3-D profiling and tracking are given in Sec. 3. Wide and active field of applications on 3-D profiling and tracking of different samples, including microparticles, bubbles, biological cells, microfibers, etc., by DHM techniques are introduced in Sec. 4, and the conclusion is provided in Sec. 5. This review has an emphasis on applications of DHM in the field of 3-D profiling and tracking and, therefore, omits some major areas, such as digital holography principle and development, theoretical studies of digital holography, special digital holography techniques, etc.


DHM Configurations

A basic DHM setup consists of an illumination source, an interferometer, a digitizing camera, and a computer with necessary programs. Most often a laser is used for illumination with the necessary coherence to produce interference. For multiwavelength techniques, two or more different lasers can be coupled into the interferometer, or a tunable laser can be employed. There are also low-coherence techniques for the purpose of reducing speckle and spurious interference noise or generating contour or tomographic images. Even an LED typically has 30 μm or so coherence length, which can be sufficient for holographic microscopy.


In-Line DHM

In this review, we refer to Gabor holography as in-line holography.57 In-line DHM is a type of microscopy without objective lenses, and as illustrated in Fig. 1, a single light beam directed onto a pinhole of a diameter of the order of a wavelength illuminates the object, typically several thousand wavelengths away from the pinhole, and the object beam is the part of the incident light that is scattered by the object and the unscattered remainder is taken as the reference beam. If the light source is coherent, the pinhole is not necessary. The object field is automatically in alignment with the reference beam, and the interference of the object and reference beams results in the holographic diffraction pattern, which is recorded by CCD camera and then transferred to a computer for numerical reconstruction.

Fig. 1

In-line digital holographic microscopy (DHM) setup. MO, microscope objectives; P, pinhole; S, sample object.


By removing background effects first and then reconstructing the object field at different planes, a 3-D image can be built up from a single two-dimensional image, and the entire hologram pixel count is utilized, which also leads to shorter minimum distances for reconstruction and higher resolution of the resultant image. With the simplicity of the apparatus and large depth of field, the in-line DHM is particularly useful for particle image analysis. In conjunction with reconstruction algorithm, the 3-D position of a particle can be determined and then applied to particle velocimetry, such as tracking small particles or swimming cells in a liquid flow.


Off-Axis DHM

In an off-axis DHM setup, a Mach-Zehnder interferometer is chosen since it offers flexibility in alignment, especially when microscopic imaging optics is used. One microscope objective lens (MO) is needed for object magnification and another one is used in the reference arm to match the curvatures of the object and reference wave fronts. We illustrated the off-axis holographic microscopy setup in Fig. 2. The object arm contains a sample stage and an MO that projects a magnified image of the object onto a CCD camera. The reference arm similarly contains another MO, and the reference and object waves are offset by an angle to avoid the overlap of the reference and the twin images, so that the holographic interference pattern contains fringes due to interference between the diffracted object field and the off-axis reference field. The captured holographic image is numerically converted into Fourier domain to obtain the angular spectrum,44,45 and a spatial filter is then applied to retain the real image peak alone. The filtered angular spectrum is propagated to appropriate distance and, by an inverse Fourier transform, reconstructed as an array of complex numbers containing the amplitude and phase images of the sample. Both amplitude and phase profiles can be further analyzed to quantitatively determine the positions of objects in a 3-D space.

Fig. 2

Off-axis DHM setup. M's, mirrors; BS’s, beam splitters; MO’s, microscope objectives; S, sample object.



Quantitative Phase Microscopy by Digital Holography

Off-axis DHM is a very effective process for achieving high-precision QPM since it allows measurement of optical thickness with nanometer-scale accuracy by single-shot, wide-field acquisition, and it yields phase profiles without some of the complications of other phase imaging methods. The phase image is immediately and directly available on calculating the two-dimensional complex array of the holographic image, and the phase profile conveys quantitative information about the physical thickness and index of refraction of cells. We have recently utilized DH-QPM to study the wrinkling of a silicone rubber film by motile fibroblasts.58 The wrinkle formation has been visualized, quantitative measures of surface deformation have been extracted, and cellular traction force has been estimated in a direct and straightforward manner. A nonwrinkling substrate, collagen-coated polyacrylamide (PAA) was also employed to make direct measurement of elastic deformations, albeit at discrete locations.59 The Young’s modulus of PAA can be adjusted by controlling the concentrations of the monomer and cross-linker. DH-QPM, with its capacity of yielding quantitative measures of deformation directly, has been employed to measure the Young’s modulus of PAA,60 which provided a very effective process for achieving high-precision quantitative phase microscopy. We have utilized DH-QPM to study motile fibroblasts deforming PAA gel,61 where the cell-substrate adhesion has been visualized and quantitative measures of surface deformation have been extracted. The substrate stiffness and quantitative measures of substrate deformation have been combined to produce estimates of the traction forces and characterize how these forces vary depending on the substrate rigidity.


Twin-Beams DHM

The twin-beams DHM setup is illustrated in Fig. 3.62 Two beams, coming from the same laser source and slightly off-axis, enter into the 20× microscope objective and result in an image on a CCD plane. Theoretically, the two beams have the same focal plane; however, the twin beams can experience aberrations resulting in some focus-shift. As shown in Fig. 3, each particle of the sample appears as two projections on the CCD array. The tracking is performed by evaluating the two out-of-focus projections of the particles due to the twin-beams onto the CCD plane, and the separation between the two projections is a function of the longitudinal position of the particle. QPM images of the micro-objects can be directly obtained from the holograms recorded by the same configuration.

Fig. 3

Twin-beams DHM setup (reprinted from Ref. 62 by permission of OSA).



Numerical Focus Methods for Three-Dimensional Profiling and Tracking by DHM


Quantification of Image Sharpness and Peak Searching

To localize an object in a 3-D volume, an automated data analysis algorithm to accurately identify depth position is required. For example, pure amplitude objects located out of focal plane appear as gray patches without sharp structures, while pure phase objects show minimum visibility when they are in focus.6364.65 Thus, focusing metrics based on the quantification of image sharpness are applied by searching on the reconstructed images for the planes that contain the objects with the peak and sharpest details to establish the all-in-focus profile and its corresponding depth map.66 At each pixel, the intensity variation in longitudinal direction is investigated to find out the peak intensity. The peak intensity value becomes a pixel of an all-in-focus profile and its location in longitudinal direction is recorded as a corresponding depth map value. The combination of resultant all-in-focus profile and depth map enables one to produce the 3-D visualization of objects, and object properties can then be extracted.

We recently utilized off-axis DHM and numerical autofocusing algorithm based on peak searching to image suspended polymer microspheres (9.6 μm in diameter) in suspension.21 Figure 4(a) presents the hologram of microspheres in suspension captured by off-axis DHM. Angular spectrum method44,45 was applied and the reconstructed amplitude image was shown in Fig. 4(b), where several in-focus and out-of-focus microspheres were seen at this reconstruction plane. The field of view was 90×90μm2 with 464×464pixels. The microsphere focused the incoming light and formed a bright spot along the optical path. Therefore, the pattern of the microsphere appeared to have a maximum intensity at the center near the in-focus plane and the surrounding was dark. Then, the autofocusing algorithm based on peak searching was applied to identify the in-focus position of each particle. The hologram was numerically reconstructed in longitudinal direction Z from 100 to 100μm, in 2-μm steps. The autofocusing algorithm was then applied by searching on the reconstructed images for the planes that contained the objects with the peak and sharpest details to establish the all-in-focus profile. The peak intensity at each pixel along the longitudinal direction was chosen and this value became a pixel of an in-focus profile; meanwhile, its location in longitudinal direction was also recorded as a corresponding depth map value. The combination of all-in-focus intensity profiles [Fig. 4(c)] and the corresponding depth maps [Fig. 4(d)] enabled to produce the 3-D visualization of microspheres and a threshold on the intensity allowed for distinguishing objects from other elements, which is illustrated in Fig. 4(e). In fact, we took every center position of brightest points in Fig. 4(c) as X-Y locations of particles, and Fig. 4(d) was used to determine the Z-location of those particles identified in Fig. 4(c). Figure 4(d) is mostly noisy and Z information of only the locations where the particles were present was actually used. Figures 4(f) to 4(h) show the XY, XZ, and YZ views of the 3-D profile, respectively. This straightforward autofocusing algorithm is able to determine focal planes for all the objects in the reconstructed volume automatically. The combination of depth map Z and axial X-Y position of objects allows for quantitative 3-D profiling. This method has been widely used by researchers, and more applications will be discussed in Sec. 4.

Fig. 4

Three-dimensional (3-D) profiling of suspended microspheres. Field of view of (a) to (d) is 90×90μm2 with 464×464pixels. (a) Hologram. (b) Amplitude image. (c) All-in-focus intensity profile. (d) Depth position profile. (e) 3-D profile (90×90×120μm3) and gray-scale representation of intensity. (f) XY view of three-dimensional profile (90×90μm2). (g) XZ view of three-dimensional profile (90×120μm2). (h) YZ view of three-dimensional profile (90×120μm2).



3-D Deconvolution Methods

Overviews of the deconvolution methods regarding the out-of-focus problem in optical microscopy can be found in Refs. 6768.69. to 70. In Ref. 71, 3-D deconvolution methods, including instant and iterative deconvolution, have been applied to restore particle distribution from holograms and the corresponding reconstructions. Instant 3-D deconvolution could be applied on separated particles spread in a certain volume and allowed for the identification of the depth position of object and removal of out-of-focus information. The way to perform the instant 3-D deconvolution was to apply it to reconstructed intensity distributions from holograms using Eq. (1):


where O is a 3-D object, |Uo|2 is the reconstructed intensity of the object wave, |UP|2 is the reconstructed intensity of a point scatter [point spread function (PSF)], β is a small constant to avoid denominator to be zero, FT is 3-D Fourier transform. A simulated particle distribution in the shape of letters α, β, γ was placed at three different planes and each of them was constituted of point scatters, shown in Fig. 5(a). The reconstructions were performed at continuous planes in the volume and the PSF was chosen as the reconstruction of the hologram of a point scatter located at a fixed plane. The results of 3-D deconvolution from Eq. (1) showed a sharp 3-D visualization of individual scatters constituting the three letters, which were seen to be distinguished from each other, as shown in Fig. 5(b).

Fig. 5

(a) and (b) Results of the instant 3-D deconvolution of the intensities reconstructed from the simulated hologram of a particle distribution. (c) and (d) Results of the iterative 3-D deconvolution of the complex fields reconstructed from the simulated hologram of continuous objects (reprinted from Ref. 71 by permission of OSA).


Iterative deconvolution could be applied on objects of continuous and extended shapes and distributions. The way to perform the iterative deconvolution was to apply it to reconstructed complex fields from holograms. The iterative loop included


where On and Uon are the complex distributions of the object and reconstructed wave front, respectively. A simulated hologram of three continuous, not-point-like letters α, β, γ shown in Fig. 5(c) was placed at three different planes. With the same parameters of reconstruction as the instant 3-D deconvolution case, the iterative deconvolution was performed by the iterative loop described above in Eq. (2). The result showed the three continuous letters, which were seen to be clearly distinguished in a 3-D volume [Fig. 5(d)]. These two methods were also demonstrated by experimental in-line holograms of polystyrene microspheres randomly locating on both sides of a thin glass. A cut-out from the original hologram of microspheres, which included only a single sphere, was taken as the PSF and the 3-D localization of microspheres was retrieved. Moreover, a simulated PSF was also tested and a comparable good result was obtained.


Rayleigh-Sommerfeld Backpropagation Method

Rayleigh-Sommerfeld backpropagation method is a general and fast volume reconstruction analysis tool for DHM.55,72,73 A well-known Gouy phase anomaly55,74 can be used to discriminate between objects lying on either side of the holographic image plane. A recorded in-line hologram was reconstructed using Rayleigh-Sommerfeld diffraction integral to refocus the scattered object at a particular position and observe the phase anomaly as a contrast inversion about the object’s position by Wilson et.al.72 If the complex field of scattering objects was obtained in the focal plane, it could be reconstructed at plane z above the focal plane as a convolution of amplitude in the focal plane with Rayleigh-Sommerfeld propagator, as in Eq. (3):55


where r=(x2+y2+z2). The intensity gradient was then extracted based on the convolution with the Sobel filter,72 and if a particle is located in positive z-direction, the field was converging when it arrived at the image sensor and had not yet passed through the geometrical focus. This allowed for determining the location of the center of the object and separating objects on different sides of the hologram plane. Particles always had twin images in the reconstructed space no matter they are located above or below the focal plane. If only the upstream of the focal plane was reconstructed, the intensity gradient was a maximum for particles with z<0 and a minimum for those with z>0. A micro-sized diameter sphere was tested by this technique and a typical volumetric reconstructed rays converging to the bright point above the focal plane was shown in Fig. 6. Another analysis of an experimental sample consisted of closely spaced particles on a glass slide showing the applicability of this method to distinguish weakly scattering particles, such as colloids or planktons locating at close lateral but separate axial positions.

Fig. 6

Example data from a single particle. Scale bars represent 2 μm in all cases. (a) Vertical slice through the center of an image stack created by physically translating the sample. (b) Image of a particle located at z9μm (downstream of the focal plane in the illumination path). (c) Optical field reconstructed from the previous panel. The hologram plane (z=0) would be located below the bottom of the image. (d) Intensitygradient<0. The dark central spot is azimuthally symmetric about the z axis and gives the particle location in all three dimensions. (e) and (f) The companion images to (b) and (c), for a particle located at z9μm (upstream in the illumination path). (g) Intensitygradient>0. The particle location is specified by a maximum of intensity gradient for those scatterers with z<0 (reprinted from Ref. 72 by permission of OSA).



Compressive Holographic Method

Compressive sensing is a technique to efficiently acquire the information of signal in the underdetermined linear system supported by the sparseness of the signal.75 Hence, the higher dimensional signal can be recovered by the lower-dimensional measurements. Recently, some research groups have been studying about the compressive holography, which applied the concept of compressive sensing to holography.7576.77.78 In holographic recording, 3-D data of the object were recorded as a two-dimensional complex hologram. Supported by the condition that the distribution of the object in 3-D space was sufficiently sparse, the compressive sensing technique could retrieve the 3-D distribution from two-dimensional complex hologram. Brady et al.75,76 had successfully shown the possibility of the compressive holography with the hologram recorded by in-line DHM. In-line DHM, because of its simplicity and efficiency, provided a potent tool for compressive sampling by recording a 3-D object onto a two-dimensional focal plan array, which was related to Fourier transform of the object field. Compressive holography collected fewer amounts of measurements than voxels in the reconstructions. The 3-D datacube was encoded and compressed into two-dimensional holographic measurement by holographic sampling process, and the encoding was then inverted using the compressive sampling theory.

Two seed parachutes were illuminated and placed away from the detector array at different distance separately [Fig. 7]. After reconstruction, the stem and petals representing the high-frequency features in the image were clearly shown and the distance between the two parachutes were also estimated. The final result presented a 3-D datacube of voxels reconstructed from a single two-dimensional hologram, and it demonstrated the applicability of compressive holography to generate multidimensional images from lower-dimensional data. A detailed description of compressive holography theory and experiment results can be found in Refs. 75 and 76.

Fig. 7

(a) Raw Gabor hologram for seed parachutes of taraxacum. (b) and (c) Photographs of the individual objects. (d) Transverse slices at various ranges of the backpropagated (numerical refocusing) field (reprinted from Ref. 75 by permission of OSA).



Twin-Beams DHM Method

As previously depicted in Fig. 3, the twin-beams DHM had two plane beams passing through one objective microscope. Theoretically, this process could be sketched as two cones superposed and transmitted through the particle sample P (Fig. 8, left). Beams B1 and B2 produced projections P1 and P2 of particle P on the CCD plane, respectively. In both simulation and experiment, the 3-D coordinates of B1 and B2 points and the radius of each beam circle on the CCD plane were first determined based on a simple calibration method.62 Then, the corresponding coordinates of P1 and P2 could be evaluated accordingly from the simple geometric schematic and the position of real object P was determined. A real situation of various motile cells floating into a microfluidic chamber was studied using this method and the 3-D paths of cells were estimated [Fig. 8(b)].62 The phase-contrast maps of the cells were also obtained by numerically reconstructing the recorded holograms.79,80

Fig. 8

Twin-beams DHM results. Sketch of interference between two beams (a) and estimated path for random motion of the three cells (b) (reprinted from Ref. 62 by permission of OSA).



Applications of Three-Dimensional Profiling and Tracking by DHM


Three-Dimensional Profiling of Suspended Particles

Micro-sized particles, such as latex spheres and ferromagnetic beads, suspended in gelatin have been imaged in a 3-D volume by in-line DHM.8182.83 A hologram of objects and another background image without the object present were both recorded. The difference of these two images eliminated the artifacts of the illumination source. Reconstruction was achieved by means of Kirchhoff-Helmholtz transform84 to generate the complex field of objects, and the combination of a stack of two-dimensional reconstructions at various planes resulted in the 3-D profile of the objects. The amplitude obtained from complex field was used to represent objects. By inspecting intensity profiles through a particle, the depth position was obtained where a sharpest intensity appeared. This allowed estimating a particle’s 3-D position coordinates within submicron. When the derivative of a second-order polynomial fitted to the intensity profiles was taken, the X, Y, and Z positions of particles could be determined within 50 nm.

In Ref. 85, the distribution of dense microparticles (3.2 μm) suspended in a liquid volume with the depth of 1 mm was measured, as shown in Fig. 9. Figure 9(a) is a section of the recorded hologram and the reconstructed images at different depths are shown in Figs. 9(b) to 9(d). It was clearly observed that the in-focus particles appeared as dark circular spots. Figure 9(e) shows the combination of all the reconstructed images. The peak-searching method was then applied and each pixel was assigned the lowest intensity obtained over the entire depth. A thresholding segmentation method based on signal-to-noise ratio [II¯]/σI was applied on the reconstructed planes,86 where I¯ is the mean intensity of the particles in the volume and σI is the standard deviation of intensity over this volume. A list of line segments after scanning through the reconstructed images were combined to two-dimensional planar blobs using a join operator, which were then united into 3-D particle traces by repeating the peak-searching procedure at different depths. The position of a particle was then determined by using the centroid of the 3-D blob and the distribution of all the particles detected within the 3-D volume was presented in Fig. 9(f).

Fig. 9

Results of 3-D profiling of particles. (a) Part of a recorded hologram using a 10× objective, containing 3.189-μm-diameter particles in a 1-mm-deep solution. (b) to (d) Reconstruction of planes located 120, 580, and 800 μm from the hologram plane. In-focus particles appear as dark spots on the bright background. (e) A combined and/or compressed image containing all the particles covered by the hologram section shown in (a). (f) Location of all the particles detected within the entire 1.5×1.5×1mm3 volume, totaling 5769 particles (reprinted from Ref. 85 by permission of OSA).



Three-Dimensional Study of Microfluidics

Kreuzer’s group has studied the microfluidic phenomena by in-line DHM with the capacity of large depth of field and numerical reconstruction.87 In a first experiment, a sphere of 150 μm diameter was attached to the top of a millimeter-wide tank, which was filled with water and seeded with small micro-sized latex beads as tracer. The water was sucked and pumped to flow by a blotting paper and the flow past a spherically shaped obstruction was recorded as a series of holograms by in-line DHM. The difference between consecutive hologram pairs was calculated pixel by pixel (h1h2,h3h4,,hn1hn) to eliminate background structure and retain only the information of moving beads, where hn is the n’th hologram of recorded frames. The resulting difference holograms were summed (h2h1+h4h3+h13h14) into a single hologram, which was then reconstructed with Kirchhoff-Helmholtz transform to obtain the focal planes of the stream lines of beads (Fig. 10). The flow information containing sequential positions of beads at successive recording times and the velocity field defined by vectors connecting two successive positions was measured and analyzed and the velocity field was proven to be in good agreement with Navier-Stokes’ solution. In this way, as indicators of flow patterns, the trajectories and velocities of latex microspheres were measured in a 3-D volume and time with subsecond and micron resolution.

Fig. 10

(a) 3-D flow around a fixed sphere. Ninety holograms were taken at intervals of 0.17 s to generate a difference hologram. (b) One of 60 reconstructions made to render the field velocity shown in (c). In (d) the solid blue lines represent the solution to the Navier–Stokes equation for our experiment and the arrows correspond to the measured velocity field (reprinted from Ref. 87 by permission of Optik).


The size and 3-D positions of particles in flow-through system, such as fast-moving bubbles in air-water mixture flows, were measured by Tian et al.88 Bubbles were generated by a motor driving a small propeller in a water tank and the density was not so high as to corrupt the illumination beam. Holograms were first recorded by in-line DHM system and the Fresnel approach was used to propagate the recorded object field to various planes in the vicinity of the object.55,56 At each propagation distance, the Fourier transform of the hologram was multiplied by the Fresnel approximated transfer function and an inverse Fourier transform was then applied to obtain the complex field of objects containing the intensity profile. A focused bubble appeared to have a minimum intensity and this was used as a focus criterion to locate bubbles. A two-dimensional projection image was generated by scanning the reconstructed images and recording the minimum intensity value for each pixel. The corresponding depth map where the minimum intensity occurred was also recorded for analysis. By applying thresholding on the minimum intensity projection image, a binary image showing only the edges of bubbles was created and then combined with the depth map to determine the depth of objects. The issue of overlapping edges in the projection image was resolved by using a Gaussian mixture model.89 Bubble sizes were determined from the areas enclosed by each edge. In this way, the 3-D visualization of the positions and sizes of bubbles was presented, as shown in Fig. 11.

Fig. 11

Data processing results from a 1024×1024pixel hologram. (a) 3-D visualization of data processing results from a single hologram (diameters not to scale). (b) Bubble size distribution (reprinted from Ref. 88 by permission of OSA).



Four-Dimensional Motility Tracking of Biological Cells

Motility is a major characteristic of certain types of cells, and it is essential to understand the motility of cells by measuring their 3-D motion. DHM offers a rapid and efficient method to monitor the 3-D motile behavior and dynamic process of various types of cells. The effectiveness of this method is quantified in terms of the maximum number and speed of the moving cells that can be tracked, as well as the direction of the motion (lateral and axial). The results from this work can be further used to implement an optical trap that can automatically track and capture cells with specified characteristics, such as speed, size, or shape.


Red blood cells

DHM has been applied in the field of hemodynamic research.90 A human red blood cells (RBC) sample was supplied to a circular microtube and the RBC motion in this microtube flow was recorded sequentially by in-line DHM. The angular spectrum method and quantification of image sharpness were applied to reconstruct and locate the focal planes of cells. The measurement accuracies in the lateral and axial directions were ±0.3 and ±1μm, respectively. In this way, the full four-dimensional trajectories of RBCs, including space and time information, and the three-dimensional velocity field were measured, and the linear movement of cells was seen in the circular microtube. The feasibilities of in-line DHM, volumetric particle tracking, and numerical reconstruction were demonstrated with RBCs.

3-D tracking of the sedimentation of RBCs has also been implemented by DH-QPM based on the phase contrast images.91 The captured holograms were reconstructed by nondiffractive spatial-phase-shifting-based reconstruction92,93 to access to the complex field of objects, which was then numerically propagated to different planes by a convolution-based algorithm.94 The phase profile was extracted and phase unwrapping algorithm was applied for future focusing evaluation. The autofocus criterion for a sharply focused phase object showed minimum visibility.64,65 The time dependence of axial displacement and depth position where the focused phase happened were combined to illustrate the 3-D process of the sedimentation of RBC due to gravitation, as illustrated in Fig. 12.

Fig. 12

Time-dependent 3-D tracking of a sedimenting human RBC by DHM. (a) to (d) Unfocused phase distributions of the RBC during the sedimentation process at t=0, t=22s, t=44s, and t=66s. (e) to (h) Quantitative phase distributions obtained of the same digital hologram as in (a) to (d) by application of digital holographic autofocus; the dotted box represents the region of interest for the x-y tracking in (j). (i) Time dependence of the axial position Δg of the sample obtained by digital holographic autofocusing. (j) x-y trajectory of the RBC obtained by determination of the coordinates with maximum phase contrast. (k) 3-D trajectories obtained by combination of the data in (i) and (j) (reprinted from Ref. 91 by permission of SPIE).



Free-swimming cells

Chilomonas are fast-moving cells, which sense and respond to their surroundings by swimming toward or away from stimuli. We have utilized off-axis DHM to track the motility of chilomonas in a 3-D volume and time.21 A time-lapse hologram movie of the movement of chilomonas (average diameter 10 μm) was first recorded at 30 fps. The field of view was 90×90μm2 with 464×464pixels. An excerpt of 14 frames was taken and the difference between consecutive hologram pairs was calculated pixel by pixel (h1h2,h3h4,hn1hn) to eliminate background structure and retain only the object information, where hn is the n’th hologram of selected frames. The resulting seven difference holograms (from a total of 14) were then summed (h2h1+h4h3+h13h14) into a single hologram [Fig. 13(a)], which contained all the information on moving cells. Reconstruction by the angular spectrum method was applied on the resultant hologram, and Fig. 13(b) was the amplitude image reconstructed at the specific plane Z=0. The 3-D trajectory of the moving chilomonas was built up by first reconstructing the two-dimensional resultant hologram from Z=80 to 120μm, in 2-μm steps. The autofocusing method based on peak searching was then employed, keeping the maximum intensity value at each pixel in the reconstructed image along Z direction to obtain the all-in-focus intensity projection [Fig. 13(c)] and combining with the depth position where in-focus intensity occurred [Fig. 13(d)] to determine the focal planes for the entire trajectory in the reconstructed volume. A threshold of intensity was applied to eliminate the unnecessary background, and the center position and the mean intensity of every point cloud representing every trajectory of chilomonas were determined. The combination of depth position Z and XY displacement of cells allowed for quantitative 3-D motility tracking. The cell was estimated to move a total path length of 93 μm at the velocity of 198μm/s. Figure 13(e) demonstrates the applicability of DHM for automatic four-dimensional tracking of living cells with temporal and spatial resolution at the subsecond and micro level. Figures 13(f) to 13(h) show the XY, XZ, and YZ views of the 3-D trajectory, respectively.21

Fig. 13

Four-dimensional tracking of a moving chilomonas cell. The field of view of (a) to (d) is 90×90μm2 with 464×464pixels. (a) Difference hologram. (b) Amplitude image. (c) All-in-focus intensity profile. (d) Depth position profile. (e) 3-D profile of the trajectory of the cell (90×90×90μm2) and gray-scale representation of intensity. (f) XY view of 3-D trajectory (90×90μm2). (g) XZ view of 3-D trajectory (90×90μm2). (h) YZ view of 3-D trajectory (90×90μm2).


In Ref. 95, a real-time hologram movie of free-swimming algae (5 to 10 μm in diameter) in seawater was recorded by in-line DHM. Subtracting the consecutive frames first and adding these difference frames then provided a final image containing the time evolution of algae trajectories free from background interference. A similar peak-searching reconstruction algorithm as described previously constructed the 3-D trajectories of algae in seawater approaching a concentrated salt solution. Free-swimming algae were seen to sense and respond to their surroundings by swimming toward or away from stimuli at an average velocity of 150μm/s.

Specimen of small-sized bacteria (<1.5μm in diameter) inside a diatom (200 to 300 μm in diameter) was tracked in three dimensions since the diatom was almost stationary in seconds and the in-line DHM and reconstruction algorithm, including peak searching, edge detection, etc., were applicable in this weaker structure.83 This demonstrated the possibility of tracking smaller and weaker scatters with high resolution by in-line DHM.


Special cells with helical trajectories

Various self-propulsive microorganisms of different shapes with various width-to-length ratios exhibit helical swimming paths. Unlike spherical particles, these structures would cause nonuniform scattering, which makes it difficult to determine the depth map due to the large depth-of-focus of the specimen. In Ref. 96, cellular behaviors and interactions among microorganisms, such as the prey-induced changes in swimming behavior of predatory dinoflagellates, were statistically studied. The captured holograms by in-line DHM were numerically reconstructed by Fresnel transformation and the intensity field of the reconstructed images can be obtained. The same algorithm as described in Refs. 85 and 86 was utilized to determine the 3-D positions of cells. Besides this, information on each cell, such as cross-section area, dimensions, aspect ratio, total reconstructed volume, and in-focus snapshots, are stored for cell tracking. Each helical trajectory was taken as a vector determined by each criterion step by step, i.e., smoothness of the trajectory segment, 3-D velocity, acceleration, similarity of cell size, in-focus cross-section, 3-D segmented volume, shape based on correlation, aspect ratio, etc. This statistical analysis provided simultaneous tracking and characterization of cellular behavior and interactions among microorganisms in dense suspensions. Figure 14 presents the trajectories of P. piscicida before introduction of prey in a 3-D volume. In Ref. 97, another type of ellipsoidal shape Cochlodinium polykrikoides cell had a size of 40μm length and 25 μm width, and it could form chains consisting of several cells. A sample of these mixed C. polykrikoides cells, also exhibiting helical swimming paths in filtered seawater, were also statistically studied by in-line DHM and similar approaches.

Fig. 14

The 3-mm-deep sample volume showing trajectories of P. piscicida before introduction of prey. Upper inset shows close-up of part of the sample volume, which is marked by red dash-dot line in the bottom image [reprinted from Ref. 96, Copyright (2007) National Academy of Sciences, U.S.A].


Moreover, helical trajectories (>1500) of human sperm (3 to 4 μm) within a large volume (>1μL) was tracked and analyzed by a dual-view and dual-color lens-free DHM setup, which consisted of two partially coherent LEDs of different wavelengths placed at 45 deg with respect to each other (Fig. 15).98 The recorded in-line holograms only reconstructed at correct combination of depth, angle, and wavelength could produce clear images, and the sperm head images projected at two different wavelengths and viewing angles could be isolated even if they were recorded in the same frame. Combining with iterative phase recovery method,99 the peak-searching algorithm, and related numerical smoothing method, the 3-D trajectories of sperms and the statistics information of the swim path, pattern, speed, etc., were revealed (Figs. 15 and 16).

Fig. 15

Dual-view lensfree 3-D tracking of human sperms. (a) The schematic diagram of the imaging system. Two partially coherent light sources (red and blue LEDs at 625 and 470 nm, respectively) are butt-coupled to multimode fibers to simultaneously illuminate the sperms at two different angles (red at 0 deg and blue at 45 deg). A CMOS sensor chip re-cords the dual-view lensfree holograms that encode the position information of each sperm. The 3-D location of each sperm is determined by the centroids of its head images reconstructed in the vertical (red) and oblique (blue) channels. (b) The reconstructed 3-D sperm trajectories. At a frame rate of 92 fps, 1575 human sperms were tracked inside a volume of 7.9 μL. The time position of each track point is encoded by its color (see the color bar) (reprinted from Ref. 98 by permission of Proceedings of National Academy of Sciences, U.S.A).


Fig. 16

Four major categories of human sperm swimming patterns. (a) The typical pattern. (b) The helical pattern. (c) The hyperactivated pattern. (d) The hyperhelical pattern. The inset in each panel represents the front view of the straightened trajectory of the sperm. The arrows indicate the directions of the sperms’ forward movement. The time position of each track point is encoded by its color (see the color bar). The helices shown in (b) and (d) are both right-handed (reprinted from Ref. 98 by permission of Proceedings of National Academy of Sciences, U.S.A).



Crawling cells in matrix gels

Ferraro’s group has presented a 3-D tracking of in vitro osteosarcoma cells in a gelified collagen matrix characterized by a high morphological change.100 DH-QPM was utilized to capture time-lapse sequences of cell migration at 12 fph for 15 h to record any movement and morphological change of cells. Numerical reconstruction allowed for evaluating both intensity and phase profiles in different planes. The cells were first detected from phase profiles by applying a simple thresholding filter, and this type of pure phase objects showed minimum visibility when they were in focus; therefore, the estimation of focal planes was performed by numerical scanning on the stack of reconstructed amplitude images instead for the sharpest details. Then, the axial position of the center of mass of each cell was determined from phase profiles at the focal plane.101 Taking into account the morphological changes of cells, the concept of minimum boundary filters100 was introduced to determine the positions of cells, and it was compared with other three methods, i.e., centroid, weighted centroid, and maximum phase value, as shown in Fig. 17.

Fig. 17

3-D tracking of in vitro cells in a gelified collagen matrix. Column (a) shows three different cells with the (X, Y) positions estimated using the four methods, and the in-focus distance is reported on the top of each (a). Column (b) shows the trajectories estimated of the four methods (reprinted from Ref. 100 by permission of OSA).


In Ref. 102, to reduce the coherent noise, a partially spatial coherent illumination was utilized as the illumination source of off-axis DHM setup and an accurate equalization between the object and reference beams must be reached. A semiautomated software package performed a four-dimensional tracking of unstained HT-1080 fibrosarcoma cells migration in matrix gels on the basis of phase profiles, and the 3-D locations of a central point in the cell as a function of time were recorded simultaneously. Quantitative measures on the displacement and velocity of cells were also obtained. By numerical autofocusing and evaluation of quantitative phase images, DH-QPM was demonstrated to be a label-free and noninvasive method for tracking living cells in a 3-D environment.


Characterization of Microfibers

Kempkes et al. have applied in-line DHM to record a hologram of a needle-shaped tilted opaque carbon fiber103 and showed the reconstruction at different planes using Fresnel-Kirchhoff diffraction integral.94 The intensity profiles were then obtained and converted to binary images by setting a threshold, which retained only the in-focus pixels. The combination of the stack of binary images produced a 3-D profile containing point clouds resulting from the focused pixels of a fiber at different reconstructed planes. A straight line was fit onto the point clouds, and the orientation and length of the fiber were estimated. A more complex sample of fibers randomly suspended in water was studied in the same way. Seven fiber segments were identified and the corresponding orientation and length were measured (Fig. 18). In-line DHM is shown to be an effective technique to characterize nonspherical particles in 3-D volume.

Fig. 18

3-D image of a volume of carbon fibers in suspension, viewed from two different directions and showing both the fitted lines (black lines), as well as the point clouds (colored dots) (reprinted from Ref. 103 by permission of OSA).


We have recently studied the nonwoven curved and randomly oriented microfibers in suspension.21 Fibers were spun onto glass coverslips from co-poly (L-gluatmicacid4, L-tyrosine1) (PLEY) dissolved in water and crosslinked as described previously.58 The diameter of fibers was 0.2 to 5.0 μm. Figure 19(a) presents a hologram of microfiber. The angular spectrum method was applied to obtain the reconstructed amplitude image in Fig. 19(b). The microfiber appeared curved and elongated, consistent with our previous work.104 The field of view was 90×90μm2 with 464×464pixels. The hologram was then reconstructed from Z=100 to 100μm, in 1-μm steps. Adding these 201 reconstructed images of amplitude yielded the axial projection shown in Fig. 19(c). The final image had the same number of pixels as each reconstructed image. The projected image was converted to binary by threshold and segmentation, as shown in Fig. 19(d). The binary image was then axially backprojected along the path of amplitude over the set of reconstructed images. We applied the same autofocusing algorithm on the resultant images to record the maximum intensity value at each pixel along the reconstruction direction and the axial XY and depth position Z of where the maximum intensity occurred at the same time. The obtained all-in-focus intensity profile [Fig. 19(e)] and the corresponding depth map [Fig. 19(f)] were then combined to generate the 3-D visualization of microfiber and a threshold on intensity was chosen to highlight the area of regional interest. However, there still existed noisy point clouds from speckle noise in the 3-D profile [Fig. 19(g)]. To address this, an average and polynomial curve fitting algorithm accounting for errors was adopted. The fitted line of the point clouds (standard deviations in the X and Z directions were 0.21 and 0.06 μm) in the 3-D coordinate was present in Fig. 19(h). Figures 19(i) to 19(k) show the XY, XZ, and YZ views of the 3-D profile, respectively. Z coordinate provides the depth information of real image of microfiber in the reconstruction volume. The length of the microfiber in 3-D volume was determined from X, Y, and Z coordinates of real image of microfiber, which was 126 μm.21

Fig. 19

3-D profiling of stationary microfibers. The field of view of (a) to (f) is 90×90μm2 with 464×464pixels. (a) Hologram. (b) Amplitude image. (c) Projection image. (d) Binary image, microfiber 1 and background 0. (e) All-in-focus intensity profile. (f) Depth position profile. (g) 3-D profile (90×90×35μm3) and gray-scale representation of intensity. (h) 3-D fit (90×90×35μm3). (i) XY view of 3-D fit (90×90μm2). (j) XZ view of 3-D fit (90×35μm2). (k) YZ view of 3-D fit (90×35μm2).



Cell–Environment Interaction

When cells swim in the local environment, they exhibit movement toward or away from the stimuli, as well as interaction with the stimuli. Knowledge of this interaction is essential for understanding the impact of cell’s behavior on the process of locating food, avoiding predators, finding mates, etc.


Three-dimensional flow field by swimming copepod

In-line DHM was utilized to measure the 3-D trajectory of a free-swimming copepod and the complex flow around it.105 A sample of Diaptomus minutus in well water was prepared and seeded with 20 μm polystyrene spheres. The in-line DHM apparatus with He-Ne laser source consisted of two inclined mirrors on the walls of the sample volume to provide two perpendicular views of the object in the same image. The series of holograms that captured the copepod with two views were reconstructed using Fresnel-Huygens principle106 with a paraxial approximation to obtain the intensity distribution. Based on the location of the traces resulting from the superposition of three successive exposures, the corresponding location of mirrored views was estimated. With each of the two views, the particle displacement was measured using cross-correlation of the intensity distribution of the three exposures,107,108 and the lateral and vertical displacements were then combined into a 3-D velocity distribution of particles, as shown in the left panel of Fig. 20. In the reconstructed images sequence movie, the copepod was seen to sink, swim upward, and then sink again. As it sinks, a recirculating flow pattern was generated due to the movement of appendages, as illustrated in the right panel of Fig. 20. Stokeslets model109,110 was employed to estimate the excess weight (7.2×109N) and excess density (6.7kg/m3) of the copepod, and also the propulsive force generated by appendages (1.8×108N).

Fig. 20

Left: measured instantaneous velocity near the copepod (a) in the ambient frame of reference and (b) in the copepod frame of reference. Right: particle streaks in the copepod reference frame obtained by combining 130 appropriately shifted reconstructed images. In all cases, the dorsal and lateral views are in focus (reprinted from Ref. 105 by permission of Journal of Experimental Biology).



Three-dimensional displacement of microfiber by swimming paramecium

A hologram movie of paramecium moving through a microfibers mesh was recorded by off-axis DHM at 17 fps.21 A cell was seen to swim through fibers and cause obvious displacement of fibers. Two hologram frames showing cell approaching to fibers [Fig. 21(a)] and swimming away after pulling fibers [Fig. 21(b)] were extracted. The field of view was 200×200μm2 with 768×768pixels. The difference hologram [Fig. 21(c)] clearly showed the displacement of the fiber due to the movement of cell and the trajectory of cell. The reconstructed amplitude images from angular spectrum method at Z=0 and Z=8μm are shown in Figs. 21(d) and 21(e), respectively. Different parts of the fiber appeared to be in focus at different reconstructed planes. Autofocusing algorithm at the reconstruction volume from Z=60 to 70μm was applied and the all-in-focus intensity and depth position information were obtained in Figs. 21(f) and 21(g).

Fig. 21

3-D displacement of microfiber by swimming paramecium. The field of view of (a) to (h) and (j) is 200×200μm2 with 768×768pixels. (a) Hologram 1. (b) Hologram 2. (c) Difference of holograms 1 and 2. (d) Amplitude reconstructed at Z=0. (e) Amplitude reconstructed at Z=8μm. (f) All-in-focus intensity profile. (g) Depth position profile. (h) All-in-focus intensity of the first track of the displacement of fiber. (i) 3-D fitted line of (h) (200×200×20μm3). (j) All-in-focus intensity of the second track of the displacement of fiber. (k) 3-D fitted line of (j) (200×200×20μm3). (l) 3-D fitted line of the displacement of fiber (200×200×20μm3). (m) XY view of (l) (200×200μm2). (n) XZ view of (l) (200×20μm2). (o) YZ view of (l) (200×20μm2). For better visualization, the scales of X-Y and Z coordinate are not in ratio.


To best visualize the displacement of the fiber in 3-D volume, we analyzed the two tracks of the displacement of the fiber independently. For brevity, most of the following descriptions referred to the first track. Numerical segmentation and thresholds were applied to Fig. 21(f), and the all-in-focus intensity image containing only one track of fiber was generated in Fig. 21(h). Combined with the obtained depth position profile [Fig. 21(g)] and fitting algorithm accounting for errors, the fitted line of the point clouds (standard deviations in X and Z directions were 1.27 and 0.27 μm) in 3-D coordinate is present in Fig. 21(i). Similar results for the second track of fiber were shown in Figs. 21(j) and 21(k), and the standard deviation of its corresponding point clouds in the X and Z directions were estimated to be 0.80 and 0.26 μm. The two tracks of the fiber were then plotted in one 3-D coordinate [Fig. 21(l)], and Figs. 21(m) to 21(o) show the XY, XZ, and YZ views of the 3-D plot, respectively. The length and the displacement of the fiber were estimated to be 210.1 and 11.1 μm. A movie of 3-D displacement of the fiber within 0.5 s period was also presented. The 3-D displacement of an individual microfiber due to interaction with paramecium cell has been measured as a function of time at subsecond and micrometer level.21

To the best of our knowledge, this is the first quantitative profiling and tracking by DHM of the curvature and displacement of individual microfiber by swimming cells in three dimensions. Displacement could vary with cell type, physiological state of the cells, microfiber preparation, etc. The prospects seem excellent for DHM for particle flow analysis and 3-D imaging of randomly oriented microfibers, in cases where quantitative measurement of characteristics (size, length, orientation, speed, displacement, etc.) is of great interest.



We present a review of a subset of DHM for 3-D profiling and tracking, focusing on techniques, reconstruction criteria, and applications in various fields. For brevity, a summary of the research progress is provided in Table 1. The table is divided into studied objects and their status, and also includes the DHM techniques used for imaging and the corresponding reconstruction methods and autofocusing algorithms. DHM is demonstrated to have the capacity of full-field, label-free, noncontact, noninvasive, real-time refocusing of an object without mechanical focus adjustment and proven to be a potent tool to perform 3-D imaging and tracking. The combination of DHM with robust autofocusing algorithms enables the remarkable study of 3-D profiling and tracking of micro-objects or nano-objects. It has merits of wide applications covering particles, living cells, and fibers with different shapes (spherical, needle shaped, and randomly oriented), size (few to hundred micrometer), and conditions (static, suspended, and flow-through). Due to the uncertainties connected with the imaging target in 3-D profiling and tracking, it is important to tailor and find strategies that can be adapted to the specific problem.

Table 1

Summary of research progress in the field of digital holographic microscopy (DHM) for three-dimensional profiling and tracking.

ObjectStatus of objectDHM configurationMethod of reconstruction & autofocusing algorithms
Microspheres/beads SuspendesusSuspended, sparseOff-axisAngular spectrum sharpness and peak searching21
Suspended, sparseIn-lineKirchhoff-Helmholtz transform sharpness and peak searching81–83
Suspended, sparseIn-lineThree-dimensional deconvolution71
Suspended, sparseIn-lineRayleigh-Sommerfeld backpropagation72
Suspended, denseIn-lineFresnel sharpness and peak searching85
Flow-throughIn-lineKirchhoff-Helmholtz transform sharpness and peak searching87
Seed parachutesSeparated positionIn-lineCompressive holographic method75
BubblesFlow-throughIn-lineFresnel sharpness and peak searching88
In vitro cellsFloatingTwin-beamsTwin-beams DHM method62
Chilomonas (10 μm)Free-swimmingOff-axisAngular spectrum sharpness and peak searching21
Algae (5 μm)Free-swimmingIn-lineKirchhoff-Helmholtz transform sharpness and peak searching95
Bacteria (<1.5  μm)Free-swimmingIn-lineKirchhoff-Helmholtz transform sharpness and peak searching83
Cochlodinium polykrikoides cell (25 to 40 μm)Free-swimmingIn-lineAngular spectrum sharpness and peak searching97
Red blood cellFlow-throughIn-lineAngular spectrum sharpness and peak searching90
SedimentationDH-QPMFresnel–Kirchhoff integral sharpness and peak searching91
Predatory dinoflagellatesHelical motionIn-lineKirchhoff-Helmholtz transform sharpness and peak searching96
Human spermsHelical motionTwo-set in-line with LEDsAngular spectrum, iterative phase recovery sharpness and peak searching98
Osteosarcoma cellsCrawling on matrix gelsDH-QPMFresnel–Kirchhoff integral sharpness and peak searching100
HT-1080 fibrosarcoma cellsCrawling on matrix gelsDH-QPM with LED illuminationFresnel–Kirchhoff integral sharpness and peak searching102
MicrofibersStraight, suspendedIn-lineFresnel–Kirchhoff integral sharpness and peak searching103
Curved, suspendedOff-axisAngular spectrum sharpness and peak searching21
Copepod, microsphereInteraction of cell and flow fieldIn-lineFresnel-Huygens sharpness and peak searching105
Microfibers, parameciumInteraction of fiber and cellOff-axisAngular spectrum sharpness and peak searching21

Note: DH-QPM, quantitative phase microscopy by digital holography.

The past years have seen DHM develop to such a state that it has become a potent tool to be considered in the various applications. Indeed, DHM revitalizes the general area of holography and is seen to be a breakthrough in microscopic imaging and application. Researchers not only in optical field, but also in microbiology, medicine, marine, metrology, etc., are realizing the new capabilities of DHM. Although it is necessary to develop and improve the performance, such as noise suppression, resolution enhancement, to name but one, we believe that with the development of optical imaging and computing techniques, DHM will be more diversified to solve different specific problems with appropriate strategies, and it deserves special attention to the future contributions in the field of energy, chemical activity, clinical behavior, remote sensing, environmental monitoring, etc.


1. R. J. Adrian, “Twenty years of particle image velocimetry,” Exp. Fluids 39(2), 159–169 (2005).EXFLDU0723-4864 http://dx.doi.org/10.1007/s00348-005-0991-7 Google Scholar

2. A. NaqwiF. DurstG. Kraft, “Sizing of submicrometer particles using a phase-Doppler system,” Appl. Opt. 30(33), 4903–4913 (1991).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.30.004903 Google Scholar

3. F. Pereiraet al., “Defocusing digital particle image velocimetry: a 3-component 3-dimensional DPIV measurement technique. Application to bubbly flows,” Exp. Fluids 29(1), S78–S84 (2000).EXFLDU0723-4864 http://dx.doi.org/10.1007/s003480070010 Google Scholar

4. K. MishimaT. Hibiki, “Development of high-frame-rate neutron radiography and quantitative measurement method for multiphase flow research,” Nucl. Eng. Des. 184(2–3), 183–201 (1998).NEDEAU0029-5493 http://dx.doi.org/10.1016/S0029-5493(98)00196-4 Google Scholar

5. C. D. MeinhartS. T. WereleyJ. G. Santiago, “Micron- resolution velocimetry techniques,” in Developments in Laser Techniques and Applications to Fluid Mechanics, R. J. Adrianet al., Eds., pp. 57–70, Springer, Berlin, Germany (1998). Google Scholar

6. Ch. Brücker, “3D scanning PIV applied to an air flow in a motored engine using digital high-speed video,” Meas. Sci. Technol. 8(12), 1480–1492 (1997).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/8/12/011 Google Scholar

7. J. P. Brodyet al., “Biotechnology at low Reynolds numbers,” Biophys. J. 71(6), 3430–3441 (1996).BIOJAU0006-3495 http://dx.doi.org/10.1016/S0006-3495(96)79538-3 Google Scholar

8. A. Hatchet al., “A rapid diffusion immunoassay in a T-Sensor,” Nat. Biotechnol. 19(5), 461–465 (2001).NABIF91087-0156 http://dx.doi.org/10.1038/88135 Google Scholar

9. J. Yanget al., “Cell separation on microfabricated electrodes using dielectrophoretic/gravitational field flow fractionation,” Anal. Chem. 71(5), 911–918 (1999).ANCHAM0003-2700 http://dx.doi.org/10.1021/ac981250p Google Scholar

10. J. G. Santiagoet al., “A particle image velocimetry system for microfluidics,” Exp. Fluids 25(4), 316–319 (1998).EXFLDU0723-4864 http://dx.doi.org/10.1007/s003480050235 Google Scholar

11. C. D. MeinhartS. T. WereleyJ. G. Santiago, “PIV measurements of a microchannel flow,” Exp. Fluids 27(5), 414–419 (1999).EXFLDU0723-4864 http://dx.doi.org/10.1007/s003480050366 Google Scholar

12. D. DeAngelisL. Gross, Individual-Based Models and Approaches in Ecology, Chapman & Hall, New York (1992). Google Scholar

13. W. FennelT. Neumann, Introduction to the Modelling of Marine Ecosystems, Vol. 72 of the Elsevier Oceanography Series, Elsevier, San Diego, California (2004). Google Scholar

14. W. FennelT. Osborn, “A unifying framework for marine ecological model comparison,” Deep Sea Research II, pp. 1334–1357, Elsevier Science B.V., Amsterdam (2005). Google Scholar

15. C. Decaesteckeret al., “Can anti-migratory drugs be screened in vitro? A review of 2D and 3D assays for the quantitative analysis of cell migration,” Med. Res. Rev. 27(2), 149–176 (2007).MRREDD1098-1128 http://dx.doi.org/10.1002/(ISSN)1098-1128 Google Scholar

16. E. Cukiermanet al., “Taking cell-matrix adhesions to the third dimension,” Science 294(5547), 1708–1712 (2001).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.1064829 Google Scholar

17. P. FriedlK. Wolf, “Tumour-cell invasion and migration: diversity and escape mechanisms,” Nat. Rev. Cancer 3, 362–374 (2003).NRCAC41474-175X http://dx.doi.org/10.1038/nrc1075 Google Scholar

18. F. Sabehet al., “Tumor cell traffic through the extracellular matrix is controlled by the membrane-anchored collagenase MT1-MMP,” J. Cell Biol. 167, 769–781 (2004).JCLBA30021-9525 http://dx.doi.org/10.1083/jcb.200408028 Google Scholar

19. K. A. BeningoM. DemboY. L. Wang, “Responses of fibroblasts to anchorage of dorsal extracellular matrix receptors,” Proc. Natl. Acad. Sci. U. S. A. 101(52), 18024–18029 (2004).PNASA60027-8424 http://dx.doi.org/10.1073/pnas.0405747102 Google Scholar

20. K. Camphausenet al., “Influence of in vivo growth on human glioma cell line gene expression: convergent profiles under orthotopic conditions,” Proc. Natl. Acad. Sci. U. S. A. 102(23), 8287–8292 (2005).PNASA60027-8424 http://dx.doi.org/10.1073/pnas.0502887102 Google Scholar

21. X. Yuet al., “Four-dimensional motility tracking of biological cells by digital holographic microscopy,” J. Biomed. Opt. (2014) (accepted). Google Scholar

22. E. Topraket al., “Three-dimensional particle tracking via bifocal imaging,” Nano Lett. 7(7), 2043–2045 (2007).NALEFD1530-6984 http://dx.doi.org/10.1021/nl0709120 Google Scholar

23. M. D. McMahonet al., “3D particle trajectories observed by orthogonal tracking microscopy,” ACS Nano 3(3), 609–614 (2009).1936-0851 http://dx.doi.org/10.1021/nn8008036 Google Scholar

24. S. Ramet al., “Improved single particle localization accuracy with dual objective multifocal plane microscopy,” Opt. Express 17(8), 6881–6898 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.006881 Google Scholar

25. M. J. Mlodzianoskiet al., “Experimental characterization of 3D localization techniques for particle-tracking and super-resolution microscopy,” Opt. Express 17(10), 8264–8277 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.008264 Google Scholar

26. H. YangN. HalliwellJ. Coupland, “Application of the digital shearing method to extract three-component velocity in holographic particle image velocimetry,” Meas. Sci. Technol. 15(4), 694–698 (2004).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/15/4/011 Google Scholar

27. H. P. KaoA. S. Verkman, “Tracking of single fluorescent particles in three dimensions: use of cylindrical optics to encode particle position,” Biophys. J. 67(3), 1291–1300 (1994).BIOJAU0006-3495 http://dx.doi.org/10.1016/S0006-3495(94)80601-0 Google Scholar

28. S. R. P. PavaniR. Piestun, “Three dimensional tracking of fluorescent microparticles using a photon-limited double-helix response system,” Opt. Express 16(26), 22048–22057 (2008).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.16.022048 Google Scholar

29. D. J. StephensV. J. Allan, “Light microscopy techniques for live cell imaging,” Science 300(5616), 82–86 (2003).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.1082160 Google Scholar

30. M. T. WeiA. Chiou, “Three-dimensional tracking of Brownian motion of a particle trapped in optical tweezers with a pair of orthogonal tracking beams and the determination of the associated optical force constants,” Opt. Express 13(15), 5798–5806 (2005).OPEXFF1094-4087 http://dx.doi.org/10.1364/OPEX.13.005798 Google Scholar

31. A. RohrbachH. KressE. H. Stelzer, “Three-dimensional tracking of small spheres in focused laser beams: influence of the detection angular aperture,” Opt. Lett. 28(6), 411–413 (2003).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.28.000411 Google Scholar

32. J. Berglundet al., “Fast, bias-free algorithm for tracking single particles with variable size and shape,” Opt. Express 16(18), 14064–14075 (2008).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.16.014064 Google Scholar

33. Y. Parket al., “Fresnel particle tracing in three dimensions using diffraction phase microscopy,” Opt. Lett. 32(7), 811–813 (2007).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.32.000811 Google Scholar

34. W. DenkJ. H. StricklerW. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science 248(4951), 73–76 (1990).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.2321027 Google Scholar

35. S. O. Isikmanet al., “Optofluidic tomography on a chip,” Appl. Phys. Lett. 98(16), 161109 (2011).APPLAB0003-6951 http://dx.doi.org/10.1063/1.3548564 Google Scholar

36. D. Huanget al., “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.1957169 Google Scholar

37. D. R. Soll, “Computer-assisted three-dimensional reconstruction and motion analysis of living, crawling cells,” Comput. Med. Imaging Graph. 23(1), 3–14 (1999).CMIGEY0895-6111 http://dx.doi.org/10.1016/S0895-6111(98)00058-5 Google Scholar

38. J. Sharpeet al., “Optical projection tomography as a tool for 3D microscopy and gene expression studies,” Science 296(5567), 541–545 (2002).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.1068206 Google Scholar

39. J. Huiskenet al., “Optical sectioning deep inside live embryos by selective plane illumination microscopy,” Science 305(5686), 1007–1009 (2004).SCIEAS0036-8075 http://dx.doi.org/10.1126/science.1100035 Google Scholar

40. T. ZhangI. Yamaguchi, “Three-dimensional microscopy with phase-shifting digital holography,” Opt. Lett. 23(15), 1221–1223 (1998).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.23.001221 Google Scholar

41. P. Massatschet al., “Time-domain optical coherence tomography with digital holographic microscopy,” Appl. Opt. 44(10), 1806–1812 (2005).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.44.001806 Google Scholar

42. S. TamanoY. HayasakiN. Nishida, “Phase-shifting digital holography with a low-coherence light source for reconstruction of a digital relief object hidden behind a light scattering medium,” Appl. Opt. 45(5), 953–959 (2006).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.45.000953 Google Scholar

43. M. Tzirakiet al., “Photorefractive holography for imaging through turbid media using low coherence light,” Appl. Phys. B 70(1), 151–154 (2000).APBOEM0946-2171 http://dx.doi.org/10.1007/s003400050023 Google Scholar

44. M. K. Kim, “Principles and techniques of digital holographic microscopy,” SPIE Rev. 1(1), 018005 (2010). http://dx.doi.org/10.1117/6.0000006 Google Scholar

45. M. K. Kim, Digital Holographic Microscopy Principles, Techniques, and Applications, Springer Series in Optical Sciences, New York (2011). Google Scholar

46. K. D. Hinsch, “Holographic particle image velocimetry,” Meas. Sci. Technol. 13(7), R61–R72 (2002).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/13/7/201 Google Scholar

47. Y. PuH. Meng, “Four-dimensional dynamic flow measurement by holographic particle image velocimetry,” Appl. Opt. 44(36), 7697–7708 (2005).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.44.007697 Google Scholar

48. D. H. BarnhartR. J. AdrianG. C. Papen, “Phase-conjugate holographic system for high-resolution particle image velocimetry,” Appl. Opt. 33(30), 7159–7170 (1994).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.33.007159 Google Scholar

49. B. TaoJ. KatzC. Meneveau, “Statistical geometry of subgridscale stresses determined from holographic particle image velocimetry measurements,” J. Fluid Mech. 457, 35–78 (2002).JFLSA70022-1120 http://dx.doi.org/10.1017/S0022112001007443 Google Scholar

50. J. ZhangB. TaoJ. Katz, “Turbulent flow measurement in a square duct with hybrid holographic PIV,” Exp. Fluids 23(5), 373–381 (1997).EXFLDU0723-4864 http://dx.doi.org/10.1007/s003480050124 Google Scholar

51. J. C. NauenG. V. Lauder, “Quantification of the wake of rainbow trout (Oncorhynchus mykiss) using three-dimensional stereoscopic digital particle image velocimetry,” J. Exp. Biol. 205(Pt 21), 3271–3279 (2002).JEBIAM0022-0949 Google Scholar

52. A. K. Prasad, “Stereoscopic particle image velocimetry,” Exp. Fluids 29(2), 103–116 (2000).EXFLDU0723-4864 http://dx.doi.org/10.1007/s003480000143 Google Scholar

53. C. Brucker, “3D scanning PIV applied to an air flow in a motored engine using digital high-speed video,” Meas. Sci. Technol. 8(12), 1480–1492 (1997).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/8/12/011 Google Scholar

54. D. Gabor, “Microscopy by reconstructed wavefronts,” Proc. R. Soc. London Ser. A 197(1051), 454–487 (1949).PRLAAZ0080-4630 http://dx.doi.org/10.1098/rspa.1949.0075 Google Scholar

55. J. W. Goodman, Introduction to Fourier Optics, 3rd ed., Roberts and Company Publishers, Englewood, Colorado (2005). Google Scholar

56. E. CucheP. MarquetC. Depeursinge, “Simultaneous amplitude contrast and quantitative phase-contrast microscopy by numerical reconstruction of Fresnel off-axis holograms,” Appl. Opt. 38(34), 6994–7001 (1999).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.38.006994 Google Scholar

57. D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948).NATUAS0028-0836 http://dx.doi.org/10.1038/161777a0 Google Scholar

58. X. Yuet al., “Measurement of the traction force of biological cells by digital holography,” Biomed. Opt. Express 3(1), 153–159 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.000153 Google Scholar

59. M. DemboY. L. Wang, “Stresses at the cell-to-substrate interface during locomotion of fibroblasts,” Biophys. J. 76(4), 2307–2316 (1999).BIOJAU0006-3495 http://dx.doi.org/10.1016/S0006-3495(99)77386-8 Google Scholar

60. X. Yuet al., “Measurement of Young’s modulus of polyacrylamide gel by digital holography,” presented at Digital Holography and Three-Dimensional Imaging, OSA Techinal Digest (CD), paper DTuC32, OSA, Washington, DC (2011). Google Scholar

61. X. Yuet al., “Quantitative imaging and measurement of cell-substrate surface deformation by digital holography,” J. Mod. Opt. 59(18), 1591–1598 (2012).JMOPEW0950-0340 http://dx.doi.org/10.1080/09500340.2012.729095 Google Scholar

62. P. Memmoloet al., “Twin-beams digital holography for 3D tracking and quantitative phase-contrast microscopy in microfluidics,” Opt. Express 19(25), 25833–25842 (2011).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.19.025833 Google Scholar

63. F. Duboiset al., “Focus plane detection criteria in digital holography microscopy by amplitude analysis,” Opt. Express 14(13), 5895–5908 (2006).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.14.005895 Google Scholar

64. P. LangehanenbergB. KemperG. von Bally, “Autofocus algorithms for digital-holographic microscopy,” Proc. SPIE 6633, 66330E (2007).PSISDG0277-786X http://dx.doi.org/10.1117/12.727784 Google Scholar

65. P. Langehanenberget al., “Autofocusing in digital holographic phase contrast microscopy on pure phase objects for life cell imaging,” Appl. Opt. 47(19), D176–D182 (2008).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.47.00D176 Google Scholar

66. J. Lükeet al., “Near real-time estimation of super-resolved depth and all-in-focus images from a plenoptic camera using graphics processing units,” Int. J. Digital Multimedia Broadcast. 2010, 942037 (2010).1687-7578 http://dx.doi.org/10.1155/2010/942037 Google Scholar

67. P. J. Shaw, “Comparison of wide-field/deconvolution and confocal microscopy for 3D imaging,” in Handbook of Biological Confocal Microscopy, J. B. Pawley, Ed., pp. 373–387, Plenum Press, New York and London (1995). Google Scholar

68. J. G. McNallyet al., “Three-dimensional imaging by deconvolution microscopy,” Methods 19(3), 373–385 (1999).MTHDE91046-2023 http://dx.doi.org/10.1006/meth.1999.0873 Google Scholar

69. P. SarderA. Nehorai, “Deconvolution methods for 3-D fluorescence microscopy images,” IEEE Signal Process. Mag. 23(3), 32–45 (2006).ISPRE61053-5888 http://dx.doi.org/10.1109/MSP.2006.1628876 Google Scholar

70. W. WallaceL. H. SchaeferJ. R. Swedlow, “A working person’s guide to deconvolution in light microscopy,” Biotechniques 31(5), 1076–1097 (2001).BTNQDO0736-6205 Google Scholar

71. T. LatychevskaiaF. GehriH. Fink, “Depth-resolved holographic reconstructions by three-dimensional deconvolution,” Opt. Express 18(21), 22527–22544 (2010).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.18.022527 Google Scholar

72. L. WilsonR. Zhang, “3D localization of weak scatterers in digital holographic microscopy using Rayleigh-Sommerfeld back-propagation,” Opt. Express 20(15), 16735–16744 (2012).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.20.016735 Google Scholar

73. S.-H. LeeD. G. Grier, “Holographic microscopy of holographically trapped three-dimensional structures,” Opt. Express 15(4), 1505–1512 (2007).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.15.001505 Google Scholar

74. M. BornE. Wolf, Principles of Optics, 6th ed., Cambridge University Press, Cambridge, UK (1998). Google Scholar

75. D. J. Bradyet al., “Compressive holography,” Opt. Express 17(15), 13040–13049 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.013040 Google Scholar

76. J. Hahnet al., “Video-rate compressive holographic microscopic tomography,” Opt. Express 19(8), 7289–7298 (2011).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.19.007289 Google Scholar

77. E. J. CandèsJ. K. RombergT. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).CPMAMV0010-3640 http://dx.doi.org/10.1002/(ISSN)1097-0312 Google Scholar

78. E. J. CandesT. Tao, “Near-optimal signal recovery from random projections: universal encoding strategies?,” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).IETTAW0018-9448 http://dx.doi.org/10.1109/TIT.2006.885507 Google Scholar

79. G. Coppolaet al., “Digital self-referencing quantitative phase microscopy by wavefront folding in holographic image reconstruction,” Opt. Lett. 35(20), 3390–3392 (2010).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.35.003390 Google Scholar

80. F. Merolaet al., “Driving and analysis of micro-objects by digital holographic microscope in microfluidics,” Opt. Lett. 36(16), 3079–3081 (2011).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.36.003079 Google Scholar

81. W. Xuet al., “Digital in-line holography for biological applications,” Proc. Natl. Acad. Sci. U. S. A. 98(20), 11,301–11,305 (2001).1091-6490 http://dx.doi.org/10.1073/pnas.191361398 Google Scholar

82. W. Xuet al., “Digital in-line holography of microspheres,” Appl. Opt. 41(25), 5367–5375 (2002).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.41.005367 Google Scholar

83. J. G. Sucerquiaet al., “Digital in-line holographic microscopy,” Appl. Opt. 45(5), 836–850 (2006).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.45.000836 Google Scholar

84. H. J. Kreuzeret al., “Theory of the point source electron microscope,” Ultramicroscopy 45(3–4), 381–403 (1992).ULTRD60304-3991 http://dx.doi.org/10.1016/0304-3991(92)90150-I Google Scholar

85. J. ShengE. MalkielJ. Katz, “Digital holographic microscope for measuring three-dimensional particle distributions and motions,” Appl. Opt. 45(16), 3893–3901 (2006).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.45.003893 Google Scholar

86. J. ShengE. MalkielJ. Katz, “Single beam two-views holographic particle image velocimetry,” Appl. Opt. 42(2), 235–250 (2003).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.42.000235 Google Scholar

87. J. Garcia-Sucerquiaet al., “4D imaging of fluid flow with digital in-line holographic microscopy,” Optik 119(9), 419–423 (2008).OTIKAJ0030-4026 http://dx.doi.org/10.1016/j.ijleo.2007.01.004 Google Scholar

88. L. Tianet al., “Quantitative measurement of size and three-dimensional position of fast-moving bubbles in air-water mixture flows using digital holography,” Appl. Opt. 49(9), 1549–1554 (2010).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.49.001549 Google Scholar

89. D. ReynoldsR. Rose, “Robust text-independent speaker identification using Gaussian mixture speaker models,” IEEE Trans. Speech Audio Process. 3(1), 72–83 (1995).IESPEJ1063-6676 http://dx.doi.org/10.1109/89.365379 Google Scholar

90. Y.-S. ChoiS.-J. Lee, “Three-dimensional volumetric measurement of red blood cell motion using digital holographic microscopy,” Appl. Opt. 48(16), 2983–2990 (2009).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.48.002983 Google Scholar

91. P. Langehanenberget al., “Automated three-dimensional tracking of living cells by digital holographic microscopy,” J. Biomed. Opt. 14(1), 014018 (2009).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3080133 Google Scholar

92. D. Carlet al., “Parameteroptimized digital holographic microscope for high-resolution living cell analysis,” Appl. Opt. 43(36), 6536–6544 (2004).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.43.006536 Google Scholar

93. B. Kemperet al., “Investigation of living pancreas tumor cells by digital holographic microscopy.” J. Biomed. Opt. 11(3), 034005 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2204609 Google Scholar

94. U. SchnarsW. Jüptner, “Digital recording and numerical reconstruction of holograms,” Meas. Sci. Technol. 13(9), R85–R101 (2002).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/13/9/201 Google Scholar

95. W. Xuet al., “Tracking particles in four dimensions with in-line holographic microscopy,” Opt. Lett. 28(3), 164–166 (2003).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.28.000164 Google Scholar

96. J. Shenget al., “Digital holographic microscopy reveals prey-induced changes in swimming behavior of predatory dinoflagellates,” Proc. Natl. Acad. Sci. U. S. A. 104(44), 17512–17517 (2007).1091-6490 http://dx.doi.org/10.1073/pnas.0704658104 Google Scholar

97. S. Leeet al., “Three-dimensional motion measurements of free-swimming microorganisms using digital holographic microscopy,” Meas. Sci. Technol. 22(6), 064004 (2011).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/22/6/064004 Google Scholar

98. T. W. SuL. XueA. Ozcan, “High-throughput lensfree 3D tracking of human sperms reveals rare statistics of helical trajectories,” Proc. Natl. Acad. Sci. U. S. A. 109(40), 16018–16022 (2012).1091-6490 http://dx.doi.org/10.1073/pnas.1212506109 Google Scholar

99. O. Mudanyaliet al., “Compact, light-weight and cost-effective microscope based on lensless incoherent holography for telemedicine applications,” Lab Chip 10(11), 1417–1428 (2010).LCAHAM1473-0197 http://dx.doi.org/10.1039/c000453g Google Scholar

100. P. Memmoloet al., “On the holographic 3D tracking of in vitro cells characterized by a highly-morphological change,” Opt. Express 20(27), 28485–28493 (2012).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.20.028485 Google Scholar

101. J. F. RestrepoJ. Garcia-Sucerquia, “Automatic three-dimensional tracking of particles with high-numerical-aperture digital lensless holographic microscopy,” Opt. Lett. 37(4), 752–754 (2012).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.37.000752 Google Scholar

102. F. Duboiset al., “Digital holographic microscopy for the three-dimensional dynamic analysis of in vitro cancer cell migration,” J. Biomed. Opt. 11(5), 054032 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2357174 Google Scholar

103. M. Kempkeset al., “Three dimensional digital holographic profiling of micro-fibers,” Opt. Express 17(4), 2938–2943 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.002938 Google Scholar

104. D. T. HaynieD. B. KhadkaM. C. Cross, “Physical properties of polypeptide electrospun nanofiber cell culture scaffolds on a wettable substrate,” Polymers 4(3), 1535–1553 (2012).2073-4360 http://dx.doi.org/10.3390/polym4031535 Google Scholar

105. E. Malkielet al., “The three-dimensional flow field generated by a feeding calanoid copepod measured using digital holography,” J. Exp. Biol. 206(Pt 20), 3657–3666 (2003).JEBIAM0022-0949 http://dx.doi.org/10.1242/jeb.00586 Google Scholar

106. E. Hecht, Optics, Addison-Wesley, Reading, MA (2002). Google Scholar

107. G. I. RothJ. Katz, “Five techniques for increasing the speed and accuracy of PIV interrogation,” Meas. Sci. Technol. 12(3), 238–245 (2001).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/12/3/302 Google Scholar

108. G. I. RothD. T. MascenikJ. Katz, “Measurements of the flow structure and turbulence within a ship bow wave,” Phys. Fluids 11(11), 3512–3523 (1999).PHFLE61070-6631 http://dx.doi.org/10.1063/1.870209 Google Scholar

109. H. S. JiangT. R. OsbornC. Meneveau, “The flow field around a freely swimming copepod in steady motion. Part I: Theoretical analysis,” J. Plankton Res. 24(3), 167–189 (2002).JPLRD90142-7873 http://dx.doi.org/10.1093/plankt/24.3.167 Google Scholar

110. C. Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow, Cambridge University Press, England, New York (1992). Google Scholar


Xiao Yu received a BS degree in physics from Nankai University, Tianjin, China, in 2008 and MS degree in applied physics from University of South Florida (USF), Tampa, Florida, in 2010. She is currently working toward her PhD degree at USF. Her main research interests are digital holographic microscopy and its applications in biological field, biomedical imaging, digital image processing, etc.

Jisoo Hong received his BS and MS degrees in electrical engineering from Seoul National University, Korea, in 2002 and 2004, respectively. From 2004 to 2008, he was a senior research engineer with LG Electronics, Korea. In August 2012, he received his PhD degree from his alma mater. Presently, he is a postdoctoral researcher in the Department of Physics at USF, Tampa, Florida. He is interested in incoherent digital holography and its application to biomedical imaging and three-dimensional photography.

Changgeng Liu obtained his BS in applied physics in 2007 and his MS in optics in 2010 from Beijing University of Technology. He is currently seeking a PhD degree in applied physics in USF. His research interests include digital holography, adaptive optics, biomedical imaging, computational optics, phase retrieval, digital imaging processing, and optical systems design.

Myung K. Kim is a professor of physics at the USF, and directs the Digital Holography and Microscopy Laboratory. He received his PhD degree from the University of California, Berkeley, in 1986. His current research interests are in the development of novel techniques and applications in digital holography, microscopy, interference imaging, optical tomography, biomedical imaging, and manipulation of particles and microbes by optical force. He has more than 200 publications and 13 patents granted or pending.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Xiao Yu, Xiao Yu, Jisoo Hong, Jisoo Hong, Changgeng Liu, Changgeng Liu, Myung K. Kim, Myung K. Kim, } "Review of digital holographic microscopy for three-dimensional profiling and tracking," Optical Engineering 53(11), 112306 (3 April 2014). https://doi.org/10.1117/1.OE.53.11.112306 . Submission:


Back to Top