Phasorial analysis of detuning error in temporal phase shifting algorithms

Phase error analysis in Temporal Phase Shifting (TPS) algorithms due to frequency detuning has been to date only performed numerically. In this paper, we show an exact analytical expression to obtain this phase error due to detuning using the spectral TPS response. The new proposed method is based on the phasorial representation of the output of the TPS quadrature filter. Doing this, the detuning problem is reduced to a ratio of two symmetrical spectral responses of the quadrature filter at the detuned frequency. Finally, some popular cases of TPS algorithms are analyzed to show the usefulness of the proposed method. 2009 Optical Society of America OCIS codes: (120.3180) Interferometry; (120.2650) Fringe analysis; (120.5050) Phase Measurement. References and links 1. J. Schwider, “Advanced evaluation techniques in interferometry,” in Progress in Optics, E. Wolf ed., (North Holland, Amsterdam, Oxford, New York, Tokyo, 1990). 2. J. H. Bruning, D. R. Herriot, J. E. Gallagher, D. P. Rosenfel, A. D. White, and D. J. Brangaccio, “Digital wavefront measuring interferometry for testing optical surfaces and lenses,” Appl. Opt. 13, 2693-2703 (1974). 3. H. Schreiber and J. H. Brunning, “Phase shifting interferometry,” in Optical Shop Testing, D. Malacara ed., (John Wiley & Sons, Inc., Hoboken, New Jersey 2007). 4. M. Servin and M. Kujawinska, “Modern fringe pattern analysis in Interferometry,” in Handbook of Optical Engineering, D. Malacara and B. J. Thompson eds., (Marcel Dekker, 2001). 5. J. Schwider, R. Burrow, K. E. Elssner, J. Grzanna, R. Spolaczyk, and K. Merkel, “Digital wave-front measuring interferometry: some systematic error sources,” Appl. Opt. 22, 3421-3432 (1983). 6. J P. Hariharan, B. Oreb, and T. Eiju, “Digital phase shifting interferometry: a simple error compensating phase calculation algorithm,” Appl. Opt. 26, 2504-2505 (1987). 7. K. Freischland and C. L. Koliopoulos, “Fourier description of digital phase measuring interferometry,” J. Opt. Soc. Am. A 7, 542-551 (1990). 8. M. Servin, D. Malacara, J. L. Marroquin and F. J. Cuevas, “Complex linear filters for phase shifting with very low detuning sensitivity,” J. Mod. Opt. 44, 1269-1278 (1997).


Introduction
Among automated interferogram analysis methods, temporal phase shifting (TPS) techniques are considered as the most accurate wave-front extraction [1][2][3][4][5][6][7][8].The measurement accuracy of the TPS depends on the estimated phase of the complex analytical signal associated with the real cosine signal given by the observed interferogram.The phase-shifting techniques utilize several TPS fringe patterns to determine this complex analytical signal whose angle is the desired signal to be estimated.All TPS algorithms give the exact fringe pattern phase (assuming very low measuring noise) when the temporal signal's carrier 0 this phase estimation error [1][2][3][4][5][6][7][8].The analytical expression herein derived for the detuning error in TPS algorithms allows us to minimize some TPS quadrature filters for detuning error.Most TPS algorithms used nowadays have no free parameters that permit us this kind of detuning optimization.However, some new techniques may be used to generate more general TPS algorithms including free parameters that allow us to improve the robustness of these new TPS filters to detuning.Even if no optimization is possible, one may always compare among the myriad of TPS algorithms available in the literature using the analytical expression given in this paper.The input of any TPS quadrature filter is a set of real cosine signals which are the intensity of a series of TPS interferograms captured by a CCD camera.The output of these TPS quadrature filters is the complex analytical signal at the temporal carrier frequency -0 ω associated with the real intensity of our interferograms.As mentioned before, when the carrier frequency of the interferograms is detuned, an additional undesired complex signal with carrier frequency + 0 ω is added to our desired signal at -0 ω .Hence the complexvectorial (or phasor) sum of both signals give rise to an erroneously (detuned) estimated phase from this set of temporal interferograms.

