Abstract

Compressive holography enables 3D reconstruction from a single 2D holographic snapshot for objects that can be sparsely represented in some basis. The snapshot mode enables tomographic imaging of microscopic moving objects. We demonstrate video-rate tomographic image acquisition of two live water cyclopses with 5.2 μm spatial resolution and 60 μm axial resolution.

© 2011 OSA

1. Introduction

Holography characterizes the optical field, including the amplitude and the phase [1,2], on a 2D surface. Quasi-3D imaging has been achieved using holography via the so-called backpropagation mechanism, which numerically focuses the optical field at several depths, as if lenses focused the field to form an image. However, as discussed in [3], the backpropagation in digital holography recovers merely 2D fields at multiple depths as opposed to the true 3D scattering density that induces the 2D fields. It is possible, however, to form truly tomographic images from holographic data using nonlinear object estimation algorithms. This paper demonstrates the utility of these algorithms for 3D imaging of microscopic moving objects. In particular, we show that tomographic images can be acquired at full frame rate at the longitudinal and transverse resolution limits of the optical system.

Our group has demonstrated that 3D tomographic reconstruction can be achieved from a single holographic projection by using the theory of compressive sensing [46] provided that the representation of the 3D scattering density in some transformation basis is sufficiently sparse. This method, called compressive holography, has a potential for fast acquisition for 3D imaging. Several papers based on a similar philosophy have been published. In particular, Denis et. al. demonstrated success application of sparse reconstruction techniques to holographic 3D imaging and field-of-view expansion [7,8]. Also, Coskun et al. applied a similar sparse estimation method for lenseless wide-field fluorescent imaging on a chip [9]. A similar philosophy used in our compressive holography work has been applied to the 2D holographic imaging as well. Marim et al. and Rivenson et al. applied decompressive inference to frequency-shifting hologram [10] and digital Fresnel holography [11]. It is also noteworthy that other 3D imaging methods such as optical sectioning [12] and Tikhonov-regularized solution [13] are available in the literature even though all these methods are not strictly 3D imaging in that the solutions from the methods do not mathematically correspond to the true 3D scattering densities.

While the resolution of a microscopic system is determined by its numerical aperture (NA), the practical resolving power and reconstruction fidelity of compressive holography largely depend upon the effective NA that is determined by object feature sizes. It is desirable that the diffraction signals fully occupy the detector area such that the NA of an optical system is fully utilized. In practice, the combined effects of these factors may be reasonably estimated from the amounts of diffraction signals produced by the object features. Increasing diffraction signals, and equivalently the effective NA, may be achieved by several ways. For example, a microscopic objective can be used to increase the numerical aperture of a system by magnifying the field entering the objective. Of course, the objective tends to increase system volume and decrease field of view.

Highly compact microscopes may be constructed, as in digital holographic lensless microscopy [14,15], by eliminating optics between the scattering object and the focal plane. This paper focuses on such designs while also obtaining detector limited NA. We illuminate the object with a spherical wave incident field to magnify the diffraction signal on the FPA. Other more sophisticated incident fields may be utilized to further increase the detectable diffraction signals to improve the performance of our methods. Our physical system is similar to the digital in-line holographic microscope [14]; our contribution is to adapt compressive holography to microscopy and to demonstrate its effectiveness in single frame 3D tomography. Building on our first report of this approach in [16], this paper presents a theoretical framework for holographic tomography and demonstrates the utility of the approach for practical imaging. We note that the theoretical framework developed in this paper can be generalized to other microscopic geometries with small modifications which may utilize more complicated incident fields.

This paper is organized as follows. Section 2 describes the mathematical forward model and reconstruction methods for compressive holographic tomography. Section 3 describes an experimental system and presents experimental results. Section 4 concludes with suggestions for future development.

2. System model

We consider a microscopic system composed of a point source for illuminating the objects and the FPA that measures the Gabor hologram. Figure 1 shows a schematic of such a holographic microscope.

 

Fig. 1 A schematic of the compressive holographic microscope.

Download Full Size | PPT Slide | PDF

The intensity of interference pattern I(x,y;z) measured using the FPA placed at a propagation distance z can be written as

