Robust 2D phase unwrapping algorithm based on the transport of intensity equation

2D phase-unwrapping algorithms (PUAs) are commonly used to obtain a continuous phase map from the sawtooth-shaped phase map. However, implementing PUAs can be time consuming, and the accuracy of those algorithms may be low if there is heavy noise. In this paper, we develop a simple and robust PUA based on the transport of intensity equation (TIE). In our method, the TIE was solved using the fast cosine transform, and a phase correction operation was introduced after the TIE was solved. Because of the phase correction operation, the proposed method can obtain a satisfactory unwrapping result even in a notably hash noise condition. The simulation and experimental results are presented to validate the effectiveness of the proposed method. The detailed software package can be found at https://ww2.mathworks.cn/matlabcentral/fileexchange/68493-robust-2d-phase-unwrapping-algorithm.

[1]  M Servin,et al.  Phase unwrapping with a regularized phase-tracking system. , 1998, Applied optics.

[2]  Zhi-Pei Liang,et al.  A model-based method for phase unwrapping , 1996, IEEE Trans. Medical Imaging.

[3]  C. Werner,et al.  Satellite radar interferometry: Two-dimensional phase unwrapping , 1988 .

[4]  Lu Zhang,et al.  2D phase unwrapping algorithm for interferometric applications based on derivative Zernike polynomial fitting technique , 2014 .

[5]  Justo Arines Least-squares modal estimation of wrapped phases: application to phase unwrapping. , 2003, Applied optics.

[6]  M. Teague Deterministic phase retrieval: a Green’s function solution , 1983 .

[7]  J A Quiroga,et al.  Phase-unwrapping algorithm based on an adaptive criterion. , 1995, Applied optics.

[8]  Tomasz Kozacki,et al.  Fast and accurate phase-unwrapping algorithm based on the transport of intensity equation. , 2017, Applied optics.

[9]  J A Quiroga,et al.  Phase-unwrapping algorithm for noisy phase-map processing. , 1994, Applied optics.

[10]  Shi Yuqing Robust phase unwrapping by spinning iteration. , 2007, Optics express.

[11]  N. Streibl Phase imaging by the transport equation of intensity , 1984 .

[12]  T. Kozacki,et al.  Computation of diffracted fields for the case of high numerical aperture using the angular spectrum method. , 2012, Applied optics.

[13]  Wei Xu,et al.  A region-growing algorithm for InSAR phase unwrapping , 1999, IEEE Trans. Geosci. Remote. Sens..

[14]  K. Khare,et al.  Two-dimensional phase unwrapping using the transport of intensity equation. , 2016, Applied optics.

[15]  Javier Vargas,et al.  Generalization of the Principal Component Analysis algorithm for interferometry , 2013 .

[16]  Louis A. Romero,et al.  Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods , 1994 .

[17]  K. Nugent,et al.  Partially coherent fields, the transport-of-intensity equation, and phase uniqueness , 1995 .

[18]  J. M. Huntley Noise-immune phase unwrapping algorithm. , 1989, Applied optics.

[19]  A Baldi Two-dimensional phase unwrapping by quad-tree decomposition. , 2001, Applied optics.