IGS-CMAES: A Two-Stage Optimization for Ground Deformation and DEM Error Estimation in Time Series InSAR Data

Interferometric synthetic aperture radar (InSAR) has become an increasingly recognized remote sensing technology for earth surface monitoring. Slow and subtle terrain displacements can be estimated using time-series InSAR (TSInSAR) data. However, a substantial increase in the availability of exclusive time series data necessitates the development of more efficient and effective algorithms. Research in these areas is usually carried out by solving complicated optimization problems, which is very computationally expensive and time-consuming. This work proposes a two-stage black-box optimization framework to jointly estimate the average ground deformation rate and terrain digital elevation model (DEM) error. The method performs an iterative grid search (IGS) to acquire coarse candidate solutions, and then a covariance matrix adaptive evolution strategy (CMAES) is adopted to obtain the final local results. The performance of our method is evaluated using both simulated and real datasets. Both quantitative and qualitative comparisons using different optimizers support the reliability and effectiveness of our work. The proposed IGS-CMAES achieves higher accuracy with a significantly fewer number of objective function evaluations than other established algorithms. It offers the possibility for wide-area monitoring, where high precision and real-time processing is essential.