Open Access
21 February 2013 Generic algorithm of phase reconstruction in phase-shifting interferometry
Vladimir Guzhov, Sergey Ilinykh, Roman Kuznetsov, Dmitry Haydukov
Author Affiliations +
Abstract
A generic reconstruction algorithm is proposed, which is obtained by using a phase-shifting technique with an arbitrary number of phase shifts between intensity measurements. The generic algorithm entirely describes the structure of known phase-shifting algorithms and permits the construction of new ones with an arbitrary number of phase shifts.

1.

Introduction

Phase-shifting interferometry methods have been widely used in optical interference systems in recent years. Widespread use of these methods has resulted from the simplicity of specifying the values of phase shifts, the low complexity of the algorithms, and the high precision that they can achieve. At the same time, the layouts of interferometers can be easily modified.

A large number of expressions for phase reconstruction with an arbitrary number of phase shifts are known. The earliest algorithms used decoding equations with three or four shifts. With an increased amount of available computational power, it became possible to use algorithms with a larger number of shifts. Thus, in Ref. 1, an algorithm that uses 15 phase shifts was introduced, and an algorithm using 101 phase shifts was described in Ref. 2.

The equations used for the reconstruction are derived by solving trigonometric equations. In Ref. 3, an analysis of formulas for phase reconstruction with an arbitrary number of the phase shifts is given, but the number of shifts should be constant. In this work, we propose a generic algorithm that allows the structure of known algorithms to be determined and the construction of new algorithms with an unlimited number of arbitrary, and not always constant, phase shifts.

2.

Synthesis of the Algorithm

Phase-shifting methods are based on the capture of several interferograms while the phase of the reference wave continues to follow the specified values. At different phase shifts, the intensity of an interferogram with phase shift δi can be represented as

Eq. (1)

Ii(x,y)=A(x,y)+B(x,y)cos(ϕ(x,y)+δi),
where A(x,y) and B(x,y) are the average brightness of the interferogram and the amplitude of the interference pattern in the point (x,y), respectively, i=0,1,2,,m1, m is the number of phase shifts, and δ0=0.

The algorithms obtained at m different values are called m-point algorithms. There are many algorithms that work for different numbers of phase shifts. Many of the possible implementations arose from an interest in determining a generic scheme for the algorithms. A generic scheme allows evaluation of the pros and cons of a specific variant and determining the methodology of the algorithm construction.

Equation (1) can be represented in vector form as

Eq. (2)

I=A+(Bcosϕ)C(Bsinϕ)S,
where I is the set of intensities for different phase shifts δi at each point of the interferogram Ii(x,y), A=A·(1,,1)T is an m-dimensional vector, C=(cosδ0,,cosδm1)T, S=(sinδ0,,sinδm1)T, and the size of the vectors is determined by the number of phase shifts.

Let us rewrite Eq. (2) as follows:

Eq. (3a)

I·C=A·C+(B·cosϕ)C·C(B·sinϕ)S·C
and

Eq. (3b)

I·S=A·S+(B·cosϕ)C·S(B·sinϕ)S·S.

To extract the quadrature components sinϕ and cosϕ, we can use a property of the dot product for orthogonal vectors (a·a)=0. Let C be the result of the cross product A and C, C=C×A, and S be the result of the cross product A and S, S=S×A. Taking into account that S is orthogonal to A and S, and the mentioned property of the dot product, we obtain

Eq. (4a)

I·C=(B·sinϕ)S·C
and

Eq. (4b)

I·S=(B·cosϕ)C·S;
then

Eq. (5)

Bsinϕ=I·CS·C
and

Eq. (6)

Bcosϕ=I·SC·S.

Considering the properties of the dot product and the cross-product c(b×a)=b(c×a) in the case of the noncyclic permutation of the vectors, we obtain

(S·C)=S(C×A)=C(S×A)=(C·S).
Then, a reconstruction formula can be represented in the vector form as

Eq. (7)

sinϕcosϕ=I·CI·S,ϕ=arctan(I·CI·S),
because in this case only the vector I is calculated. For the case of three shifts,

Eq. (8)

