摘要
磁共振成像得到的信号是复数形式,包含幅度和相位。常规诊断中,往往用到的是磁共振的幅度图像。然而,其相位也包含了大量信息,例如自旋原子核的移动速度、磁场的不均匀度和磁化率变化等。因此磁共振相位可以用来估计主磁场均匀性和获取临床相关的生理参数。
但从复数信号中提取真实相位时,相位值却会被限制在(-π,π)弧度区间内,位于该区间外的真实相位被缠绕到这一区间内。此现象称为相位缠绕,得到的相位称为缠绕相位。从缠绕相位恢复真实相位的过程就叫做相位解缠。噪声、欠采样和物体不连续的存在使相位解缠变得困难。
本文提出了三种新的二维相位解缠方法。模拟数据和实际磁共振相位数据被用来评估这些方法的表现。
第一种方法是基于离散粒子群优化算法的枝切线法。这种方法先将整幅图像的残差分成几组;在每组内使用离散粒子群优化算法对正负极性残差进行配对;用枝切线连接每组内配好对的正负极性残差;最后绕过这些枝切线进行相位解缠。
与最新的基于人工智能的枝切法对比,这种方法能在较短时间内得到合理的枝切线连接,且通过对残差分组进一步降低了计算时间。
第二种相位解缠方法是基于直接求解法的加权最小Lp范数法。它将整个相位图像的解缠相位梯度与缠绕相位梯度之间差值的加权Lp范数作为优化目标函数;将这个目标函数转化成一个方程组,其系数矩阵采用稀疏结构储存和表达;最后使用直接求解法求解方程组。由于方程组的系数矩阵与解缠相位有关,因此采取迭代方式得到最终的解缠结果。
与一些常用方法相比,这种方法能有效减少计算时间且解缠结果更准确。
第三种方法是基于掩码的区域增长法。这种方法采用一种新的掩码提取方式将残差合理地连接起来作为掩码中的零点;将掩码与相位导数方差结合构成最终的质量图,这样连接残差经过的点均被当成零质量(也就是质量最差)的点,会被滞留到最后才被相位解缠;接着根据质量图将整幅图像分成多个区域,在每个区域内单独进行相位解缠,其中质量最差的那个区域从多个方向进行相位加权平均;最后将多个区域融合在一起。
与最新的区域增长相位解缠方法(PHUN)相比,这种方法运算速度快并能够限制误差的传播。
In magnetic resonance imaging (MRI), the complex signal contains both the magnitude and phase parts. Usually the magnitude of the MRI signal has been mainly considered. However, the phase of MRI signal offers very important information on the velocity of the moving spins, the main Bo field inhomogeneity, the magnetic susceptibility variations, etc. So the phase can be used to estimate the main B0field inhomogeneity and obtain clinically relevant physiological parameters.
when extracting the phase from a measured complex MR dataset through some mathematical operation, the result is typically wrapped into the principal interval of (-π,π] radians, producing the wrapped phase. The process of estimating the true phase from the wrapped phase is called phase unwrapping. Because of the presence of the noise, undersampling and/or object discontinuities, phase unwrapping becomes intractable and nontrivial. In the literature, there are quite a few existing phase unwrapping algorithms.
In this thesis, three two-dimensional (2D) phase unwrapping methods were proposed. Both simulated and MR data were used to evaluate these algorithms' performances.
(1) A new branch-cut method based on discrete particle swarm optimization (dPSO) algorithm was proposed to solve the phase unwrapping problem of MR data. In this method, all the residues were first grouped by dividing the phase image into sub-regions. Then dPSO was performed region by region to match the opposite polarity residues which were connected by branch cuts afterward. Finally, flood-fill method was used to unwrap phases avoiding these branch cuts. Compared with conventionally used branch-cut phase unwrapping algorithms, the dPSO algorithm is rather robust and effective.
(2) A direct-solver-based weighted minimum Lp-norm algorithm was proposed for MRI phase unwrapping. First, the algorithm converted the weighted minimum-Lp-norm objective function for phase unwrapping into a linear system of equations whose system (coefficient) matrix was a large, symmetric one. Then, the coefficient-matrix was represented in the sparse structure. Finally, standard direct solvers were employed to solve this linear system. The results demonstrate that the proposed algorithm is reliable and robust.
(3) A mask-based region growing approach was proposed for MRI phase unwrapping. The residues were first connected by zero-value pixels that were then served as mask. And the mask and the original quality map were combined into a new quality map. Guided by the new quality map, unwrapping was carried out within multiple regions. At last, the regions were merged by adjusting the offset between one another. The results show the advantage of the proposed approach over the recent region-growing phase unwrapping method (PHUN) in accuracy.
引文
[1]Markl M. Velocity Encoding and Flow Imaging. http://ee-classes.usc.edu/ee591/library/Markl-FlowImaging.pdf; 2005.
[2]Chavez S, Xiang QS, An L. Understanding phase maps in MRI:A new outline phase unwrapping method. IEEE Transactions on Medical Imaging 2002;21(8):966-977.
[3]Young IR, Bydder GM. Phase imaging. In:Stark DD, Bradley WG, editors. Magnetic Resonance Imaging.2nd ed. St. Louis, MO.:Mosby Year Book; 1992.
[4]Liu W, Tang X, Ma Y, Gao J-H.3D phase unwrapping using global expected phase as a reference:Application to MRI global shimming. Magnetic Resonance in Medicine 2013;70(1):160-168.
[5]Maudsley A, Simon H, Hilal S. Magnetic field measurement by NMR imaging. Journal of Physics E:Scientific Instruments 1984;17(3):216-220.
[6]Jezzard P, Balaban RS. Correction for geometric distortion in echo planar images from BO field variations. Magnetic Resonance in Medicine 1995;34(1):65-73.
[7]Cusack R, Papadakis N. New robust 3-D phase unwrapping algorithms: application to magnetic field mapping and undistorting echoplanar images. Neuroimage 2002; 16(3 Pt 1):754-764.
[8]Robinson S, Jovicich J. B0 mapping with multi-channel RF coils at high field. Magnetic Resonance in Medicine 2011;66(4):976-988.
[9]Wharton S, Schafer A, Bowtell R. Susceptibility mapping in the human brain using threshold - based k - space division. Magnetic Resonance in Medicine 2010;63(5):1292-1304.
[10]Schweser F, Deistung A, Lehr BW, Reichenbach JR. Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase:an approach to in vivo brain iron metabolism? Neuroimage 2011;54(4):2789-2807.
[11]Liu T, Spincemaille P, de Rochefort L, Kressler B, Wang Y. Calculation of susceptibility through multiple orientation sampling (COSMOS):a method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in MRI. Magnetic Resonance in Medicine 2009;61(1):196-204.
[12]Reichenbach JR, Venkatesan R, Schillinger DJ, Kido DK, Haacke EM. Small vessels in the human brain:MR venography with deoxyhemoglobin as an intrinsic contrast agent. Radiology 1997;204(1):272-277.
[13]Bagher-Ebadian H, Jiang Q, Ewing JR. A modified fourier-based phase unwrapping algorithm with an application to MRI venography. Journal of Magnetic Resonance Imaging 2008;27(3):649-652.
[14]Reichenbach JR, Essig M, Haacke EM, et al. High-resolution venography of the brain using magnetic resonance imaging. Magnetic Resonance Materials in Physics, Biology and Medicine 1998;6(1):62-69.
[15]Rydell J, Knutsson H, Pettersson J, et al. Phase sensitive reconstruction for water/fat separation in MR imaging using inverse gradient. Lecture Notes in Computer Science 2007;4791:210-218.
[16]Haacke EM, Xu Y, Cheng Y-CN, Reichenbach JR. Susceptibility weighted imaging (SWI). Magnetic Resonance in Medicine 2004;52(3):612-618.
[17]Gilbert G, Savard G, Bard C, Beaudoin G. Quantitative comparison between a multiecho sequence and a single-echo sequence for susceptibility-weighted phase imaging. Magnetic Resonance Imaging 2012;30(5):722-730.
[18]Chiang P, Cai Y, Mak KH, Zheng J. A B-spline approach to phase unwrapping in tagged cardiac MRI for motion tracking. Magnetic Resonance in Medicine 2013;69(5):1297-1309.
[19]Huifang W, John BW, Irina IP, Marvin MD, Keith DP. A three-dimensional quality-guided phase unwrapping method for MR elastography. Physics in Medicine and Biology 2011;56(13):3935.
[20]Nayler G, Firmin D, Longmore D. Blood flow imaging by cine magnetic resonance. Journal of computer assisted tomography 1986; 10(5):715-722.
[21]Szumowski J, Coshow WR, Li F, Quinn SF. Phase unwrapping in the three-point Dixon method for fat suppression MR imaging. Radiology 1994;192(2):555-561.
[22]Skinner TE, Glover GH. An extended two - point dixon algorithm for calculating separate water, fat, and BO images. Magnetic Resonance in Medicine 1997;37(4):628-630.
[23]刘亚涛,俎栋林,包尚联.水、脂分离磁共振成像Dixon方法.中国医学物理学杂志2012;29(6):3759-3768.
[24]Moon-Ho Song S, Napel S, Pelc NJ, Glover GH. Phase unwrapping of MR phase images using Poisson equation. Image Processing, IEEE Transactions on 1995;4(5):667-676.
[25]Glover G, Schneider E. Three - point dixon technique for true water/fat decomposition with BO inhomogeneity correction. Magnetic Resonance in Medicine 1991;18(2):371-383.
[26]Yang GZ, Burger P, Kilner PJ, Karwatowski SP, Firmin DN. Dynamic range extension of cine velocity measurements using motion registered spatiotemporal phase unwrapping. Journal of Magnetic Resonance Imaging 1996;6(3):495-502.
[27]Xiang Q-S. Temporal phase unwrapping for cine velocity imaging. Journal of Magnetic Resonance Imaging 1995;5(5):529-534.
[28]Andris P, Frollo I. Simple and accurate unwrapping phase of MR data. Measurement 2009;42(5):737-741.
[29]Robinson S, Grabner G, Witoszynskyj S, Trattnig S. Combining phase images from multi-channel RF coils using 3D phase offset maps derived from a dual-echo scan. Magnetic Resonance in Medicine 2011;65(6):1638-1648.
[30]Feng W, Neelavalli J, Haacke EM. Catalytic multiecho phase unwrapping scheme (CAMPUS) in multiecho gradient echo imaging:Removing phase wraps on a voxel-by-voxel basis. Magnetic Resonance in Medicine 2013;70(1):117-126.
[31]Robinson S, Schodl H, Trattnig S. A method for unwrapping highly wrapped multi-echo phase images at very high field:UMPIRE. Magnetic Resonance in Medicine 2013.
[32]Dagher J, Reese T, Bilgin A. High-resolution, large dynamic range field map estimation. Magnetic Resonance in Medicine 2014;71(1):105-117.
[33]Haacke EM, Reichenbach JR. Susceptibility weighted imaging in MRI:basic concepts and clinical applications:John Wiley & Sons; 2011.
[34]Xu B, Liu T, Spincemaille P, Prince M, Wang Y. Flow compensated quantitative susceptibility mapping for venous oxygenation imaging. Magnetic Resonance in Medicine 2013.
[35]Zebker HA, Lu Y. Phase unwrapping algorithms for radar interferometry: residue-cut, least-squares, and synthesis algorithms. J Opt Soc Am A 1998;15(3):586-598.
[36]Spagnolini U.2-D phase unwrapping and phase aliasing. Geophysics 1993;58(9):1324-1334.
[37]Cho B-L, Kong Y-K, Kim Y-S. Quality map extraction for radar interferometry using weighted window. Electronics Letters 2004;40(8):472-473.
[38]Hooper A, Zebker HA. Phase unwrapping in three dimensions with application to InSAR time series. J Opt Soc Am A 2007;24(9):2737-2747.
[39]Huntley JM, Saldner H. Temporal phase-unwrapping algorithm for automated interferogram analysis. Appl Opt 1993;32(17):3047-3052.
[40]Qianqian F, Meaney PM, Paulsen KD. The Multidimensional Phase Unwrapping Integral and Applications to Microwave Tomographical Image Reconstruction. Image Processing, IEEE Transactions on 2006; 15(11):3311-3324.
[41]Schofield MA, Zhu Y. Fast phase unwrapping algorithm for interferometric applications. Opt Lett 2003;28(14):1194-1196.
[42]Ahmad FH, Helms GL, Castellane RM, Durst BP. A discrete phase-unwrapping technique for interferograms. Microwave and Optical Technology Letters 2002;32(2):101-104.
[43]Fried DL. Least-square fitting a wave-front distortion estimate to an array of phase-difference measurements. J Opt Soc Am 1977;67(3):370-375.
[44]Meneses J, Gharbi T, Humbert P. Phase-unwrapping algorithm for images with high noise content based on a local histogram. Appl Opt 2005;44(7):1207-1215.
[45]Parkhurst J, Price G, Sharrock P, Moore C. Phase unwrapping algorithms for use in a true real-time optical body sensor system for use during radiotherapy. Appl Opt 2011;50(35):6430-6439.
[46]Su X, Chen W. Reliability-guided phase unwrapping algorithm:A review. Optics and Lasers in Engineering 2004;42(3):245-261.
[47]Karout S. Two-Dimensional Phase Unwrapping. Liverpool John Moores University; 2007.
[48]Volkov W, Zhu Y. Deterministic phase unwrapping in the presence of noise. Opt Lett 2003;28(22):2156-2158.
[49]Saldner HO, Huntley JM. Temporal phase unwrapping:application to surface profiling of discontinuous objects. Appl Opt 1997;36(13):2770-2775.
[50]Ghiglia DC, Romero LA. Minimum Lp-norm two-dimensional phase unwrapping. J Opt Soc Am A 1996;13(10):1999-2013.
[51]Ghiglia DC, Pritt MD. Two-dimensional phase unwrapping:Theory, algorithms, and software. New York:Wiley; 1998.493 p.
[52]Jenkinson M. Fast, automated, N-dimensional phase-unwrapping algorithm. Magnetic Resonance in Medicine 2003;49(1):193-197.
[53]Herraez MA, Gdeisat MA, Burton DR, Lalor MJ. Robust, Fast, and Effective Two-Dimensional Automatic Phase Unwrapping Algorithm Based on Image Decomposition. Appl Opt 2002;41(35):7445-7455.
[54]Abdul-Rahman HS, Gdeisat MA, Burton DR, Lalor MJ. Three-dimensional phase unwrapping algorithms:a comparison. In:Photon06 Conference, Manchester, UK, 2006:4-7.
[55]Huntley JM. Three-Dimensional Noise-Immune Phase Unwrapping Algorithm. Appl Opt 2001;40(23):3901-3908.
[56]Marklund O, Huntley J, Cusack R. Robust unwrapping algorithm for 3-D phase volumes of arbitrary shape containing knotted phase singularity loops. Research report, Lulea University of Technology, Sweden; 2005.
[57]Gdeisat M, Arevalillo-Herraez M, Burton D, Lilley F. Three-dimensional phase unwrapping using the Hungarian algorithm. Optics letters 2009;34(19):2994-2996.
[58]Salfity M, Huntley J, Graves M, Marklund O, Cusack R, Beauregard D. Extending the dynamic range of phase contrast magnetic resonance velocity imaging using advanced higher-dimensional phase unwrapping algorithms. Journal of The Royal Society Interface 2006;3(8):415-427.
[59]Kuhn HW. The Hungarian method for the assignment problem. Naval research logistics quarterly 1955;2(1-2):83-97.
[60]Brakke KA. The surface evolver. Experimental mathematics 1992; 1(2):141-165.
[61]Abdul-Rahman HS, Gdeisat MA, Burton DR, Lalor MJ, Lilley F, Moore CJ. Fast and robust three-dimensional best path phase unwrapping algorithm. Appl Opt 2007;46(26):6623-6635.
[62]Abdul-Rahman H, Gdeisat M, Burton D, Lalor M. Fast three-dimensional phase-unwrapping algorithm based on sorting by reliability following a non-continuous path.2005;(5856):33.
[63]Herraez MA, Burton DR, Lalor MJ, Gdeisat MA. Fast two-dimensional phase-unwrapping algorithm based on sorting by reliability following a noncontinuous path. Applied optics 2002;41(35):7437-7444.
[64]Abdul-Rahman H, Arevalillo-Herraez M, Gdeisat M, et al. Robust three-dimensional best-path phase-unwrapping algorithm that avoids singularity loops. Appl Opt 2009;48(23):4582-4596.
[65]Arevalillo-Herraez M, Gdeisat MA, Burton DR. Hybrid robust and fast algorithm for three-dimensional phase unwrapping. Applied optics 2009;48(32):6313-6323.
[66]Arevalillo-Herraez M, Burton DR, Lalor MJ. Clustering-based robust three-dimensional phase unwrapping algorithm. Applied optics 2010;49(10):1780-1788.
[67]Langley J, Zhao Q. A model-based 3D phase unwrapping algorithm using Gegenbauer polynomials. Physics in Medicine and Biology 2009;54(17):5237.
[68]Su X, Chen W, Zhang Q, Chao Y. Dynamic 3-D shape measurement method based on FTP. Optics and Lasers in Engineering 2001;36(1):49-64.
[69]Venkatesh BA, Schiros CG, Gupta H, Lloyd SG, Dell'Italia L, Denney TS. Three-dimensional plus time biventricular strain from tagged MR images by phase-unwrapped harmonic phase. Journal of Magnetic Resonance Imaging 2011;34(4):799-810.
[70]Liang ZP. A model-based method for phase unwrapping. IEEE Transactions on Medical Imaging 1996;15(6):893-897.
[71]Lei Y, Zhi-Pei L, Munson DC, Jr., Koetter R, Frey BJ. Unwrapping of MR phase images using a Markov random field model. IEEE Transactions on Medical Imaging 2006;25(1):128-136.
[72]Nico G, Palubinskas G, Datcu M. Bayesian approaches to phase unwrapping: Theoretical study. IEEE Transactions on Signal Processing 2000;48(9):2545-2556.
[73]Langley J, Zhao Q. Unwrapping magnetic resonance phase maps with Chebyshev polynomials. Magnetic Resonance Imaging 2009;27(9):1293-1301.
[74]Bioucas-Dias JM, Valadao G. Phase unwrapping via graph cuts. IEEE Transactions on Image Processing 2007;16(3):698-709.
[75]Flynn TJ. Two-dimensional phase unwrapping with minimum weighted discontinuity. J Opt Soc Am A 1997; 14(10):2692-2701.
[76]Langley JA, Brice RG, Zhao Q. Recursive approach to the moment-based phase unwrapping method. Appl Opt 2010;49(16):3096-3101.
[77]Bioucas-Dias JM, Valadao G. Phase unwrapping:a new max-flow/min-cut based approach. In:Image Processing,2005 ICIP 2005 IEEE International Conference on, 2005;(2):Ⅱ-77-80.
[78]Itoh K. Analysis of the phase unwrapping algorithm. Applied optics 1982;21(14):2470-2470.
[79]Kreyszig E. Advanced engineering mathematics:John Wiley & Sons; 2007.
[80]Goldstein R, Zebker H, Werner C. Satellite radar interferometry:Two-dimensional phase unwrapping. Radio Science 1988;23(4):713-720.
[81]Bone DJ. Fourier fringe analysis:the two-dimensional phase unwrapping problem. Appl Opt 1991;30(25):3627-3632.
[82]Hosseini MS. Derivative Compressive Sampling with Application to Inverse Problems and Imaging. Volume Master. Waterloo:University of Waterloo; 2010.
[83]Gutmann B, Weber H. Phase Unwrapping with the Branch-Cut Method: Clustering of Discontinuity Sources and Reverse Simulated Annealing. Appl Opt 1999;38(26):5577-5593.
[84]Buckland JR, Huntley JM, Turner SRE. Unwrapping noisy phase maps by use of a minimum-cost-matching algorithm. Appl Opt 1995;34(23):5100-5108.
[85]Gens R. Two-dimensional phase unwrapping for radar interferometry: Developments and new challenges. International Journal of Remote Sensing 2003;24(4):703-710.
[86]Wilkinson AJ. Techniques for 3-D surface reconstruction using synthetic aperture radar interferometry. Volume Doctor. London, United Kingdom:University College London; 1997.
[87]Pritt MD. Phase unwrapping by means of multigrid techniques for interferometric SAR. IEEE Transactions on Geoscience and Remote Sensing 1996;34(3):728-738.
[88]Frost VS, Stiles JA, Shanmugan KS, Holtzman J. A Model for Radar Images and Its Application to Adaptive Digital Filtering of Multiplicative Noise. Pattern Analysis and Machine Intelligence, IEEE Transactions on 1982;PAMI-4(2):157-166.
[89]Eineder M, Hubig M, Milcke B. Unwrapping large interferograms using the minimum cost flow algorithm. In:IGARSS'98, Seattle, Washington,1998. IEEE;(1):83-87.
[90]Refice A, Satalino G, Stramaglia S, Chiaradia MT, Veneziani N. Weights determination for minimum cost flow InSAR phase unwrapping. In:Geoscience and Remote Sensing Symposium,1999 IGARSS '99 Proceedings IEEE 1999 International, 1999;(2):1342-1344 vol.1342.
[91]CarbaIlo GF, Fieguth PW. Probabilistic cost functions for network flow phase; unwrapping. IEEE Transactions on Geoscience and Remote Sensing 2000;38(5):2192-2201.
[92]Witoszynskyj S, Rauscher A, Reichenbach JR, Barth M. Phase unwrapping of MR images using ΦUN-A fast and robust region growing algorithm. Medical Image Analysis 2009;13(2):257-268.
[93]Flynn TJ. Consistent 2-D phase unwrapping guided by a quality map.1996. IEEE;(4):2057-2059 vol.2054.
[94]Salfity MF, Ruiz PD, Huntley JM, Graves MJ, Cusack R, Beauregard DA. Branch cut surface placement for unwrapping of undersampled three-dimensional phase data: application to magnetic resonance imaging arterial flow mapping. Appl Opt 2006;45(12):2711-2722.
[95]Fang S, Meng L, Wang L, Yang P, Komori M. Quality-guided phase unwrapping algorithm based on reliability evaluation. Applied optics 2011;50(28):5446-5452.
[96]Zhou K, Zaitsev M, Bao S. Reliable two-dimensional phase unwrapping method using region growing and local linear estimation. Magnetic Resonance in Medicine 2009;62(4):1085-1090.
[97]An L, Xiang QS, Chavez S. A fast implementation of the minimum spanning tree method for phase unwrapping. Medical Imaging, IEEE Transactions on 2000;19(8):805-808.
[98]Xu W, Cumming I. A region-growing algorithm for InSAR phase unwrapping. IEEE Transactions on Geoscience and Remote Sensing 1999;37(1):124-134.
[99]Chen CW, Zebker HA. Network approaches to two-dimensional phase unwrapping: Intractability and two new algorithms. J Opt Soc Am A 2000;17(3):401-414.
[100]Ahuja RK, Magnanti TL, Orlin JB. Network flows:theory, algorithms, and applications. Englewood Cliffs, NJ:Prentice hall; 1993.
[101]Huntley JM. Noise-immune phase unwrapping algorithm. Appl Opt 1989;28(16):3268-3270.
[102]Cusack R, Huntley JM, Goldrein HT. Improved noise-immune phase-unwrapping algorithm. Appl Opt 1995;34(5):781-789.
[103]Wilson DR. Minimax theorems in graph theory. In:Beineke LW, Wilson RJ, editors. Selected Topics in Graph Theory. New York:Academic Press; 1978. p. 250-252.
[104]Karout SA, Gdeisat MA, Burton DR, Lalor MJ. Two-dimensional phase unwrapping using a hybrid genetic algorithm. Applied optics 2007;46(5):730-743.
[105]Huntley JM, Buckland JR. Characterization of sources of 2π phase discontinuity in speckle interferograms. J Opt Soc Am A 1995; 12(9):1990-1996.
[106]Ma J. Dixon techniques for water and fat imaging. Journal of Magnetic Resonance Imaging 2008;28(3):543-558.
[107]Costantini M. A novel phase unwrapping method based on network programming. IEEE Transactions on Geoscience and Remote Sensing 1998;36(3):813-821.
[108]Chen CW, Zebker HA. Phase unwrapping for large SAR interferograms: statistical segmentation and generalized network models. IEEE Transactions on Geoscience and Remote Sensing 2002;40(8):1709-1719.
[109]Chen CW. Statistical-cost network-flow approaches to two-dimansional phase unwrapping for radar interferometry. Volume Doctor of Philosophy:Stanford University; 2001.
[110]Kramer R, Loffeld O. A novel procedure for cutline detection. International Journal Of Electronics and Communications 1996;50(2):112-116.
[111]Karout SA, Gdeisat MA, Burton DR, Lalor MJ. Residue vector, an approach to branch-cut placement in phase unwrapping:theoretical study. Applied optics 2007;46(21):4712-4727.
[112]Spottiswoode B.2D phase unwrapping algorithms. http://www.mathworks.com/matlabcentral/fileexchange/22504-2d-phase-unwrapping-a lgorithms:2008.
[113]Ching NH, Rosenfeld D, Braun M. Two-dimensional phase unwrapping using a minimum spanning tree algorithm. Image Processing, IEEE Transactions on 1992;l(3):355-365.
[114]Lilley F, Lalor MJ, Burton DR. Robust fringe analysis system for human body shape measurement. Optical Engineering 2000;39(1):187-195.
[115]Herraez MA, Boticario JG, Lalor MJ, Burton DR. Agglomerative clustering-based approach for two-dimensional phase unwrapping. Applied optics 2005;44(7):1129-1140.
[116]Baldi A. Two-dimensional phase unwrapping by quad-tree decomposition. Applied optics 2001;40(8):1187-1194.
[117]Herraez MA, Burton DR, Lalor MJ, Clegg DB. Robust, simple, and fast algorithm for phase unwrapping. Applied optics 1996;35(29):5847-5852.
[118]Stephenson P, Burton DR, Lalor MJ. Data validation techniques in a tiled phase unwrapping algorithm. Optical Engineering 1994;33(11):3703-3708.
[119]Gierloff JJ. Phase unwrapping by regions. In:Proc SPIE,1987;(818):2-9.
[120]Hung KM, Yamada T. Phase unwrapping by regions using least-squares approach. Optical Engineering 1998;37(11):2965-2970.
[121]Schneider E, Glover G. Rapid in vivo proton shimming. Magnetic Resonance in Medicine 1991;18(2):335-347.
[122]Kemao Q, Gao W, Wang H. Windowed Fourier filtered and quality guided phase unwrapping algorithm:on locally high-order polynomial phase. Applied optics 2010;49(7):1075-1079.
[123]Chen CW, Zebker HA. Two-dimensional phase unwrapping with use of statistical models for cost functions in nonlinear optimization. Journal of the Optical Society of America A 2001;18(2):338-351.
[124]Ghiglia DC, Romero LA. Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods. J Opt Soc Am A 1994;11(1):107-117.
[125]Pritt MD, Shipman JS. Least-squares two-dimensional phase unwrapping using FFT's. Geoscience and Remote Sensing, IEEE Transactions on 1994;32(3):706-708.
[126]Hunt BR. Matrix formulation of the reconstruction of phase values from phase differences. J Opt Soc Am 1979;69(3):393-399.
[127]Kim S, Kim Y-S. Least squares phase unwrapping in wavelet domain. IEE Proceedings-Vision, Image and Signal Processing 2005;152(3):261-267.
[128]Marroquin JL, Rivera M. Quadratic regularization functionals for phase unwrapping. JOSAA 1995;12(11):2393-2400.
[129]Huang HYH, Tian L, Zhang Z, Liu Y, Chen Z, Barbastathis G. Path-independent phase unwrapping using phase gradient and total-variation (TV) denoising. Opt Express 2012;20(13):14075-14089.
[130]Fornaro G, Franceschetti G, Lanari R. Interferometric SAR phase unwrapping using Green's formulation. Geoscience and Remote Sensing, IEEE Transactions on 1996;34(3):720-727.
[131]Fornaro G, Franceschetti G, Lanari R, Sansosti E. Robust phase-unwrapping techniques:a comparison. JOSAA 1996;13(12):2355-2366.
[132]Marklund 0. An anisotropic evolution formulation applied in 2-D unwrapping of discontinuous phase surfaces. Image Processing, IEEE Transactions on 2001;10(11):1700-1711.
[133]Wilkinson AJ, Datcu M. Interferometric SAR topographic mapping, a Bayesian model-based approach. In:Communications and Signal Processing,1997 COMSIG '97, Proceedings of the 1997 South African Symposium on,1997:83-88.
[134]Boykov Y, Kolmogorov V. An experimental comparison of min-cut/max-flow algorithms for energy minimization in vision. Pattern Analysis and Machine Intelligence, IEEE Transactions on 2004;26(9):1124-1137.
[135]Tipper DJ, Burton DR, Lalor MJ. A Neural Network Approach To The Phase Unwrapping Problem In Fringe Analysis. Nondestructive Testing and Evaluation 1996; 12(6):391-400.
[136]Hamzah S, Pearson JD, Lisboa PJ, Hobson CA. Phase unwrapping in 3-D shape measurement using artificial neural networks. In:Sixth International Conference on Image Processing and Its Applications 1997;(2):680-683.
[137]Schwartzkopf W, Milner TE, Ghosh J, Evans BL, Bovik AC. Two-Dimensional Phase Unwrapping Using Neural Networks. Proceedings of the 4th IEEE Southwest Symposium on Image Analysis and Interpretation:IEEE Computer Society; 2000. p.274-277.
[138]Stramaglia S, Guerriero L, Pasquariello G, Veneziani N. Mean-Field Annealing for Phase Unwrapping. Applied optics 1999;38(8):1377-1383.
[139]Collaro A, Franceschetti G, Palmieri F, Ferreiro MS. Phase unwrapping by means of genetic algorithms. Journal of the Optical Society of America A 1998;15(2):407-418.
[140]Hui PTH, Wang Yin C. Towards hybrid 2D phase unwrapping using fuzzy clustering and neuro-fuzzy learning for SAR images:a case study on IFSAR phase image. In:Proceedings of International Symposium on Intelligent Multimedia, Video and Speech Processing,2001:271-274.
[141]Guerriero L, Nico G, Pasquariello G, Stramaglia S. New Regularization Scheme for Phase Unwrapping. Applied optics 1998;37(14):3053-3058.
[142]Liang JJ, Qin AK, Suganthan PN, Baskar S. Comprehensive learning particle swarm optimizer for global optimization of multimodal functions. Evolutionary Computation, IEEE Transactions on 2006;10(3):281-295.
[143]Juang C-F. A hybrid of genetic algorithm and particle swarm optimization for recurrent network design. IEEE Transactions on Systems, Man, and Cybernetics--Part B:Cybernetics 2004;34(2):997-1006.
[144]杨维,李歧强.粒子群优化算法综述.中国工程科学2004;6(5):87-94.
[145]Kennedy J, Eberhart R. Particle swarm optimization. In:IEEE international conference on neural networks,1995;(4):1942-1948.
[146]Shi Y, Eberhart R. A modified particle swarm optimizer. In:IEEE World Congress on Computational Intelligence, The 1998 IEEE International Conference on Evolutionary Computation Proceedings,1998. IEEE:69-73.
[147]OTSU N. A Threshold Selection Method from Gray-Level Histograms. IEEE Transactions on Systems, Man and Cybernetics 1979;9(1):62-66.
[148]Cuiru W, Jiangwei Z, Jing Y, Chaoju H, Jun L. A Modified Particle Swarm Optimization Algorithm and its Application For Solving Traveling Salesman Problem. In:International Conference on Neural Networks and Brain (ICNN&B'05), 2005;(2):689-694.
[149]Wang K-P, Huang L, Zhou C-G, Pang W. Particle swarm optimization for traveling salesman problem. In:International Conference on Machine Learning and Cybernetics,2003. IEEE;(3):1583-1585.
[150]Wikipedia(?). Flood fill. http://en.wikipedia.org/wiki/Flood fill.
[151]Chen K, Xi J, Yu Y, Chicharo JF. Fast quality-guided flood-fill phase unwrapping algorithm for three-dimensional fringe pattern profilometry. In:Photonics Asia 2010, Beijing,2010. International Society for Optics and Photonics:78550X.
[152]Dias JMB, Leitao JMN. The ZπM algorithm:A method for interferometric image reconstruction in SAR/SAS. IEEE Transactions on Image Processing 2002;11(4):408-422.
[153]Dardyk G, Yavneh I. A multigrid approach to two-dimensional phase unwrapping. Numerical Linear Algebra with Applications 2004;11(2-3):241-259.
[154]Shalem I, Yavneh I. A multilevel graph algorithm for two dimensional phase unwrapping. Computing and Visualization in Science 2008;11(2):89-100.
[155]Tewarson RP. Sparse Matrices. New York:Academic Press Inc.; 1973.160 p.
[156]Pooch UW, Nieder A. A Survey of Indexing Techniques for Sparse Matrices. ACM Comput Surv 1973;5(2):109-133.
[157]Scott JA, Hu Y. Experiences of sparse direct symmetric solvers. ACM Trans Math Soft 2007;33(3):18.
[158]Watkins DS. Fundamentals of matrix computations,2nd edition. New York: John Wiley & Sons; 2002.618 p.
[159]Krishnamoorthy A, Menon D. Matrix Inversion Using Cholesky Decomposition. arXiv preprint arXiv:11114144 2011.
[160]Duff IS, Erisman AM, Reid JK. Direct methods for sparse matrices. New York: Oxford University Press, Inc.; 1986.341 p.
[161]Gould NI, Scott JA, Hu Y. A numerical evaluation of sparse direct solvers for the solution of large sparse symmetric linear systems of equations. ACM Trans Math Softw 2007;33(2):10.
[162]Chen Y, Davis TA, Hager WW, Rajamanickam S. Algorithm 887:CHOLMOD, Supernodal Sparse Cholesky Factorization and Update/Downdate. ACM Trans Math Softw 2008;35(3):1-14.
[163]Duff IS. MA57---a code for the solution of sparse symmetric definite and indefinite systems. ACM Trans Math Softw 2004;30(2):118-144.
[164]He W, Cheng Y, Xia L, Liu F. A New Particle Swarm Optimization-Based Method for Phase Unwrapping of MRI Data. Computational and Mathematical Methods in Medicine 2012;2012:9.
[165]Rauscher A, Barth M, Reichenbach JR, Stollberger R, Moser E. Automated unwrapping of MR phase images applied to BOLD MR-venography at 3 Tesla. Journal of Magnetic Resonance Imaging 2003; 18(2):175-180.