Compressed SAR Interferometry in the Big Data Era

Modern Synthetic Aperture Radar (SAR) missions provide an unprecedented massive interferometric SAR (InSAR) time series. The processing of the Big InSAR Data is challenging for long-term monitoring. Indeed, as most deformation phenomena develop slowly, a strategy of a processing scheme can be worked on reduced volume data sets. This paper introduces a novel ComSAR algorithm based on a compression technique for reducing computational efforts while maintaining the performance robustly. The algorithm divides the massive data into many mini-stacks and then compresses them. The compressed estimator is close to the theoretical Cramer–Rao lower bound under a realistic C-band Sentinel-1 decorrelation scenario. Both persistent and distributed scatterers (PSDS) are exploited in the ComSAR algorithm. The ComSAR performance is validated via simulation and application to Sentinel-1 data to map land subsidence of the salt mine Vauvert area, France. The proposed ComSAR yields consistently better performance when compared with the state-of-the-art PSDS technique. We make our PSDS and ComSAR algorithms as an open-source TomoSAR package. To make it more practical, we exploit other open-source projects so that people can apply our PSDS and ComSAR methods for an end-to-end processing chain. To our knowledge, TomoSAR is the first public domain tool available to jointly handle PS and DS targets.

Fusion of SAR Interferometry and Polarimetry Methods for Landslide Reactivation Study, the Bureya River (Russia) Event Case Study

In this paper, we demonstrate the estimation capabilities of landslide reactivation based on various SAR (Synthetic Aperture Radar) methods: Cloude-Pottier decomposition of Sentinel-1 dual polarimetry data, MT-InSAR (Multi-temporal Interferometric Synthetic Aperture Radar) techniques, and cloud computing of backscattering time series. The object of the study is the landslide in the east of Russia that took place on 11 December 2018 on the Bureya River. H-α-A polarimetric decomposition of C-band radar images not detected significant transformations of scattering mechanisms for the surface of the rupture, whereas L-band radar data show changes in scattering mechanisms before and after the main landslide. The assessment of ground displacements along the surface of the rupture in the 2019–2021 snowless periods was carried out using MT-InSAR methods. These displacements were 40 mm/year along the line of sight. The SBAS-InSAR results have allowed us to reveal displacements of great area in 2020 and 2021 snowless periods that were 30–40 mm/year along the line-of-sight. In general, the results obtained by MT-InSAR methods showed, on the one hand, the continuation of displacements along the surface of the rupture and on the other hand, some stabilization of the rate of landslide processes.

Parameterized Modeling and Calibration for Orbital Error in TanDEM-X Bistatic SAR Interferometry over Complex Terrain Areas

The TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X) bistatic system provides high-resolution and high-quality interferometric data for global topographic measurement. Since the twin TanDEM-X satellites fly in a close helix formation, they can acquire approximately simultaneous synthetic aperture radar (SAR) images, so that temporal decorrelation and atmospheric delay can be ignored. Consequently, the orbital error becomes the most significant error limiting high-resolution SAR interferometry (InSAR) applications, such as the high-precision digital elevation model (DEM) reconstruction, subway and highway deformation monitoring, landslide monitoring and sub-canopy topography inversion. For rugged mountainous areas, in particular, it is difficult to estimate and correct the orbital phase error in TanDEM-X bistatic InSAR. Based on the rigorous InSAR geometric relationship, the orbital phase error can be attributed to the baseline errors (BEs) after fixing the positions of the master SAR sensor and the targets on the ground surface. For the constraint of the targets at a study scene, the freely released TanDEM-X DEM can be used, due to its consistency with the TanDEM-X bistatic InSAR-measured height. As a result, a parameterized model for the orbital phase error estimation is proposed in this paper. In high-resolution and high-precision TanDEM-X bistatic InSAR processing, due to the limited precision of the navigation systems and the uneven baseline changes caused by the helix formation, the BEs are time-varying in most cases. The parameterized model is thus built and estimated along each range line. To validate the proposed method, two mountainous test sites located in China (i.e., Fuping in Shanxi province and Hetang in Hunan province) were selected. The obtained results show that the orbital phase errors of the bistatic interferograms over the two test sites are well estimated. Compared with the widely applied polynomial model, the residual phase corrected by the proposed method contains little undesirable topography-dependent phase error, and avoids unexpected height errors ranging about from −6 m to 3 m for the Fuping test site and from −10 m to 8 m for the Hetang test site. Furthermore, some fine details, such as ridges and valleys, can be clearly identified after the correction. In addition, the two components of the orbital phase error, i.e., the residual flat-earth phase error and the topographic phase error caused by orbital error, are separated and quantified based on the parameterized expression. These demonstrate that the proposed method can be used to accurately estimate and mitigate the orbital phase error in TanDEM-X bistatic InSAR data, which increases the feasibility of reconstructing high-resolution and high-precision DEM. The rigorous geometric constraint, the refinement of the initial baseline parameters, and the assessment for height errors based on the estimated BEs are investigated in the discussion section of this paper.