I=I×[111]=M·I=[011101110]·I=[I2I3I3I1I1I2],
where M is a matrix, which calculates the cross-product of the vectors.4 Then, we obtain the equation described in Ref. 5:

Eq. (9)

ϕ=arctan(I2I3)cosδ1+(I3I1)cosδ2+(I1I2)cosδ3(I3I2)sinδ1+(I1I3)sinδ2+(I2I1)sinδ3.

If δ1=π/4, δ2=3π/4, and δ3=5π/4, we obtain the expression described in Ref. 6,

Eq. (10)

ϕ=arctanI3I2I1I2.

The m-dimensional matrix M (m3) can be presented as

Eq. (11)

M=[0100110100010000000110010]m×m.
and for the case of four shifts,

Eq. (12)

I=[0101101001011010]·[I1I2I3I4]=[I2I4I1+I3I2+I4I1I3].

If δ1=0, δ2=π/2, δ3=π, and δ4=3π/2, we obtain the expression found in Ref. 7:

Eq. (13)

ϕ=arctanI4I2I1I3.

Obtaining Eqs. (9), (10), and (13) is not necessary if a generic algorithmic procedure that implements Eq. (7) is used:

Eq. (14)

ϕ=arctani=0m1[(Imod(i+1,m)Imod(m+i1,m))·(δi)]i=0m1[(Imod(i+1,m)Imod(m+i1,m))·sin(δi)],
where mod(i, m), the remainder when i is divided by m.

3.

Conclusions

We have proposed an algorithm that can be used to obtain reconstruction formulas without solving systems of trigonometric equations. Our procedure can be used for any set of arbitrary phase shifts. Note that the proposed method allows corresponding real, but not assumed, values of phase shifts to be used in the reconstruction formulas. This property can undoubtedly help to increase the precision of phase calculation in phase-shifting interferometry.

References

1. 

P. de Groot, “Phase-shift calibration errors in interferometers with spherical Fizeau cavities,” Appl. Opt., 34 (16), 2856 –2863 (1995). http://dx.doi.org/10.1364/AO.34.002856 APOPAI 0003-6935 Google Scholar

2. 

P. de Groot, “101-frame algorithm for phase shifting interferometry,” Proc. SPIE, 3098 283 –292 (1997). http://dx.doi.org/10.1117/12.281171 PSISDG 0277-786X Google Scholar

3. 

P. A. A. Magalhãeset al., “New equations for phase evaluation in measurements with an arbitrary but constant phase shift between captured intensity signs,” Opt. Eng., 48 (11), 113602 (2009). http://dx.doi.org/10.1117/1.3265438 OPEGAR 0091-3286 Google Scholar

4. 

S. LiuG. Trenkler, “Hadamard, Khatri-Rao, Kronecker and other matrix products,” Int. J. Inform. Syst. Sci., 4 (1), 160 –177 (2008). 1708-296X Google Scholar

5. 

P. HariharanB. F. OrebN. Brown, “Digital phase-measurement system for real-time holographic interferometry,” Opt. Commun., 41 (6), 393 –398 (1982). http://dx.doi.org/10.1016/0030-4018(82)90163-8 OPCOB8 0030-4018 Google Scholar

6. 

J. C. WyantK. Creath, “Recent advances in interferometric optical testing,” Laser Focus, 21 (11), 118 –132 (1985). LAFOAK 0740-2511 Google Scholar

7. 

J. C. Wyant, “Interferometric optical metrology: basic system and principles,” Laser Focus, 18 (5), 65 –67 (1982). LAFOAK 0740-2511 Google Scholar
© 2013 Society of Photo-Optical Instrumentation Engineers (SPIE) 0091-3286/2013/$25.00 © 2013 SPIE
Vladimir Guzhov, Sergey Ilinykh, Roman Kuznetsov, and Dmitry Haydukov "Generic algorithm of phase reconstruction in phase-shifting interferometry," Optical Engineering 52(3), 030501 (21 February 2013). https://doi.org/10.1117/1.OE.52.3.030501
Published: 21 February 2013
Lens.org Logo
CITATIONS
Cited by 17 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Reconstruction algorithms

Phase shifts

Phase interferometry

Interferometry

Interferometers

Optical engineering

RELATED CONTENT


Back to Top