I(x,y;zFPA)=|U(x,y;zFPA)|2=|U0(x,y;zFPA)|2+|Us(x,y;zFPA)|2+2Re{U0*(x,y;zFPA)Us(x,y;zFPA)},
where U0 represents the incident field, and Us represents the complex field at distance zFPAscattered off of the objects. The incident field U0 is generated by a point source and thus forms a spherical wave. Thus,
Us(x,y; zFPA)=U(x,y;zFPA)U0(x,y;zFPA)=dxdydzU0(x,y;z)β(x,y,z)h(xx,yy;zFPAz).
Here, (x,y)and (x,y,z)represent the coordinates in the FPA plane and in the object space, respectively. β(x,y,z) is the scattering density of a 3D object and h is the Huygens-Fresnel point-spread function [17]. Since the NA of the system is relatively large, the Fresnel approximation may be inaccurate. Thus, the use of the angular spectrum transfer function for propagation would be desirable. However, the effective NA is smaller than the numerical aperture of the system, which allows us to use the Fresnel approximation without suffering from the loss of numerical accuracy. The concept of the effective NA is discussed in details in Sec. 3.

The incident field U0(x,y;z) is a spherical field originating from z=0 and thus can be expressed usinghF(x,y;z). The scattered wave from the objects is represented as

Us(x,y;zFPA)=dxdydzhF(x,y;z)β(x,y,z)hF(xx,yy;zFPAz)=hF(x,y;zFPA)dxdydzC(z)β(xzzFPA,yzzFPA,z)hF(xx,yy;zFPAzz/zFPA)=hF(x,y;zFPA)dzC(z)F2D1{F2D{β(xzzFPA,yzzFPA,z)}HF(kx,ky;zFPAzz/zFPA)}.
Here, HF(kx,ky;z) is the Fresnel transfer function and F2D and F2D1 denote the 2D Fourier transform and its inverse, respectively. The coordinates (x,y) is defined as x=xzFPA/z and y=yzFPA/z. A uniform phase delay according to z is given by
C(z)=exp[jkzFPA(1zFPA/z)].
Considering that the FPA samples the scattered field, the last term in Eq. (1) can be accordingly discretized as
U0*(n1Δ,n2Δ,zFPA)Us(n1Δ,n2Δ,zFPA)=hF*(n1Δ,n2Δ,zFPA)Us(n1Δ,n2Δ,zFPA)=1N2lC(lΔzzFPA)m1m2n1n2β(n1Δ,n2Δ,lΔz)exp(2πjm1n11+m2n2N)×HF(m1Δk,m2Δk,zFPAlΔzlΔz/zFPA)exp(2πjn1m1+n2m2N),
where Δ denotes the sampling spacing in the FPA plane and Δ denotes the sampling spacing in the object space. Δz is the sampling spacing in the axial direction in the object space. The sample spacings in the Fourier domain Δk satisfy different relations: ΔΔk=2π and ΔΔk=2πz/zFPA.

Equation (5) can be rewritten as

U0,n1n2*Us,n1n2=lClF2D1{F2D{βn1n2l}m1m2HF,m1m2l}n1n2.
Here, βn1n2l=β(n1Δ,n2Δ,lΔz)and HF,m1m2l=HF(m1Δk,m2Δk,zFPA(zFPAlΔz)/lΔz). Thus, considering Eq. (1) as in [3], Eq. (6) can be algebraically expressed as
g=2Re{Hf}+e+n,
where fand g denote vectorized versions ofβn1n2land U0,n1n2*Us,n1n2, and H denotes a measurement matrix whose element is given by Hij=[CF2D1HFF2D]ij. Here, HF denotes a diagonal matrix representing the Fresnel transfer function and C is a matrix representing C(z) in Eq. (4). The term e represents the measurement error resulting from the autocorrelation |U0(x,y;zFPA)|2+|Us(x,y;zFPA)|2, and n denotes additive noise.

Equation (7) is solved by finding f that minimizes the total variation (TV) [18]. Minimizing the total variation (TV) is equivalent to enforcing the sparsity of f in the variational domain. The minimum TV estimate f can be obtained by solving

