Systematic spectral analysis of Phase Shifting Interferometry (PSI) algorithms was first proposed in 1990 by Freischlad and Koliopoulos (F&K). This analysis was proposed with the intention that “in a glance” the main properties of the PSI algorithms would be highlighted. However a major drawback of the F&K spectral analysis is that it changes when the PSI algorithm is rotated or its reference signal is time-shifted. In other words, the F&K spectral plot is different when the PSI algorithm is rotated or its reference is time-shifted. However, it is well known that these simple operations do not alter the basic phase demodulation properties of PSI algorithms, except for an unimportant piston. Here we propose a new way to analyze the spectra of PSI algorithms which is invariant to rotation and/or reference time-shift among other advantages over the nowadays standard PSI spectral analysis by F&K.

©2009 Optical Society of America

Full Article  |  PDF Article
OSA Recommended Articles
The general theory of phase shifting algorithms

M. Servin, J. C. Estrada, and J. A. Quiroga
Opt. Express 17(24) 21867-21881 (2009)

Design of phase-shifting algorithms by fine-tuning spectral shaping

A. Gonzalez, M. Servin, J. C. Estrada, and J. A. Quiroga
Opt. Express 19(11) 10692-10697 (2011)


  • View by:
  • |
  • |
  • |

  1. K. Freischlad and C. L. Koliopoulos, “Fourier description of digital phase-measuring interferometry,” J. Opt. Soc. Am. A 7(4), 542–551 (1990).
  2. D. W. Phillion, “General methods for generating phase-shifting interferometry algorithms,” Appl. Opt. 36(31), 8098–8115 (1997).
  3. Y. Surrel, “Design of algorithms for phase measurements by the use of phase stepping,” Appl. Opt. 35(1), 51–60 (1996).
    [Crossref] [PubMed]
  4. D. Malacara, M. Servin, and Z. Malacara, Interferogram analysis for Optical Testing, 2th ed., (Marcel Deker, 2003).
  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. J. G. Proakis, and D. G. Manolakis, Digital Signal Processing, 4th-ed., (Prentice Hall, 2007).
  7. M. Servin, J. C. Estrada, J. A. Quiroga, J. F. Mosiño, and M. Cywiak, “Noise in phase shifting interferometry,” Opt. Express 17(11), 8789–8794 (2009).
    [Crossref] [PubMed]

2009 (1)

1997 (1)

1996 (1)

1995 (1)

1990 (1)

Creath, K.

Cywiak, M.

Estrada, J. C.

Freischlad, K.

Koliopoulos, C. L.

Mosiño, J. F.

Phillion, D. W.

Quiroga, J. A.

Schmit, J.

Servin, M.

Surrel, Y.

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.

Figures (3)

Fig. 1
Fig. 1 Spectral magnitude |H(ω)| of two PSI algorithms. We show the input interferogram as I(t), and the output complex signal as Ic(t). We can clearly see that these two filters pass only the left side analytical signal at frequency -ω 0, while rejecting the background a(x,y) at ω = 0 and the right side analytical signal at + ω 0.
Fig. 2
Fig. 2 Spectral magnitude of the two 4-step PSI algorithms in Eq. (17). We see that they are actually proportional and as a consequence, have identical phase demodulation properties. These PSI filters let pass the left side analytical signal at ω = −1.0 while rejecting the background at ω = 0 and the complex signal at ω = 1.0.
Fig. 3
Fig. 3 Spectral magnitude |H(ω)| of the 5-step Schwider-Hariharan algorithm and its rotated twin. Both are of course identical. Only the left side analytical signal at ω = −1.0 passes through, while rejecting the background and the complex signal at ω = + 1.0.

Equations (28)

Equations on this page are rendered with MathJax. Learn more.

H(ω0)0,    and      H(ω0)=H(0)=0.
{F[hr(t)]iF[hi(t)]   }|ω=ω0=0,      and      F[hr(t)]|ω=0=F[hi(t)]|ω=0=0.
h(t)exp[iΔ0],    and    h(tt0).
h(tt0)=exp[iω0(tt0)]{   nanδ(tnT)}     =h(t)exp[iω0t0].
tan[φ1(x,y)]=I1+I2I3I4I1I2I3+I4,  and  tan[φ2(x,y)]=I2I4I1I3.
h1(t)=δ(t+3π/4)δ(t+π/4)δ(tπ/4)+δ(t3π/4)+                   i[δ(t+3π/4)+δ(t+π/4)δ(tπ/4)δ(t3π/4)]h2(t)=δ(tπ/2)δ(t+π/2)+i   [δ(t)δ(t+π)].
Hr1(ω)=2cos(πω/4)2cos(3πω/4),Hi1(ω)=2i[sin(πω/4)+sin(3πω/4)]Hr2(ω)=2isin(πω/2),             Hi2(ω)=1exp(iπω).
2   h(t)exp[iπ/4]=δ(tπ)+2δ(tπ/2)2δ(t)2δ(t+π/2)+δ(t+π)+                     i[δ(tπ)2δ(tπ/2)2δ(t)+2δ(t+π/2)+δ(t+π)].