## Abstract

A two-step self-tuning phase-shifting method is presented. The phase-step between the two interferograms is not known when the experiment is performed. Our demodulating method finds, in a robust way, this unknown phase-step. Once the phase-step is estimated we proceed to phase demodulate the interferograms. Moreover our method only requires the fringe patterns to have a constant unknown phase-shift between them. As a consequence, this technique can be used to demodulate open and closed-fringed patterns without phase-sign ambiguity. The method may be regarded as a self-tuning quadrature filter, which determines the phase-shift between the two fringe patterns and finally estimates the demodulated phase map. The proposed technique has been tested with simulated and real interferograms obtaining satisfactory results.

©2011 Optical Society of America

## 1. Introduction

Phase-Shifting Interferometry (PSI) is a useful technique in optical metrology for measuring the modulated phase of interferograms [1]. In standard PSI, a sequence of *N* interferograms is obtained having a known temporal carrier. Usually a minimum of three phase-shifted interferograms are need to retrieve the phase. Other phase-shifting methods have been reported as the four-step, five-step or the least-squares algorithm, for example. A good general overview of the phase-shifting algorithms can be found in [2]. The use of fewer interferograms in PSI simplifies the computation process and the processing speed. Moreover, it increases the robustness against uncontrolled mechanical vibrations, air turbulence or temperature changes in the measuring system. In the case of open-fringe interferograms, there have been proposed many single-frame demodulating methods but in the case of closed-fringes, it is needed at least-two frames to solve the local phase-sign ambiguity. Note that a global sign ambiguity in the retrieved phase is irrelevant.

In the past there have been reported works about phase reconstruction with only two-frames. Ref [3]. proposes a method for demodulating two interferograms with a known phase-shifting of $\frac{\pi}{2}$ between them. In Refs [4,5] are proposed two-step phase-shifting algorithms that don’t need knowing the phase step between interferograms. The main drawbacks of [4,5] are that the methods don’t work if the phase-shift is close to π radians and in the case of interferograms with low spatial frequency fringes. Ref [6]. shows an alternative approach that uses the least-squares method to extract the unknown phase-shift between two interferograms and then obtain the phase distribution. This method requires many iterations to converge (typically 17 iterations) and is sensitive to noise as only two frames are used. More recently, Ref [7]. shows a method to extract the arbitrary unknown phase-shift between two interferograms and then reconstruct the complex object wave without iteration process. This technique only can be used in amplitude division interferometry as it needs to measure, additionally with the two phase-shifted frames, the intensity of the reference beam. On the other hand, this method requires the reference intensity to be constant along the full aperture that can be a non realistic case when we deal with real interferograms. In addition, the method assumes that the intensity of the reference beam is greater than the maximum of the object beam that can be problematic in some cases.

In this work, we propose a self-tuning two-step phase-shifting algorithm that doesn’t need to known the phase-step between interferograms. This phase-step can be any value between [-π, π]. The proposed method doesn’t need to iterate, it can be used to any kind of fringe pattern as it doesn’t need to measure the intensity of the reference beam and no suppositions are made to the fringe spatial frequency, background and modulation maps.

The structure of the paper is as follow: in the next section we introduce our tunable two-step phase-shifting algorithm to determine the temporal carrier. After this, the modulated phase is obtained. In sections 3 and 4 the simulation and experimental results are given. In the simulation results section we have applied the proposed method to two simulations, one with open and other with closed fringes. Finally, in section 5 some conclusions are drawn.

## 2. Proposed method

In PSI, a discrete temporal interferogram can be described using the following expression,

*t*denotes the

*t*

^{th}interferogram obtained at time

*t*. In our case

*t*= [0, 1]. For simplicity in the following we will drop the spatial dependence. Expression (1) can be rewritten as follow,

*δ*is the Dirac delta function. In order to obtain the desired phase map, we can define a complex quadrature filter with frequency response $H\left(\omega \right)$ and with the following requirements$H\left(0\right)=0$, $H\left({\omega}_{0}\right)=0$ and $H\left(-{\omega}_{0}\right)\ne 0$ [2]. The filtered complex function, $\tilde{G}\left(\omega \right)$, is obtained by multiplication of $G\left(\omega \right)$ and $H\left(\omega \right)$ as,

*t*= 0 obtaining the angle and the norm of $\tilde{g}\left(0\right)$ respectively as,

The mathematical analysis shown above assumes that the temporal frequency ${\omega}_{0}$ is known. On the other hand, as we said, in order to built a practical quadrature filter that achieves the requirement $H\left(0\right)=0$ and $H\left({\omega}_{0}\right)=0$, it is necessary at least three temporal interferograms [1].

