About this Journal Submit a Manuscript Table of Contents
ISRN Optics
Volume 2012 (2012), Article ID 396960, 6 pages
http://dx.doi.org/10.5402/2012/396960
Research Article

Single Frame Fringe Pattern Analysis for Phase Recovery with Analytic Signal

Laboratoire Instrumentation de Mesure et de Contrôle, Equipe Metrologie Optique et Traitement Image, Université Chouaib Doukkali, Faculté des Sciences, BP 20, El Jadida 24000, Morocco

Received 25 March 2012; Accepted 9 May 2012

Academic Editors: O. Frazao and D. Poitras

Copyright © 2012 Khalid Assid et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We consider a new application of the normalized Hilbert-Huang transform to extract directly the phase from a single fringe pattern. We present a technique to provide, with good accuracy, the phase distribution from a single interferogram without unwrapping step and this by a new exploitation of the analytic signal corresponding to each intrinsic mode function, resulting from one-dimensional empirical mode decomposition of the fringe pattern. A theoretical analysis was carried out for this technique, followed by computer simulations and a real experimental fringe pattern for verification.

1. Introduction

In optical techniques, several physical magnitudes in a variety of scientific and engineering fields are encoded as the phase of a periodic intensity profile. The developments of more sophisticated phase evaluation algorithms are continuously needed [1, 2].

The use of a single fringe pattern for phase retrieval is more suitable in dynamic processes. The unwrapping step, which provides, a continuous phase distribution over its definition domain implies several difficulties and sophisticated algorithms are then needed [3].

Recently, the normalized Hilbert-Huang transform (NHHT) has emerged as a promising and powerful tool for processing nonlinear and nonstationary signals. It consists of three parts: empirical mode decomposition (EMD), intrinsic mode functions (IMF) normalization, and Hilbert spectral analysis (HSA). The NHHT was developed by Huang et al. [47] to overcome the limitations set by the Bedrosian and Nuttall theorems when applying the regular Hilbert transform [812].

In our work, we attempt to retrieve the unwrapped phase from a single image. Our fringe analysis technique consists of decomposing the fringe pattern into IMFs, applying the HSA to each normalized IMF, then constructing individual analytic signals and computing the phase gradient. This leads directly to the phase distribution avoiding the complex step of phase unwrapping.

Computer simulation results obtained using our method are in very good agreement with those produced using phase-shifting methods [13]. To validate the performance of the proposed technique, it has been further tested on experimental fringe patterns. It was shown that the result totally agrees with that retrieved by phase-shifting techniques.

2. Empirical Mode Decomposition for Fringe Pattern Analysis

The fringe pattern, derived from a two-beam interferometer is characterized by the sinusoidal dependence of the intensity on the spatial coordinates (𝑥,𝑦) of the image plane [1, 2]: 𝐼(𝑥,𝑦)=𝑎(𝑥,𝑦)+𝑏(𝑥,𝑦)cos𝜙(𝑥,𝑦),(1) where 𝑎(𝑥,𝑦) is the background. The terms 𝑏(𝑥,𝑦) and cos𝜙(𝑥,𝑦) are the modulation intensity and the optical phase.

The empirical mode decomposition (EMD) analysis is an adaptive and fully driven data method to decompose any signal into a set of components called intrinsic mode functions (IMF) and a residue [5]. Adding all the IMFs together with the residue reconstructs the original signal without information loss or distortion [14]: 𝐼(𝑥)=𝑛1𝑘=1IMF𝑘(𝑥)+Res(𝑥).(2) Adding all the IMFs by just skipping the DC component that is, Res(𝑥), we remove the background. The filtered fringe pattern is 𝐽(𝑥)=𝑏(𝑥)cos𝜙(𝑥)=𝑛1𝑘=1IMF𝑘(𝑥).(3) When decomposing the intensity by EMD, the resulting IMFs are not restricted to narrow band signals and can be both amplitude and frequency modulated [8]. Each IMF should be representable in the form IMF𝑘(𝑥)=𝑏𝑘𝜙(𝑥)cos𝑘(𝑥),(4) where the amplitude 𝑏𝑘(𝑥) and phase 𝜙𝑘 are both physically and mathematically meaningful. Once a suitable polar parameterization is determined, it is possible to analyse 𝐽(𝑥) by processing these individual components. The most common procedure to determine a polar representation is the analytic method of signal processing [15].

