A Parallel InSAR Phase Unwrapping Method Based on Separated Continuous Regions

Phase unwrapping is an imperative step in interferometry processing that has a significant influence on the quality of subsequent products. Many existing phase unwrapping algorithms have been designed to solve for the unwrapped phase under the assumption that noisy areas with discontinuities are small or that reliable continuity can be recovered there. They attempt to restore the unwrapped phase by using continuity and data quality measures, such as residues. However, when the observing field is divided into separate zones of continuous phase due to a large range of noise, such as those caused by rivers or mountains, it is difficult to use traditional phase unwrapping techniques to recover global continuity in these noisy areas. To address this challenge, we present a two-dimensional parallel phase unwrapping method that is designed to handle cases where the continuity of the phase is separated by closed noisy loops. Based on continuity distances, this method aims to identify continuous regions that are free of hidden phase discontinuities and restore phase continuity between the separated regions. A heterogeneous residual diffusion scheme is used to restore the unwrapped phase outside continuous regions. The parallel algorithm for extracting continuous regions, restoring continuity between the regions, and diffusing residuals was implemented on a GPU device to increase the processing efficiency. We applied our method to typical TanDEM-X data covering rivers, islands, and mountains and demonstrated that it is a promising solution for large-scale, heavily noisy phase unwrapping problems.