A Parallel InSAR Phase Unwrapping Method Based on Separated Continuous Regions
Abstract
:1. Introduction
2. Principle and Method
2.1. Defining and Extracting Continuous Regions
2.2. Region Residue Pairs’ Detection and Region Phase Alignment
2.2.1. Definitions of Region Phase Difference and Region Residue
2.2.2. Balancing Region Residues and Aligning Region Phases
2.3. Region Growing Based on Heterogeneous Residual Diffusion
2.3.1. Initializing Unwrapped Phase Outside Regions
2.3.2. Heterogeneous Residual Diffusion
2.3.3. Continuous Region Growing
3. Numerical Implementation
- Construction of the phase continuity distance grid: the first step in constructing the phase continuity distance grid is the detection and marking of phase residues in the wrapped phase. This is followed by the calculation of the quality grid, which is defined as the inverse of the wrapped gradient variance, for exampleAs an optional operation, the procedure constructs a preliminary minimum balanced tree (MBT) for coupling dual-phase residues and marking simple discontinuity boundaries. This step can reduce the number of unbalanced residues and simplify the distribution of the following continuous regions. The parallel calculation of the distance grid and implementation of the dual residue pair coupling in the sampling phase can be found in our previous work [31]:
- Unwrapping within the continuous regions: With the given unwrapped phase references, local direct unwrapping is performed independently in each region. The unwrapping process is based on the relations between and its unwrapped neighbors such thatIf there are any neighboring phase pixels and in the continuous regions, away from known discontinuity boundaries, that meet the condition:
- Continuous region boundary expansion: The local unwrapped phase in the continuous regions is expanded outwards as follows with known , if is greater than or undefined:
- Corner continuity distance calculation: In order to improve the efficiency of corner continuity distance calculation, we utilized a parallel implementation of the fast march method (FMM) [34,35]. The positive and negative continuity distance of each corner point is determined through the examination of the discontinuity cost of boundaries and the distribution of region residues. The calculation is performed using known , in accordance with the following equation:
- Balancing region residues: We utilized the procedure defined in Equations (10)–(12) to identify pairs of detectable residues and assigned them the appropriate continuity distance values. Subsequently, the residues and boundary indexes are adjusted in accordance with the process outlined in Algorithm 1. Here, the boundary couple index, C, takes a value of zero in the case of an uncoupled pair, and the direction index of the pair, b, is defined as , depending on whether the direction of the pair and the direction of the boundary are the same.
Algorithm 1 Balancing region residues. - 1:
- repeat
- 2:
- for all unbalanced residue ri in parallel do
- 3:
- test coupling conditions in Equations (10)–(12)
- 4:
- if found a pair then
- 5:
- 6:
- for Equations (10) and (12);
- 7:
- trace pair boundaries ;
- 8:
- for in do
- 9:
- ;
- 10:
- end for
- 11:
- end if
- 12:
- end for
- 13:
- update distance and sources ;
- 14:
- until residues are all balanced.
- Aligning region phase: A reference region was chosen, and the unwrapped phase of all other regions was adjusted accordingly, using Equation (13), on boundaries without a couple for C. This resulted in the establishment of a global phase continuity framework.
- Initializing solution outside regions: The initial unwrapped phase outside the regions is obtained after detecting phase discontinuities in the horizontal and vertical directions. This was accomplished by solving the Poisson equations according to the linear discontinuities defined in Equation (14). The local Poisson equation in the horizontal or vertical direction outside the continuous regions can be written asIn our GPU implementation, we employed the parallel cyclic reduction (PCR) method [36] to solve the Poisson equation.
- Heterogeneous diffusion: Heterogeneous diffusion is applied to ensure that the unwrapped phase remains congruent with the wrapped phase, in accordance with Equations (17) and (18). This process is repeated to expand the regions until the phase stabilizes and converges. The final unwrapped result is obtained once convergence is achieved. It is important to note that the time interval should not break the stability of the diffusion process, and we set in our implementation.
4. Experiments and Analysis
5. Discussion
6. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Ghiglia, D.C.; Pritt, M.D. Two-Dimensional Phase Unwrapping: Theory, Algorithms, and Software; Wiley: Hoboken, NJ, USA, 1998. [Google Scholar]
- Yu, H.; Lan, Y.; Yuan, Z.; Xu, J.; Lee, H. Phase Unwrapping in InSAR: A Review. IEEE Geosci. Remote Sens. Mag. 2019, 7, 40–58. [Google Scholar] [CrossRef]
- Xu, W.; Cumming, I. A region-growing algorithm for InSAR phase unwrapping. IEEE Trans. Geosci. Remote Sens. 1999, 37, 124–134. [Google Scholar] [CrossRef] [Green Version]
- Karsa, A.; Shmueli, K. SEGUE: A Speedy rEgion-Growing Algorithm for Unwrapping Estimated Phase. IEEE Trans. Med Imaging 2019, 38, 1347–1357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Ojha, C.; Manunta, M.; Pepe, A.; Paglia, L.; Lanari, R. An innovative region growing algorithm based on Minimum Cost Flow approach for Phase Unwrapping of full-resolution differential interferograms. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 22–27 July 2012; pp. 5582–5585. [Google Scholar] [CrossRef]
- Herráez, M.A.; Burton, D.R.; Lalor, M.J.; Gdeisat, M.A. Fast two-dimensional phase-unwrapping algorithm based on sorting by reliability following a noncontinuous path. Appl. Opt. 2002, 41, 7437–7444. [Google Scholar] [CrossRef] [PubMed]
- Xie, X.; Zeng, Q. Multi-baseline extended particle filtering phase unwrapping algorithm based on amended matrix pencil model and quantized path-following strategy. J. Syst. Eng. Electron. 2019, 30, 78–84. [Google Scholar] [CrossRef]
- Dudczyk, J.; Kawalec, A. Optimizing the Minimum Cost Flow Algorithm for the Phase Unwrapping Process in SAR Radar. Bull. Pol. Acad. Sci. Tech. Sci. 2014, 62, 511–516. [Google Scholar] [CrossRef] [Green Version]
- Huang, Q.; Zhou, H.; Dong, S.; Xu, S. Parallel Branch-Cut Algorithm Based on Simulated Annealing for Large-Scale Phase Unwrapping. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3833–3846. [Google Scholar] [CrossRef]
- Chen, C.; Zebker, H. Phase unwrapping for large SAR interferograms: Statistical segmentation and generalized network models. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1709–1719. [Google Scholar] [CrossRef] [Green Version]
- Sawaf, F.; Tatam, R.P. Finding minimum spanning trees more efficiently for tile-based phase unwrapping. Meas. Sci. Technol. 2006, 17, 1428. [Google Scholar] [CrossRef]
- Dong, J.; Zhuo, Z.; Li, J.; He, Y. A new phase image reconstruction method using Markov random fields. In Proceedings of the 2017 IEEE/ACIS 16th International Conference on Computer and Information Science (ICIS), Wuhan, China, 24–26 May 2017; pp. 1–5. [Google Scholar] [CrossRef]
- Xu, J.; An, D.; Huang, X.; Wang, G. Phase Unwrapping for Large-Scale P-Band UWB SAR Interferometry. IEEE Geosci. Remote Sens. Lett. 2015, 12, 2120–2124. [Google Scholar] [CrossRef]
- Zhang, Y.; Xing, Z. A Hybrid Phase Unwrapping Algorithm Based on Quality-Guided and Surface-Fitting. In Proceedings of the 2018 IEEE International Workshop on Electromagnetics: Applications and Student Innovation Competition (iWEM), Nagoya, Japan, 29–31 August 2018; pp. 1–2. [Google Scholar] [CrossRef]
- Zhang, Y.; Xing, Z. A Phase Unwrapping Algorithm Based on Branch-cut and B-Spline Fitting in InSAR. In Proceedings of the 2018 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Hong Kong, China, 14–16 September 2018; pp. 1–3. [Google Scholar] [CrossRef]
- Yu, H.; Lee, H. A convex hull algorithm based fast large-scale two-dimensional phase unwrapping method. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 3824–3827. [Google Scholar] [CrossRef]
- Yu, H.; Lan, Y.; Xu, J.; An, D.; Lee, H. Large-Scale L0-Norm and L1-Norm 2-D Phase Unwrapping. IEEE Trans. Geosci. Remote Sens. 2017, 55, 4712–4728. [Google Scholar] [CrossRef]
- Zhang, Y.; Zhang, S.; Gao, Y.; Li, S.; Jia, Y.; Li, M. Adaptive Square-Root Unscented Kalman Filter Phase Unwrapping with Modified Phase Gradient Estimation. Remote Sens. 2022, 14, 1229. [Google Scholar] [CrossRef]
- Gao, Y.; Tang, X.; Li, T.; Lu, J.; Li, S.; Chen, Q.; Zhang, X. A Phase Slicing 2-D Phase Unwrapping Method Using the L 1 -Norm. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
- Li, H.; Zhong, H.; Ning, M.; Zhang, P.; Tang, J. Using neural networks to create a reliable phase quality map for phase unwrapping. Appl. Opt. 2023, 62, 1206. [Google Scholar] [CrossRef] [PubMed]
- Zhou, L.; Yu, H.; Lan, Y. Deep Convolutional Neural Network-Based Robust Phase Gradient Estimation for Two-Dimensional Phase Unwrapping Using SAR Interferograms. IEEE Trans. Geosci. Remote Sens. 2020, 58, 4653–4665. [Google Scholar] [CrossRef]
- Li, L.; Zhang, H.; Tang, Y.; Wang, C.; Gu, F. InSAR Phase Unwrapping by Deep Learning Based on Gradient Information Fusion. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
- Zhao, J.; Liu, L.; Wang, T.; Wang, X.; Du, X.; Hao, R.; Liu, J.; Liu, Y.; Zhang, J. VDE-Net: A two-stage deep learning method for phase unwrapping. Opt. Express 2022, 30, 39794. [Google Scholar] [CrossRef] [PubMed]
- Pu, L.; Zhang, X.; Zhou, Z.; Li, L.; Zhou, L.; Shi, J.; Wei, S. A Robust InSAR Phase Unwrapping Method via Phase Gradient Estimation Network. Remote Sens. 2021, 13, 4564. [Google Scholar] [CrossRef]
- Sawaf, F.; Groves, R.M. Phase discontinuity predictions using a machine-learning trained kernel. Appl. Opt. 2014, 53, 5439–5447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wu, Z.; Wang, T.; Wang, Y.; Wang, R.; Ge, D. Deep-Learning-Based Phase Discontinuity Prediction for 2-D Phase Unwrapping of SAR Interferograms. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–16. [Google Scholar] [CrossRef]
- Spoorthi, G.E.; Gorthi, S.; Gorthi, R.K.S.S. PhaseNet: A Deep Convolutional Neural Network for Two-Dimensional Phase Unwrapping. IEEE Signal Process. Lett. 2019, 26, 54–58. [Google Scholar] [CrossRef]
- Sica, F.; Calvanese, F.; Scarpa, G.; Rizzoli, P. A CNN-Based Coherence-Driven Approach for InSAR Phase Unwrapping. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
- Zhou, L.; Yu, H.; Lan, Y.; Xing, M. Deep Learning-Based Branch-Cut Method for InSAR Two-Dimensional Phase Unwrapping. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–15. [Google Scholar] [CrossRef]
- Li, H.; Zhong, H.; Zhang, P.; Ning, M.; Tang, J. Two-Dimensional Phase Unwrapping Based on Residual Prediction Neural Network. IEEE Geosci. Remote Sens. Lett. 2023, 20, 1–5. [Google Scholar] [CrossRef]
- Gao, J.; Sun, Z. Phase unwrapping method based on parallel local minimum reliability dual expanding for large-scale data. J. Appl. Remote Sens. 2019, 13, 038506. [Google Scholar] [CrossRef] [Green Version]
- Chen, C.W.; Zebker, H.A. Two-dimensional phase unwrapping with use of statistical models for cost functions in nonlinear optimization. JOSA A 2001, 18, 338–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Lyuboshenko, I.V.; Maitre, H.; Maruani, A. Least-mean-squares phase unwrapping by use of an incomplete set of residue branch cuts. Appl. Opt. 2002, 41, 2129–2148. [Google Scholar] [CrossRef]
- Capozzoli, A.; Curcio, C.; Liseno, A.; Savarese, S. A comparison of Fast Marching, Fast Sweeping and Fast Iterative Methods for the solution of the eikonal equation. In Proceedings of the 2013 21st Telecommunications Forum Telfor (TELFOR), Belgrade, Serbia, 26–28 November 2013; pp. 685–688. [Google Scholar]
- Jeong, W.K.; Whitaker, R.T. A Fast Iterative Method for Eikonal Equations. SIAM J. Sci. Comput. 2008, 30, 2512–2534. [Google Scholar] [CrossRef] [Green Version]
- Swarztrauber, P.N. The Methods of Cyclic Reduction, Fourier Analysis and the FACR Algorithm for the Discrete Solution of Poisson’s Equation on a Rectangle. SIAM Rev. 1977, 19, 490–501. [Google Scholar] [CrossRef]
- Tophu. Available online: https://github.com/opera-adt/tophu (accessed on 16 January 2023).
- isce-framework/isce3: InSAR Scientific Computing Environment. Available online: https://github.com/isce-framework/isce3 (accessed on 16 January 2023).
Data | (%) | (%) | (%) | Regions | Coverage (%) | |
---|---|---|---|---|---|---|
Wuhan | 192,475 | 0.1 | 0.1 | 0.8 | 7504 | 84.10 |
Singapore | 3,013,342 | 0.1 | 0.1 | 0.3 | 552 | 87.64 |
Songshan | 128,245 | 1.0 | 0.1 | 1.9 | 7246 | 53.87 |
Data | Bound | Bound | Corner | Pair | |||||
---|---|---|---|---|---|---|---|---|---|
Wuhan | 19,253 | 526 | 12,890 | 188 | 182 | 2 | 5 | 0 | 192 |
Singapore | 1402 | 35 | 973 | 96 | 95 | 2 | 1 | 0 | 105 |
Songshan | 20,298 | 200 | 13,762 | 1004 | 997 | 9 | 11 | 0 | 1054 |
Data | SNAPHU | Region Growing | GAMMMA | PHASS | Proposed |
---|---|---|---|---|---|
Wuhan | 98.072 | 85.919 | 73.510 | 67.049 | 70.695 |
Singapore | 86.992 | 27.809 | 64.580 | 23.046 | 23.324 |
Songshan | 93.014 | 56.511 | 161.610 | 43.068 | 35.997 |
Data | Constructing Map | Extracting Regions | Aligning Regions | Solving Outside |
---|---|---|---|---|
Wuhan | 15.8% | 15.4% | 3.2% | 65.6% |
Singapore | 16.3% | 3.6% | 2.4% | 77.7% |
Songshan | 21.8% | 17.9% | 4.1% | 56.2% |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Gao, J.; Jiang, H.; Sun, Z.; Wang, R.; Han, Y. A Parallel InSAR Phase Unwrapping Method Based on Separated Continuous Regions. Remote Sens. 2023, 15, 1370. https://doi.org/10.3390/rs15051370
Gao J, Jiang H, Sun Z, Wang R, Han Y. A Parallel InSAR Phase Unwrapping Method Based on Separated Continuous Regions. Remote Sensing. 2023; 15(5):1370. https://doi.org/10.3390/rs15051370
Chicago/Turabian StyleGao, Jian, Houjun Jiang, Zhongchang Sun, Ruisheng Wang, and Youmei Han. 2023. "A Parallel InSAR Phase Unwrapping Method Based on Separated Continuous Regions" Remote Sensing 15, no. 5: 1370. https://doi.org/10.3390/rs15051370