Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review

Interferometric Synthetic Aperture Radar (InSAR) is a new measurement technology, making use of the phase information contained in the Synthetic Aperture Radar (SAR) images. InSAR has been recognized as a potential tool for the generation of digital elevation models (DEMs) and the measurement of ground surface deformations. However, many critical factors affect the quality of InSAR data and limit its applications. One of the factors is InSAR data processing, which consists of image co-registration, interferogram generation, phase unwrapping and geocoding. The co-registration of InSAR images is the first step and dramatically influences the accuracy of InSAR products. In this paper, the principle and processing procedures of InSAR techniques are reviewed. One of important factors, tie points, to be considered in the improvement of the accuracy of InSAR image co-registration are emphatically reviewed, such as interval of tie points, extraction of feature points, window size for tie point matching and the measurement for the quality of an interferogram

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Monitoring von Hangbewegungen mit InSAR Techniken im Gebiet Ciloto, Indonesien

In this doctoral thesis, the InSAR techniques are applied to detect the ground movement phenomenon and to assess the InSAR result geometrically in the Ciloto area, Indonesia. Mainly, one of those techniques, the SB-SDFP algorithm, overcomes the limitations of conventional InSAR in monitoring rural and agricultural areas and can observe extremely slow landslides. The InSAR strategy is positively known as a promising option to detect and quantify the kinematics of active landslides on a large areal scale. To minimize the bias of the InSAR displacement result, the correction of the tropospheric phase delay was carried out in a first step. This procedure is demonstrated in experiments both in the small study area in Ciloto and in a larger area. The latter is an area located in Northern Baja California, Mexico and is dominated by tectonic activity as well as groundwater-induced subsidence. A detailed investigation of the slope movement's behavior in the Ciloto district was conducted utilizing multi-temporal and multi-band SAR data from ERS1/2 (1996-1999), ALOS PALSAR (2007-2009) and Sentinel-1 (2014-2018) satellites. The region was successfully identified as a permanent active landslide prone area, especially in the vicinity of the Puncak Pass and Puncak Highway. The full 3D velocity field and the displacement time series were estimated using the inversion model. The velocity rate was classified from extremely slow to slow movement. To comprehend the landslide's behavior, a further examination of the relationship between InSAR results and physical characteristics of the area was carried out. For the long period of a slow-moving landslide, the relationship between precipitation and displacement trend shows a weak correlation. It is concluded that the extremely slow to slow deformation is not directly influenced by the rainfall intensity, yet it effectuates the subsurface and the groundwater flow. The run-off process with rainfall exceeding a soil's infiltration capacity was suspected as the main driver of the slow ground movement phenomenon. However, when analyzing rapid and extremely rapid landslide events at Puncak Pass, a significant increase in the correlation coefficient between precipitation and displacement rate could be observed.In dieser Doktorarbeit wird die Anwendung von erweiterten Verarbeitungsstrategien von InSAR Daten zur Erkennung und geometrischen Bewertung der Bodenbewegungen im Ciloto - Indonesien dargestellt. Dieser Ansatz überwindet die Beschränkungen konventioneller SAR-Interferometrie und ermöglicht sowohl ein kontinuierliches Monitoring dieses landwirtschaftich geprägten Gebietes als auch die Erfassung extrem langsamer Hangrutschungen. Um eine Verzerrung der InSAR Deformationsergebnisse zu minimieren, wurde zunächst eine Korrektur der troposphärischen Phase durchgeführt. Diese neuartige Strategie wird sowohl im Forschungsgebiet Ciloto als auch an einem größeren Gebiet demonstriert. Bei letzterem handelt es sich um einen Küstenstreifen im nördlichen Niederkalifornien, Mexiko, welcher durch hohe tektonische Aktivität und grundwasserinduzierte Landsetzungen charakterisiert ist. Die detaillierte Untersuchung des Verhaltens von Hangrutschungen im Ciloto erfolgte durch die Verarbeitung multi-temporaler SAR-Daten unter Nutzung verschiedener Frequenzbänder, darunter ESR1/2 (1996-1999), ALOS PALSAR (2007-2009) und Sentinel-1 (2014-2018) Daten. Die Region konnte erfolgreich als permanent aktives Hangrutschungsgebiet identifiziert werden, wobei der Puncak Pass und der Puncak Highway ein erhöhtes Gefahrenpotential aufweisen. Ein 3D- Geschwindig-keitsfeld der Deformation und die zugehörigen Zeitreihen wurden mit dem Inversionsmodell berechnet. Die Geschwindigkeitsrate wurde als langsam bis extrem langsam klassifiziert. Um das dynamische Verhalten der Hangrutschung zu verstehen wurde, in einer weiteren Untersuchung die Beziehung zwischen dem InSAR-Ergebnis und den physikalischen Begebenheiten im Forschungsgebiet analysiert. Es wird der Schluss gezogen, dass die langsame bis extrem langsame Verformung nicht direkt von der Niederschlagsintensität beeinflusst wird, diese sich aber auf den Untergrund und die Grundwasserströmung auswirkt. Es wird vermutet, dass der Oberflächenablauf, welcher die Infiltrationskapazität des Bodens übersteigt, ausschlaggebend für das Phänomen der langsamen Bodenbewegung ist. Für die schnellen und extrem schnellen Hangrutschungen jedoch konnte eine signifikante Erhöhung des Korrelationskoeffizienten zwischen Niederschlag und Verschiebungsrate bei Untersuchungen der Hangrutschung am Puncak-Pass nachgewiesen werden

