Creating an extended focus image of a tilted object in Fourier digital holography

We present a new method to numerically reconstruct images on a tilted plane by digital holography in Fourier configuration. The proposed technique is based on a quadratic deformation of spatial coordinates of the digital hologram. By this approach we demonstrate that it is possible to recover the extended focus image (EFI) of a tilted object in a single reconstruction step from the deformed hologram. ©2009 Optical Society of America OCIS codes: (090.1995) Holography: Digital holography; (180.3170) Microscopy: Interference microscopy; (110.0180) Imaging systems: Microscopy; (180.6900) Microscopy: Threedimensional microscopy; (100.3010) Image processing: Image reconstruction techniques. References and links 1. W. E. Ortyn, D. J. Perry, V. Venkatachalam, L. Liang, B. E. Hall, K. Frost, and D. A. Basiji, “Extended depth of field imaging for high speed cell analysis,” Cytometry A 71A(4), 215–231 (2007), http://dx.doi.org/10.1002/cyto.a.20370. 2. E. R. Dowski, Jr., and W. T. Cathey, “Extended depth of field through wavefront coding,” Appl. Opt. 34(11), 1859–1866 (1995). 3. Z. Zalevsky, and S. Ben-Yaish, “Extended depth of focus imaging with birefringent plate,” Opt. Express 15(12), 7202–7210 (2007), http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-12-7202. 4. G. Mikula, Z. Jaroszewicz, A. Kolodziejczyk, K. Petelczyc, and M. Sypek, “Imaging with extended focal depth by means of lenses with radial and angular modulation,” Opt. Express 15(15), 9184–9193 (2007), http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-15-9184. 5. H. Zhao, Q. Li, and H. Feng, “Improved logarithmic phase mask to extend the depth of field of an incoherent imaging system,” Opt. Lett. 33(11), 1171–1173 (2008), http://www.opticsinfobase.org/abstract.cfm?URI=ol-3311-1171. 6. S. Bagheri, and B. Javidi, “Extension of depth of field using amplitude and phase modulation of the pupil function,” Opt. Lett. 33(7), 757–759 (2008), http://www.opticsinfobase.org/abstract.cfm?URI=ol-33-7-757. 7. M. Kanka, R. Riesenberg, and H. J. Kreuzer, “Reconstruction of high-resolution holographic microscopic images,” Opt. Lett. 34(8), 1162–1164 (2009). 8. D. Wang, J. Zhao, F. Zhang, G. Pedrini, and W. Osten, “High-fidelity numerical realization of multiple-step Fresnel propagation for the reconstruction of digital holograms,” Appl. Opt. 47(19), D12–D20 (2008). 9. Y. S. Hwang, S.-H. Hong, and B. Javidi, “Free View 3-D Visualization of Occluded Objects by using Computational Synthetic Aperture Integral Imaging,” J. Disp. Technol. 3(1), 64–70 (2007). 10. Y. Takaki, and H. Ohzu, “Hybrid holographic microscopy: visualization of three-dimensional object information by use of viewing angles,” Appl. Opt. 39(29), 5302–5308 (2000). 11. E. Malkiel, J. N. Abras, and J. Katz, “Automated scanning and measurements of particle distributions within a holographic reconstructed volume,” Meas. Sci. Technol. 15(4), 601–612 (2004). 12. P. Ferraro, G. Coppola, S. De Nicola, A. Finizio, and G. Pierattini, “Digital holographic microscope with automatic focus tracking by detecting sample displacement in real time,” Opt. Lett. 28(14), 1257–1259 (2003). 13. M. Antkowiak, N. Callens, C. Yourassowsky, and F. Dubois, “Extended focused imaging of a microparticle field with digital holographic microscopy,” Opt. Lett. 33(14), 1626–1628 (2008). 14. F. Dubois, C. Schockaert, N. Callens, and C. Yourassowsky, “Focus plane detection criteria in digital holography microscopy by amplitude analysis,” Opt. Express 14(13), 5895–5908 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-13-5895. 15. M. Liebling, and M. Unser, “Autofocus for digital Fresnel holograms by use of a Fresnelet sparsity criterion,” J. Opt. Soc. Am. A 21(12), 2424–2430 (2004). 16. R. J. Pieper, and A. Korpel, “Image processing for extended depth of field,” Appl. Opt. 22(10), 1449–1453 (1983). 17. P. Ferraro, S. Grilli, D. Alfieri, S. De Nicola, A. Finizio, G. Pierattini, B. Javidi, G. Coppola, and V. Striano, “Extended focused image in microscopy by digital Holography,” Opt. Express 13(18), 6738–6749 (2005), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-13-18-6738. #114254 $15.00 USD Received 14 Jul 2009; revised 9 Sep 2009; accepted 11 Sep 2009; published 23 Oct 2009 (C) 2009 OSA 26 October 2009 / Vol. 17, No. 22 / OPTICS EXPRESS 20546 18. S. De Nicola, A. Finizio, G. Pierattini, P. Ferraro, and D. Alfieri, “Angular spectrum method with correction of anamorphism for numerical reconstruction of digital holograms on tilted planes,” Opt. Express 13(24), 9935– 9940 (2005). 19. S. J. Jeong, and C. K. Hong, “Pixel-size-maintained image reconstruction of digital holograms on arbitrarily tilted planes by the angular spectrum method,” Appl. Opt. 47(16), 3064–3071 (2008), http://www.opticsinfobase.org/abstract.cfm?URI=ao-47-16-3064. 20. K. Matsushima, “Formulation of the rotational transformation of wave fields and their application to digital holography,” Appl. Opt. 47(19), D110–D116 (2008). 21. C. P. McElhinney, B. M. Hennelly, and T. J. Naughton, “Extended focused imaging for digital holograms of macroscopic three-dimensional objects,” Appl. Opt. 47(19), D71–D79 (2008). 22. D. Leseberg and C. Frère, Detlef Leseberg and Christian Frère, “Computer-generated holograms of 3-D objects composed of tilted planar segments,” Appl. Opt. 27(14), 3020–3024 (1988), http://www.opticsinfobase.org/ao/abstract.cfm?URI=ao-27-14-3020. 23. P. Ferraro, M. Paturzo, P. Memmolo and A. Finizio, “Controlling depth-of-focus in 3D image reconstructions by flexible and adaptive deformation of digital holograms,” Opt. Lett. 34, 2787-2789 (2009).