f^=argminf12g2Re(Hf)2+τfTV.
In this equation, fTV denotes the total variation of f and is defined as
fTV=zxy|(fz)x,y|,
where fz denotes a 2D transverse slice of the 3D object datacube, and |(fz)x,y| represents the magnitude of the gradient at location (x,y)in the z-th transverse slice. We solve the optimization problem in Eq. (9) by adapting the two-step iterative shrinkage/thresholding (TwIST) algorithm [19].

Separation of holographic signal from the background intensity terms (e in Eq. (7)) and from the conjugate image is traditionally challenging for in-line holography. As previously reported in [3], we note that the estimator described in Eq. (8) achieves this separation by localizing the conjugate image appropriately outside the reconstruction volume and by localizing the background terms in the z=0 plane. One expects this separation because the total variation of the propagated background terms is minimized in the z=0 plane where these terms are real (thus eliminating variation in the complex component).

3. Experimental results and discussion

We constructed a holographic microscope using a He-Ne laser with 632.8nm wavelength as a light source. A microscope objective with 0.65NA manufactured by LOMO is used to generate a spherical wave for illuminating the sample. While the illumination source is not compact, we note that alternative sources using fiber or holographic components could achieve similar NA. The microscope objective is chosen to have a larger NA (0.65) than the microscope system NA (0.27) to ensure uniform illumination intensity on the FPA. A Lumenera CMOS sensor records the hologram. The sensor has 1280×1024 resolution, 5.2μmpixel pitch, 10bit digitization, and the maximum frame rate 15fps. Figure 2(a) shows a photograph of the microscope instrumentation. The relative positions of the microscope objective and the sample are adjusted by two 3-axis stages. Figure 2(b) shows a photograph of the container in which two live water cyclopses of roughly the same size are floating in water. The depth of the container is 3 mm. The distance from the FPA to the microscope objective for illumination is 10mm.

 

Fig. 2 Photographs of (a) the holographic microscope and (b) a water container in which water cyclopses are swimming.

Download Full Size | PPT Slide | PDF

Figure 3 shows a raw Gabor hologram. From this raw measurement, the spherical-wave incident field intensity |U0(x,y;zFPA)|2 is removed by simply deleting DC term in Fourier transform. Approximately 20 pixels along each dimension around the boundaries are set to zero, and only 1240×984 data are used to avoid wraparound effects. This preprocessed data is defined as g in the algorithm in Eq. (8) to reconstruct the object scattering density f(or equivalently, β).

 

Fig. 3 Raw image of a Gabor hologram.

Download Full Size | PPT Slide | PDF

Figure 4 compares the backpropagation reconstructions with compressive holographic reconstructions obtained with a single 2D measurement of the two live water cyclopses. Figures 4(a) and 4(b) show transverse slices of the backpropagation reconstructions at z=3.31mm and z=1.87mm , respectively. Figures 4(c) and 4(d) show transverse slices at the same axial positions as those in Figs. 4(a) and 4(b). It is clear that the compressive holography reconstructions show significantly better localization (or sectioning) capability. Also, the compressive holography reconstructions suffer less from the undesired background “noise” resulting in better image contrast. This reflects as better reconstruction fidelity. For example, the tails and antennae of both water cyclopses are remarkably sharper and more discernible in Fig. 4(c) and 4(d) compared to those in Fig. 4(a) and 4(b).

 

Fig. 4 A comparison of reconstructions at chosen axial positions: (a-b) the backpropagation reconstructions and by (c-d) the compressive holographic reconstructions using the data shown in Fig. 3. All the transverse slices of the reconstructions are sequentially shown in (e) obtained by the backpropagation method (Media 1) and in (f) obtained by the compressive holography method (Media 2) for the full range of water container.

Download Full Size | PPT Slide | PDF

The theoretical limits of the resolution of the proposed system may be expressed as [20]

