## Abstract

We present a new method for numerically reconstructing digital holograms on tilted planes. The method is based on the angular spectrum of plane waves. Fast Fourier transform algorithm is used twice and coordinate rotation in the Fourier domain enables to reconstruct the object field on the tilted planes. Correction of the anamorphism resulting from the coordinate transformation is performed by suitable interpolation of the spectral data. Experimental results are presented to demonstrate the method for a single-axis rotation. The algorithm is especially useful for tomographic image reconstruction.

©2005 Optical Society of America

## 1. Introduction

The basic advantage of digital holography is that it makes possible to extract quantitatively the three-dimensional (3D) information of the object from the numerical reconstruction of a single digitally recorded hologram [1–3]. This open the field to a variety applications, such as 3D microscopic investigations of biological specimens where wavelength-scanning digital interference holography [4] is used to reconstruct the 3D volume from a set of scanning tomographic images, hybrid holographic microscopy [5], three-dimensional microscopy by optical scanning holography [6], phase shifting digital holography, particle holography, just to quote few of them. Accurate reconstruction of the 3D image volume requires reconstructing the object field at several axial distances from the recording plane and this is a demanding task from the numerical point of view since, for example, if we use the Fresnel transform method (FTM) [2,3] for hologram reconstructions, reconstructed images at different distances need to be resized to account of the dependence of the reconstruction pixel on the distance [7].

Although various methods have been devised for controlling the lateral resolution in the reconstruction plane [7,8], the task is still demanding when many reconstructions along the axial distance have to be stuck together and the axial resolution differs from the lateral one. However, if one is interested in inspecting characteristics of the object on a plane or a set of planes tilted with respect to the recording hologram one, such as in tomographic applications, it is more efficient to develop a method capable of reconstructing the hologram at arbitrarily tilted planes. Basically, this means simulating light propagation through calculation of the diffraction between arbitrarily oriented planes. Leseberg and Frére were the first [9] that addressed the problem of using the Fresnel approximation to find the diffraction pattern of a tilted plane. Later on, a general-purpose numerical method for analyzing optical systems by use of full scalar diffraction theory was proposed [10]. FTM for numerical reconstruction of digital holograms with changing viewing angles was described [11]. In this letter we present a different method for digital holography on tilted plane (DHT) that makes use of the angular spectrum of plane waves [12] for fast calculation of optical diffraction. Fast Fourier transformation is used twice and coordinate rotation of the spectrum enables to reconstruct the hologram on the tilted plane. Interpolation of the spectral data is shown to be effective for correcting the anhamorphism of the reconstructed image.

## 2. Experimental approach

The experimental set-up for DHT shown in Fig. 1 is a Mach-Zender interferometer designed for reflection imaging of microscopic objects. A linearly polarised collimated beam from a diode-pumped, frequency-doubled Nd:YVO_{4} laser with *λ*=*532 nm* is divided by the combination of half-wave plate *λ*/2 and polarizing beam splitter into two beams.

In each arm, spatial filters and beam expanders are introduced to produce plane waves. The object beam is reflected by the test object consisting of the lithographic patterned logo “MEMS” covering an area of (644 × 1568) μm^{2} on a reflective silicon substrate. As shown in Fig. 1, the object is tilted with respect to the hologram plane (*ξ* - *η*). The interfering object beam and the reference beam are tilted at small angle with respect to each other to produce off-axis holograms. The digital holograms are recorded by a standard black and white CCD camera and digitized into a square array *NxN* =*(1024 x 1024) pixel ^{2}*, Δ

*ξ*×Δ

*η*=

*(6.7 × 6.7)μm*pixel width.

^{2}The microscope produces a magnified image of the object and the hologram plane is located between the MO and the image plane *I* at a distance *d* from the recording hologram plane *H*. In digital holographic microscopy we can consider the object wave emerging from the magnified image and not from the object itself [13].

The geometry of the reconstruction algorithm is illustrated in Fig. 2. The image plane *I* has coordinates (*x*̂ - *y*̂) and its origin located along the optical *z* axis at a distance *d* from the hologram plane (*ξ* - *η*). The plane *I* is tilted at angle *θ* with respect to the plane (*x* - *y*), that is (*x* - *y*) is parallel to the hologram plane and share the same origin of the tilted plane *I*. Reconstructing by digital holography the complex wave distribution *o*(*x*̂,*y*̂) on the tilted image plane means simulating light propagation from the hologram plane to the image plane through calculation of the diffraction of the object wave distribution *o*(*x*,*y*) on the plane *H* (at *z*=0) to the tilted plane.

The angular spectrum-based algorithm (ASA), for reconstructing the wave field *o*(*x*̂, *y*̂), consists basically of two steps. In the first stage we reconstruct the wave distribution on the intermediate plane (*x*̂ - *y*̂) at distance *d* and calculate the spectrum *O*̂(*d*;*u*,*v*)where *u* = *ξ*/*λd* and *v* = *η*/*λd* denotes the Fourier frequencies in the intermediate plane. A convenient way of performing the spectrum calculation is through the application of the so called convolution method (CM) for reconstructing the wave field at the translational distance *d* [1]. The formulation makes use of the convolution theorem to calculate the spectrum as follows:

