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

Landslide monitoring using multi-temporal SAR interferometry with advanced persistent scatterers identification methods and super high-spatial resolution terraSAR-X Images

Landslides are one of the most common and dangerous threats in the world that generate considerable damage and economic losses. An efficient landslide monitoring tool is the Differential Synthetic Aperture Radar Interferometry (DInSAR) or Persistent Scatter Interferometry (PSI). However, landslides are usually located in mountainous areas and the area of interest can be partially or even heavily vegetated. The inherent temporal decorrelation that dramatically reduces the number of Persistent Scatters (PSs) of the scene limits in practice the application of this technique. Thus, it is crucial to be able to detect as much PSs as possible that can be usually embedded in decorrelated areas. High resolution imagery combined with efficient pixel selection methods can make possible the application of DInSAR techniques in landslide monitoring. In this paper, different strategies to identify PS Candidates (PSCs) have been employed together with 32 super high-spatial resolution (SHR) TerraSAR-X (TSX) images, staring-spotlight mode, to monitor the Canillo landslide (Andorra). The results show that advanced PSI strategies (i.e., the temporal sub-look coherence (TSC) and temporal phase coherence (TPC) methods) are able to obtain much more valid PSs than the classical amplitude dispersion (DA) method. In addition, the TPC method presents the best performance among all three full-resolution strategies employed. The SHR TSX data allows for obtaining much higher densities of PSs compared with a lower-spatial resolution SAR data set (Sentinel-1A in this study). Thanks to the huge amount of valid PSs obtained by the TPC method with SHR TSX images, the complexity of the structure of the Canillo landslide has been highlighted and three different slide units have been identified. The results of this study indicate that the TPC approach together with SHR SAR images can be a powerful tool to characterize displacement rates and extension of complex landslides in challenging areasPeer ReviewedPostprint (published version

Urban Deformation Monitoring using Persistent Scatterer Interferometry and SAR tomography

This book focuses on remote sensing for urban deformation monitoring. In particular, it highlights how deformation monitoring in urban areas can be carried out using Persistent Scatterer Interferometry (PSI) and Synthetic Aperture Radar (SAR) Tomography (TomoSAR). Several contributions show the capabilities of Interferometric SAR (InSAR) and PSI techniques for urban deformation monitoring. Some of them show the advantages of TomoSAR in un-mixing multiple scatterers for urban mapping and monitoring. This book is dedicated to the technical and scientific community interested in urban applications. It is useful for choosing the appropriate technique and gaining an assessment of the expected performance. The book will also be useful to researchers, as it provides information on the state-of-the-art and new trends in this fiel

Multi-temporal landslide activity investigation by spaceborne SAR interferometry: The case study of the Polish Carpathians

The main goal of this research is to verify the activity state of landslides provided by an existing landslide inventory map using Persistent Scatterers (PS) Interferometry (PSInSAR). The study was conducted in the Małopolskie municipality, a rural setting with sparse urbanization in the Polish Flysch Carpathians. PSInSAR has been applied using Synthetic Aperture Radar (SAR) data from ALOS PALSAR and Sentinel 1A/B with different acquisition geometries (ascending and descending orbit) to increase PS coverage and mitigate the geometric effects due to layover and shadowing. The Line-Of-Sight PSInSAR measurements were projected to the steepest slope, which allowed to homogenize the results from diverse acquisition modes and to compare the displacement velocities with different slope orientations. Additionally, landslide intensity (motion rate) and expected damage maps were generated and verified during field investigations. A high correlation between PSInSAR results and in-situ damage observations was confirmed. The activity state and landslide-related expected damage maps have been confirmed for 43 out of a total of 50 landslides investigated in the field. The short temporal baseline provided by both Sentinel satellites (1A/B data) increases the PS density significantly. The study substantiates the usefulness of SAR based landslide activity monitoring for land use and land development, even in rural areas

Advanced exploitation of Sentinel-1 data for supporting landslide risk analysis

Tesi en modalitat de compendi de publicacionsSatellite Synthetic Aperture Radar Interferometry (InSAR) and Persistent Scatterer Interferometry (PSI) are now consolidated tools for ground movement detection and monitoring. Sentinel-1 (S1) is the first satellite providing free data access and ensuring a regular acquisition worldwide, every 6 days, increasing its potential for long-term monitoring applications. Several regional and national ground motion services are already active, providing products based on S1 data. Soon in 2022 the first European Ground Motion Service (EGMS) will be available and freely provide a displacement map over the whole Europe, with annual updates. This implies a strong expansion of availability of PSI-based displacement maps and an easy access for anyone, with an increasing interest among a wider range of users, including public or governmental institutions, academia, industry, and citizens. The analysis and interpretation of this amount of data is difficult and time consuming, mostly for non-expert InSAR users. The objective of this work is developing methodologies to simplify the operational use of PSI displacement maps, generating derived products with a clear message, easy-to-interpret, and fast to read. We propose a method to be applied over regional scale PSI displacement maps, to fast detect the most significant Active Deformation Areas (ADAs). The ADA map is a first product that allows a fast focusing on the active areas, to prioritize further analysis and investigation. Starting from the ADAs, the potential phenomena are attributed to each area through a preliminary interpretation based on auxiliary data, to derive the Geohazard Activity Map. In this work, a methodology to include the ADA information in the Civil Protection Activities is proposed, with the main output called Vulnerable Elements Activity Maps (VEAM). An application of the VEAM is illustrated in the Canary Islands. Furthermore, the ADA map is used in the Valle d'Aosta Region (Northern of Italy) to generate vulnerability and potential loss maps. Finally, a methodology to derive potential damage maps of the exposed buildings, based on the spatial gradients of movement, is proposed, and applied in a coastal area of the Province of Granada (Spain). A pack of software tools has been developed based on the proposed methods to extract ADA and then classify them to generate a Geohazard Activity Map. The set of tools is called ADATools, it is open-access, easy to use and fast, improving the operational exploitation of PSI regional-scale displacement maps. All the methodologies have been developed in the frame of several European projects (Safety, U-Geohaz, MOMIT and RISKCOAST), and are aimed at supporting the multi-scale territorial management and risk analysis activities, with a specific focus on landslides.La interferometría satelital radar (InSAR) y la interferometría de dispersores persistentes (PSI) son herramientas consolidadas para la detección y el monitoreo de movimientos de la superficie de la Tierra. Sentinel-1 (S1) es el primer satélite que proporciona acceso gratuito a los datos y garantiza una adquisición regular en todo el mundo, cada 6 días, aumentando su potencial para aplicaciones de monitoreo a largo plazo. Varios Ground Motion Services regionales y nacionales ya están activos, proporcionando productos basados en datos S1. Pronto, en 2022, el primer servicio europeo (European Ground Motion Service - EGMS) estará disponible y facilitará libremente un mapa de movimientos de toda Europa, con actualizaciones anuales. Esto implica un aumento de la disponibilidad de mapas de movimientos basados en PSI y un fácil acceso para cualquier persona, con un interés creciente entre una amplia gama de usuarios, incluyendo instituciones públicas o gubernamentales, academias, industrias y ciudadanos. El análisis e interpretación de esta cantidad de datos es difícil y consume mucho tiempo, mayormente para usuarios no expertos en la técnica. El objetivo de este trabajo es desarrollar metodologías para simplificar el uso operativo de los mapas de desplazamiento PSI, generando productos derivados con un mensaje claro, fácil de interpretar, y rápido de leer. Se propone un método para detectar rápidamente las Áreas de Deformación Activas (ADAs) más significativas, a partir de mapas de desplazamiento PSI de escala regional. El mapa de las ADAs es un primer producto que permite un enfoque rápido en las áreas activas, útil para priorizar el análisis y las investigaciones adicionales. A partir de las ADAs, se propone una interpretación preliminar basada en datos auxiliares, que atribuye a cada área el fenómeno que está detrás del movimiento, generando el Geohazard Activity Map (GAM). Después, se propone una metodología para incluir la información de las ADAs en las actividades de protección civil, generando los Vulnerable Element Activity Maps (VEAM), a través de su aplicación en las Islas Canarias. Además, el mapa de las ADAs se utiliza en la región de Valle D'Aosta (norte de Italia) para generar mapas de vulnerabilidad y posibles pérdidas económicas. Finalmente, se propone una metodología para obtener mapas de daños potenciales de los edificios expuestos, basados en los gradientes espaciales de movimiento, y se aplica en un área costera de la provincia de Granada (España). A partir de los métodos propuestos para extraer y clasificar las ADAs, y de otros métodos de análisis existentes, se ha desarrollado un paquete de herramientas, los ADAtools, de acceso abierto, fáciles de usar y rápidas, que optimizan la explotación operativa de los mapas de desplazamiento de escala regional. Todas las metodologías se han desarrollado en el marco de varios proyectos europeos (Safety, U-Geohaz, MOMIT y RISKCOAST), y están dirigidos a apoyar las actividades de gestión territorial y análisis de riesgos, con un enfoque específico a los deslizamientos de tierra.Postprint (published version

Advanced exploitation of Sentinel-1 data for supporting landslide risk analysis

Tesi en modalitat de compendi de publicacionsSatellite Synthetic Aperture Radar Interferometry (InSAR) and Persistent Scatterer Interferometry (PSI) are now consolidated tools for ground movement detection and monitoring. Sentinel-1 (S1) is the first satellite providing free data access and ensuring a regular acquisition worldwide, every 6 days, increasing its potential for long-term monitoring applications. Several regional and national ground motion services are already active, providing products based on S1 data. Soon in 2022 the first European Ground Motion Service (EGMS) will be available and freely provide a displacement map over the whole Europe, with annual updates. This implies a strong expansion of availability of PSI-based displacement maps and an easy access for anyone, with an increasing interest among a wider range of users, including public or governmental institutions, academia, industry, and citizens. The analysis and interpretation of this amount of data is difficult and time consuming, mostly for non-expert InSAR users. The objective of this work is developing methodologies to simplify the operational use of PSI displacement maps, generating derived products with a clear message, easy-to-interpret, and fast to read. We propose a method to be applied over regional scale PSI displacement maps, to fast detect the most significant Active Deformation Areas (ADAs). The ADA map is a first product that allows a fast focusing on the active areas, to prioritize further analysis and investigation. Starting from the ADAs, the potential phenomena are attributed to each area through a preliminary interpretation based on auxiliary data, to derive the Geohazard Activity Map. In this work, a methodology to include the ADA information in the Civil Protection Activities is proposed, with the main output called Vulnerable Elements Activity Maps (VEAM). An application of the VEAM is illustrated in the Canary Islands. Furthermore, the ADA map is used in the Valle d'Aosta Region (Northern of Italy) to generate vulnerability and potential loss maps. Finally, a methodology to derive potential damage maps of the exposed buildings, based on the spatial gradients of movement, is proposed, and applied in a coastal area of the Province of Granada (Spain). A pack of software tools has been developed based on the proposed methods to extract ADA and then classify them to generate a Geohazard Activity Map. The set of tools is called ADATools, it is open-access, easy to use and fast, improving the operational exploitation of PSI regional-scale displacement maps. All the methodologies have been developed in the frame of several European projects (Safety, U-Geohaz, MOMIT and RISKCOAST), and are aimed at supporting the multi-scale territorial management and risk analysis activities, with a specific focus on landslides.La interferometría satelital radar (InSAR) y la interferometría de dispersores persistentes (PSI) son herramientas consolidadas para la detección y el monitoreo de movimientos de la superficie de la Tierra. Sentinel-1 (S1) es el primer satélite que proporciona acceso gratuito a los datos y garantiza una adquisición regular en todo el mundo, cada 6 días, aumentando su potencial para aplicaciones de monitoreo a largo plazo. Varios Ground Motion Services regionales y nacionales ya están activos, proporcionando productos basados en datos S1. Pronto, en 2022, el primer servicio europeo (European Ground Motion Service - EGMS) estará disponible y facilitará libremente un mapa de movimientos de toda Europa, con actualizaciones anuales. Esto implica un aumento de la disponibilidad de mapas de movimientos basados en PSI y un fácil acceso para cualquier persona, con un interés creciente entre una amplia gama de usuarios, incluyendo instituciones públicas o gubernamentales, academias, industrias y ciudadanos. El análisis e interpretación de esta cantidad de datos es difícil y consume mucho tiempo, mayormente para usuarios no expertos en la técnica. El objetivo de este trabajo es desarrollar metodologías para simplificar el uso operativo de los mapas de desplazamiento PSI, generando productos derivados con un mensaje claro, fácil de interpretar, y rápido de leer. Se propone un método para detectar rápidamente las Áreas de Deformación Activas (ADAs) más significativas, a partir de mapas de desplazamiento PSI de escala regional. El mapa de las ADAs es un primer producto que permite un enfoque rápido en las áreas activas, útil para priorizar el análisis y las investigaciones adicionales. A partir de las ADAs, se propone una interpretación preliminar basada en datos auxiliares, que atribuye a cada área el fenómeno que está detrás del movimiento, generando el Geohazard Activity Map (GAM). Después, se propone una metodología para incluir la información de las ADAs en las actividades de protección civil, generando los Vulnerable Element Activity Maps (VEAM), a través de su aplicación en las Islas Canarias. Además, el mapa de las ADAs se utiliza en la región de Valle D'Aosta (norte de Italia) para generar mapas de vulnerabilidad y posibles pérdidas económicas. Finalmente, se propone una metodología para obtener mapas de daños potenciales de los edificios expuestos, basados en los gradientes espaciales de movimiento, y se aplica en un área costera de la provincia de Granada (España). A partir de los métodos propuestos para extraer y clasificar las ADAs, y de otros métodos de análisis existentes, se ha desarrollado un paquete de herramientas, los ADAtools, de acceso abierto, fáciles de usar y rápidas, que optimizan la explotación operativa de los mapas de desplazamiento de escala regional. Todas las metodologías se han desarrollado en el marco de varios proyectos europeos (Safety, U-Geohaz, MOMIT y RISKCOAST), y están dirigidos a apoyar las actividades de gestión territorial y análisis de riesgos, con un enfoque específico a los deslizamientos de tierra.Enginyeria del terren

An investigation of ongoing displacements of active faults in the Gobi desert using persistent scatterer interferometric synthetic aperture radar technique to support the permanent disposal of high-level waste in Beishan, China

This research demonstrated the application of PSInSAR method in identifying and characterising the micro-displacements along active faults in Beishan to support the selection of GDF host rock. This research first distinguishes and separates the tectonic induced and non-tectonic induced deformation within three study areas at Suanjingzi, Jiujing and Xinchang. Through the application of coherence change detection, it found the granite outcrop areas characterised by high coherence provide more robust results of tectonic activity. The Quaternary sediments covered areas which are characterised by low coherence usually show higher deformation rates due to the impacts of erosion and deposition. The tectonic induced displacements generally range from -0.4 to 0.4 mma-1 and are dominated by fault bound tectonic movements. As a part of wrench faut zone, Beishan is impacted by a NE-SW trended maximum in situ compressive stress field (σ1). To correlate the visible valleys, gullies, or cracks in Google Earth imagery with the SAR image deformation discontinuities, this study mapped and characterised more than 40 active faults in the three study areas, these include (1) the NE-SW trended sinistral strike-slip faults triggered by extension and (2) the NW-SE/W-E trended reverse faults triggered by maximum compression. The fault activity is characterised by subtle (minor) displacement rate value difference between the two sides of the fault plane. This research successfully improved the understanding of local structural geology and provided moderate guidance for the selection of HLW disposal sites in China. It was indicated that Xinchang has the highest tectonic stability, and this is then followed by Jiujing and Suanjingzi. This kind of displacement rate difference is possible due to the angle difference towards the Sanweishan Fault Zone. To trace and characterise the undiscovered active fault planes, the PSInSAR approach also benefits the prediction of earthquake by improving the positioning of the potential epicentres.Open Acces