3. Analytic Signal Construction

To construct the analytic signal corresponding to each IMF, the simplest way is by employing the Hilbert transform [6]. The real part of an analytic signal is the original input, whereas its imaginary component is the Hilbert transform HT of that signal.

For each IMF𝑘(𝑥), the corresponding analytic signal is 𝑍𝑘(𝑥)=IMF𝑘(𝑥)+𝑖HTIMF𝑘(𝑥).(5) Using the exponential notation 𝑍𝑘(𝑥)=𝑏𝑘(𝑥)𝑒𝑖𝜙𝑘(𝑥)(6) the following instantaneous attributes of 𝑍𝑘(𝑥) can be defined: 𝑏𝑘(𝑥)=IMF𝑘2(𝑥)+HTIMF𝑘(𝑥)21/2,(7)IMFN𝑘(𝑥)=IMF𝑘(𝑥)𝑏𝑘𝜙(𝑥)(8)𝑘𝑏(𝑥)=artctan𝑘HTIMFN𝑘(𝑥)IMF𝑘𝑓(𝑥)(9)𝑘1(𝑥)=2𝜋𝑑𝜙𝑘(𝑥),𝑑𝑥(10) where 𝑏𝑘(𝑥) is the instantaneous amplitude, IMFN𝑘(𝑥) is the normalized IMF𝑘(𝑥), 𝜙𝑘(𝑥) is the instantaneous phase, and 𝑓𝑘(𝑥) is the instantaneous frequency [16].

In (9), we chose to conduct Hilbert transform on IMFN𝑘(𝑥) rather than on IMF𝑘(𝑥), and this to avoid the Bedrosian and Nuttall restrictions [912].

The expansion of the signal in terms of the functions given in (6) is 𝑍(𝑥)=𝑛1𝑘=1𝑍𝑘(𝑥)=𝑛1𝑘=1𝑏𝑘(𝑥)𝑒𝑖𝜙𝑘(𝑥),𝑏(11)(𝑥)𝑒𝑖𝜙(𝑥)=𝑛1𝑘=1𝑏𝑘(𝑥)𝑒𝑖𝜙𝑘(𝑥).(12) Typically, the modulation term 𝑏(𝑥) is a slowly and smoothly varying function compared to 𝜙(𝑥), The gradient of (12) leads to: 𝑖𝑏(𝑥)𝜕𝜙(𝑥)𝑒𝜕𝑥𝑖𝜙(𝑥)=𝑛1𝑘=1𝜕𝑏𝑘(𝑥)𝜕𝑥+𝑖2𝜋𝑓𝑘(𝑥)𝑏𝑘𝑒(𝑥)𝑖𝜙𝑘(𝑥).(13) Dividing (13) by (12) gives then the phase gradient by the real part of the complex expression in 𝜕𝜙(𝑥)||||||𝜕𝑥=𝑛1𝑘=12𝜋𝑓𝑘(𝑥)𝑏𝑘(𝑥)𝑖𝜕𝑏𝑘𝑒(𝑥)/𝜕𝑥𝑖𝜙𝑘(𝑥)𝑛1𝑘=1𝑏𝑘(𝑥)𝑒𝑖𝜙𝑘(𝑥)||||||(14) which leads to the phase by integration of the gradient.

As we can see, the IMF𝑘 contribution to the phase gradient is ||||||2𝜋𝑓𝑘(𝑥)𝑏𝑘(𝑥)𝑖𝜕𝑏𝑘𝑒(𝑥)/𝜕𝑥𝑖𝜙𝑘(𝑥)𝑛1𝑘=1𝑏𝑘(𝑥)𝑒𝑖𝜙𝑘(𝑥)||||||.(15)

4. Numerical Simulations

The numerical simulation consists in generating digitally fringe patterns to verify the ability of the method to determine the phase distribution. The test phase function we used is 𝜙(𝑥,𝑦)=0.0009(𝑥128)2+(𝑦128)2.(16) The intensity distribution of fringe pattern shown in Figure 1 is 𝐼(𝑥,𝑦)=1+cos(𝜙(𝑥,𝑦)).(17) We illustrated in Figure 2 the IMF is extracted by EMD from a data row corresponding to 𝑦=50.