where *O*(*u*,*v*) = *ℑ*{*o*(*x*,*y*)} is the spectrum of the object wave distribution *o*(*x*,*y*) at *z* = 0 is defined by the Fourier transform

In CM the field distribution at distance *d* is derived by taking the inverse Fourier transform of the propagated angular spectrum (1), but in our case we have to take into account of the tilting of the reconstruction plane (*x*̂ - *y*̂) with respect to the intermediate plane. Standard transformation matrix is used to rotate the wave vector coordinates **k** = 2*π*[*u*,*v*,*w*] where *w* = (*λ*
^{-2} - *u*
^{2} -*v*
^{2})^{1/2} , on the *y* axis with the angle of *θ* to determine the spatial frequencies *u*̂ and *v*̂ associated to the components of the rotated wave vector components **k** = 2*π*[*u*̂,*v*̂,*w*̂], namely *u*̂ = *u*cos*θ*-*w* sin *θ* and *v*̂ = *v*. According to the spatial frequencies transformation, the rotated angular spectrum (RAS) can be written as *O*̂(*d*;*u*̂,*v*̂) = *O*̂(*d*; *u*̂ cos θ + *w*̂ sin *θ*, *v*̂).

In the second stage of ASA, the rotated spectrum is inverse Fourier transformed to calculate the reconstructed wave field *o*(*x*̂, *y*̂) on the tilted plane, namely

It should be remarked that reconstructing the field according to Eq. (3) is valid if we stay within the paraxial approximation, as it is applies to our experimental conditions. However, Eq. (3) can be generalized to include frequency dependent terms of the Jacobian associated to the RAS. However, the inclusion of these terms does not affect the reconstruction, since they are small in the paraxial case. Equation (3) is a good approximation for the numerical calculation of the wave field on tilted plane from a digitally recorded hologram as it will be shown in the discussion of the experimental results.

## 3. Experimental results

Figure 3(a) shows the hologram of the reflective target recorded in off-axis configuration. The plane reference wave interferes with the object waves at small angle (≤ 0.5), as required by the sampling theorem. The recording distance *d* is set to be 265 mm Figs. 3(b), 3(c) and 3(d) show a sequence of amplitude images reconstructed at different *z* planes located at distances 240 mm, 265 mm and 290 mm, respectively, from the hologram plane. The real and unsharp virtual images of the object are separated because of the off-axis geometry. The real amplitude reconstructions, calculated by the FTM with the zero-order term filtered out, can be seen in the bottom-left part of the reconstructed area.

Note the three reconstructions by FTM cannot be compared each other, their size being different in each case, since the reconstruction pixel given by *RP* = *λd*/*N*Δ*ξ* depends on the reconstruction distance *d*.

Looking at the reconstructions shown in the Fig. 3, different features of the amplitude images are more or less highlighted or darkened as the reconstruction distance is changed. At 240 mm the left border of the letter “S” is in- focus; at 265 mm the “E” is in-focus in the middle of the logo script and at 290 mm the “M” letter on the left is in-focus. The result of application of ASA is shown in Fig. 3(e). The target is reconstructed on a plane tilted at angle *θ* = 45°. Now the features of the patterned letters appear to be highlighted uniformly but the reconstructed image appears clearly shrunk in the horizontal *x* direction, since the pixel size Δ*x* × Δ*y* is different along the *x* and *y* directions of the tilted plane.

In fact, the spatial frequency resolution Δ*u*̂ in the horizontal *x* direction is related to the pixel size Δ*x* by the relation Δ*u*̂ = 1/*N*Δ*x* and, neglecting the *w* dependent term in the transformation law of *u*̂, Δ*u*̂ is given to a first approximation by Δ*u*̂ = Δ*u* cos *θ*. Since Δ*u* = Δ*ξ*/*λd* we readly obtain Δ*x* = *λd*/(*N*Δ*ξ* cos *θ*). However, the spatial frequency along the *y* direction is unchanged, i.e. Δ*v*̂ = Δ*v* = Δ*η*/*λd*, and the pixel size along the *y* direction is Δ*y* = *λd*/*N*Δ*η* = Δ*x*cos *θ*, where in the last step we have assumed Δ*ξ* = Δ*η*.

Clearly, the greater the tilt angle, the larger the size of the reconstruction pixel in the tilted plane and the anamorphism of the reconstructed image.

Let us discuss a simple method for compensating this anamorphism which works straightforwardly in our case. In ASA we are in general confronted with the numerical problem of applying FFT algorithm to the RAS (see. Eq. (3)). The spacing of the sampling points of the RAS is not equal, due to the transformation law of the spatial frequencies and a method of interpolation must be unavoidable introduced into the sampled values of the RAS, given that FFT algorithm works properly for an equidistant sampled grid.