In a self-tuning phase-shifting algorithm the temporal frequency ${\omega}_{0}$ is unknown and it is introduced as a free parameter of the algorithm [8]. The algorithm has to determine this temporal frequency by its own. As we only have available two interferograms, it is not possible to obtain a quadrature filter that achieves $H\left(0\right)=0$ and therefore the filter will not subtract the interferograms background. Consequently, a necessary previous step is to eliminate the background component in expression (1). One possibility is to normalize the interferograms. If the interferograms are normalized, the background is filtered-out and the zero-frequency filter condition, $H\left(0\right)=0$ is relaxed. In this case, the quadrature filter requirements are only $H\left({\omega}_{0}\right)=0$ and $H\left(-{\omega}_{0}\right)\ne 0$. The normalized version of expression (1) is

The normalizing method used in this work has been reported in [9]. The method shown in [9] consists in employing an isotropic operator. This method is robust to wide-band fringe spatial frequencies, obtaining good solutions.

The proposed self-tuning phase-shifting algorithm is based on analyzing ${\left|\tilde{g}\left(0\right)\right|}^{2}$ parameter when a detuning error exists. This detuning error appears when the quadrature filter $H\left(\omega \right)$ is tuned in a temporal frequency ${\omega}_{f}$ different to the assumed one ${\omega}_{0}$ [2]. Figure 1 shows a scheme of a $H\left(\omega \right)$ filter when there is non detuning error (Fig. 1 (a)) and when there is (Fig. 1 (b)). In Fig. 1(a) it is shown the case where ${\omega}_{f}$ coincides with ${\omega}_{0}$ and there is non detuning error. In Fig. 1(b) it is shown the case where ${\omega}_{f}\ne {\omega}_{0}$ and therefore, $H\left({\omega}_{0}\right)\ne 0$ and $H\left(-{\omega}_{0}\right)\ne 0$. In this case, a detuning error does exist and expression (4) is not longer valid.

This expression is rewritten in the case of ${\omega}_{f}\ne {\omega}_{0}$ as,

*Φ*is determined using expressions (4) to (6).

In this work, in order to accurately obtaining the temporal carrier ${\omega}_{0}$ by the proposed self-tuning method explained above, it is very important to use a quadrature filter $H\left(\omega \right)$ with a very sharp frequency response near ${\omega}_{f}$. This sharp response allows one to obtain a clear and narrow minimum in $\sigma \left[{\left|\tilde{g}\left(0\right)\right|}^{2}\right]$ and therefore determining accurately${\omega}_{0}$. The quadrature filter that we have used is [2]

And its corresponding temporal impulse response is given by,*Φ*is obtained.

This method uses exhaustive search, thus, it is as fast as the size of the discrete searching tuning frequency set. This is not a serious concern given that the computational cost of exploring a sampled one-dimensional space $\left[-\pi ,\pi \right]$ exhaustively with, say 600 samples, is not high.

## 3. Simulations

In order to show the performance of the proposed method, we have tested it with a simulation using closed-fringe interferograms. The utility of this method has special interest in the case of closed-fringe interferograms as two-frames are the minimum number of interferograms to reconstruct the phase without local phase-sign ambiguity. In Fig. 2 are shown the two fringe patterns of the simulation. As can be seen from Fig. 2, the interferograms have a size of 411x411 pixels2. The background and modulation maps of the fringe patterns are equal to zero and 125 a.u (arbitrary units) in every pixel. The noise is Gaussian with zero mean and a variance of 13 a.u. The temporal frequency ${\omega}_{0}$ is equals to 1.2 rad.

In Fig. 3
, it is shown the plot between the standard deviation of the modulation map $\sigma \left[{\left|\tilde{g}\left(0\right)\right|}^{2}\right]$ and the different temporal frequencies${\omega}_{f}$. The discrete searching tuning-frequency set is composed by 600 frequencies equally spaced between [0, 2*π*] and the processing time is of 5 s using a 1.6 GHz laptop and MATLAB. The detected minimum corresponds to ${\stackrel{\u2322}{\omega}}_{0}$ = 1.203 rad and it is marked in the figure with a red square.

Observe in Fig. 3 that are obtained two minimums, one at ${\stackrel{\u2322}{\omega}}_{0}$ and other at $2\pi -{\stackrel{\u2322}{\omega}}_{0}$. Note that it is irrelevant to select one or the other as the only consequence will be a global phase-sign change in the reconstructed phase *Φ*. In Fig. 4
, it is shown the obtained wrapped phase. Finally, in Fig. 5
are shown the reconstructed phase and the computed error between actual and obtained phases. The *rms* (root-mean-square) error of the difference between the actual and measured phases is 0.08 rad.

