21 February 2013 Generic algorithm of phase reconstruction in phase-shifting interferometry
Author Affiliations +
Optical Engineering, 52(3), 030501 (2013). doi:10.1117/1.OE.52.3.030501
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.
Guzhov, Ilinykh, Kuznetsov, and Haydukov: Generic algorithm of phase reconstruction in phase-shifting interferometry



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.


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


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


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:





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









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

Then, a reconstruction formula can be represented in the vector form as


because in this case only the vector I is calculated. For the case of three shifts,


where M is a matrix, which calculates the cross-product of the vectors.4 Then, we obtain the equation described in Ref. 5:



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



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


and for the case of four shifts,



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



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


where mod(i, m), the remainder when i is divided by m.



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.



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


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


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).OPEGAR0091-3286http://dx.doi.org/10.1117/1.3265438Google Scholar


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


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


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


J. C. Wyant, “Interferometric optical metrology: basic system and principles,” Laser Focus 18(5), 65–67 (1982).LAFOAK0740-2511Google Scholar

Vladimir Guzhov, Sergey Ilinykh, Roman Kuznetsov, Dmitry Haydukov, "Generic algorithm of phase reconstruction in phase-shifting interferometry," Optical Engineering 52(3), 030501 (21 February 2013). http://dx.doi.org/10.1117/1.OE.52.3.030501
Submission: Received ; Accepted

Reconstruction algorithms

Phase shifts

Phase interferometry


Optical engineering



Back to Top