Abstract
The main purpose of this paper is to present a method to design tunable quadrature filters in phase shifting interferometry. The algorithm is obtained from a generalized Fourier transform of a symmetrical quadrature filter. This formalism allows us to represent the detuning phase shift error and bias modulation as geometrical conditions. Therefore, the design of the filter becomes a set of solvable linear equations. Hence, to prove our method, several general tunable filters, like three and four frame algorithms, are obtained. Finally, from our results we reproduce particular symmetrical four frame algorithms reported in literature.
©2009 Optical Society of America
1. Introduction
Among regular methods to analyze an interferogram, temporal phase shifting (TPS) techniques are considered as one of the most accurate wave-front extraction methods [1–8]. However, the accuracy of these measurements is limited by the presence of several kinds of systematic errors. Two examples are the detuning error due to the miscalibration of the phase shifter and the bias modulation introduced by the use of a source light like a laser diode. In both cases, the errors cannot be eliminated totally, and they may be only minimized. The usual technique to minimize systematic errors is the use of a quadrature filter insensitive to such errors [1–8]. There are several methods to design quadrature filters and a myriad of algorithms exist in the literature [1–8]; however, most of them only give us an algorithm designed for a specific phase step. Hence, in this paper we introduce a simple and practical method to obtain tunable filters. First, the filter is obtained from the general form of the Fourier transform of a quadrature filter. The quadrature conditions of the filter and errors like the detuning error and bias modulation are represented as geometrical conditions to be satisfied by the filter. That is, from the quadrature conditions and the desired errors to minimize, we obtain a set of solvable linear equations. In other words, from the Fourier transform of the general filter, we apply the geometrical conditions to recover a desired filter. Furthermore, we apply the method to obtain some examples like three and four frame algorithms which are tunable in all cases. Finally, we show that our results for tunable four frame algorithms coincide with particular cases reported in literature, and we also report a table containing the best known ones.
2. Phase-Shifting Interferometry (PSI)
Assuming very low noise, the measured intensity of an interferogram obtained from a CCD detector is expressed as [1–8]
Where, denotes the unknown phase, is the background illumination, and is the contrast of interference fringes. The temporal carrier is a linear phase shift among the set of interferograms. Meanwhile, t corresponds with the temporal sampling, which is taken as a natural number in this paper. To recover the phase in time domain we convolve a set of several temporal phase shifted interferograms with a quadrature filter that has a temporal impulse response [7–9], and it yields an output signal given by where (*) denotes one dimensional temporal convolution. Therefore, in the Fourier domain we have that this temporal convolution becomes , where and are given byNotice that the spatial dependence (x,y) of the functions I, a, b and φ has been dropped.
Now, considering a symmetric quadrature filter, the Fourier transform H(ω) becomes a real function [9]. Then, assuming that the quadrature filter is tuned onto the right side of the frequency axis, as depicted in Fig. (1) by the solid curve, at exactly the carrier frequency , we have the quadrature conditions as and , while . Therefore, from Eq. (3) output signal becomes
That is, the estimated phase for a filter tuned onto the right side is given by or .
Another possibility to recover the desired phase is to find a quadrature filter tuned onto the left side of the frequency axis. Then; as shown in Fig. (1) by the dashed curve, at exactly the carrier frequency with , the quadrature conditions are given by and . Therefore, output signal becomes,
Then, the estimated phase is φ. Thus, we can conclude that the design of a quadrature filter can be viewed as a geometrical design problem, where the interception with the frequency axis gives us the quadrature conditions or the cancellation of a specific harmonic frequency. In the same way, the derivate that has a value of zero at a specific frequency gives us the condition for a quadrature filter insensitive to linear phase shift error at that frequency. Finally, when the derivate equals zero at the origin, we have the conditions for a filter which is insensitive to bias variation error [2–5]. Therefore, assuming a symmetrical algorithm with order N and tuned onto the left side of the frequency axis, the estimated phase is
where the coefficients and for symmetrical algorithms satisfy the relationships and . Notice that for odd N, the middle term is [7,8]. Then, the temporal impulse response [2–9] of this filter is given by
where is chosen to obtain a symmetrical filter. Since Eq. (7) corresponds with a quadrature filter, the real part must be an even function, whereas the imaginary part is an odd function. Then, the Fourier transform of h(t) generates the real function H(ω), given by
From Eq. (6) we can observe that to design a desired quadrature filter with N frames as support, only N-1 parameters are independent. Then, taking a step α, the geometrical conditions to satisfy a quadrature filter tuned onto the left side are given by
That is, the filter cuts-off both ω = 0 and ω = −α frequencies. Then, the number of independent parameters to design the filter is reduced to three, like in the three frame algorithm case. However, if filter order N becomes greater than three, we can use these other conditions to eliminate or compensate some kinds of errors, like detuning phase shift or bias variation errors. That is, to design a filter tuned onto the left side and insensitive to the nth order detuning phase shift error, we just apply the following geometrical conditions,
Where, n gives us the order of insensitivity to this phase shift error, and is the nth derivate of with regard to ω.
On the other hand, to design a filter that is also insensitive to the mth order bias variation error, we can add the following conditions
Hence, to obtain a high order of insensitivity to one error, like detuning or bias variation, we simply increase the order of derivate conditions. That is, the maximum number of independent conditions to design a quadrature filter becomes N-3. Then, we can decide the conditions to apply and compensate for some kind of errors like the detuning phase shift error or bias modulation, or simply cancel a specific frequency. These conditions always give us a set of solvable linear equations whose solution produces the coefficients needed for the filter.
In conclusion, the design of a tunable quadrature filter is a geometrical problem where the design conditions are interpreted as the interception over frequencies ω = 0 and ω = −α of H(ω), while the insensitivity to some errors, like phase shift and bias modulation, is interpreted as the n derivates equal to zero. Finally, the problem is reduced to solving a set of linear equations and recovering from Eq. (6) the desired TPS algorithm.
3. Designing some examples of a tunable quadrature filter
The above method was applied to two known algorithms, one odd and another even.
3.1 The tunable three frame algorithm
From Eq. (6), the general symmetrical form of a three frame algorithm is,
As mentioned above, in Eq. (13) only two parameters are independent, Therefore, by simplicity for N odd we choose . Then, from Eq. (9) Fourier transform of Eq. (13) for a step α given by Eq. (9) is,
Tuning this filter onto the left side, it must satisfy the quadrature conditions given by Eq. (10). Then, H(ω) satisfies , and from Eq. (14) we obtain the linear system
Therefore, the from Eq. (15) and from Eq. (16)
Then, from Eq. (13) we obtain the tunable three frame algorithm as,
It should be noticed that choosing a quadrature filter tuned onto the right side of the frequency axis we obtain the same result. Notice that using another value non zero for the same result Eq. (18) is obtained.
3.2 The tunable four frame algorithms
To design a four fame (TPS) algorithm we have only one independent parameter used to minimize one kind of error. Therefore, we have only three options; the first one to obtain a filter equal to the least squares method, or class A filter. The second one corresponds to a filter that is insensitive to the linear detuning phase shift error. Finally, the last option is for a filter that is insensitive to the linear bias variation error.
3.2.3 Tunable four frame filter class A
From Eq. (6) we have that the general form for a four frame algorithm tuned onto right is
Then, from Eq. (8) and taking , the Fourier transform of the filter becomes
Then, to obtain an algorithm through the least squares method, geometrical conditions are
That is, the two conditions given by Eq. (10) are the quadrature conditions that any quadrature filter must satisfy, and is the particular condition of this filter to be insensitive to linear variations of the value H(α). Then, applying Eq. (20) into Eq. (21) we have
And the solution for this filter tuned onto the right of the frequency axis is given by,
Then, by substituting into Eq. (19), we obtain our tunable four frame algorithm as
According to the authors in [5], this filter can be called a tunable four frame filter class A. In a previous work [9], we present a formalism to obtain the detuning error of a symmetrical filter. Then, for a small detuning phase Δ, the detuning error Δφ of this quadrature filter is
As we expect, this algorithm is very sensitive to errors like detuning or bias modulation.
3.2.2 Tunable four frame algorithm insensitive to linear detuning phase shift error
From Eq. (11), a filter that is insensitive to the linear detuning error satisfies the condition
Then, by using the quadrature conditions Eq. (22) and Eq. (23) with Eq. (28), our solution is
Therefore, the tunable algorithm that is insensitive to linear detuning phase shift error is
This filter, according to [5], can be called a tunable filter class B. Then from [9], for small phase detuning, the detuning error calculated is . That is, as we expect, the algorithm is insensitive to the linear detuning phase shift error.
3.2.3 Tunable four frame algorithm insensitive to linear bias variation error
From Eq. (12) to obtain a filter that is insensitive to the linear bias variation error, H(ω) must satisfy the following geometrical condition
Then, combining the quadrature conditions Eq. (22) and Eq. (23) with Eq. (31), we obtain
And from Eq. (20), the tunable four frame filter insensitive to linear bias variation error is
Therefore, this algorithm can be called a tunable filter class C. Then, from [9], the detuning error is Δφ = tan(Δ/2). Finally, in Table 1 , from our solutions showed above, we reproduce particular cases reported in literature [5–8].
4. Conclusions
An exact, analytical algorithm to design quadrature filters as a geometrical problem based on the Fourier description of the quadrature filter is presented. The quadrature conditions and the error are interpreted geometrically, and the problem is reduced to solving a set of linear equations. Then, from such solution we can obtain tunable filters with any combination of error, instead of using other known procedures [1–8]; foremost, we apply the method to obtain some tunable filters for three and four frame cases. From our solutions, we generalize all four frame symmetrical algorithm cases reported in literature, and a table with some particular cases is reported.
References and links
1. K. Freischland and C. L. Koliopoulos, “Fourier description of digital phase measuring interferometry,” J. Opt. Soc. Am. A 7(4), 542–551 (1990). [CrossRef]
2. Y. Surrel, “Design of phase-detection algorithms insensitive to bias modulation,” Appl. Opt. 36(4), 805–807 (1997). [CrossRef] [PubMed]
3. D. W. Phillion, “General methods for generating phase-shifting interferometry algorithms,” Appl. Opt. 31(36), 8098–8115 (1997). [CrossRef]
4. Y. Surrel, “Design of algorithms for phase measurements by the use of phase stepping,” Appl. Opt. 35(1), 51–60 (1996). [CrossRef] [PubMed]
5. J. Schmit and K. Creath, “Extended averaging technique for derivation of error-compensating algorithms in phase-shifting interferometry,” Appl. Opt. 34(19), 3610–3619 (1995). [CrossRef] [PubMed]
6. H. Schreiber, J. H. Brunning, and J. E. Greivenkamp, “Phase shifting interferometry,” in Optical Shop Testing, D. Malacara ed., (John Wiley & Sons, Inc., Hoboken, New Jersey 2007).
7. 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).
8. D. Malacara, M. Servin, and Z. Malacara, “Phase Detection Algorithms,” in Interferogram Analysis for Optical Testing, D. Malacara ed., (Taylor & Francis Group, 2005).
9. J. F. Mosiño, M. Servin, J. C. Estrada, and J. A. Quiroga, “Phasorial analysis of detuning error in temporal phase shifting algorithms,” Opt. Express 17(7), 5618–5623 (2009). [CrossRef] [PubMed]