In order to show the performance of the two-step self-tuning algorithm to noisy interferograms, we have obtained for different signal to noise ratios the computed temporal frequency${\widehat{\omega}}_{0}$, the absolute difference between the actual and computed temporal frequencies $\left|{\omega}_{0}-{\widehat{\omega}}_{0}\right|$, the *rms* error between the actual and computed phase maps and the processing times. The computed phase map is retrieved by interferograms without noise in order to show the detuning error caused by the mismatch between the actual and obtained temporal frequencies. In Table 1
, the first row corresponds to the different signal to noise ratios that quantify how much the phase has been corrupted by noise. As can be seen from Table 1, the proposed two-step self-tuning algorithm gives appropriate results even in the case of interferograms with high noise.

## 4. Experimental results

We have applied the proposed algorithm to real interferograms. In order to check the obtained results, we have compared the retrieved phase from our method with the phase obtained from 10 phase-shifted interferograms, and using the least-squares phase demodulation method [10]. In Fig. 6 are shown the two interferograms used in our approach and the normalized version of them.

In Fig. 7(a) and 7(b) are shown the acquired wrapped phase obtained from our method and with the least-squares method using 10 interferograms. The processing times using the proposed method and the least-squares method are 12 s and 30 s respectively.

As can be seen from Figs. 7(a) and 7(b) both wrapped phases are similar. In Figs. 8(a) and 8(b) are shown the reconstructed phases with both methods and finally in Fig. 9 it is shown the difference between the reconstructed phases.

As can be seen from Fig. 9 there is a good agreement between both measures. The *rms* (root-mean-square) error of the difference between both reconstructed phases is 0.3 rad. On the other hand, the temporal frequency between interferograms obtained by the proposed and the least-squares methods is ${\stackrel{\u2322}{\omega}}_{0}$ = 0.7 rad and ${\stackrel{\u2322}{\omega}}_{0}$ = 0.85 rad respectively.

## 5. Conclusions

In this work, we have proposed a self-tuning two-step phase-shifting algorithm that doesn’t need to known the phase-step between interferograms. This phase-step can be any value between [-π, π]. The proposed method doesn’t need to iterate, however our technique uses an exhaustive search strategy. This method can be used with any kind of fringe pattern as it doesn’t need to measure the intensity of the reference beam and no suppositions are made to the fringe spatial frequencies, background and modulation maps. We have tested the proposed method to simulated and real interferograms, and we have shown a very good performance of the proposed method.

## References and links

**1. **D. Malacara, M. Servín, and Z. Malacara, “Interferogram analisis for optical testing”, Cambridge University Press, (2004), Marcel Dekker, Inc, (1998)

**2. **M. Servin, J. C. Estrada, and J. A. Quiroga, “The general theory of phase shifting algorithms,” Opt. Express **17**(24), 21867–21881 (2009). [CrossRef] [PubMed]

**3. **F. Mendoza-Santoyo, D. Kerr, and J. R. Tyrer, “Interferometric fringe analysis using a single phase step technique,” Appl. Opt. **27**, 4362–4364 (1988).

**4. **S. Almazán-Cuéllar and D. Malacara-Hernandez, “Two-step phase-shifting algorithm,” Opt. Eng. **42**(12), 3524–3531 (2003). [CrossRef]

**5. **Y. Zhu, L. Liu, Z. Luana, and J. Sun, “Discussions about FFT-based two-step phase-shifting algorithm,” Optik (Stuttg.) **119**(9), 424–428 (2008). [CrossRef]

**6. **X. F. Xu, L. Z. Cai, Y. R. Wanga, X. F. Meng, H. Zhang, G. Y. Dong, and X. X. Shen, “Blind phase shift extraction and wavefront retrieval by two-frame phase-shifting interferometry with an unknown phase shift,” Opt. Commun. **273**(1), 54–59 (2007). [CrossRef]

**7. **X. F. Xu, L. Z. Cai, Y. R. Wang, X. F. Meng, W. J. Sun, H. Zhang, X. C. Cheng, G. Y. Dong, and X. X. Shen, “Simple direct extraction of unknown phase shift and wavefront reconstruction in generalized phase-shifting interferometry: algorithm and experiments,” Opt. Lett. **33**(8), 776–778 (2008). [CrossRef] [PubMed]

**8. **M. Servin, J. C. Estrada, and J. A. Quiroga, “Spectral analysis of phase shifting algorithms,” Opt. Express **17**(19), 16423–16428 (2009). [CrossRef] [PubMed]

**9. **J. A. Quiroga and M. Servín, “Isotropic n-dimensional fringe pattern normalization,” Opt. Commun. **224**(4-6), 221–227 (2003). [CrossRef]

**10. **Z. Wang and B. Han, “Advanced iterative algorithm for phase extraction of randomly phase-shifted interferograms,” Opt. Lett. **29**(14), 1671–1673 (2004). [CrossRef] [PubMed]