396960.fig.001
Figure 1: The fringe pattern I.
396960.fig.002
Figure 2: The tree IMFs component resulting from EMD decomposition of the data row 𝐽(𝑦=50).

As we can see in Figure 3, the physical significance of each IMF is proven since we can see their individual contributions to the global phase phaset.

396960.fig.003
Figure 3: The individual contribution of each IMF to construct the global phase phaset.

For determining the phase distribution of two dimensional fringe pattern, the previous procedure is generalized to all rows of the image.

The difference map between the simulated phase in Figure 4(a) and the retrieved phase in Figure 4(b) is given in Figure 4(c) where we can see that the error is small in general. Even over the large error location, the error values are smaller than 5%, except for the end region, where the end effect of the Hilbert transform causes additional problems. The occasionally large error values offer an indication where the method failed. All such locations occur at the minimum amplitude. Overall, the NHHT method gives a more consistent, stable quadrature to construct the analytic signal.

fig4
Figure 4: (a) Phase gradient in 𝑥 direction retrieved by (14), (b) phase distribution of (16), and (c) the difference map between the simulated and the retrieved phase.

5. Experimental Results

To illustrate the use of our method for real applications, we tested its performance on a rough ground finish, aluminium surface from a hard disk drive assembly. Changes were induced by tilting the surface and recording interferograms before and after. Figure 5 shows a fringe pattern with 9 fringes of tilt.

396960.fig.005
Figure 5: Fringe pattern for an aluminium transducer suspension assembly from a hard disk drive that was tilted between measurements.

The 1D empirical mode decomposition method is also applied successively to reduce speckle noise in the correlation fringes produced in digital speckle pattern interferometry (DSPI) [1719]. When decomposing a data row of the fringe pattern, it is clearly shown in Figure 6 that the first and second IMF contain only speckle noise and must be omitted to construct a smooth fringe pattern.

396960.fig.006
Figure 6: The four IMFs resulting from EMD decomposition of the data row 𝐽(𝑦=110).

The results are of the same order of accuracy to DSPI measurements using the Carré algorithm [20] as it is shown in Figures 7(a) and 7(b). It is worth noting that our method gives the phase distribution with a single image and without unwrapping, while most phase evaluation methods require unwrapping phase over its definition domain, in the manner to provide a continuous phase distribution.

fig7
Figure 7: (a) Phase retrieved using our method, and (b) phase retrieved with Carré algorithm.

6. Conclusion

In this paper, we have presented and tested a new method of phase extraction from a single fringe pattern. This was achieved by a new exploitation of the analytic signal corresponding to each normalized intrinsic mode function. Moreover, we gave a physical sense of IMFS derived from fringe pattern decomposition by EMD. We used the NHHT since it has helped to overcome many of the difficulties of the regular Hilbert transform, it gives consistent and stable quadrature when constructing analytic signals. The method applied to simulated and experimental fringe patterns is compared to phase-shifting methods. It shows great robustness and a very close agreement in the computation of the phase.

In summary, the advantages of this technique over other methods include the following.(1) This method is adequate for real-time acquisition since it provides the optical phase distribution from only one single fringe pattern with a high accuracy. (2) It is possible to perform an accurate phase demodulation on the fringe pattern even when the image is noisy.(3) This technique requires less processing time in computation and is easy to implement.(4) This technique leads directly to the phase distribution avoiding the complex step of phase unwrapping.(5) The NHHT applied to fringe pattern analysis has defined the qualitative and quantitative contribution of IMFs to recover the phase. This way of fringe pattern analysis is as far as we know has not yet been reported.

Acknowledgments

The authors want to thank Dr. Joseph W. Goodman for his constructive and helpful comments. they also want to thank Dr. Neal Brock and Dr. J. C. Wyant from 4D Technology for providing them the experimental shifted fringe patterns.