Error Detuning in Phase-Shifting Interferometry (PSI)
The measured intensity of an interferogram on a CCD detector can be expressed as [3,4]: Where, ( ) where the spatial dependence (x,y) of the functions I, a, b and φ has been dropped.Now, we can observe in Fig. 1, that at (exactly) the carrier frequency Therefore, the spectrum of the analytical signal at the output of the TPS algorithm Then the searched phase φ is recovered from the complex analytical signal ) , , ( t y x g which is its inverse Fourier transform.Now let us consider that our phase stepping interferograms are not sampled at the expected frequency rate 0 ω , but at an erroneous temporal frequency given by The two (complex) output signals given by Eq. ( 5) are also graphically depicted in Fig. 1 in the frequency domain.We may also use a graphical-phasor representation of these two complex signals given in Eq. ( 5).This is shown in Fig. 2, where we have denoted the desired analytical signal as c and the undesired or spurious signal as ε, ( ) then Eq. ( 5) becomes, ) Equation ( 7) may be graphically represented as phasors as in Fig. 2, where we show the desired signal c with an angle ) , ( y x φ and the erroneous spurious signal ε with angle ) , ( y x φ − .( ) The detuning error may be defined as the difference between the desired phase φ and the undesired phase φ′ .This expression for the detuning error φ ∆ is widely used to evaluate the robustness to detuning in TPS and it is has been evaluated only numerically [1][2][3][4][5][6][7][8].
From of Eq. ( 8) we obtain φ′ , then by using this into φ φ φ , we find that the detuning error is expressed as, We can see in this expression that when 0 ε → no detuning error is present and the erroneous phase φ′ becomes the desired phase φ φ →

'
. The main objective in this paper is to find an easy, useful and analytical detuning robustness equation in terms of the frequency response ) (ω H of any TPS algorithm.To this end, we must substitute Eq. ( 8) into Eq.( 9) and after using some trigonometric relations, the searched expression for the detuning error φ ∆ in terms of the frequency response of the TPS algorithm is, ( ) ( ) ( ) This analytical expression is a novel, very simple and an exact way to evaluate the detuning error.This expression may be further simplified by taking into account that for all practical purposes 0 . 1 ) / ( << c ε , as a consequence of this the finally Eq. (10) becomes, This shows the well known fact (without proof [1][2][3][4][5][6][7][8]) that the detuned phase always has a component which depends on the interferogram's fringes having twice the interference fringes which is ) 2 sin( φ .The Eq. ( 11) is the most important result in this paper.A short-cut way to obtain the last result given by Eq. (11) from Eq. ( 8), is to consider that the detuning error is small.Therefore the following approximations apply , as a consequence the detuning error becomes which is equal to Eq. ( 11).

Some examples of error detuning in Phase-Shifting Interferometry.
In this section we analyze some popular TPS algorithms just for illustrative purposes to show how to apply the general techniques developed so far.Let us start with the Schwider-Hariharan five step algorithm [6],

The five frames algorithm:
The most used five steps TPS algorithm (Schwider-Hariharan) is given by [1,3, is the interferogram's phase for phase step / π α = .We first obtain the time response of the quadrature filter from the TPS formula Eq. ( 13) with as, ; then, we obtain the frequency response ( ) This filter is tuned at frequency 1 = ω .Now assuming a detuning error in the phase step of ∆ + = 2 / π α and according to Eq. ( 8) we obtain, Finally using this result into Eq.( 10) we have the exact detuned phase error for the five frames TPS algorithm as, This equation is the exact detuning error for any detuned step ∆ + = 2 / π α .This result may be approximated for small ∆ (as reported in [3,4,6]) to obtain, . It can be seen that the maximum detuning is when the sine function equals one so 4 / 2 max ∆ = ∆φ renders the well known quadratic detuning error behavior of the five step Schwider-Hariharan TPS algorithm [5,7].

Seven frames algorithm
The seven steps TPS algorithm given by the Servin's algorithm [8] is, Then, the frequency response of this TPS algorithm is, Using the first order approximation given by cos(n∆)=1.0and sin(n∆)=n∆ for any n in Eq. (20) the maximum detuning error is We can observe the extremely low sensitivity to the detuning reported only numerically in [4,8].

Eleven frames algorithm
Finally, just for illustrative purposes we briefly analyze the eleven steps TPS algorithm (De Groot in [3]) which is, ( Applying the same easy procedure as before we obtain, ( ) For small ∆ we may approximate this detuning error as 16 / 4 max ∆ ≈ ∆φ .As a final point, we graphically compare in Fig. 3 all the TPS algorithms that were analyzed.

Conclusions
An exact, analytical algorithm to evaluate the error detuning in phase shifting algorithms was obtained from a phasor representation of the detuned output signal of any TPS algorithm.The derived expression was applied to evaluate useful approximations to the error detuning from some well known quadrature filters with five, seven, and eleven steps.

Fig. 2 .
Fig. 2. Phasor representation of G(x,y,ω).Now we proceed to the analysis of the phasor diagram given in Fig.2, by applying the sinus law to the triangle formed by c, ε and R which stand for,

Fig. 3 .
Fig. 3. Error detuning for TPS algorithms for five, seven and eleven steps.