Δx=λ2NA,
Δz=2λNA2,
where the NA is defined as the half width of the FPA over the distance of objects to the FPA plane. The NA is 0.27 in our experiments. From the equations, the theoretical resolution limits are estimated asΔx=1.17μmand Δz=17.3μm. In practice, however, the resolutions are limited by other factors such as the quantization of FPA, various sources of noise, and object feature sizes. These factors reflect on the hologram as discernible diffraction signals. Practical resolutions may thus be estimated by considering amounts of diffraction signals [21]. Note that such resolutions can depend on objects because diffractions depend on object feature sizes. For example, one tail of the water cyclops on the left-hand side in Fig. 3 produced diffractions over approximately 400 pixels each of which pitch is 5.2μm. This results in the effective aperture size De=2.1mm. Since the propagation distance is approximately 7.0mm, the estimated effective numerical aperture NAeis approximately 0.15. This NAe results in Δx2.2μm and Δz59μm. Therefore, in the reconstruction, the axial sample spacing, namely the distance between two adjacent transverse slices, is set to 60μm.

Figure 5 compares the magnified tails, marked by rectangles in Figs. 4(b) and 4(d), of (a) the backpropagation reconstruction and (b) the compressive holography reconstruction. The detailed comparison shows clear difference between the two reconstructions and their changes in the axial direction.

 

Fig. 5 A comparison of the magnified tails, marked by rectangles in Figs. 4(b) and 4(d), of (a) the backpropagation reconstruction and (b) the compressive holography reconstruction.

Download Full Size | PPT Slide | PDF

A 3D visualization of the compressive holography reconstruction is shown in Fig. 6 . The size of this 3D data cube is 1280×1024×35 pixels where the number of transverse slices is 35 obtained by dividing 2mm depth of a container by 60μmresolution. The reconstruction takes about one and a half hours with the code written in Matlab 7.0.4 when the data is processed on Intel Xeon CPU X5650 at 2.67GHz and 24GB of RAM. Readers may notice that two water cyclopses are located in two different depths as shown in Fig. 6.

 

Fig. 6 A 3D visualization of the compressive holography reconstruction (Media 3).

Download Full Size | PPT Slide | PDF

Note that the Born approximation is valid when the product of the index contrast and object size is less than one-quarter wavelength since the field inside the object is approximated using the incident field [22,23]. For example, the platelets of Sapphirinidae that belongs to the same subclass as water cyclopses have a high refractive index as 1.8 [24]. So the bodies of cyclopses are not in the range that the Born approximation is valid and its inner structure are not correctly reconstructed by our method based on the Born approximation. The reconstruction along the propagation direction is more inaccurate as in Fig. 6 since the holographic microscope is modeled as a linear system along the propagation direction.

Figure 7 presents an alternative visualization of the same 3D datacube. Figure 7(a) shows an image of the maximum intensity values along the propagation directions in 3D. Figure 7(b) represents a map of the axial positions corresponding to the maximum values in Fig. 7(a). From these two images, a range map is constructed as shown in Fig. 7(c). The map represents the HSV space. The hue (h) represents Fig. 7(b), and the value (v) represents Fig. 7(a). The saturation (s) is set to 1.

 

Fig. 7 Images of (a) the maximum intensity values of the reconstructed density (f) along the propagation directions and (b) a map of the axial positions corresponding to the maximum values in Fig. 7(a). (c) A range colormap represents the HSV space.

Download Full Size | PPT Slide | PDF

Figure 8 shows the videos of the raw measurements and range maps associated with the compressive holography reconstructions. The raw measurements capture the movement of three water cyclopses in a 3mmdepth container. In Fig. 8(a), the frame rate for the 2D hologram recording was set to 15 fps whose frames are played at the rate of 3 frames per second for readers’ visual convenience. Figure 8(b) shows a video of the range maps of the compressive holography reconstructions whose sizes are 1024×1024×35. This reconstructed video has 109frames and the data processing was performed parallel on a computer cluster of the Scalable Computing Support Center in Duke University. The axial resolution for these reconstructions is set to Δz90μm. A water cyclops around the center is swimming down over the range from 1.6mm to 4.5mm . The movement is illustrated by its color changes according to its depths.

 

Fig. 8 Videos of (a) raw measurements (Media 4) and (b) range maps of the compressive holography reconstructions (Media 5).

Download Full Size | PPT Slide | PDF

4. Conclusion

We have demonstrated a video-rate compressive holographic microscopic tomography. The method has been applied to microscopy where an illumination wave is generated with a microscope objective to produce a spherical wave incident field. This framework for compressive holographic microscopy can be generalized to any computational holographic microscopy with moderate modifications. We have also analyzed both the theoretical and practical resolutions of the proposed microscope.

