Phase unwrapping of MR phase images using Poisson equation

The authors have developed a technique based on a solution of the Poisson equation to unwrap the phase in magnetic resonance (MR) phase images. The method is based on the assumption that the magnitude of the inter-pixel phase change is less than pi per pixel. Therefore, the authors obtain an estimate of the phase gradient by "wrapping" the gradient of the original phase image. The problem is then to obtain the absolute phase given the estimate of the phase gradient. The least-squares (LS) solution to this problem is shown to be a solution of the Poisson equation allowing the use of fast Poisson solvers. The absolute phase is then obtained by mapping the LS phase to the nearest multiple of 2 K from the measured phase. The proposed technique is evaluated using MR phase images and is proven to be robust in the presence of noise. An application of the proposed method to the 3-point Dixon technique for water and fat separation is demonstrated.

[1]  R. Kuc,et al.  A direct relation between a signal time series and its unwrapped phase , 1982 .

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

[3]  G. Glover,et al.  Encoding strategies for three‐direction phase‐contrast MR imaging of flow , 1991, Journal of magnetic resonance imaging : JMRI.

[4]  D J Bone,et al.  Fourier fringe analysis: the two-dimensional phase unwrapping problem. , 1991, Applied optics.

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

[6]  Thomas S. Huang,et al.  Application of phase unwrapping to image restoration , 1981 .

[7]  Magnus R. Hestenes,et al.  Conjugate Direction Methods in Optimization , 1980 .

[8]  N J Pelc,et al.  Noise reduction in three‐dimensional phase‐contrast MR velocity measurementsl , 1993, Journal of magnetic resonance imaging : JMRI.

[9]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[10]  R. J. Green,et al.  Phase Unwrapping Using A-Priori Knowledge About The Band Limits Of A Function. , 1989, Other Conferences.

[11]  Hiroaki Takajo,et al.  Least-squares phase recovery from the bispectrum phase: an algorithm for a two-dimensional object , 1991 .

[12]  T E Conturo,et al.  Signal‐to‐noise in phase angle reconstruction: Dynamic range extension using phase reference offsets , 1990, Magnetic resonance in medicine.

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

[14]  Richard M. Leahy,et al.  Determining cardiac velocity fields and intraventricular pressure distribution from a sequence of ultrafast CT cardiac images , 1994, IEEE Trans. Medical Imaging.

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

[16]  L. Mertz,et al.  Speckle imaging, photon by photon. , 1979, Applied optics.

[17]  L Axel,et al.  Correction of phase wrapping in magnetic resonance imaging. , 1989, Medical physics.

[18]  C.-C. Jay Kuo,et al.  Discretization and solution of elliptic PDEs-a digital signal processing approach , 1990, Proc. IEEE.

[19]  B. L. Buzbee,et al.  The direct solution of the discrete Poisson equation on irregular regions , 1970 .

[20]  D. Luenberger.  Hyperbolic Pairs in the Method of Conjugate Gradients , 1969 .

[21]  O. Axelsson,et al.  Finite element solution of boundary value problemes - theory and computation , 2001, Classics in applied mathematics.

[22]  A. A.M. Maas,et al.  New Image Processing Algorithms For The Analysis Of Speckle Interference Patterns , 1989, Optics & Photonics.

[23]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[24]  Louis A. Romero,et al.  A Cellular Automata Method for Phase Unwrapping , 1986, Topical Meeting On Signal Recovery and Synthesis II.

[25]  Trevor J. Hall,et al.  Phase Unwrapping Using the Complex Zeros of a Band-limited Function and the Presence of Ambiguities in Two Dimensions , 1984 .

[26]  H. Takajo,et al.  Least-squares phase estimation from the phase difference , 1988 .

[27]  G H Glover,et al.  Three‐point dixon technique for true water/fat decomposition with B0 inhomogeneity correction , 1991, Magnetic resonance in medicine.

[28]  Hamid Al-Nashi.  Phase unwrapping of digital signals , 1989, IEEE Trans. Acoust. Speech Signal Process..

[29]  Bruno O. Shubert,et al.  Random variables and stochastic processes , 1979 .

[30]  José M. F. Moura,et al.  Phase unwrapping of signals propagated under the Arctic ice crust: a statistical approach , 1988, IEEE Trans. Acoust. Speech Signal Process..

[31]  Michael Braun,et al.  Two-dimensional phase unwrapping using a minimum spanning tree algorithm , 1992, IEEE Trans. Image Process..

[32]  D C Ghiglia,et al.  Direct phase estimation from phase differences using fast elliptic partial differential equation solvers. , 1989, Optics letters.

[33]  Mark Hedley,et al.  A new two‐dimensional phase unwrapping algorithm for MRI images , 1992, Magnetic resonance in medicine.

[34]  D. Luenberger.  Optimization by Vector Space Methods , 1968 .

[35]  Henry Stark,et al.  Image recovery: Theory and application , 1987 .

[36]  José Tribolet,et al.  A new phase unwrapping algorithm , 1977 .