Markov random field based phase demodulation of interferometric images

We present a novel method to solve the sign ambiguity for phase demodulation from a single interferometric image that possibly contains closed fringes. The problem is formulated in a Markov random field (MRF) energy minimization framework with the assumption of phase gradient orientation continuity. The binary pairwise objective function is non-submodular and therefore its minimization is an NP-hard problem, for which we devise a multigrid hierarchy of quadratic pseudoboolean optimization problems that can be improved iteratively to approximate the optimal solution. We name the method MSARI algorithm, for Markov based sign ambiguity resolution in interferometry. Compared with traditional path-following phase demodulation methods, the new approach does not require any heuristic scanning strategy, is not subject to the propagation of error, and the extension to three dimensional fringe patterns is straightforward. A set of experiments with synthetic data and real prelens tear film interferometric images of the human eye demonstrate the effectiveness and robustness of the proposed algorithm as compared with existing state-of-the-art phase demodulation methods.

[1]  Mariano Rivera,et al.  Fast half-quadratic regularized phase tracking for nonnormalized fringe patterns. , 2006, Journal of the Optical Society of America. A, Optics, image science, and vision.

[2]  Vladimir Kolmogorov,et al.  Minimizing Nonsubmodular Functions with Graph Cuts-A Review , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[3]  Manuel Servin,et al.  Adaptive quadrature filters and the recovery of phase from fringe pattern images , 1997 .

[4]  M. Servin,et al.  Fringe-follower regularized phase tracker for demodulation of closed-fringe interferograms , 2001 .

[5]  K Itoh,et al.  Analysis of the phase unwrapping algorithm. , 1982, Applied optics.

[6]  Kenneth H. Womack,et al.  Interferometric Phase Measurement Using Spatial Synchronous Detection , 1983 .

[7]  Pierre Hansen,et al.  Roof duality, complementation and persistency in quadratic 0–1 optimization , 1984, Math. Program..

[8]  Cho Jui Tay,et al.  A new method for phase extraction from a single fringe pattern , 2004 .

[9]  Manuel Servin,et al.  Local phase from local orientation by solution of a sequence of linear systems , 1998 .

[10]  José M. N. Leitão,et al.  The ZπM algorithm: a method for interferometric image reconstruction in SAR/SAS , 2002, IEEE Trans. Image Process..

[11]  T. Kreis Holographic Interferometry: Principles and Methods , 1996 .

[12]  P. King-Smith,et al.  Interferometric imaging of the full thickness of the precorneal tear film. , 2006, Journal of the Optical Society of America. A, Optics, image science, and vision.

[13]  Qian Kemao,et al.  Frequency guided methods for demodulation of a single fringe pattern. , 2009, Optics express.

[14]  Jesús Villa,et al.  Phase recovery from a single fringe pattern using an orientational vector-field-regularized estimator. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[15]  Daniel Malacara Handbook of Optical Engineering , 2001 .

[16]  Vladimir Kolmogorov,et al.  Optimizing Binary MRFs via Extended Roof Duality , 2007, 2007 IEEE Conference on Computer Vision and Pattern Recognition.

[17]  M. Servin,et al.  Isotropic n-dimensional fringe pattern normalization , 2003 .

[18]  K. Larkin Natural demodulation of two-dimensional fringe patterns. II. Stationary phase analysis of the spiral phase quadrature transform. , 2001, Journal of the Optical Society of America. A, Optics, image science, and vision.

[19]  Vladimir Kolmogorov,et al.  Convergent Tree-Reweighted Message Passing for Energy Minimization , 2006, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[20]  José M. Bioucas-Dias,et al.  Phase Unwrapping via Graph Cuts , 2007, IEEE Trans. Image Process..

[21]  Richard Szeliski,et al.  A Comparative Study of Energy Minimization Methods for Markov Random Fields with Smoothness-Based Priors , 2008, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[22]  M Rivera,et al.  Multigrid algorithms for processing fringe-pattern images. , 1998, Applied optics.

[23]  Daniel P. Huttenlocher,et al.  Efficient Belief Propagation for Early Vision , 2004, Proceedings of the 2004 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2004. CVPR 2004..

[24]  Manuel Servin,et al.  Regularized quadrature and phase tracking from a single closed-fringe interferogram. , 2004, Journal of the Optical Society of America. A, Optics, image science, and vision.

[25]  Manuel Servin,et al.  Two-dimensional Phase Locked Loop Demodulation of Interferograms , 1993 .

[26]  M. Doane,et al.  An instrument for in vivo tear film interferometry. , 1989, Optometry and vision science : official publication of the American Academy of Optometry.

[27]  D. J. Brangaccio,et al.  Digital wavefront measuring interferometer for testing optical surfaces and lenses. , 1974, Applied optics.

[28]  Andrew Blake,et al.  Digital tapestry [automatic image synthesis] , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[29]  Olga Veksler,et al.  Fast Approximate Energy Minimization via Graph Cuts , 2001, IEEE Trans. Pattern Anal. Mach. Intell..

[30]  Mariano Rivera,et al.  Fast phase recovery from a single closed-fringe pattern. , 2008, Journal of the Optical Society of America. A, Optics, image science, and vision.

[31]  Jerry L. Prince,et al.  Snakes, shapes, and gradient vector flow , 1998, IEEE Trans. Image Process..

[32]  M. Takeda,et al.  Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry , 1982 .

[33]  Manuel Servin,et al.  Modulo 2pi fringe orientation angle estimation by phase unwrapping with a regularized phase tracking algorithm. , 2002, Journal of the Optical Society of America. A, Optics, image science, and vision.