Our use of the term “video-rate” refers to the acquisition time of the image, processing to reconstruct the tomographic data cube is far from real-time in our current implementation. Future work may focus on the use of application specific software and hardware to speed the reconstruction and visualization process. Alternative sampling strategies using multiple illumination wavelengths or more complex illumination patterns may be effective improving the longitudinal resolution and increasing the range of observable objects.

References and links

1. D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948). [CrossRef]   [PubMed]  

2. E. N. Leith and J. Upatnieks, “Reconstructed wavefronts and communication theory,” J. Opt. Soc. Am. 52(10), 1123–1130 (1962). [CrossRef]  

3. D. J. Brady, K. Choi, D. L. Marks, R. Horisaki, and S. Lim, “Compressive holography,” Opt. Express 17(15), 13040–13049 (2009). [CrossRef]   [PubMed]  

4. E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006). [CrossRef]  

5. E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006). [CrossRef]  

6. D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Theory 52(4), 1289–1306 (2006). [CrossRef]  

7. L. Denis, D. Lorenz, E. Thiébaut, C. Fournier, and D. Trede, “Inline hologram reconstruction with sparsity constraints,” Opt. Lett. 34(22), 3475–3477 (2009). [CrossRef]   [PubMed]  

8. C. Fournier, L. Denis, and T. Fournel, “On the single point resolution of on-axis digital holography,” J. Opt. Soc. Am. A 27(8), 1856–1862 (2010). [CrossRef]  

9. A. F. Coskun, I. Sencan, T.-W. Su, and A. Ozcan, “Lensless wide-field fluorescent imaging on a chip using compressive decoding of sparse objects,” Opt. Express 18(10), 10510–10523 (2010). [CrossRef]   [PubMed]  

10. M. M. Marim, M. Atlan, E. Angelini, and J.-C. Olivo-Marin, “Compressed sensing with off-axis frequency-shifting holography,” Opt. Lett. 35(6), 871–873 (2010). [CrossRef]   [PubMed]  

11. Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010). [CrossRef]  

12. H. Kim, S.-W. Min, B. Lee, and T.-C. Poon, “Optical sectioning for optical scanning holography using phase-space filtering with Wigner distribution functions,” Appl. Opt. 47(19), D164–D175 (2008). [CrossRef]   [PubMed]  

13. E. Y. Lam, X. Zhang, H. Vo, T.-C. Poon, and G. Indebetouw, “Three-dimensional microscopy and sectional image reconstruction using optical scanning holography,” Appl. Opt. 48(34), H113–H119 (2009). [CrossRef]   [PubMed]  

14. J. Garcia-Sucerquia, W. Xu, S. K. Jericho, P. Klages, M. H. Jericho, and H. J. Kreuzer, “Digital in-line holographic microscopy,” Appl. Opt. 45(5), 836–850 (2006). [CrossRef]   [PubMed]  

15. S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009). [CrossRef]   [PubMed]  

16. J. Hahn, S. Lim, K. Choi, R. Horisaki, D. L. Marks, and D. J. Brady, “Compressive Holographic Microscopy,” in Biomedical Optics, OSA Technical Digest (CD) (Optical Society of America, 2010), paper JMA1.

17. R. E. Blahurt, Theory of Remote Image Formation (Cambridge University Press, 2005).

18. L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992). [CrossRef]  

19. J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007). [CrossRef]   [PubMed]  

20. D. J. Brady, Optical Imaging and Spectroscopy (Wiley, 2009).

21. D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009). [CrossRef]  

22. A. C. Kak, and M. Slaney, Principles of Computerized Tomographic Imaging (SIAM, 2001).

23. D. L. Marks, “A family of approximations spanning the Born and Rytov scattering series,” Opt. Express 14(19), 8837–8848 (2006). [CrossRef]   [PubMed]  