In order to obtain an evenly spaced grid of *N*×*N* interpolated values of the RAS, we use Δ*v*̂ as a frequency sampling interval for both the *u*̂ and *v*̂ ranges of the RAS. Since the *u*̂ range is less than *v*̂ range, a zero is introduced into the spectral array when the corresponding sampled *u*̂ frequency does not fit into the *u*̂ range. Taking the inverse Fourier transform we obtain a reconstructed image with a square pixel of size Δ*x*×Δ*y* =*λd*/*N*Δ*η*×*λd*/*N*Δ*η* with corrected anamorphism as shown in Fig. 3(f). This correction procedure works in the Fourier space essentially the same the zero-padding method [14] does in the hologram space variables, to recover image resolution by FTM.

## 4. Conclusion

In conclusion, a two stage reconstruction algorithm for DHT has been proposed and experimentally demonstrated. The method is implemented by using FFT twice and coordinate rotation of the angular spectrum in the Fourier domain. The correction of the anamorphism resulting from the coordinate transformation is performed by suitable interpolation of the spectral data. It is worth to point out that although the described algorithm consider single-axis rotation of the reference coordinate it can be easily extended to more general case. The proposed method can be fruitfully exploited as a computationally fast method in tomographic applications by reconstructing the object field along a set of selected planes or for observing special features of an assembly of objects along selected directions, as in particle holography.

## References and links

**
1
. **
L.
Yaroslavsky
and
M.
Eden
, Fundamentals of Digital Optics, (
BirkhŠuser, Boston
,
1996
).

**
2
. **
O.
Schnars
and
W.
Juptner
, “
Direct recording of holograms by a CCD target and numerical reconstruction
,”
Appl. Opt.
**
33
**
,
179
–
181
(
1994
). [CrossRef] [PubMed]

**
3
. **
T. M.
Kreis
, “
Frequency analysis of digital holography
,”
Opt. Eng.
**
41
**
,
771
–
778
(
2002
). [CrossRef]

**
4
. **
M.
Kim
, “
Tomographic three-dimensional imaging of a biological specimen using wavelength-scanning digital interference holography
,”
Opt. Express
**
7
**
,
305
–
310
(
2000
).
http://www.opticsexpress.org/abstract.cfm?URI=OPEX-7-9-305
[CrossRef] [PubMed]

**
5
. **
Y.
Takaki
and
H.
Ohzu
, “
Hybrid holographic microscopy: visualization of three-dimensional object information by use of viewing angles
,”
Appl. Opt.
**
39
**
,
5302
–
5308
(
2000
). [CrossRef]

**
6
. **
T.C.
Poon
,
K.B.
Doh
, and
B.W.
Schilling
, “
Three-dimensional microscopy by optical scanning holography
,”
Opt. Eng.
**
34
**
,
1338
–
44
(
1995
). [CrossRef]

**
7
. **
P.
Ferraro
,
S. De
Nicola
,
G.
Coppola
,
A.
Finizio
,
D.
Alfieri
, and
G.
Pierattini
, “
Controlling image size as a function of distance and wavelength in Fresnel-transform reconstruction of digital hologram
,”
Opt. Lett.
**
29
**
,
854
–
856
(
2004
). [CrossRef] [PubMed]

**
8
. **
T.
Zhang
and
I.
Yamaguchi
, “
Three-dimensional microscopy with phase-shifting digital holography
”
Opt. Lett.
**
23
**
,
1221
–
3
(
1998
). [CrossRef]

**
9
. **
D.
Leserberg
and
C.
Frére
, “
Computer generated holograms of 3-D objects composed of tilted planar segments
,”
Appl. Opt.
**
27
**
,
3020
–
2024
(
1988
) [CrossRef]

**
10
. **
N.
Delen
and
B.
Hooker
, “
Free-space beam propagation between arbitrarily oriented planes based on full diffraction theory: a fast Fourier approach
,”
J. Opt. Soc. A.
**
15
**
,
857
–
867
(
1998
). [CrossRef]

**
11
. **
L.
Yu
,
Y.
An
, and
L.
Cai
, “
Numerical reconstruction of digital holograms with variable viewing angles
,”
Opt. Express
**
10
**
,
1250
–
1257
(
2002
).
http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-22-1250
[PubMed]

**
12
. **
K.
Matsushima
,
H.
Schimmel
, and
F.
Wyrowski
,“
Fast calculation method for optical diffraction on tilted planes by use of the angular spectrum of plane waves
,”
J. Opt. Soc. A.
**
20
**
,
1755
–
1762
(
2003
). [CrossRef]

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

**
14
. **
P.
Ferraro
,
S. De
Nicola
,
A.
Finizio
,
G.
Pierattini
, and
G.
Coppola
, “
Recovering image resolution in reconstructing digital off-axis holograms by Fresnel-transform method
,”
Appl. Phys. Lett.
**
85
**
,
2709
–
2711
(
2004
). [CrossRef]