Table of Contents Author Guidelines Submit a Manuscript
Volume 2017 (2017), Article ID 7293905, 7 pages
Research Article

Manipulations of Wavefront Propagation: Useful Methods and Applications for Interferometric Measurements and Scanning

1Department of Applied Physics/Electro-Optics Engineering, Faculty of Engineering, Lev Academic Center, 9116001 Jerusalem, Israel
2Shlonsky St. 31, 69400 Tel Aviv, Israel

Correspondence should be addressed to Avi Karsenty

Received 5 April 2017; Revised 29 June 2017; Accepted 20 July 2017; Published 21 August 2017

Academic Editor: Richard Arinero

Copyright © 2017 Avi Karsenty et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Phase measurements obtained by high-coherence interferometry are restricted by the 2π ambiguity, to height differences smaller than λ/2. A further restriction in most interferometric systems is for focusing the system on the measured object. We present two methods that overcome these restrictions. In the first method, different segments of a measured wavefront are digitally propagated and focused locally after measurement. The divergent distances, by which the diverse segments of the wavefront are propagated in order to achieve a focused image, provide enough information so as to resolve the 2π ambiguity. The second method employs an interferogram obtained by a spectrum constituting a small number of wavelengths. The magnitude of the interferogram’s modulations is utilized to resolve the 2π ambiguity. Such methods of wavefront propagation enable several applications such as focusing and resolving the 2π ambiguity, as described in the article.

1. Introduction

Highly accurate phase measurements are very important in applications of surface measurements [13]. Nevertheless, phase measurements obtained by high-coherence interferometry are restricted by the 2π ambiguity to height differences smaller than λ/2. In order to increase the dynamic range of 3D measurements for continuous objects, known unwrapping algorithms can be used [4, 5]. However, for step-like objects, several different measurements, each consisting of a different wavelength, must be applied [6]. By comparing the phases obtained by these measurements at each point of the object, the 2π ambiguity is resolved and a considerably larger dynamic range can be measured. In this article, we suggest to overcome this restriction regarding a step-like object by digitally propagating the measured wavefront reflected by the object and refocusing it. By focusing each segment of the wavefront locally and by retaining the distances by which each part of the wavefront was propagated and focused, the information as to resolve the 2π ambiguity is provided. Moreover, this procedure obviates the need for autofocusing mechanism simply by taking a measurement and propagating it to the image plane.

2. Model

2.1. Gibbs Phenomenon

The complex amplitude of a wavefront obtained by a uniformly collimated monochromatic beam of light reflected from an object is given bywhere is the reflectance of the object and is the phase map of the reflected wavefront which is related to the height map of the object. By approximation, the phase map is given bywhere λ is the wavelength and mod is the phase modulus of . Accordingly, the height map of the object is restricted by the 2π ambiguity, and it can be reconstructed unambiguously from the phase map , provided that the pixel-to-pixel height differences are less than λ/2. Moreover, when focusing the interferometer on a certain plane ringing artifacts occur as a result of unfocused levels which influence focused planes. This effect, which is created due to the propagation of the wavefront and which belongs to the unfocused sections of the measured surface to the same level as the focused sections, is known as the Gibbs phenomenon. The ringing artifacts for a piston object with height differences of 108 μm are illustrated in Figure 1. The figure shows a two-dimensional grayscale representation of the intensities of interferometric measurement of a wavefront of light with μm reflected from the two height levels. Figure 2 is a horizontal intersection of the reconstructed heights.

Figure 1: 2D grayscale representation of the height levels.
Figure 2: Horizontal intersection of that same paradigm.

The gray levels and the reconstructed heights, of both the upper and lower segments of the surface, are identical due to the 2π ambiguity. Moreover, due to the focusing problem, a distortion of ±15 nm is perceived on the border of the adjoining height levels. It is important to note that unwrapping is inefficient in this case. The variety of potential levels, complying with the condition of , makes it impossible to deduce the real deviation.

2.2. The Focusing Criteria Package