24. J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948).
    [CrossRef] [PubMed]
  2. E. N. Leith and J. Upatnieks, “Reconstructed wavefronts and communication theory,” J. Opt. Soc. Am. 52(10), 1123–1130 (1962).
    [CrossRef]
  3. D. J. Brady, K. Choi, D. L. Marks, R. Horisaki, and S. Lim, “Compressive holography,” Opt. Express 17(15), 13040–13049 (2009).
    [CrossRef] [PubMed]
  4. E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
    [CrossRef]
  5. E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).
    [CrossRef]
  6. D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Theory 52(4), 1289–1306 (2006).
    [CrossRef]
  7. L. Denis, D. Lorenz, E. Thiébaut, C. Fournier, and D. Trede, “Inline hologram reconstruction with sparsity constraints,” Opt. Lett. 34(22), 3475–3477 (2009).
    [CrossRef] [PubMed]
  8. C. Fournier, L. Denis, and T. Fournel, “On the single point resolution of on-axis digital holography,” J. Opt. Soc. Am. A 27(8), 1856–1862 (2010).
    [CrossRef]
  9. A. F. Coskun, I. Sencan, T.-W. Su, and A. Ozcan, “Lensless wide-field fluorescent imaging on a chip using compressive decoding of sparse objects,” Opt. Express 18(10), 10510–10523 (2010).
    [CrossRef] [PubMed]
  10. M. M. Marim, M. Atlan, E. Angelini, and J.-C. Olivo-Marin, “Compressed sensing with off-axis frequency-shifting holography,” Opt. Lett. 35(6), 871–873 (2010).
    [CrossRef] [PubMed]
  11. Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010).
    [CrossRef]
  12. H. Kim, S.-W. Min, B. Lee, and T.-C. Poon, “Optical sectioning for optical scanning holography using phase-space filtering with Wigner distribution functions,” Appl. Opt. 47(19), D164–D175 (2008).
    [CrossRef] [PubMed]
  13. E. Y. Lam, X. Zhang, H. Vo, T.-C. Poon, and G. Indebetouw, “Three-dimensional microscopy and sectional image reconstruction using optical scanning holography,” Appl. Opt. 48(34), H113–H119 (2009).
    [CrossRef] [PubMed]
  14. J. Garcia-Sucerquia, W. Xu, S. K. Jericho, P. Klages, M. H. Jericho, and H. J. Kreuzer, “Digital in-line holographic microscopy,” Appl. Opt. 45(5), 836–850 (2006).
    [CrossRef] [PubMed]
  15. S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
    [CrossRef] [PubMed]
  16. J. Hahn, S. Lim, K. Choi, R. Horisaki, D. L. Marks, and D. J. Brady, “Compressive Holographic Microscopy,” in Biomedical Optics, OSA Technical Digest (CD) (Optical Society of America, 2010), paper JMA1.
  17. R. E. Blahurt, Theory of Remote Image Formation (Cambridge University Press, 2005).
  18. L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
    [CrossRef]
  19. J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007).
    [CrossRef] [PubMed]
  20. D. J. Brady, Optical Imaging and Spectroscopy (Wiley, 2009).
  21. D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
    [CrossRef]
  22. A. C. Kak, and M. Slaney, Principles of Computerized Tomographic Imaging (SIAM, 2001).
  23. D. L. Marks, “A family of approximations spanning the Born and Rytov scattering series,” Opt. Express 14(19), 8837–8848 (2006).
    [CrossRef] [PubMed]
  24. J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994).
    [CrossRef]

2010

2009

2008

2007

J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007).
[CrossRef] [PubMed]

2006

J. Garcia-Sucerquia, W. Xu, S. K. Jericho, P. Klages, M. H. Jericho, and H. J. Kreuzer, “Digital in-line holographic microscopy,” Appl. Opt. 45(5), 836–850 (2006).
[CrossRef] [PubMed]

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
[CrossRef]

E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).
[CrossRef]

D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Theory 52(4), 1289–1306 (2006).
[CrossRef]

D. L. Marks, “A family of approximations spanning the Born and Rytov scattering series,” Opt. Express 14(19), 8837–8848 (2006).
[CrossRef] [PubMed]

1994

J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994).
[CrossRef]

1992

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
[CrossRef]

1962

1948

D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948).
[CrossRef] [PubMed]

Angelini, E.

Atlan, M.

Bioucas-Dias, J. M.

J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007).
[CrossRef] [PubMed]

Brady, D. J.

Candes, E. J.

