Advanced pixel selection and optimization algorithms for Persistent Scatterer Interferometry (PSI)

Tesi amb diferents seccions retallades per dret de l'editorPremi Extraordinari de Doctorat, promoció 2018-2019. Àmbit de les TICGround deformation measurements can provide valuable information for minimization of associated loss and damage caused by natural and environmental hazards. As a kind of remote sensing technique, Persistent Scatterer Interferometry (PSI) SAR is able to measure ground deformation with high spatial resolution, efficiently. Moreover, the ground deformation monitoring accuracy of PSI techniques can reach up to millimeter level. However, low coherence could hinderthe exploitation of SAR data, and high-accuracy deformation monitoring can only be achieved by PSI for high quality pixels. Therefore, pixel optimization and identification of coherent pixels are crucial for PSI techniques. In this thesis, advanced pixel selection and optimization algorithms have been investigated. Firstly, a full-resolution pixel selection method based on the Temporal Phase Coherence (TPC) has been proposed. This method first estimates noise phase term of each pixel at interferogram level. Then, for each pixel, its noise phase terms of all interferograms are used to assess this pixel’s temporal phase quality (i.e., TPC). In the next, based on the relationship between TPC and phase Standard Deviation (STD), a threshold can be posed on TPC to identify high phase quality pixels. This pixel selection method can work with both Deterministic Scatterers (PSs) and Distributed Scatterers (DSs). To valid the effectiveness of the developed method, it has been used to monitor the Canillo (Andorra) landslide. The results show that the TPC method can obtained highest density of valid pixels among the employed three approaches in this challenging area with X-band SAR data. Second, to balance the polarimetric DInSAR phase optimization effect and the computation cost, a new PolPSI algorithm is developed. This proposed PolPSI algorithm is based on the Coherency Matrix Decomposition result to determine the optimal scattering mechanism of each pixel, thus it is named as CMD-PolPSI. CMDPolPSI need not to search for solution within the full space of solution, it is therefore much computationally faster than the classical Equal Scattering Mechanism (ESM) method, but with lower optimization performance. On the other hand, its optimization performance outperforms the less computational costly BEST method. Third, an adaptive algorithm SMF-POLOPT has been proposed to adaptive filtering and optimizing PolSAR pixels for PolPSI applications. This proposed algorithm is based on PolSAR classification results to firstly identify Polarimetric Homogeneous Pixels (PHPs) for each pixel, and at the same time classify PS and DS pixels. After that, DS pixels are filtered by their associated PHPs, and then optimized based on the coherence stability phase quality metric; PS pixels are unfiltered and directly optimized based on the DA phase quality metric. SMF-POLOPT can simultaneously reduce speckle noise and retain structures’ details. Meanwhile, SMF-POLOPT is able to obtain much higher density of valid pixels for deformation monitoring than the ESM method. To conclude, one pixel selection method has been developed and tested, two PolPSI algorithms have been proposed in this thesis. This work make contributions to the research of “Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer InterferometryLes mesures de deformació del sòl poden proporcionar informació valuosa per minimitzar les pèrdues i els danys associats causats pels riscos naturals i ambientals. Com a tècnica de teledetecció, la interferometria de dispersors persistents (Persistent Scatter Interferometry, PSI) SAR és capaç de mesurar de forma eficient la deformació del terreny amb una alta resolució espacial. A més, la precisió de monitorització de la deformació del sòl de les tècniques PSI pot arribar a arribar a nivells del mil·límetre. No obstant això, una baixa coherència pot dificultar l’explotació de dades SAR i el control de deformació d’alta precisió només es pot aconseguir mitjançant PSI per a píxels d’alta qualitat. Per tant, l’optimització de píxels i la identificació de píxels coherents són crucials en les tècniques PSI. En aquesta tesi s¿han investigat algorismes avançats de selecció i optimització de píxels. En primer lloc, s'ha proposat un mètode de selecció de píxels de resolució completa basat en la coherència temporal de fase (Temporal Phase Coherence, TPC). Aquest mètode estima per primera vegada el terme de fase de soroll de cada píxel a nivell d’interferograma. A continuació, per a cada píxel, s'utilitzen els termes de la fase de soroll de tots els interferogrames per avaluar la qualitat de fase temporal d'aquest píxel (és a dir, TPC). A la següent, basant-se en la relació entre el TPC i la desviació estàndard de fase (STD), es pot plantejar un llindar de TPC per identificar píxels de qualitat de fase alta. Aquest mètode de selecció de píxels es capaç de detectar tant els dispersors deterministes (PS) com els distribuïts (DS). Per validar l’eficàcia del mètode desenvolupat, s’ha utilitzat per controlar l’esllavissada de Canillo (Andorra). Els resultats mostren que el mètode TPC pot obtenir la major densitat de píxels vàlids, comparat amb els mètodes clàssics de selecció, en aquesta àrea difícil amb dades de SAR de banda X. En segon lloc, per equilibrar l’efecte d’optimització de fase DInSAR polarimètrica i el cost de càlcul, es desenvolupa un nou algorisme de PolPSI. Aquest algorisme proposat de PolPSI es basa en el resultat de la descomposició de la matriu de coherència per determinar el mecanisme de dispersió òptim de cada píxel, de manera que es denomina CMD-PolPSI. CMDPolPSI no necessita buscar solucions dins de l’espai complet de la solució, per tant, és molt més eficient computacionalment que el mètode clàssic de mecanismes d’igualtat de dispersió (Equal Scattering Mechanism, ESM), però amb un efecte d’optimització no tant òptim. D'altra banda, el seu efecte d'optimització supera el mètode BEST, el que te un menor cost computacional. En tercer lloc, s'ha proposat un algoritme adaptatiu SMF-POLOPT per al filtratge adaptatiu i l'optimització de píxels PolSAR per a aplicacions PolPSI. Aquest algorisme proposat es basa en els resultats de classificació PolSAR per identificar primer els píxels homogenis polarimètrics (PHP) per a cada píxel i, alhora, classificar els píxels PS i DS. Després d'això, els píxels DS es filtren pels seus PHP associats i, a continuació, s'optimitzen en funció de la mètrica de qualitat de la fase d'estabilitat de coherència; els píxels classificats com PS no es filtren i s'optimitzen directament en funció de la mètrica de qualitat de la fase DA. SMF-POLOPT pot reduir simultàniament el soroll de la fase interferomètrica i conservar els detalls de les estructures. Mentrestant, SMF-POLOPT aconsegueix obtenir una densitat molt més alta de píxels vàlids per al seguiment de la deformació que el mètode ESM. Per concloure, en aquesta tesi s’ha desenvolupat i provat un mètode de selecció de píxels, i s’han proposat dos algoritmes PolPSI. Aquest treball contribueix a la recerca en "Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer Interferometry"Postprint (published version

Basin scale assessment of landslides geomorphological setting by advanced InSAR analysis

An extensive investigation of more than 90 landslides affecting a small river basin in Central Italy was performed by combining field surveys and remote sensing techniques. We thus defined the geomorphological setting of slope instability processes. Basic information, such as landslides mapping and landslides type definition, have been acquired thanks to geomorphological field investigations and multi-temporal aerial photos interpretation, while satellite SAR archive data (acquired by ERS and Envisat from 1992 to 2010) have been analyzed by means of A-DInSAR (Advanced Differential Interferometric Synthetic Aperture Radar) techniques to evaluate landslides past displacements patterns. Multi-temporal assessment of landslides state of activity has been performed basing on geomorphological evidence criteria and past ground displacement measurements obtained by A-DInSAR. This step has been performed by means of an activity matrix derived from information achieved thanks to double orbital geometry. Thanks to this approach we also achieved more detailed knowledge about the landslides kinematics in time and space

Building deformation assessment by means of Persistent Scatterer Interferometry analysis on a landslide-affected area: the Volterra (Italy) case study

In recent years, space-borne InSAR (interferometric synthetic aperture radar) techniques have shown their capabilities to provide precise measurements of Earth surface displacements for monitoring natural processes. Landslides threaten human lives and structures, especially in urbanized areas, where the density of elements at risk sensitive to ground movements is high. The methodology described in this paper aims at detecting terrain motions and building deformations at the local scale, by means of satellite radar data combined with in situ validation campaigns. The proposed approach consists of deriving maximum settlement directions of the investigated buildings from displacement data revealed by radar measurements and then in the cross-comparison of these values with background geological data, constructive features and on-field evidence. This validation permits better understanding whether or not the detected movements correspond to visible and effective damages to buildings. The method has been applied to the southwestern sector of Volterra (Tuscany region, Italy), which is a landslide-affected and partially urbanized area, through the use of COSMO-SkyMed satellite images as input data. Moreover, we discuss issues and possible misinterpretations when dealing with PSI (Persistent Scatterer Interferometry) data referring to single manufactures and the consequent difficulty of attributing the motion rate to ground displacements, rather than to structural failures

Railways’ Stability Observation by Satellite Radar Images

Remote sensing has many vital civilian applications. Space-borne Interferometric Synthetic Aperture Radar has been used to measure the Earth’s surface deformation widely. In particular, Persistent Scatterer Interferometry (PSI) is designed to estimate the temporal characteristics of the Earth’s deformation rates from multiple InSAR images acquired over time. This chapter reviews the space-borne Differential Interferometric Synthetic Aperture Radar techniques that have shown their capabilities in monitoring of railways displacements. After description of the current state of the art and potentials of the available radar remote sensing techniques, one case study is examined, pertaining to a railway bridge in the Campania region, Italy

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

Leveraging power of deep learning for fast and efficient elite pixel selection in time series SAR interferometry

This work proposes an improved convolutional long short-term memory (ConvLSTM) based architecture for selection of elite pixels (i.e., less noisy) in time series interferometric synthetic aperture radar (TS-InSAR). Compared to previous version, the model can process InSAR stacks of variable time steps and select both persistent (PS) and distributed scatterers (DS). We trained the model on ~20,000 training images (interferograms), each of size 100 by 100 pixels, extracted from InSAR time series interferograms containing both artificial features (buildings and infrastructure) and objects of natural environment (vegetation, forests, barren or agricultural land, water bodies). Based on such categorization, we developed two deep learning models, primarily focusing on urban and coastal sites. Training labels were generated from elite pixel selection outputs generated from the wavelet-based InSAR (WabInSAR) software developed by Shirzaei (2013) and improved in Lee and Shirzaei (2023). With 4 urban and 7 coastal sites used for training and validation, the predicted elite pixel selection maps reveal that the proposed models efficiently learn from WabInSAR-generated labels, reaching a validation accuracy of 94%. The models accurately discard pixels affected by geometric and temporal decorrelation while selecting pixels corresponding to urban objects and those with stable phase history unaffected by temporal and geometric decorrelation. The density of pixels in urban areas is comparable to and higher for coastal areas compared to WabInSAR outputs. With significantly reduced time computation (order of minutes) and improved selection of elite pixels, the proposed models can efficiently process long InSAR time series stacks and generate rapid deformation maps

Satellite Monitoring of Railways using Interferometric Synthetic Aperture Radar (InSAR)

There is over 15,600 km of track in the Swedish railroad network. This network is vital for the transportation of people and goods across the country. It is important that this network is monitored and maintained to ensure good function and safety. A tool for monitoring and measuring ground deformation over a large area remotely with high frequency and accuracy was developed in recent decades. This tool is known as Interferometric Synthetic Aperture Radar (InSAR), and is used by researchers, geo-technicians, and engineers. The purpose of this study has been to evaluate the use and feasibility of the InSAR technique for track condition monitoring and compare it to conventional track condition monitoring techniques. Malmbanan, which is primarily used to transport iron-ore from mines in Sweden to the ports of Luleå, Sweden and Narvik, Norway, is used as a case study for this project; specifically, the section between Kiruna and Riksgränsen. Coordinate matching of measurements from the provided Persistent Scatterer Interferometry (PSI) InSAR data and Optram data from survey trains were performed. Then measured changes over different time spans within the two systems were overlapped and classified with different thresholds to see if there is correlation between the two systems. An extensive literature review was also conducted in order to gain an understanding of InSAR technologies and uses.The literature review showed that there is a large potential and a quickly growing number of applications of InSAR to monitor railways and other types of infrastructure, and that the tools and algorithms for this are being improved. The case study, on the other hand, shows that it can be difficult to directly compare measurement series from different tools, each working on different resolutions in terms of both time and space. InSAR is thus not about to replace techniques such as those behind Optram (using measurement trains). Instead, the approaches offer complementary perspectives, each highlighting different types of issues. We find that InSAR offers a good way to identify locations with settlements or other types of ground motions. Especially transition zones between settlements and more stable ground can be challenging from a maintenance point of view and can clearly be identified and monitored using InSAR. With the rollout of national InSAR-data, and the large increase in data accessibility, we see a considerable potential for future studies that apply the technique to the railway area

Application of Differential and Polarimetric Synthetic Aperture Radar (SAR) Interferometry for Studying Natural Hazards

In the following work, I address the problem of coherence loss in standard Differential Interferometric SAR (DInSAR) processing, which can result in incomplete or poor quality deformation measurements in some areas. I incorporate polarimetric information with DInSAR in a technique called Polarimetric SAR Interferometry (PolInSAR) in order to acquire more accurate and detailed maps of surface deformation. In Chapter 2, I present a standard DInSAR study of the Ahar double earthquakes (Mw=6.4 and 6.2) which occurred in northwest Iran, August 11, 2012. The DInSAR coseismic deformation map was affected by decorrelation noise. Despite this, I employed an advanced inversion technique, in combination with a Coulomb stress analysis, to find the geometry and the slip distribution on the ruptured fault plane. The analysis shows that the two earthquakes most likely occurred on a single fault, not on conjugate fault planes. This further implies that the minor strike-slip faults play more significant role in accommodating convergence stress accumulation in the northwest part of Iran. Chapter 3 presents results from the application of PolInSAR coherence optimization on quad-pol RADARSAT-2 images. The optimized solution results in the identification of a larger number of reliable measurement points, which otherwise are not recognized by the standard DInSAR technique. I further assess the quality of the optimized interferometric phase, which demonstrates an increased phase quality with respect to those phases recovered by applying standard DInSAR alone. Chapter 4 discusses results from the application of PolInSAR coherence optimization from different geometries to the study of creep on the Hayward fault and landslide motions near Berkeley, CA. The results show that the deformation rates resolved by PolInSAR are in agreement with those of standard DInSAR. I also infer that there is potential motion on a secondary fault, northeast and parallel to the Hayward fault, which may be creeping with a lower velocity

Satellite SAR Interferometry for Earth’s Crust Deformation Monitoring and Geological Phenomena Analysis

Synthetic aperture radar interferometry (InSAR) and the related processing techniques provide a unique tool for the quantitative measurement of the Earth’s surface deformation associated with certain geophysical processes (such as volcanic eruptions, landslides and earthquakes), thus making possible long-term monitoring of surface deformation and analysis of relevant geodynamic phenomena. This chapter provides an application-oriented perspective on the spaceborne InSAR technology with emphasis on subsequent geophysical investigations. First, the fundamentals of radar interferometry and differential interferometry, as well as error sources, are briefly introduced. Emphasis is then placed on the realistic simulation of the underlying geophysics processes, thus offering an unfolded perspective on both analytical and numerical approaches for modeling deformation sources. Finally, various experimental investigations conducted by acquiring SAR multitemporal observations on areas subject to deformation processes of particular geological interest are presented and discussed