ABSTRACT
The phase unwrapping is the process of attempting to reconstruct the true phase from modulo 2π phase values. This procedure requires that we have an important congruent constraint i.e. rewrapping unwrapped image should be identical to the original wrapped image. This constraint condition causes a discrete optimization problem. However, many methods have ignored this constraint condition to solve a continuous minimization problem directly, making it difficult to solve the solution correctly. We recently presented new continuous minimum norm method that is based on the Lovász extension. Our method can reach the optimal solution with the congruent constraint condition. Note that in our work, we have taken the subgradient method for minimization, but it is time consuming. In this paper, we introduce the incremental subgradient method to the minimization procedure, which is faster than the subgradient method. To solve the phase unwrapping, first, we also use the Lovász extension to transform the phase unwrapping problem to equivalent continuous minimization problem which consists of the sum of large numbers of component functions. Then we take the incremental subgradient method to solve the minimization problem in which operate on a single component at each iteration, rather than on the entire cost function. On the other hand, we also introduce new minimal function to deal with high lever noise. Several simulations show, compared with the subgradient method, the new algorithm only takes one quarter of (or less) the numbers of iteration for the convergence.
- A. Nedić and D.P. Bertsekas, Incremental subgradient methods for nondifferentiable optimization, SIAM Journal on Optimization 56 (2001), no. 1, 109--138.Google Scholar
- Bach, F. 2011. Learning with submodular functions: A convex optimization perspective. arXiv: 1111.6453 (2011).Google Scholar
- Bioucas-Dias, J. M., and Valadão, G. 2007. Phase unwrapping via graph cuts. IEEE Trans. Image Process. 16(2007), 698--709.Google ScholarDigital Library
- Campbell, N. D. F., Vogiatzis, G., Hernández, C., and Cipolla, R. 2010. Automatic 3D object segmentation in multiple views using volumetric graph-cuts. Image and Vision Computing, 28(2010), 14--25.Google ScholarDigital Library
- Dong, J., Chen, F., Zhou, D., Liu, T., Yu, Z., and Y. Wang. 2017. Phase unwrapping with graph cuts optimization and dual decomposition acceleration for 3D high-resolution MRI data. Magn Reson Med 77(2017), 1353--1358.Google ScholarCross Ref
- Flynn, T. 1997. Two-dimensional phase unwrapping with minimum weighted discontinuity. J. Opt. Soc. Amer. A. 14(1997), 2692--2701.Google ScholarCross Ref
- Ghiglia, D., and Pritt, M. 1998. Two-Dimensional Phase Unwrapping. Theory, Algorithms, and Software (New York: Wiley, 1998).Google Scholar
- Huang, Howard Y. H., Tian, L., Zhang, Z., Liu, Y., Chen, Z., and Barbastathis, G. 2012. Path-independent phase unwrapping using phase gradient and total-variation(TV) denoising. Opt. Express. 20(2012), 14075--14089.Google ScholarCross Ref
- Itoh, K. 1982. Analysis of the phase unwrapping problem. Appl. Opt. 21(1982), 2470--2470.Google ScholarCross Ref
- Jezzard, P., and Balaban, R. S. 1995. Correction for geometric distortion in echo planar images from B_0 field variations. Magnetic Resonance in medicine. 34(1995) 65--73.Google Scholar
- Kudo, H., Lian, S., Wada, K., Two-Dimensional Phase Unwrapping with Continuous Submodular Minimization. In Proc. the 2018 International Conference on Digital Medicine and Image Processing, 6--13.Google Scholar
- Loecher, M., Schrauben, E., Johnson, K., and Wieben, O. 2016. Phase unwrapping in 4D MR Flow with a 4D single-step Laolacian Algorithm. Journal of Magneloechertic Resonance Imajing. 43(2016), 833--842.Google ScholarCross Ref
- Lovász, L. 1983. Submodular functions and convexity. In Mathematical programming: The state of the art, A. Bachem, M. Grötschel, and B. Korte, eds., (Springer-Verlag, Berlin, 1983), 235--257.Google Scholar
- Martin, D., Fowlkes, C., Tal. D., and Malik, J. 2001. A database of human segmented natural images and its application to evaluating segmentation algorithms and measuring ecological statistics. In Proc. 8th Int. Conf. Comput. Vis. 2(2001), 416--423.Google ScholarCross Ref
- Massonnet, D., and Feigl, K. L. 1998. Radar interferometry and its application to changes in the earth's surface. Rev. Geophys. 36(1998), 441--500.Google ScholarCross Ref
- Pfeiffer, F., Weitkamp, T., Oliver, B., Christian, D. 2006. Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources. Nature Phys. 2(2006), 258--261.Google ScholarCross Ref
- Rosen, P. A., Hensley, S., Joughin, I. R., Li, F. K., Madsen, S. N., Rodriguez, E., and Goldstein, R. 2000. Synthetic aperture radar interferometry. Proc. IEEE. 88(2000), 333--382.Google ScholarCross Ref
- Weitkamp, T., Diaz, A., David, C., Christian, D., Pfeiffer, F., Stampanoni, M., Cloetens, P., Ziegler, E. 2005. Quantitative X-ray phase imaging with a grating interferometer. Opt. Express. 13(2005), 6296--6304.Google ScholarCross Ref
- Ying L. Phase unwrapping. In: M Akay, editor. Wiley encyclopedia of biomedical engineering. 1st Ed. Hoboken, NJ: John Wiley & Sons; (2006).Google Scholar
Index Terms
- Improved algorithm for Phase Unwrapping with Continuous Submodular Minimization
Recommendations
Two-Dimensional Phase Unwrapping with Continuous Submodular Minimization
DMIP '18: Proceedings of the 2018 International Conference on Digital Medicine and Image ProcessingThe phase unwrapping is recovering true phase from its 2π modulo observations which are related to some discrete optimization problems. The challenge is to exactly solve the discrete optimization problem arising from noisy data. In this paper, we ...
Two-dimensional Phase Unwrapping with Priori Boundary Condition
ICVISP 2019: Proceedings of the 3rd International Conference on Vision, Image and Signal ProcessingPhase unwrapping is essential for many applications such as Interferometric synthetic aperture radar, MRI, and X-ray phase imaging etc. In these applications, the phase is determined only in the principal value range of [-π,π). The unwrapping that ...
Phase Unwrapping Using Energy Minimization Methods for MRI Phase Image
Phase Unwrapping PU is reconstruction of absolute phase data from its wrapped phase. The absolute phase cannot be extracted from the wrapped phase data directly. Without phase noise, singularity, and aliasing problems, the phase information can be ...
Comments