E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).
[CrossRef]

Candès, E. J.

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
[CrossRef]

Chae, J.

J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994).
[CrossRef]

Choi, K.

Coskun, A. F.

Denis, L.

Donoho, D. L.

D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Theory 52(4), 1289–1306 (2006).
[CrossRef]

Erlinger, A.

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Fatemi, E.

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
[CrossRef]

Figueiredo, M. A. T.

J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007).
[CrossRef] [PubMed]

Fournel, T.

Fournier, C.

Gabor, D.

D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948).
[CrossRef] [PubMed]

Garcia-Sucerquia, J.

Hennelly, B. M.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Horisaki, R.

Indebetouw, G.

Javidi, B.

Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010).
[CrossRef]

Jericho, M. H.

Jericho, S. K.

Kelly, D. P.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Kim, H.

Klages, P.

Kreuzer, H. J.

Lam, E. Y.

Lee, B.

Leith, E. N.

Lim, S.

Lorenz, D.

Marim, M. M.

Marks, D. L.

Min, S.-W.

Naughton, T. J.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Nishida, S.

J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994).
[CrossRef]

Olivo-Marin, J.-C.

Osher, S.

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
[CrossRef]

Ozcan, A.

A. F. Coskun, I. Sencan, T.-W. Su, and A. Ozcan, “Lensless wide-field fluorescent imaging on a chip using compressive decoding of sparse objects,” Opt. Express 18(10), 10510–10523 (2010).
[CrossRef] [PubMed]

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Pandey, N.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Poon, T.-C.

Rhodes, W. T.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Rivenson, Y.

Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010).
[CrossRef]

Romberg, J. K.

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
[CrossRef]

Rudin, L. I.

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
[CrossRef]

Sencan, I.

Seo, S.

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Stern, A.

Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010).
[CrossRef]

Su, T. W.

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Su, T.-W.

Tao, T.

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
[CrossRef]

E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).
[CrossRef]

Thiébaut, E.

Trede, D.

Tseng, D. K.

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Upatnieks, J.

Vo, H.

Xu, W.

Zhang, X.

Appl. Opt.

Commun. Pure Appl. Math.

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure Appl. Math. 59(8), 1207–1223 (2006).
[CrossRef]

IEEE Trans. Image Process.

J. M. Bioucas-Dias and M. A. T. Figueiredo, “A new twIst: two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007).
[CrossRef] [PubMed]

IEEE Trans. Inf. Theory

E. J. Candes and T. Tao, “Near-optimal signal recovery from random projections: Universal encoding strategies?” IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006).
[CrossRef]

D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Theory 52(4), 1289–1306 (2006).
[CrossRef]

J. Disp. Technol.

Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Disp. Technol. 6(10), 506–509 (2010).
[CrossRef]

J. Opt. Soc. Am.

J. Opt. Soc. Am. A

Lab Chip

S. Seo, T. W. Su, D. K. Tseng, A. Erlinger, and A. Ozcan, “Lensfree holographic imaging for on-chip cytometry and diagnostics,” Lab Chip 9(6), 777–787 (2009).
[CrossRef] [PubMed]

Mar. Biol.

J. Chae and S. Nishida, “Integumental ultrastructure and colour patterns in the iridescent copepods of the family Sapphirinidae (Copepoda: Poecilostomatoida),” Mar. Biol. 119(2), 205–210 (1994).
[CrossRef]

Nature

D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948).
[CrossRef] [PubMed]

Opt. Eng.

D. P. Kelly, B. M. Hennelly, N. Pandey, T. J. Naughton, and W. T. Rhodes, “Resolution limits in practical digital holographic systems,” Opt. Eng. 48(9), 095801 (2009).
[CrossRef]

Opt. Express

Opt. Lett.

Physica D

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60(1–4), 259–268 (1992).
[CrossRef]

Other

A. C. Kak, and M. Slaney, Principles of Computerized Tomographic Imaging (SIAM, 2001).

J. Hahn, S. Lim, K. Choi, R. Horisaki, D. L. Marks, and D. J. Brady, “Compressive Holographic Microscopy,” in Biomedical Optics, OSA Technical Digest (CD) (Optical Society of America, 2010), paper JMA1.