Split-Band Interferometric SAR Processing Using TanDEM-X Data

Most recent SAR sensors use wide band signals to achieve metric range resolution. One can also take advantage of wide band to split it into sub-bands and generate several lower-resolution images, centered on slightly different frequencies, from a single acquisition. This process, named Multi Chromatic Analysis (MCA) corresponds to performing a spectral analysis of SAR images. Split-Band SAR interferometry (SBInSAR) is based on spectral analysis performed on each image of an InSAR pair, yielding a stack of sub-band interferograms. Scatterers keeping a coherent behaviour in each subband interferogram show a phase that varies linearly with the carrier frequency, the slope being proportional to the absolute optical path difference. This potentially solves the problems of phase unwrapping on a pixelper-pixel basis. In this paper, we present an SBInSAR processor and its application using TanDEM-X data over the Nyiragongo volcano.Fil: Derauw, Dominique Maurice. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Patagonia Norte. Instituto de Investigación en Paleobiología y Geología; Argentina. Centre Spatial de Liège; BélgicaFil: Kervyn, François. Royal Museum Of Central Africa; BélgicaFil: d'Oreye, Nicolas. European Centre For Geodynamics And Seismology; Luxemburgo. National Museum of Natural History; LuxemburgoFil: Smets, Benoit. Royal Museum Of Central Africa; Bélgica. European Centre For Geodynamics And Seismology; Luxemburgo. Vrije Unviversiteit Brussel; BélgicaFil: Albino, Fabien. Royal Museum Of Central Africa; BélgicaFil: Barbier, Christian. Centre Spatial de Liège; BélgicaAdvances in the Science and Applications of SAR Interferometry and Sentinel-1 InSAR WorkshopFrascatiItaliaEuropean Space Agenc

Application of Synthetic Aperture Radar Interferometry (InSAR) in defining groundwater-withdrawal-related subsidence, Diamond Valley, Nevada

Interferometric Synthetic Aperture Rader (InSAR) technique has been recently used in detecting and monitoring ground displacements such as volcanic activities, earthquakes, landslides and surface deformations caused by fluid extraction. Ground subsidence related to groundwater withdrawal has occurred in many places such as urban areas and large scale agricultural areas. This study utilizes interferometry to detect ground subsidence in an agricultural field, where a large volume of groundwater has been pumped for decades, in Diamond Valley, Nevada. InSAR has proven to have great potential to detect and quantify ground subsidence caused by aquifer system compaction. It mapped ground deformation signals with high spatial detail and resolution of displacement, developed in a groundwater basin in the area, using radar data collected from the ERS-1/ERS-2 and Envisat satellites. The subsidence signal at the south part of the valley, where irrigation wells exist, shows a minimum of 37.6 cm of cumulative subsidence between July 17, 1992 and November 27, 1999 and a 17.5 cm of cumulative subsidence between October 16, 2004 and December 15, 2007.The profile views of the subsidence signals assist in visualizing the deformation geometry which indicates that the subsurface lithology can increase or decrease the deformation. The subsurface model estimated from the history of water table decline and the subsurface lithology distribution approximately correlate with the subsidence signals. Poor correlation occurred where limited availability of good subsurface data and limited spatial coverage of well logs existed

Ground-based synthetic aperture radar (GBSAR) interferometry for deformation monitoring