There are many existing autofocus methods, which are used for various types of optical devices. Most of them are based on sharpness function, which is a real-valued estimate of the image’s sharpness [7]. In most of these methods, several images are taken in different planes, and the sharpness function is calculated where the focal plane is defined as the plane for which it reaches a single optimum. As reviewed in [7], the existing sharpness functions are based on several techniques and criteria. Among others, we can find the image derivative [8] and the Fourier Transform [9].

As part of the above presented methods, several “focusing criteria” have been adapted and developed by our team. These criteria define the wavefront of a focused image and are based on the changes in diffraction and Gibbs phenomena when a wavefront is propagated from a plane to another one.

(i) The Fourier Transform Criterion. In this criterion [10], the spatial behavior of the wavefront is checked according to the integral of the absolute value of the Fourier Transform of the squared wavefront [11]:The rationale behind this method is as follows: since in focused image the amplitude and phase changes in the image, as a consequence of the diffraction and Gibbs phenomena, are the smallest, at least locally, the value of the integral will be also the smallest for a focused image.

(ii) The Standard Deviation Criterion. This criterion checks the shape of the spatial wavefront behavior as a function of the local standard deviation of the squared wavefront [10]:where represents the standard deviation operator. The assumption, on which this criterion is based, is as follows: since, for nonfocused image, the changes in the amplitude in the image are much stronger than for a focused image, the value of the overall standard deviation of the squared wavefront will be the smallest at the focused image.

(iii) The Gradient Criterion. Based on this criterion, the focusing level of the image is checked as a function of the local integral of the modulus of the local gradient of the square wavefront [10]:Since the changes in amplitude and phase in a focused image, as a function of diffraction and Gibbs phenomena [11, 12], are the smallest, the value of the integral will locally remain the smallest for a focused image.

(iv) The Entropy Criterion. At the end, the last criterion is based on the previous one and is the local integral of the modulus of the gradient of the complex wavefront. This criterion takes also into account not only the intensity function as the abovementioned methods [8, 9] but also the phase function of the complex amplitude of the image [10]:

2.3. Beam Propagation

There are two mathematically equivalent formulas which enable the propagation of a wavefront, depending on whether it is preferred to decompose the wavefront as a sum of spherical waves or as a sum of plane waves. For the spherical wave decomposition, the formula of wavefront propagation to a plane located at a distance from the plane positioned at is [10]where is the wavelength and is the sign of , and the symbol stands for convolution. For the plane waves’ decomposition, the formula of wavefront propagation, known also as the “angular spectrum,” is given by [11]where and are the spatial frequencies and stands for the Fourier Transform. For both (7) and (9), the Fast-Fourier Transform (FFT) algorithm can be implemented and the diffracted wave can therefore be computed very efficiently.

3. Results

3.1. Refocusing

Figures 3(a), 3(b), and 3(c) show an image of a measured wavefront of an object illuminated by coherent light where in 3(a) the object is located at image plane and in 3(b) and 3(c) the object is, respectively, in defocus of 400 μm and −400 μm. As a result the defocused image is blurred and possesses a greater amount of ripples which occur around edges, relative to the focused image.

Figure 3: An image of an object illuminated by coherent light. (a) At Image plane, (b) defocus of 400 μm, and (c) defocus of −400 μm.
3.2. Virtual Autofocus

The defocused wavefront was propagated by means of BPM program, towards the image plane while calculating the focus criteria values. In Figure 4 we present the calculated focus criteria values for continuous propagating distances. In this figure, we can clearly observe that the criteria values are changing as a function of the wavefront propagation and that the minimum value is near the focal plan which is 400 μm. The added value of such an approach (Figure 5) is that there is no need for any additional mechanical equipment, since it is realizable through the algorithms of the wavefront propagation. Such a technique can be called Virtual Autofocus.

Figure 4: Focus criteria values measured for defocus of 390 μm. (a) Std, (b) entropy, (c) FFT, and (d) gradient.
Figure 5: Flow chart of the first method for resolving the 2π ambiguity.
3.3. 2π Ambiguity