References

  1. S. Sirohi and S. Chau, Optical Methods of Measurement, Marcel Dekker, New york, NY, USA, 1999.
  2. B. V. Dorrío and J. L. Fernández, “Phase-evaluation methods in whole-field optical measurement techniques,” Measurement Science and Technology, vol. 10, no. 3, pp. R33–R55, 1999. View at Publisher · View at Google Scholar · View at Scopus
  3. D. Ghiglia and M. Pritt, Two-Dimensional Phase Unwrapping Theory, Algorithms and Applications, John Wiley & Sons, New York, NY, USA, 1998.
  4. N. E. Huang, Z. Shen, and S. R. Long, “A new view of nonlinear water waves: the Hilbert spectrum,” Annual Review of Fluid Mechanics, vol. 31, pp. 417–457, 1999. View at Scopus
  5. N. E. Huang, Z. Shen, S. R. Long et al., “The empirical mode decomposition and the Hubert spectrum for nonlinear and non-stationary time series analysis,” Proceedings of the Royal Society, vol. 454, no. 1971, pp. 903–995, 1998. View at Scopus
  6. N. E. Huang and Z. Wu, “A review on Hilbert-Huang transform: method and its applications to geophysical studies,” Reviews of Geophysics, vol. 46, no. 2, Article ID RG2006, 2008. View at Publisher · View at Google Scholar · View at Scopus
  7. N. E. Huang and S. R. Long, A generalized Zero Crossing for Local Frequency Determination, United States Patent pending, 2003.
  8. N. E. Huang, Z. Wu, S. R. Long, K. C. Arnold, X. Chen, and K. Blank, “On instantaneous frequency,” Advances in Adaptive Data Analysis, vol. 1, no. 2, pp. 177–229, 2009. View at Publisher · View at Google Scholar · View at Scopus
  9. E. Bedrosian, “A product theorem for Hilbert transform,” Proceedings of the IEEE, vol. 51, no. 2, pp. 868–869, 1972. View at Scopus
  10. A. H. Nuttall, “On the quadrature approximation to the Hilbert transform of modulated signals,” Proceeding of IEEE, vol. 54, pp. 1458–1459, 1966.
  11. N. E. Huang and M. Bethesda, “Computing instantaneous frequency by normalizing Hilbert transform,” 2005, US Patent 6901353.
  12. K. Assid, V. Dembele, F. alaoui, and A. Nassim, “Normalized Hilbert Huang transform-NHHT applied to phase extraction in wavelet domain,” The Open Optics Journal, vol. 6, pp. 9–13, 2012.
  13. K. Creath, “V phase-measurement interferometry techniques,” Progress in Optics, vol. 26, no. C, pp. 349–393, 1988. View at Publisher · View at Google Scholar · View at Scopus
  14. G. Rilling, P. Flandrin, and P. Goncalves, “On empirical mode decomposition and its algorithms,” in Proceedings of IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing, 2003.
  15. R. C. Sharpley and V. Vatchev, “Analysis of the intrinsic mode functions,” Journal of Constructive Approximation, vol. 24, no. 1, pp. 17–47, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. B. Boashash, “Estimating and interpreting the instantaneous frequency of a signal,” Proceeding of IEEE, vol. 80, pp. 520–583, 1992.
  17. M. B. Bernini, G. E. Galizzi, A. Federico, and G. H. Kaufmann, “Evaluation of the 1D empirical mode decomposition method to smooth digital speckle pattern interferometry fringes,” Optics and Lasers in Engineering, vol. 45, no. 6, pp. 723–729, 2007. View at Publisher · View at Google Scholar · View at Scopus
  18. K. Assid, V. Dembele, F. alaoui, and A. Nassim, “Bidimensional empirical mode decomposition BEMD applied to speckle denoising for wavelet phase evaluation,” Physical Chemical News, vol. 61, pp. 17–23, 2011.
  19. P. Flandrin, G. Rilling, and P. Gonçalvés, “Empirical mode decomposition as a filter bank,” IEEE Signal Processing Letters, vol. 11, no. 2, pp. 112–114, 2004. View at Publisher · View at Google Scholar · View at Scopus
  20. P. Carré, “Installation et utilisation du comparateur photoélectrique et interférentiel du Bureau International des Poids et Mesures,” Metrologia, vol. 2, no. 1, article 005, pp. 13–23, 1966. View at Publisher · View at Google Scholar · View at Scopus