Off-axis holography [tutorial]


Off-axis holography is a popular technique to numerically reconstruct a hologram. It allows to retrieve amplitude and phase of a wavefront by measuring only one image with a digital camera. It relies on an intereferometric setup with a non-zero angle with the reference beam and the signal beam and requires a numerical spatial frequency filtering. 


This technique was initially presented in [1]. The schematic of the experimental setup is shown in Fig. 1. The idea is to make interfere the unknown "signal" wavefront of complex field with a plane wave of amplitude  with a given angle  between the direction of propagation of the two beams. We assume that  and  are in the same state of polarization.


Figure 1. Schematic of the setup.


The intensity measured by the camera reads:



The three terms correspond to signals centered around the spatial frequencies 0,  and . If the spatial frequency  is larger than the maximal spatial frequency in the signal to reconstruct, the Fourier transform of the measured intensity pattern shows three components corresponding orders 0, +1 and -1 that do not overlap. We note that, since  is constant, the order -1 is directly proportional to the complex field  we want to retrieve but simply shifted by the spatial frequency . The idea is then to numerically select the order -1 in the Fourier domain and center the signal around the zero frequency. The inverse Fourier transform then directly gives access to the complex image . The steps of this process are illustrated in Fig. 2 and Fig. 3.

The following images were obtained with the Matlab code attached to this article and used to simulate an off-axis holography measurement (the file is available at the bottom of this page, the code includes comments). For the sake of simplicity, in this example, we try to recover an image with a flat amplitude and with a random phase pattern, similar to the phase of speckle pattern [Fig 2.a.]. After interference with the tilted plane wave [Fig 2.b.], we calculate the intensity pattern as it would be measured on the camera [Fig 2.c.].



Figure 2. Phase wavefront to measure (a), tilted plane wave reference (b) and resulting interference pattern.


 We then calculate the Fourier transform of this intensity pattern [Fig 3.a.]. We note that the angle of the tilted plane wave has to be high enough so that order 0 and -1 do not overlap. We then select only the -1 order [Fig 3.b.] and shift it around the zero frequency [Fig 3.b.].


Figure 3. Numerical spatial frequency filtering. (a) Complete Fourier transform of the intensity pattern. (b) Fourier transform filtered to keep only the -1 order. (c) Signal shifted around the zero spatial frequency.


The reconstructed complex field is obtained by calculating the inverse Fourier transform of the last step. We compare in Fig 4. the phase of the in initial image and the phase of the recovered one.


 Figure 4. Comparison between the actual phase pattern (a) and the reconstructed one (b).

This technique allows to perform a one shot measurement of the complex field for one state of polarization at the cost of a loss of resolution due to the need to resolve the fringes inside one speckle grain. Such a technique can be extended to obtain a full image for both polarizations by using two reference beams with orthogonal polarizations with tilt angles and  in two different directions.


[1] E. Cuche et al., Opt. Express, 39, (2000)


Submit to FacebookSubmit to Google BookmarksSubmit to TwitterSubmit to LinkedIn
Download this file (simOffAxis.m)simOffAxis.m[Simple Matlab code to simulate an off-axis holography measurement]3 kB
Download this file (simOffAxis.p)[Simple Python code to simulate an off-axis holography measurement]3 kB


#3 Davood Khodadad 2014-09-23 11:43
Thank you very much, it is very useful to understand whatever is going on in Holography
#2 Yongzhao Du 2014-04-23 23:00
Good,Thank you so much!
#1 Lei Gong 2013-10-12 09:06
Nice!Thanks for your sharing!

Additional information