PSI deformation map retrieval by means of temporal sublook coherence on reduced sets of SAR images

Prior to the application of any persistent scatterer interferometry (PSI) technique for the monitoring of terrain displacement phenomena, an adequate pixel selection must be carried out in order to prevent the inclusion of noisy pixels in the processing. The rationale is to detect the so-called persistent scatterers, which are characterized by preserving their phase quality along the multi-temporal set of synthetic aperture radar (SAR) images available. Two criteria are mainly available for the estimation of pixels' phase quality, i.e., the coherence stability and the amplitude dispersion or permanent scatterers (PS) approach. The coherence stability method allows an accurate estimation of the phase statistics, even when a reduced number of SAR acquisitions is available. Unfortunately, it requires the multi-looking of data during the coherence estimation, leading to a spatial resolution loss in the final results. In contrast, the PS approach works at full-resolution, but it demands a larger number of SAR images to be reliable, typically more than 20. There is hence a clear limitation when a full-resolution PSI processing is to be carried out and the number of acquisitions available is small. In this context, a novel pixel selection method based on exploiting the spectral properties of point-like scatterers, referred to as temporal sublook coherence (TSC), has been recently proposed. This paper seeks to demonstrate the advantages of employing PSI techniques by means of TSC on both orbital and ground-based SAR (GB-SAR) data when the number of images available is small (10 images in the work presented). The displacement maps retrieved through the proposed technique are compared, in terms of pixel density and phase quality, with traditional criteria. Two X-band datasets composed of 10 sliding spotlight TerraSAR-X images and 10 GB-SAR images, respectively, over the landslide of El Forn de Canillo (Andorran Pyrenees), are employed for this study. For both datasets, the TSC technique has showed an excellent performance compared with traditional techniques, achieving up to a four-fold increase in the number of persistent scatters detected, compared with the coherence stability approach, and a similar density compared with the PS approach, but free of outliers.Peer ReviewedPostprint (published version

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

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft

Coherency Matrix Decomposition-Based Polarimetric Persistent Scatterer Interferometry

© 2019 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.The rationale of polarimetric optimization techniques is to enhance the phase quality of the interferograms by combining adequately the different polarization channels available to produce an improved one. Different approaches have been proposed for polarimetric persistent scatterer interferometry (PolPSI). They range from the simple and computationally efficient BEST, where, for each pixel, the polarimetric channel with the best response in terms of phase quality is selected, to those with high-computational burden like the equal scattering mechanism (ESM) and the suboptimum scattering mechanism (SOM). BEST is fast and simple, but it does not fully exploit the potentials of polarimetry. On the other side, ESM explores all the space of solutions and finds the optimal one but with a very high-computational burden. A new PolPSI algorithm, named coherency matrix decomposition-based PolPSI (CMD-PolPSI), is proposed to achieve a compromise between phase optimization and computational cost. Its core idea is utilizing the polarimetric synthetic aperture radar (PolSAR) coherency matrix decomposition to determine the optimal polarization channel for each pixel. Three different PolSAR image sets of both full- (Barcelona) and dual-polarization (Murcia and Mexico City) are used to evaluate the performance of CMD-PolPSI. The results show that CMD-PolPSI presents better optimization results than the BEST method by using either $D_{\mathrm{ A}}$ or temporal mean coherence as phase quality metrics. Compared with the ESM algorithm, CMD-PolPSI is 255 times faster but its performance is not optimal. The influence of the number of available polarization channels and pixel's resolutions on the CMD-PolPSI performance is also discussed.Peer ReviewedPostprint (author's final draft

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

SMF-POLOPT: an adaptive multitemporal pol(DIn)SAR filtering and phase optimization algorithm for PSI applications