Ph. D ThesisGround-based synthetic aperture radar (GBSAR), together with interferometry, represents a powerful tool for deformation monitoring. GBSAR has inherent flexibility, allowing data to be collected with adjustable temporal resolutions through either continuous or discontinuous mode. The goal of this research is to develop a framework to effectively utilise GBSAR for deformation monitoring in both modes, with the emphasis on accuracy, robustness, and real-time capability. To achieve this goal, advanced Interferometric SAR (InSAR) processing algorithms have been proposed to address existing issues in conventional interferometry for GBSAR deformation monitoring. The proposed interferometric algorithms include a new non-local method for the accurate estimation of coherence and interferometric phase, a new approach to selecting coherent pixels with the aim of maximising the density of selected pixels and optimizing the reliability of time series analysis, and a rigorous model for the correction of atmospheric and repositioning errors. On the basis of these algorithms, two complete interferometric processing chains have been developed: one for continuous and the other for discontinuous GBSAR deformation monitoring. The continuous chain is able to process infinite incoming images in real time and extract the evolution of surface movements through temporally coherent pixels. The discontinuous chain integrates additional automatic coregistration of images and correction of repositioning errors between different campaigns. Successful deformation monitoring applications have been completed, including three continuous (a dune, a bridge, and a coastal cliff) and one discontinuous (a hillside), which have demonstrated the feasibility and effectiveness of the presented algorithms and chains for high-accuracy GBSAR interferometric measurement. Significant deformation signals were detected from the three continuous applications and no deformation from the discontinuous. The achieved results are justified quantitatively via a defined precision indicator for the time series estimation and validated qualitatively via a priori knowledge of these observing sites.China Scholarship Council (CSC), Newcastle Universit

Ionospheric correction of interferometric SAR data with application to the cryospheric sciences

Thesis (Ph.D.) University of Alaska Fairbanks, 2018The ionosphere has been identified as an important error source for spaceborne Synthetic Aperture Radar (SAR) data and SAR Interferometry (InSAR), especially for low frequency SAR missions, operating, e.g., at L-band or P-band. Developing effective algorithms for the correction of ionospheric effects is still a developing and active topic of remote sensing research. The focus of this thesis is to develop robust and accurate techniques for ionospheric correction of SAR and InSAR data and evaluate the benefit of these techniques for cryospheric research fields such as glacier ice velocity tracking and permafrost deformation monitoring. As both topics are mostly concerned with high latitude areas where the ionosphere is often active and characterized by turbulence, ionospheric correction is particularly relevant for these applications. After an introduction to the research topic in Chapter 1, Chapter 2 will discuss open issues in ionospheric correction including processing issues related to baseline-induced spectrum shifts. The effect of large baseline on split spectrum InSAR technique has been thoroughly evaluated and effective solutions for compensating this effect are proposed. In addition, a multiple sub-band approach is proposed for increasing the algorithm robustness and accuracy. Selected case studies are shown with the purpose of demonstrating the performance of the developed algorithm. In Chapter 3, the developed ionospheric correction technology is applied to optimize InSAR-based ice velocity measurements over the big ice sheets in Greenland and the Antarctic. Selected case studies are presented to demonstrate and validate the effectiveness of the proposed correction algorithms for ice velocity applications. It is shown that the ionosphere signal can be larger than the actual glacier motion signal in the interior of Greenland and Antarctic, emphasizing the necessity for operational ionospheric correction. The case studies also show that the accuracy of ice velocity estimates was significantly improved once the developed ionospheric correction techniques were integrated into the data processing flow. We demonstrate that the proposed ionosphere correction outperforms the traditionally-used approaches such as the averaging of multi-temporal data and the removal of obviously affected data sets. For instance, it is shown that about one hundred multi-temporal ice velocity estimates would need to be averaged to achieve the estimation accuracy of a single ionosphere-corrected measurement. In Chapter 4, we evaluate the necessity and benefit of ionospheric-correction for L-band InSAR-based permafrost research. In permafrost zones, InSAR-based surface deformation measurements are used together with geophysical models to estimate permafrost parameters such as active layer thickness, soil ice content, and permafrost degradation. Accurate error correction is needed to avoid biases in the estimated parameters and their co-variance properties. Through statistical analyses of a large number of L-band InSAR data sets over Alaska, we show that ionospheric signal distortions, at different levels of magnitude, are present in almost every InSAR dataset acquired in permafrost-affected regions. We analyze the ionospheric correction performance that can be achieved in permafrost zones by statistically analyzing correction results for large number of InSAR data. We also investigate the impact of ionospheric correction on the performance of the two main InSAR approaches that are used in permafrost zones: (1) we show the importance of ionospheric correction for permafrost deformation estimation from discrete InSAR observations; (2) we demonstrate that ionospheric correction leads to significant improvements in the accuracy of time-series InSAR-based permafrost products. Chapter 5 summarizes the work conducted in this dissertation and proposes next steps in this field of research

InSAR Modeling of Geophysics Measurements

In the present work, the geometry and basic parameters of interferometric synthetic aperture radar (InSAR) geophysics system are addressed. Equations of pixel height and displacement evaluation are derived. Synthetic aperture radar (SAR) signal model based on linear frequency modulation (LFM) waveform and image reconstruction procedure are suggested. The concept of pseudo InSAR measurements, interferogram, and differential interferogram generation is considered. Interferogram and differential interferogram are generated based on a surface model and InSAR measurements. Results of numerical experiments are provided