This technique can also be utilized for resolving the 2π ambiguity when refocusing each segment of the wavefront locally. Once the wavefront has been measured for a certain plane , each segment of the wavefront can be propagated to a certain plane at which the image has a minimum amount of ripples around edges. This plane can be defined as the image plane of that segment. The distance by which each part of the wavefront has been propagated provides the additional information for resolving the 2π ambiguity. We have measured the wavefront reflected from an object containing two segments, a substrate and a step height as shown in Figure 6. The object is a chrome-plated certified step-height standard of 929 nm height (by VLSI-standards) which is well beyond the 2π ambiguity range. The measurements were taken with an imaging interferometer [3, 11] mounted on a standard Olympus white light microscope, using a magnification of ×100, a NA of 0.9, and a 680 nm wavelength illumination. The different segments of the measured wavefront, the substrate and the step height, were propagated separately to divergent planes. As before, the image of each segment which appeared to include minimum distortions was defined as the focused image of that segment. The difference between the propagation distances to the image plane of the step height and the image plane of the substrate is the height of the step relative to the substrate. The height step that was calculated by the global approach was 896 nm while the height step calculated by the local approach was 912 nm. The variation between these results is within 1.7% from each other and within 1.5% from the nominal height of the object. However, as long as the accuracy of the wavefront propagation result is better than λ/2, it can be combined with the interferometric measurements to resolve the 2π ambiguity. Accordingly, the interferometric measurement gives the most accurate result, while the digital wavefront propagation provides additional information for the solution of 2π ambiguities. It should be noted that this method requires some a priori knowledge about the configuration of the step-like object and a minimal number of pixels within each object’s segment in order to obtain an unambiguous result.

Figure 6: A certified step-height standard of 929 nm height (by VLSI-standards).
3.4. A Suggested Method for Resolving the 2π Ambiguity

Here we suggest another method (Figure 7) for resolving the 2π ambiguity by measuring the object using two or more discreet wavelengths.

Figure 7: Flow chart of the second method for resolving the 2π ambiguity.

The object’s heights are obtained in the same fashion as in noncoherent interferometry, by adding the interference patterns obtained by the different discreet wavelengths to obtain the interferogram and calculating the position of the centered maximum of the interferogram. The main advantage here is that, as opposed to noncoherent interferometry, there is no need for moving mirror to obtain the interferogram. Since the light for each wavelength is coherent, the measured wavefronts can be propagated electronically. It should be noted that using discreet number of wavelengths the interferogram is periodic and the resolution is limited. However, since this technique may be used only for resolving the 2π ambiguity, the resolution is not so important.

We have simulated this technique by carefully selecting three wavelengths that will yield satisfactory results: = 720 nm, = 740 nm, and = 760 nm. An interferogram generated by propagating the three wavefronts is shown in Figure 8. On the other hand, Figure 9 shows a sequence of the discrete spectrum interferograms at divergent levels.

Figure 8: Discrete spectrum interferogram.
Figure 9: Sequence of interferograms.

At differing heights of the sample, the same interferogram initializes at distinct positions, allowing the retrieval of sample’s shape. As illustrated above, the centered maximum value relocates itself linearly as a function of depth changes. Nevertheless, as perceived from the figure, amplitudes can coincide for two different levels, and, as a consequence, several multiples of the periodic cycle are exactly lying, one on each other. Truly, such confinement can be neglected. As for ranges of 15 to 55 microns it can be detected, depending on the opted wavelengths. For longer wavelengths this range increases, and therefore usage of infrared illumination is indicated. Another criterion influencing the periodic segment is . A small is preferable since closer wavelengths generate large periodic segments. We have checked the immunity of this method to noise by adding random noise to the interferograms as shown in Figure 10.

Figure 10: Maxima for both surface sections considering a noised interferogram.