© 2019 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Speckle noise and decorrelation can hamper the application and interpretation of PolSAR images. In this paper, a new adaptive multitemporal Pol(DIn)SAR filtering and phase optimization algorithm is proposed to address these limitations. This algorithm first categorizes and adaptively filters permanent scatterer (PS) and distributed scatterer (DS) pixels according to their polarimetric scattering mechanisms [i.e., the scattering-mechanism-based filtering (SMF)]. Then, two different polarimetric DInSAR (POLDInSAR) phase OPTimization methods are applied separately on the filtered PS and DS pixels (i.e., POLOPT). Finally, an inclusive pixel selection approach is used to identify high-quality pixels for ground deformation estimation. Thirty-one full-polarization Radarsat-2 SAR images over Barcelona (Spain) and 31 dual-polarization TerraSAR-X images over Murcia (Spain) have been used to evaluate the performance of the proposed algorithm. The PolSAR filtering results show that the speckle of PolSAR images has been well reduced with the preservation of details by the proposed SMF. The obtained ground deformation monitoring results have shown significant improvements, about ×7.2 (the full-polarization case) and ×3.8 (the dual-polarization case) with respect to the classical full-resolution single-pol amplitude dispersion method, on the valid pixels' densities. The excellent PolSAR filtering and ground deformation monitoring results achieved by the adaptive Pol(DIn)SAR filtering and phase optimization algorithm (i.e., the SMF-POLOPT) have validated the effectiveness of this proposed scheme.Peer ReviewedPostprint (author's final draft

Orbital differential SAR interferometry with polarimetric data

Space-borne Synthetic Aperture Radar (SAR) systems have demonstrated great effectiveness to survey wide areas, independently of the weather conditions and the day/night cycle. One of the most exploited applications is Differential SAR Interferometry (DInSAR), and by extension, advanced multi-temporal techniques addressed as Persistent Scatterers Interferometry (PSI). These techniques provide the means to retrieve the displacement on the observed terrain surface. Taking advantage of the high coverage provided by this technology, natural and human induced phenomena may be monitored and detected in order to avoid and prevent hazardous or even catastrophic situations. A good characterization of such events is often related to the density and quality of the information delivered. In many cases, the success in the interpretation of localized surface motion phenomena relies on having as much measurements as possible. Conventional insitu techniques, when used properly, are able to provide reliable measurements. However, they are limited to their localization, and deploy a high number of instruments is not cost effective. The good characterization provided by PSI measurements is limited by the acquisition parameters and the actual scene under monitoring. SAR images sensitivity over the scene depends strongly on the geometric shape and distribution of targets. Historically, satellite SAR sensors had single polarimetric capabilities, so sensitivity over the scene was restricted to a certain antenna configuration. On the contrary, polarimetric acquisitions emit an receive with different antenna configurations, providing sensitivity to different geometries in the scenario. The launch in these later years of several satellites with polarimetric capabilities has triggered its use in several fields of application. The objective of this thesis is to evaluate the use of Polarimetric SAR (PolSAR) data for its application in the PSI framework. The extra information contained in PolSAR datasets is bound to improve the effectiveness of PSI techniques. Obtaining a higher density and quality of motion measurements will ease the detection and interpretation of terrain displacement phenomena.Els Radars d'Obertura Sintètica (SAR) embarcats en satèl·lits han demostrat ser molt efectius per monitoritzar grans àrees, independentment de les condicions atmosfèriques i del cicle dia/nit. L'aplicació més explotada es l'Interferometria Diferencial SAR (DInSAR), i per extensió, les tècniques avançades de processament multi-temporal anomenades Interferometria de Reflectors Persistents (PSI). Aquestes tècniques són capaces de mesurar el desplaçament en superfície del terreny observat. Aprofitant la gran cobertura que s'aconsegueix amb aquesta tecnologia, fenòmens d'origen tant natural com humà poden ser detectats i controlats per evitar i prevenir situacions de perill o fins i tot catastròfiques. La bona caracterització d'aquestes situacions sovint es relaciona tant amb la densitat com en la qualitat de la informació obtinguda. En molts casos, l’èxit en la interpretació de moviments de terreny localitzats depèn de tenir tants punts de mesura com sigui possible. Les tècniques convencionals de mesura in-situ, quan s'usen de forma adequada, són capaces de donar mesures fiables. No obstant, la seva cobertura està limitada a la localització de l'instrument, i la instal·lació d’un gran nombre d'instruments no és una solució econòmica. La bona caracterització de l'àrea mesurada amb tècniques PSI està limitada pels paràmetres d'adquisició i les característiques de la pròpia escena observada. La sensitivitat de les imatges SAR depèn en gran mesura de la distribució i la forma geomètrica dels blancs mesurats. Històricament, els sensors SAR emetien en polarització simple, i per tant la sensitivitat a la escena estava restringida a aquesta configuració específica de la antena. Per altra banda, durant les adquisicions polarimètriques s'emet i es rep en diferents configuracions de l'antena. Això permet tenir sensitivitat a geometries diferents dins la mateixa escena. El llançament aquest últims anys de diferents satèl·lits amb capacitats polarimètriques ha permès el seu ús en diversos camps d'aplicació. L'objectiu d’aquesta tesi és avaluar l'ús de dades polarimètriques SAR (PolSAR) per la seva aplicació en el marc de les tècniques PSI. La informació addicional que contenen els conjunts d’adquisicions PolSAR ha de millorar l'efectivitat de les tècniques PSI. L'obtenció d'una millor densitat i qualitat de les mesures permet una millor detecció i interpretació dels fenòmens de desplaçament superficial del terreny.Postprint (published version

