2024/01/18 by Benjamin Dubois-Taine, Dubois-Taine, Benjamin, Roland Akiki +3
Computer Science · Earth and Planetary Sciences · Engineering · #3D Surveying and Cultural Heritage #Advanced Measurement and Metrology Techniques #FOS: Mathematics #Optical measurement and interference techniques #Optimization and Control (math.OC)
paper · pdf · doi:10.48550/arxiv.2401.09961
openalex publication_date 2024/01/18 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
The 2D phase unwrapping problem seeks to recover a phase image from its observation modulo 2π, and is a crucial step in a variety of imaging applications. In particular, it is one of the most time-consuming steps in the interferometric synthetic aperture radar (InSAR) pipeline. In this work we tackle the L1-norm phase unwrapping problem. In optimization terms, this is a simple sparsity-inducing problem, albeit in very large dimension. To solve this high-dimensional problem, we iteratively solve a series of numerically simpler weighted least squares problems, which are themselves solved using a preconditioned conjugate gradient method. Our algorithm guarantees a sublinear rate of convergence in function values, is simple to implement and can easily be ported to GPUs, where it significantly outperforms state of the art phase unwrapping methods.