According to our results, even in presence of noise artifacts it is feasible to distinguish the precise position of each height according to the amplitude of intensity modulations. Moreover, it is well known that the usage of multiple wavelengths is sensitive to noise, especially for close wavelengths. Our method is not influenced by such noise since it seeks for a graph’s maximum obtained as a result of multiple propagations and for a maxima obtained from the global propagation of each wavefront in each field of vision and not only for local results. So, even if there is any local concern, the global propagation enables its repair.

4. Conclusions

In conclusion, we have introduced two new methods in order to resolve the 2π ambiguity in interferometric measurements. The first method consists in beam propagation of the measured wavefront. This method enables the extension of the dynamic range of single wavelength interferometry beyond the limitation of λ/2 even for step-like objects. The second technique utilizes several wavelengths in order to acquire a discrete wavelengths interferometry measurements as well as obtaining the interferograms simply by electronic propagating and adding the wavefronts without the need for moving the interferometer’s mirror. Divergent sections of the interferogram represent differing heights, while each section undergoes different intensity modulations. We have demonstrated and implemented these methods, proving their promising usefulness. In addition to resolving the 2π ambiguity, these procedures obviate the need for focusing the interferometer on the measured object.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


  1. J. E. Greivenkamp and J. H. Bruning, “Phase shifting interferometry,” in Optical Shop Testing, D. Malacara, Ed., Wiley, New York, NY, USA, 1992, Chapter 14. View at Google Scholar
  2. K. Creath, Phase Measurements Interferometry Techniques, Elsevier Science Publishers, Amsterdam, 1988. View at Publisher · View at Google Scholar · View at Scopus
  3. S. Wolfling, D. Banitt, N. Ben-Yosef, and Y. Arieli, “Innovative metrology method for the 3-dimensional measurement of MEMS structures,” in Reliability, Testing and Characterization of MEMS/MOEMS III, Proceedings of SPIE, D. E. Tanner and R. Ramesham, Eds., vol. 5343, pp. 255–263, SPIE, Bellingham, Washington DC, Wash, USA.
  4. G. Fornaro, G. Franceschetti, R. Lanari, and E. Sansosti, “Robust phase-unwrapping techniques: A comparison,” Journal of the Optical Society of America A: Optics and Image Science, and Vision, vol. 13, no. 12, pp. 2355–2366, 1996. View at Publisher · View at Google Scholar · View at Scopus
  5. A. Collaro, G. Franceschetti, G. Franceschetti, F. Palmieri, and M. S. Ferreiro, “Phase unwrapping by means of genetic algorithms,” Journal of the Optical Society of America A: Optics and Image Science, and Vision, vol. 15, no. 2, pp. 407–418, 1998. View at Publisher · View at Google Scholar · View at Scopus
  6. Y.-Y. Cheng and J. C. Wyant, “Two-wavelength phase shifting interferometry,” Applied Optics, vol. 23, no. 24, pp. 4539–4543, 1984. View at Publisher · View at Google Scholar · View at Scopus
  7. M. E. Rudnaya, H. G. ter Morsche, J. M. Maubach, and R. M. Mattheij, “A derivative-based fast autofocus method in electron microscopy,” Journal of Mathematical Imaging and Vision, vol. 44, no. 1, pp. 38–51, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  8. G. Yang and B. J. Nelson, “Wavelet-based autofocusing and unsupervised segmentation of microscopic images,” in Procceeding of the Conference on Intelligent Robots and Systems, Canada, 2003.
  9. M. E. Rudnaya, R. M. Mattheij, and J. M. Maubach, “Evaluating sharpness functions for automated scanning electron microscopy,” Journal of Microscopy, vol. 240, no. 1, pp. 38–49, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. E. Novoselski, A. Yifrach, and Y. Arieli, Wavefront Focusing of Coherent Light, Jerusalem College of Technology, 2006.
  11. J. W. Goodman, Introduction to Fourier Optics, McGraw-Hill, 2nd edition, 1996.
  12. E. G. Steward, Fourier Optics - An Introduction, Ellis Horwood Limited, 2nd edition, 1987.