Introduction
One of main problems in optical imaging is the limited depth of focus.All optical systems are affected by this limitation that hinder to get in focus in a single image plane objects that are located at different distances, but falling inside the same field of view.This is especially true in optical systems like microscopes where the depth of field is squeezed due to the request of an high magnification.This problem is an important issue in biological applications when a fast analysis for dynamic processes is necessary [1].In fact, in classical optical microscopes the Extended Focus Image (EFI) is built-up by the optical longitudinal mechanical scanning and by the acquisition of a sequence of images during the scanning, called "image stack".To create an EFI, in each image of the "stack" the in-focus parts are selected and then stitched together in a single image (i.e.EFI).However, the long acquisition time required for digitizing many images prevents it to be used in case of dynamic events.To overcome this problem, some approaches have been developed for incoherent imaging system.One approach is based on the so called "wavefront coding" using an opportune phase plate that introduces a wavefront distortion giving the possibility to have a EFI in a single image acquisition [2].Recently, the problem to extend the depth of field has been successfully addressed for incoherent imaging system by using a birefringent plate [3], by means of lenses with radial and angular modulation [4], by logarithmic phase mask [5] or by amplitude and phase modulation of the pupil function [6].
The need of an EFI image is also important in 3D optical imaging systems both for incoherent and coherent light sources [7][8][9][10][11][12][13][14][15][16].Of course holographic microscopes allows 3D imaging capability but suffer at the same time of a limited depth-of-focus.Recently various solutions have been proposed to obtain EFI in DH systems [13,[16][17][18][19][20][21].Some of these approaches use the angular spectrum transformation, others exploit the DH capability to reconstruct image planes at different distances, or image processing during the reconstruction process.The diffraction of a tilted planar object is also studied in ref. 22, where the authors generate a computer-generated hologram of a 3D object composed of planar segment tilted in respect to the optical axis.
Here we show that by means of an opportune coordinates transformation of an hologram recorded in Fourier configuration, it is possible to obtain an EFI image of a tilted object.In our recent work we have shown that in Fresnel holography a linear transformation allows to extend the depth of focus [23].However, such approach does not work properly for Fouriertype holograms.We extend here the theory about the hologram deformation to Fourier holograms, showing that only a quadratic geometrical deformation allows to obtain a tilted object with all points in focus by a single reconstruction step.Moreover, we report the results of some specific experiments with the aim to demonstrate the validity of the method.

Deformation of Fourier-type holograms: theoretical model
In our previous paper we demonstrated how it is possible to change the in-focus distance by stretching holograms acquired in Fresnel configuration [23].
For holograms recorded in Fourier configuration, the simple linear stretching gives a different result as can be demonstrated by simple calculations reported below.In fact, in the Fourier configuration the numerical reconstruction of the hologram ( ) , f ξ η is obtained simply calculating its spatial FT ( ) ) Because of the Fourier transform properties we have that ( ) where ( ) ˆ, h x y corresponds to the new numerical reconstruction and it results to be equal to the reconstruction of the initial hologram with scaled dimensions and a scaled intensity.
Therefore, there is no change in reconstruction distance.Consequently by application of a linear transformation to a Fourier hologram no changes can be observed about image focusing properties during reconstruction.
Nevertheless, if, instead of a linear transformation described by Eq. ( 1), we apply a quadratic deformation, it results that the in-focus distance changes, as will be shown later in this paper, with a linear law given by where l' is the horizontal coordinate in the reconstruction plane, α is the deformation parameter and L is the reconstruction distance of the recorded (not deformed) hologram.Therefore, by this kind of hologram deformations we can put correctly in-focus the numerical reconstruction of a tilted object recovering its EFI.It important to note that we are considering the possibility to deform the hologram only along the direction of the tilt, that means in our case the x-axis.
To demonstrate how the change of the in-focus distance depends on the quadratic deformation we need to remind some Fourier Transform properties for a composite function If we consider the following composite function where ( ) G l is the Fourier Transform of ( ) g y , we have that the FT of ( ) .For a quadratic coordinate transformation of the hologram g(x) we have that ( ) Then, it results that ( ) .Therefore the final reconstruction, i.e. the reconstruction of the deformed hologram, can be interpreted as the reconstruction of the initial hologram at distance L, and further propagated to a distance d' that depends linearly on the coordinate l' (x-axis in the reconstruction plane) and it is a function of the deformation parameter α and of L, that is the reconstruction distance of the not-deformed hologram (see Fig. 1).
Summarizing, we can claim that, if we apply a quadratic deformation to the hologram of a tilted object along the tilting direction, we obtain an all in-focus reconstructed image of the object.In fact the reconstruction distance changes pixel by pixel along the tilting direction l' according to Eq. ( 3) as will be experimentally demonstrated in the next section.
The quadratic deformation is applied to the original recorded hologram through the coordinates transformation [ ]  where the operator is .

