Orbit error removal in InSAR/MTInSAR with a patch-based polynomial model

Yanan Du,Haiqiang Fu,Lin Liu,Guangcai Feng,Xing Peng,Debao Wen

DOI: https://doi.org/10.1016/j.jag.2021.102438

IF: 7.5

2021-10-01

International Journal of Applied Earth Observation and Geoinformation

Abstract:<p>The orbit error caused by the inaccuracy of the orbit state vector can lead to fringes in differential interferograms, which can impede the estimation of deformation in differential SAR interferometry (DInSAR) applications. Usually, a set of polynomial coefficients for an entire SAR image is obtained for orbit error removal. However, the orbit error plane is influenced by overfitting in the case that the SAR satellites do not have a precise orbit. In this paper, a patch-based polynomial method is proposed to fit the orbit error plane. The new method divides an SAR image into several overlapping patches in the azimuth and range directions. Every patch obtains its own polynomial coefficients, and an iterative least-square method is used to mosaic the orbit plane. This method is tested and validated via a simulated dataset and then applied to ALOS1/2 PALSAR and Sentinel-1A datasets. The accuracy of deformation is evaluated by <em>in situ</em> GPS datasets. The results show that the patch-based method can fit the orbit phase plane more accurately than the traditional polynomial model with millimeter-level displacement improvement, especially in the margin areas of ALOS1/2 and for the wide-coverage Sentinel-1A datasets. Moreover, in the MTInSAR parameter calculations, the new method improves the accuracy of mean velocity calculations for ALOS1 time series, with a reduction of RMSE from 4.47 mm/yr to 3.17 mm/yr. Additionally, the new method reduces the spatial correlation of the residual topographic phase, with a mean value reduction from 0.32 m to 0.13 m.</p>

remote sensing

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the orbital error problem in Differential Synthetic Aperture Radar Interferometry (DInSAR) caused by the inaccuracy of satellite orbit state vectors. This kind of orbital error will generate fringes in the differential interferogram, affecting the accuracy of deformation estimation. The traditional solution is to obtain a set of polynomial coefficients for the entire SAR image to remove the orbital error, but this method is prone to over - fitting problems when the satellite orbit is not accurate enough, especially performing poorly on edge regions and wide - coverage data sets.

To overcome these problems, the paper proposes a patch - based polynomial model method. This method divides the SAR image into multiple overlapping patches, obtains polynomial coefficients independently on each patch, and uses the iterative least - squares method to splice the orbital plane. In this way, the orbital phase plane can be fitted more accurately, especially on edge regions and wide - coverage data sets such as ALOS1/2 and Sentinel - 1A data sets. Experimental results show that the new method is superior to the traditional polynomial model in terms of millimeter - level displacement improvement, especially on the edge regions of ALOS1/2 and the wide - coverage data set of Sentinel - 1A. In addition, the new method also improves the accuracy of multi - temporal InSAR (MTInSAR) parameter calculation and reduces the spatial correlation of the residual topographic phase.

Specifically, the main contributions of the new method are as follows:
1. **Improve the fitting accuracy of the orbital phase plane**: By dividing the image into multiple patches and obtaining polynomial coefficients independently for each patch, the over - fitting problem in the traditional method is avoided.
2. **Improve the processing effect in the edge region**: In the edge region of the image, the new method can remove the orbital error more accurately and reduce the appearance of fringes.
3. **Improve the accuracy of MTInSAR parameter calculation**: The new method shows better performance in time - series analysis, especially in terms of the RMSE (root - mean - square error) of the average velocity calculation.
4. **Reduce the spatial correlation of the residual topographic phase**: The new method can separate the orbital error and deformation more effectively and reduce the influence of the residual topographic phase.

These improvements make the new method more reliable and accurate when dealing with complex terrains and large - range - coverage data sets.