R. E. Blahurt, Theory of Remote Image Formation (Cambridge University Press, 2005).

D. J. Brady, Optical Imaging and Spectroscopy (Wiley, 2009).

Supplementary Material (5)

» Media 1: MPEG (2668 KB)     
» Media 2: MPEG (1792 KB)     
» Media 3: MPEG (1774 KB)     
» Media 4: MPEG (3884 KB)     
» Media 5: MPEG (3846 KB)     

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1

A schematic of the compressive holographic microscope.

Fig. 2
Fig. 2

Photographs of (a) the holographic microscope and (b) a water container in which water cyclopses are swimming.

Fig. 3
Fig. 3

Raw image of a Gabor hologram.

Fig. 4
Fig. 4

A comparison of reconstructions at chosen axial positions: (a-b) the backpropagation reconstructions and by (c-d) the compressive holographic reconstructions using the data shown in Fig. 3. All the transverse slices of the reconstructions are sequentially shown in (e) obtained by the backpropagation method (Media 1) and in (f) obtained by the compressive holography method (Media 2) for the full range of water container.

Fig. 5
Fig. 5

A comparison of the magnified tails, marked by rectangles in Figs. 4(b) and 4(d), of (a) the backpropagation reconstruction and (b) the compressive holography reconstruction.

Fig. 6
Fig. 6

A 3D visualization of the compressive holography reconstruction (Media 3).

Fig. 7
Fig. 7

Images of (a) the maximum intensity values of the reconstructed density (f) along the propagation directions and (b) a map of the axial positions corresponding to the maximum values in Fig. 7(a). (c) A range colormap represents the HSV space.

Fig. 8
Fig. 8

Videos of (a) raw measurements (Media 4) and (b) range maps of the compressive holography reconstructions (Media 5).

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

I ( x , y ; z F P A ) = | U ( x , y ; z F P A ) | 2 = | U 0 ( x , y ; z F P A ) | 2 + | U s ( x , y ; z F P A ) | 2 + 2 Re { U 0 * ( x , y ; z F P A ) U s ( x , y ; z F P A ) } ,
U s ( x , y ;   z F P A ) = U ( x , y ; z F P A ) U 0 ( x , y ; z F P A ) = d x d y d z U 0 ( x , y ; z ) β ( x , y , z ) h ( x x , y y ; z F P A z ) .
U s ( x , y ; z F P A ) = d x d y d z h F ( x , y ; z ) β ( x , y , z ) h F ( x x , y y ; z F P A z ) = h F ( x , y ; z F P A ) d x d y d z C ( z ) β ( x z z F P A , y z z F P A , z ) h F ( x x , y y ; z F P A z z / z F P A ) = h F ( x , y ; z F P A ) d z C ( z ) F 2 D 1 { F 2 D { β ( x z z F P A , y z z F P A , z ) } H F ( k x , k y ; z F P A z z / z F P A ) } .
C ( z ) = exp [ j k z F P A ( 1 z F P A / z ) ] .
U 0 * ( n 1 Δ , n 2 Δ , z F P A ) U s ( n 1 Δ , n 2 Δ , z F P A ) = h F * ( n 1 Δ , n 2 Δ , z F P A ) U s ( n 1 Δ , n 2 Δ , z F P A ) = 1 N 2 l C ( l Δ z z F P A ) m 1 m 2 n 1 n 2 β ( n 1 Δ , n 2 Δ , l Δ z ) exp ( 2 π j m 1 n 1 1 + m 2 n 2 N ) × H F ( m 1 Δ k , m 2 Δ k , z F P A l Δ z l Δ z / z F P A ) exp ( 2 π j n 1 m 1 + n 2 m 2 N ) ,
U 0 , n 1 n 2 * U s , n 1 n 2 = l C l F 2 D 1 { F 2 D { β n 1 n 2 l } m 1 m 2 H F , m 1 m 2 l } n 1 n 2 .
g = 2 Re { H f } + e + n ,
f ^ = arg min f 1 2 g 2 Re ( H f ) 2 + τ f T V .
f T V = z x y | ( f z ) x , y | ,
Δ x = λ 2 N A ,
Δ z = 2 λ N A 2 ,

Metrics