Experimental validation
For the experimental validation, we made an experiment with a Mach-Zehnder interferometer in Fourier configuration as depicted in Fig. 1.The object was an USAF target positioned in a tilted way in respect to the laser light illumination direction.The Fourier holographic configuration is such that the curvature of the reference beam matches that of the light scattered by the left side of the object (where is the number "0"), i.e. the part of the target closer to the CCD (at a distance L).Consequently this region results to be in-focus in the numerical reconstruction, that is the hologram Fourier Transform.On contrary the right side of the target (where is the number "1") results to be completely out-of-focus in the numerical reconstruction.In Fig. 2 we show the numerical reconstruction of the acquired hologram for an object tilted with an angle φ = 55°.It is clear that the left part of the object, where there is the number "0", results to be in-focus, while the right part, with the number "1", is out of focus and that the focusing gradually worsen going from left to right side.Finally we performed another experiment for a very high tilted angle of φ = 75°.In Figs.4(a) and 4(b) are shown the numerical reconstructions of the original and deformed holograms of the same object tilted with the angle of 75°.In the reconstructions of the deformed objects, all details of the target result to be in-focus.We encircled number the number "2" in all the reconstructions to better guide the reader to visualize the comparison between the corrected and uncorrected images.In all the images some replicas of the target elements appear.They are due to the multiple reflections of the laser light on the target surfaces, in fact they are more evident in the image of the object tilted with an angle φ = 75°.Moreover, it is clear that the replicas are not due to the applied deformation as they are visible in the reconstructions of the not-deformed holograms (see Fig. 2 and Fig. 4a).

Conclusion
We have demonstrated that by an opportune deformation (numerical transformation) of a digital hologram it is possible to obtain an all in-focus image (i.e.EFI) for tilted objects in case of the Fourier-type holograms.The experimental results are in agreement with the theoretical model.It has been demonstrate that the EFI can be recovered for very high tilted angles.The approach is completely new in respect to the previous methods developed to obtain EFI both for incoherent and coherent imaging systems.In addition, this new approach can open the possibility to make more complex numerical manipulation of digital holograms to correct aberrations.The method is simple because it require just a single recorded hologram dispensing from the need to records multiple images and avoiding mechanical scanning.
The method is easily applied to digital holograms where the coordinate deformation is performed numerically.Moreover for holograms recorded on a specific material (physically recorded on a photosensitive media) it is possible to foresight to apply an adaptive mechanical deformation to extend the depth of focus even in that case.

(
$15.00 USD Received 14 Jul 2009; revised 9 Sep 2009; accepted 11 Sep 2009; published 23 Oct 2009 change the coordinates from the spatial frequency domain to the space domain, then in Eq. (

Fig. 2 .
Fig. 2. Numerical reconstruction of the holograms for an object tilted with an angle of 55 degrees as acquired by the CCD (Movie_1).

Fig. 3 .
Fig. 3. Numerical reconstruction of the holograms for an object tilted with an angle of 55 degrees after the quadratic deformation.If we apply the quadratic deformation before performing the reconstruction (with an opportune choice of the deformation parameter α), we obtain the image shown in Fig.3.In this case the α value is 2.1* 10 −5 .Moreover, we show in the Movie_1 how the focus changes varying the parameter α between 0 and the value 2.1* 10 −5 with a step of 2.1* 10 −6 /frame.

Fig. 4 .
Fig. 4. Numerical reconstruction of the holograms for an object tilted with an angle of 75 degrees as acquired by the CCD (a) and after the quadratic deformation (b).