Enhancing Coherence Images for Coherent Change Detection: An Example on Vehicle Tracks in Airborne SAR Images

In Synthetic Aperture Radar (SAR) interferometry, one of the most widely used measures for the quality of the interferometric phase is coherence. However, in favorable conditions coherence can also be used to detect subtle changes on the ground, which are not visible in the amplitude images. For such applications, i.e., coherent change detection, it is important to have a good contrast between the unchanged (high-coherence) parts of the scene and the changed (low-coherence) parts. In this paper, an algorithm is introduced that aims at enhancing this contrast. The enhancement is achieved by a combination of careful filtering of the amplitude images and the interferometric phase image. The algorithm is applied to an airborne interferometric SAR image pair recorded by the SmartRadar experimental sensor of Hensoldt Sensors GmbH. The data were recorded during a measurement campaign over the Bann B installations of POLYGONE Range in southern Rhineland-Palatinate (Germany), with a time gap of approximately four hours between the overflights. In-between the overflights, several vehicles were moved on the site and the goal of this work is to enhance the coherence image such that the tracks of these vehicles can be detected as completely as possible in an automated way. Several coherence estimation schemes found in the literature are explored for the enhancement, as well as several commonly used speckle filters. The results of these filtering steps are evaluated visually and quantitatively, showing that the mean gray-level difference between the low-coherence tracks and their high-coherence surroundings could be enhanced by at least 28%. Line extraction is then applied to the best enhancement. The results show that the tracks can be detected much more completely using the coherence contrast enhancement scheme proposed in this paper.

Application of InSAR satellite method for mapping of active landslides in Bulgaria – opportunities and perspectives

Landslides are geological phenomena that are spread on Bulgarian territory mainly along the northern Black Sea coast and on the right banks of the Danube in the western part of the country. Mitigation of the negative effects of these destructive geological phenomena is the compilation of inventory maps of their distribution and registers with the main characteristics of the individual landslides. Conventional methods for making such maps are time-consuming and resource-intensive. Modern satellite, air and ground-based remote sensing technologies facilitate the production of landslide maps, reducing the time and resources required to compile and systematically update them. In this paper, we demonstrate the applicability of Differential Sentinel-1A satellite SAR interferometry (DInSAR) to assess the movement activity and use the information for further updating the national landslide inventories in Bulgaria. We perform several analyses based on multi-temporal InSAR techniques of Sentinel-1A data over selected areas prone to landslides. The use of new opportunities for free access to satellite images, which can be applied in conjunction with other methods, greatly facilitates the processes of inventory, mapping and study of landslides.

Antarctic snow-covered sea ice topography derivation from TanDEM-X using polarimetric SAR interferometry

Single-pass interferometric synthetic aperture radar (InSAR) enables the possibility for sea ice topographic retrieval despite the inherent dynamics of sea ice. InSAR digital elevation models (DEMs) are measuring the radar scattering center height. The height bias induced by the penetration of electromagnetic waves into snow and ice leads to inaccuracies of the InSAR DEM, especially for thick and deformed sea ice with snow cover. In this study, an elevation difference between the satellite-measured InSAR DEM and the airborne-measured optical DEM is observed from a coordinated campaign over the western Weddell Sea in Antarctica. The objective is to correct the penetration bias and generate a precise sea ice topographic map from the single-pass InSAR data. With the potential of retrieving sea ice geophysical information by the polarimetric-interferometry (Pol-InSAR) technique, a two-layer-plus-volume model is proposed to represent the sea ice vertical structure and its scattering mechanisms. Furthermore, a simplified version of the model is derived, to allow its inversion with limited a priori knowledge, which is then applied to a topographic retrieval scheme. The experiments are performed across four polarizations: HH, VV, Pauli 1 (HH + VV), and Pauli 2 (HH − VV). The model-retrieved performance is validated with the optically derived DEM of the sea ice topography, showing an excellent performance with root-mean-square error as low as 0.26 m in Pauli-1 (HH + VV) polarization.