We introduce depth-filtered digital holography (DFDH) as a method for quantitative tomographic phase imaging of buried layers in multilayer samples. The procedure is based on the acquisition of multiple holograms for different wavelengths. Analyzing the intensity over wavelength pixel wise and using an inverse Fourier transform leads to a depth-profile of the multilayered sample. Applying a windowed Fourier transform with a narrow window, we choose a depth-of interest (DOI) which is used to synthesize filtered interference patterns that just contain information of this limited depth. We use the angular spectrum method to introduce an additional spatial filtering and to reconstruct the corresponding holograms. After a short theoretical framework we show experimental proof-of-principle results for the method.
© 2012 OSA
Digital holography is a powerful technique for the acquisition of the morphology and shape of samples [1–3]. The numerical reconstruction allows for filtering procedures and numerical propagation of the wavefield to any plane of interest where both amplitude and phase information can be calculated. Especially the phase information is widely used in digital holography as it yields impressive depth resolutions of a few nanometer with a relatively simple setup .
On the other hand, the performance of digital holography is limited when the sample consists of multiple semi-transparent layers. The reflections at the different depths will interfere with each other and with the reference beam leading to a superposition of interference patterns on the digital camera. Phase analysis in such cases is a great challenge, especially when the depth of interest (DOI) is a buried layer and not at the surface of the sample. Nevertheless, tomographic measurements of multiple layer samples using digital holography have been performed with a similar approach as is used in low-coherence interferometry and time domain optical coherence tomography (OCT) [4–6]. In this approach the interferences created by the different layers can be separated exploiting the low coherence of broad-band light sources. Therefore multiple scattered photons outside the DOI can be suppressed by the coherence window [7–9]. The depth-resolution in such cases is limited by the coherence length of the light source which is typically a few micrometers. In combination with spatiotemporal and phase-shifting methods this led to first tomographic phase images with nm-resolution [10, 11] behind scattering layers.
Another approach is wavelength-scanning digital holography where holograms are recorded for multiple wavelengths and processed such that tomographic measurements are possible [12–15]. Most of these techniques just analyze the amplitude of the backscattered light. Tomographic phase reconstruction in general is rare and is mainly applied in transmission measurements e.g. through cells to obtain refractive index variations [16–19]. In OCT the phase analysis has mainly been used in topographic measurements [20–22] while tomographic phase imaging has just been reported for differential measurements .
In this paper we introduce an approach for wavelength scanning holography that combines depth sectioning by means of low-coherence interferometry in the spectral domain with digital holography. In this approach a Fourier analysis of multiple digital holograms is followed by a synthesis of interference patterns that are created at specific depths. This approach allows analyzing the depth-filtered phase and amplitude of selected depths in multilayer samples and thus yields tomographic nanometer resolution.
2. Depth-filtered digital holography
In digital holography a digital camera commonly acquires interference patterns utilizing monochromatic light. We extend this method and use a similar approach and data set as in wavelength-scanning digital holography. To do so we record a set of holograms with several distinct wavelengths. The data we obtain is thus the intensity over two lateral dimensions and wavelength. The interference pattern for each wavelength can be analyzed by means of digital holography, but as described above this will only allow for surface analysis. Samples that consist of semitransparent layers cannot be imaged with this approach because the reflections from the different depths overlap on the camera and lead to a superposition of interference patterns. Thus, standard digital holographic reconstruction techniques will fail.
In the following we introduce a technique that filters parasitic reflections and preserves the required information of the depth-of interest (DOI). This information is used to synthesize interference patterns which thus contain information for just the limited depth range. The procedure is depicted schematically in Fig. 1 . As a first step, multiple holograms are recorded for different wavelengths. Then the intensity is analyzed as a function of wavelength for each pixel of the digital holograms. Each path length mismatch between sample and reference wavefield leads to a spectral modulation which can be described by the law of spectral interference expressed as Eq. (1) 
For simplicity we now consider the sample to be just one reflecting surface at distance. Applying the inverse Fourier transform leads to Eq. (2):24].
In this paper we extend this method by using an approach which is usually performed in spectroscopic optical coherence tomography to calculate depth resolved spectroscopic data . The depth of interest (DOI) in the depth profile is separated by means of a windowed Fourier transform, i.e. by multiplying the complex data with a window-function.
This window-function is only non-zero in the depth of interest and zero for other depths. The simplest form of such a window function is a rectangular window. More sophisticated approaches, with elaborate shapes of the window function can show a better performance because the Fourier transformed signal is convoluted with the Fourier Transform of the window function.
In the experiments described in this paper we use a Hann-window to filter a certain depth of the profile. The Hann-window can be expressed as Eq. (3):
It follows directly that the interference pattern at the filtered depth can be expressed as. After the Fourier transform we obtain which describes the filtered spectral modulations and can be used to form a 3D data-cube. The whole procedure leads to synthesized interference patterns, which just contain the information of the DOI. The filtered interference patterns are now reconstructed using the angular spectrum method [26, 27]. To do so, are analyzed at a single wave-number .
The complex object wave-field in the hologram plane which is the plane of the digital sensor placed at distance c = 0 is. The interference pattern we obtain by this procedure is complex. We drop the imaginary part and just use the real part for the reconstructions.
We use an off-axis configuration as introduced by Leith and Upatnieks  meaning that reference and object beams enclose a small angle φ. The spatial frequencies that occur in the spectrum depend on this angle. Similar to the Fourier transform of the spectral modulation, a very strong DC term is present, positioned in the middle of the angular spectrum and surrounded by the real and the virtual image. A part of these spatial frequencies is filtered by a numerical band pass , in our case a Gaussian filter. By appropriate selection of the center frequency of this filter unwanted interferences of light reflected by the sample and the reference can be filtered out. The filtered field can be rewritten as the inverse Fourier transform of its filtered angular spectrum (Eq. (5)):Eq. (6):27].
The reconstructed complex wavefield at any plane at the distance c = d perpendicular to the propagation axis is given by Eq. (7):
Both the depth-filtering procedure and the spatial filtering of the angular spectrum method have basically a comparable approach as just a small part of the spectrum is filtered and processed. A combination of both filtering-techniques leads to a very powerful tool that has the potential to filter object light travelling from a specific depth into a specific direction.
In the next section we present experimental proof-of-principle results for this technique and compare the results with those obtained without depth-filtering.
3. Experimental results
In order to analyze the performance of our technique we employ the sample configuration depicted in Fig. 2 . The metallic island we used in former experiments  is the object to be imaged. An exemplary reconstruction of the pure uncovered island is shown in Fig. 2(a). It has lateral dimensions of 1.0 mm x 1.3 mm and a height of 200 nm, as verified in  using a profilometer-scan.
In the following we place four 150 µm thick cover-glasses in front of the island that are arranged as stairs (Figs. 2(b) and 2(c)) to simulate a multiple layer sample as they create reflections at different depths. The cover-glasses are clamped together, so that they are separated by intermediate air layers. The distance from the last cover-glass to the metallic island is approximately 600 µm. This sample configuration is placed into a Mach-Zehnder like setup as sketched in Fig. 3 . The two lenses have a focal length of f = 150 mm and form a 4f-geometry. The CCD camera consists of 1376 x 1040 pixels with a pixel size of 6.45 µm x 6.45 µm. The laser source is a tuneable diode laser in Littman configuration as we already used in . The whole setup is in off-axis geometry with an angle of approximately 1° between object and reference arm.
We record a digital hologram of the sample and use the angular spectrum method for reconstruction. The wavelength is set to be λ = 820 nm.
For quantitative measurements we acquire a reference phase of a plane mirror in order to remove the system errors introduced to the phase by the optical components [3,26]. The reconstructions of this hologram are depicted in Fig. 4 . The amplitude (Fig. 4(a)) and the wrapped phase (Fig. 4(b)) clearly show that the cover glasses strongly affect the reconstructions. Just the part of the island that is covered by only one glass is reconstructed correctly, as is shown by the line-scan of the unwrapped phase-profile (Fig. 4(c)) across the line marked in Fig. 4(a) with a red dotted line. The 200 nm height of the island is correctly reconstructed. The other parts of the island that are covered with two or three glasses are severely disturbed and cannot be reconstructed.
Figure 5(a) shows the angular spectrum of the interference pattern used. The component in the middle is the DC term which is the centre of point-symmetry for the other parts that correspond to different spatial frequencies formed by different angles between the reference beam and the object beam. The frequency component used for reconstruction of the results in Fig. 4 is marked with the green arrow. As can be seen the cover glasses introduce many parasitic components to the interference pattern, which are caused by multiple reflections.
These multiple reflections have different angles of incidence onto the CCD camera and thus lead to different positions in the angular spectrum. The angular spectrum method allows for selecting and reconstructing each component independently. This is shown in Fig. 5. The parasitic components marked with red arrows are each reconstructed with the angular spectrum method leading to the amplitude reconstructions depicted in the figure. The reconstruction in Fig. 5(a) shows a very smooth amplitude which is connected to the top-reflection on the surface of the cover glasses. The reconstruction in Fig. 5(b) clearly contains information of the sample and is created by multiple reflection of the sample light whereas 5(c) belongs to multiple reflections on the cover glasses. The reconstructions in Fig. 4 have been performed with the spectral components marked with the green arrow.
In the following we use the DFDH-approach to generate synthesized interference patterns that only contain information of the metallic island. We use our tunable laser to record 80 holograms at equidistantly spaced wavelengths between 820.0 nm and 827.9 nm with a wavelength spacing of Δλ = 0.1 nm. The holograms are recorded sequentially, i.e. we check the exact wavelength setting to have accurate values. We do not use a direct wavelength sweeping. The acquisition time for each hologram is set to 50 ms. The recording procedure is comparable to the approach described in . The reconstruction of each individual hologram leads to results comparable to those shown in Figs. 4 and 5.
In order to increase the quality of the depth-profile calculated by the inverse Fourier transform of the modulated spectra we use zero-padding by a factor of ten. Thus we have 800 data points, from which only 400 can be used, due to the symmetry of the Fourier transform. An inverse Fourier transform of one spectral modulation leads to a peak in the depth of the corresponding layer. Repeating the procedure for all pixels, we obtain the depth profile of the whole sample as depicted in Fig. 6 . The obtained depth is proportional to the spectral modulation frequency and thus to the path-length mismatch between object and reference arms . This is in contrast to other techniques, like e.g. standing wave microscopy, where the used modulation frequency is constant .
The cover glasses appear in a depth between 0.8 and 1.3 mm. Due to multiple overlapping reflections and autocorrelations the signal is blurred. The island is placed between 0.1 and 0.4 mm. The performance of DFDH strongly depends on the quality of the depth-profile. The better the depth-separation in the profile the better the result should be. Since we use a narrow spectral bandwidth, which gives a low depth separation, the profile we obtain is of limited quality. In the experiment presented here the depth-separation is approximately 70 µm. As is shown by the angular spectrum in Fig. 5 the sample information is just a small fraction of the interference patterns that is superimposed by the DC part and many parasitic reflections at the different layers.
As we pointed out in  the procedure has inaccuracies as just 80 different wavelengths are used to scan the spectral modulation. Using a broader spectral bandwidth should lead to much better results. Typical depth-resolutions using this approach are around 10 µm and reach 2 µm in sophisticated setups [32, 33]. But the usage of broadband-sources will also increase the influence of dispersion, which will affect the calculation of the depth-profile. This will have an impact on the synthesized interference patterns and consequently on the digital holographic reconstruction.
Besides, numerical errors of the Fourier-transform and the signal processing could also degrade the results. But as we will show, for this proof of principle experiment the set up and signal processing is sufficient to analyse the phase and to obtain nanometer resolution in multilayer samples. As already mentioned we use a Hann-window to filter a specific depth of the profile. We combine the filtered modulations and synthesize interference patterns as described in section 2. Then we use the angular spectrum method for digital holographic reconstruction of the synthesized patterns.
In the following example we apply the depth-filtering for 8 different DOIs which we name window “a” to window “h” as shown in Fig. 7 . The depth-profile in Fig. 6 consists of 400 data-points into axial-direction. The windows we use have a width of 10 data-points except window “d” which has a width of 50 data-points. The different windows in Fig. 7 are marked with colored lines and are spread in different depths. Their center is given in the right column.
The corresponding amplitude reconstruction from the synthesized interference pattern belonging to each named depth is plotted in Fig. 8 . Depending on the DOI, the reconstruction strongly changes. For windows “e” to “h” the amplitude reconstruction is strongly disturbed and the island is not reconstructed. This is different for windows “a”, “b” and “c”. Here the metallic island is clearly visible and the regions where the other reconstructions have their main distortions are clearly reconstructed. Moreover, due to the changing optical path length depending on the number of cover glasses in front of the corresponding part of the island, different parts of the island can be addressed using windows in different depths. The addressed parts are clearly highlighted in the reconstructions. Window “d” ranges over windows “a”, “b” and “c”. Here the whole sample is reconstructed. These results prove that the depth-filtering has advantageous impact on the digital holographic amplitude reconstruction. The whole island is clearly reconstructed using the DFDH-approach.
In the following we analyze the reconstruction of the more important phase information. The results are depicted in Fig. 9 . The columns of Fig. 9 correspond to the amplitude, wrapped and unwrapped phase reconstructions, respectively. While the first row shows unfiltered data the second row corresponds to filtering with the broadest window “d” These results prove that the filtering procedure enables the phase reconstruction of the island.
The 3D view of the unwrapped phases is shown in Fig. 10 . Using DFDH the quality of the phase reconstruction is clearly enhanced. It is evident that the reconstructed phase at the covered parts of the island is not destroyed by the cover-glasses. But due to the change of the optical path-length introduced by the cover-glasses phase jumps occur that lead to ambiguous phase values in the borders of the glasses. This has to be taken into account for quantitative discussions. Additionally as can be seen, at the part of the island covered by one glass in Fig. 10, the phase has a small curvature which is introduced by the cover-glasses. This is especially a significant challenge for tomographic phase analysis of unknown samples. Phase distortions created by layers that are located above the DOI have a strong impact on the reconstruction. In such a case it is possible to use DFDH to analyse the phase distortions iteratively. Starting at the top surface of the sample and scanning the depth with small depth-steps, allows tracing the phase change into the depth dimension. This opens up the possibility to correct the phase in the DOI. But to do so, much better depth-profiles are needed that allow for error-free depth-selection.
In the example shown here the phase distortion can be reconstructed using a two-wavelength interferometry approach that applies reconstruction with a synthetic wavelength . The wavelengths used in this example are λ1 = 821.0 nm and λ2 = 822.0 nm, which leads to a synthetic wavelength of Λ = 675 µm. The results for the unfiltered case are shown in Fig. 11(a) . The 3D view and the unwrapped 2D phase profile of the filtered island are shown in Figs. 11(b) and 11(c). As can be seen the 150 µm path-length change introduced by each cover-glass is correctly reconstructed (Fig. 11(d)). The contour of the island can be clearly seen in the unwrapped phase, while the 200 nm height is naturally not resolvable due to the much bigger synthetic wavelength. This shows that the usage of DFDH allows for correct reconstruction of the path-length changes and that using synthesized interference patterns and synthetic wavelength leads to successful reconstructions.
To reconstruct the height of the island quantitatively a different approach has to be used. We reconstruct the phase across the cover-glasses which means that we circumvent the optical path-length change. We use the line-scans marked in the amplitude reconstruction shown in Fig. 12(a) . The corresponding height-profile is plotted in Fig. 12(b). The green line across the region covered by one glass shows the 200 nm height of the island. The regions where two and three cover glasses prohibit the reconstruction with the common DH-approach show good quantitative results for the height profile using DFDH. The line-scans across the two and three glasses correctly reconstruct the 200 nm height, but are superimposed by a modulation which can be attributed to the limited quality of the depth-profile and thus not ideally windowed and synthesized interference patterns. This modulation degrades the phase resolution which in this example is in the region of around 50 nm. Nevertheless, the reconstruction is good enough to demonstrate that the usage of DFDH enables the tomographic phase analysis of our sample. Improved experimental conditions or advanced processing promise the resolution to be in the nanometer range.
Depth-filtered digital holography is a powerful technique that enables quantitative tomographic phase analysis of covered interfaces. This is possible due to band pass filtering of modulation frequencies of the recorded spectra and spatial frequencies of the interference patterns. The combination of the two filtering techniques enables the analysis of light coming from a certain depth and at a certain angle. The ultimate limits of this technique have still to be explored, yet first results are very promising. One main point will be the improvement of the depth-profiles which will allow for better filtering and will lead to higher quality of the synthesized interference patterns. However, the usage of broadband light-sources will lead to problems with dispersion, which will have to be addressed in detail. The application of our technique through weakly scattering layers should also be possible as the first processing steps are based on a similar approach as frequency domain OCT. However, the performance and limitations of our technique in scattering environment have to be analysed in detail. A priori information of the sample is not needed. The structure of the sample can be obtained iteratively by scanning the DOI through the sample and tracking the phase changes with depth. Depth-filtered digital holography has the potential to enable quantitative measurements with nm resolution in multilayer samples. We see enormous potential for applications in the semiconductor technology for analysis of buried samples and for biomedical applications.
References and links
1. G. Pedrini, P. Fröning, H. J. Tiziani, and F. Mendoza Santoyo, “Shape measurement of microscopic structures using digital holograms,” Opt. Commun. 164(4-6), 257–268 (1999). [CrossRef]
2. F. Charrière, J. Kühn, T. Colomb, F. Montfort, E. Cuche, Y. Emery, K. Weible, P. Marquet, and C. D. Depeursinge, “Characterization of microlenses by digital holographic microscopy,” Appl. Opt. 45(5), 829–835 (2006). [CrossRef] [PubMed]
3. Y. W. Lai, N. Koukourakis, N. C. Gerhardt, M. R. Hofmann, R. Meyer, S. Hamann, M. Ehmann, K. Hackl, E. Darakis, and A. Ludwig, “Integrity of micro-hotplates during high-temperature operation monitored by digital holographic microscopy,” IEEE J. Microelectromechan. Syst. 19, 1–5 (2010).
4. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991). [CrossRef] [PubMed]
7. P. Massatsch, F. Charrière, E. Cuche, P. Marquet, and C. D. Depeursinge, “Time-domain optical coherence tomography with digital holographic microscopy,” Appl. Opt. 44(10), 1806–1812 (2005). [CrossRef] [PubMed]
10. G. Indebetouw and P. Klysubun, “Imaging through scattering media with depth resolution by use of low-coherence gating in spatiotemporal digital holography,” Opt. Lett. 25(4), 212–214 (2000). [CrossRef] [PubMed]
11. S. Tamano, Y. Hayasaki, and N. 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). [CrossRef] [PubMed]
13. L. Yu and M. K. Kim, “Wavelength-scanning digital interference holography for tomographic three-dimensional imaging by use of the angular spectrum method,” Opt. Lett. 30(16), 2092–2094 (2005). [CrossRef] [PubMed]
14. F. Montfort, T. Colomb, F. Charrière, J. Kühn, P. Marquet, E. Cuche, S. Herminjard, and C. D. Depeursinge, “Submicrometer optical tomography by multiple-wavelength digital holographic microscopy,” Appl. Opt. 45(32), 8209–8217 (2006). [CrossRef] [PubMed]
17. C. Joo, T. Akkin, B. Cense, B. H. Park, and J. F. de Boer, “Spectral-domain optical coherence phase microscopy for quantitative phase-contrast imaging,” Opt. Lett. 30(16), 2131–2133 (2005). [CrossRef] [PubMed]
18. B. Rappaz, P. Marquet, E. Cuche, Y. Emery, C. Depeursinge, and P. J. Magistretti, “Measurement of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy,” Opt. Express 13(23), 9361–9373 (2005). [CrossRef] [PubMed]
21. T. Anna, C. Shakher, and D. Singh Mehta, “Simultaneous tomography and topography of silicon integrated circuits using full-field swept source optical coherence tomography,” J. Opt. A, Pure Appl. Opt. 11(4), 045501 (2009). [CrossRef]
22. M. Sticker, C. K. Hitzenberger, R. Leitgeb, and A. F. Fercher, “Quantitative differential phase measurement and imaging in transparent and turbid media by optical coherence tomography,” Opt. Lett. 26(8), 518–520 (2001). [CrossRef] [PubMed]
23. M. Born and E. Wolf, Principles of Optics, 7th ed. (Cambridge University Press, 2003)
24. N. Koukourakis, C. Kasseck, D. Rytz, N. C. Gerhardt, and M. R. Hofmann, “Single-shot holography for depth resolved three dimensional imaging,” Opt. Express 17(23), 21015–21029 (2009). [CrossRef] [PubMed]
25. C. Kasseck, V. Jaedicke, N. C. Gerhardt, H. Welp, and M. R. Hofmann, “Substance identification by depth-resolved spectroscopic pattern reconstruction in frequency domain optical coherence tomography,” Opt. Commun. 283(23), 4816–4822 (2010). [CrossRef]
26. N. Koukourakis, T. Abdelwahab, M. Y. Li, H. Höpfner, Y. W. Lai, E. Darakis, C. Brenner, N. C. Gerhardt, and M. R. Hofmann, “Photorefractive two-wave mixing for image amplification in digital holography,” Opt. Express 19(22), 22004–22023 (2011). [CrossRef] [PubMed]
27. J. W. Goodman, Introduction to Fourier-Optics 2nd Edition, (McGraw Hill, 1996)
28. E. N. Leith, J. Upatnieks, and K. A. Haines, “Microscopy by wavefront reconstruction,” J. Opt. Soc. Am. 55(5), 981–986 (1965). [CrossRef]
30. N. Koukourakis, M. Breede, N. C. Gerhardt, M. R. Hofmann, S. Köber, M. Salvador, and K. Meerholz, “Depth-resolved holographic imaging with variable depth-resolution using spectrally tunable diode laser,” Electron. Lett. 45(1), 46–48 (2009). [CrossRef]
32. M. Wojtkowski, V. J. Srinivasan, T. H. Ko, J. G. Fujimoto, A. Kowalczyk, and J. S. Duker, “Ultrahigh-resolution, high-speed, Fourier domain optical coherence tomography and methods for dispersion compensation,” Opt. Express 12(11), 2404–2422 (2004). [CrossRef] [PubMed]
33. S. H. Yun, G. J. Tearney, B. E. Bouma, B. H. Park, and J. F. de Boer, “High-speed spectral-domain optical coherence tomography at 1.3 mum wavelength,” Opt. Express 11(26), 3598–3604 (2003). [CrossRef] [PubMed]