25 years of satellite InSAR monitoring of ground instability and coastal geohazards in the archaeological site of Capo Colonna, Italy

For centuries the promontory of Capo Colonna in Calabria region, southern Italy, experienced land subsidence and coastline retreat to an extent that the archaeological ruins of the ancient Greek sanctuary are currently under threat of cliff failure, toppling and irreversible loss. Gas extraction in nearby wells is a further anthropogenic element to account for at the regional scale. Exploiting an unprecedented satellite Synthetic Aperture Radar (SAR) time series including ERS-1/2, ENVISAT, TerraSAR-X, COSMO-SkyMed and Sentinel-1A data stacks acquired between 1992 and 2016, this paper presents the first and most complete Interferometric SAR (InSAR) baseline assessment of land subsidence and coastal processes affecting Capo Colonna. We analyse the regional displacement trends, the correlation between vertical displacements with gas extraction volumes, the impact on stability of the archaeological heritage, and the coastal geohazard susceptibility. In the last 25 years, the land has subsided uninterruptedly, with highest annual line-of-sight deformation rates ranging between -15 and -20 mm/year in 2011-2014. The installation of 40 pairs of corner reflectors along the northern coastline and within the archaeological park resulted in an improved imaging capability and higher density of measurement points. This proved to be beneficial for the ground stability assessment of recent archaeological excavations, in an area where field surveying in November 2015 highlighted new events of cliff failure. The conceptual model developed suggests that combining InSAR results, geomorphological assessments and inventorying of wave-storms will contribute to unveil the complexity of coastal geohazards in Capo Colonna. © (2016) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only

Monitoring of remedial works performance on landslide-affected areas through ground- and satellite-based techniques

Differential Interferometric Synthetic Aperture Radar (DInSAR) techniques have repeatedly proved to be an effective tool for built environments monitoring in areas affected by geological hazards. This paper describes how the Coherent Pixel Technique (CPT) approach has been successfully applied to assess the response of an unstable slope to the different phases of remedial works following a landslide event. The CPT technique was performed on 59 COSMO-SkyMed images obtained between May 2011 and August 2016 and centred on the Quercianella settlement (a small hamlet of Livorno municipality, Tuscany, Italy), where the reactivation of a dormant shallow slide had occurred in March 2011 and, hereafter, a geotechnical intervention, designed with the aim of mitigating the risks, has been conducted from August 2013, lasting thirteen months. The time series of CPT results show a deformation pattern with sudden accelerations (up to 21 mm in few months) corresponding to the beginning of the interventions, during which the area has been excavated to install a drainage well, followed by mild decelerations resulting from the stabilization of the area after the conclusion of the works. In particular, the integration of ground-based subsurface monitoring (inclinometers and piezometers) and DInSAR superficial data has provided consistent results for landslide characterization and helped defining the state of activity and the areal distribution of the sliding surface. Moreover, the performance of remedial works in the landslide-affected area has been observed, showing stabilization in the upper part of the hamlet and the ongoing movement in the lower part. The combined monitoring system also led the geotechnical company in charge of remedial works to design further stabilization works in order to preserve buildings and roads in the moving area. Therefore, the integration of remote sensing techniques and in situ instruments represents a timely and cost-efficient solution for intervention works monitoring, opening new perspectives on designing engineering solutions for the stabilization of unstable slopes

Biblioteca del Campus del Baix Llobregat

Postprint (author’s final draft