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

Sentinel-1 data exploitation for terrain deformation monitoring

Persistent Scatterer interferometry (PSI) is a group of advanced differential interferometric Synthetic Aperture Radar (SAR) techniques used to measure and monitor terrain deformation. Sentinel-1 has improved the data acquisition throughout and, compared to previous sensors, increased considerably the Differential Interferometric SAR (DInSAR) and PSI deformation monitoring potential. The effect of the refractive atmosphere on the interferometric phase and phase unwrapping ambiguity are two critical issues of InSAR. The low density of Persistent Scatterer (PS) in non-urban areas, another critical issue, has inspired the development of alternative approaches and refinement of the PS chains. Along with the efforts to develop methods to mitigate the three above-mentioned problems, the work presented in this thesis also deals with the presence of a new signal in multilooked interferograms which cannot be explained by noise, atmospheric or earth surface topography changes. This paper describes a method for atmospheric phase screen estimation using rain station weather data and three different data driven procedures to obtain terrain deformation maps. These approaches aim to exploit Sentinel-1 highly coherent interferograms and their short revisit time. The first method called the splitting makes uses of the power spectrum of the interferograms to split the signals into high and low frequency, and following a mutually exclusive consecutive processing chain for the two sets. This approach has resulted in greater density of PSs with decreased phase unwrapping errors. The second approach, called Direct Integration (DI), aims at providing a very fast and straightforward approach to screen wide areas and easily detect active areas. This approach fully exploits the coherent interferograms from the consecutive images provided by Sentinel-1 resulting in a very high sampling density. However, it lacks robustness and its usability lays on the operator experience. The third method, called PSIG (Persistent Scatterer Interferometry Geomatics) short temporal baseline, provides a constrained application of the PSIG chain, the CTTC approach to the PSI. It uses short temporal baseline interferograms and do not assume any deformation model for point selection. It is also quite a straightforward approach and a perfect complement to the direct integration approach. It improves the performances of the standard PSIG approach, increasing the PS density and providing robust measurements. The effectiveness of the approaches is illustrated through analyses performed on different test sites.La técnica Persistent Scatterer Interferometry (PSI) es un grupo de técnicas avanzadas de radar de apertura sintética interferométrica diferencial (SAR) que se utiliza para medir y monitorear losmovimientos del terreno. Sentinel-1 ha mejorado sensiblemente la adquisición de datos y, en comparación con los sensores SAR anteriores, ha aumentado considerablemente el potencial uso de la interferometría diferencial SAR y del PSI para medir y monitorizar desplazamientos del terreno. El efecto de la atmósfera sobre la fase interferométrica y la naturaleza ambigua de esta son dos cuestiones críticas de InSAR. Además, la baja densidad de Persistent Scatterer (PSs) en áreas no urbanas, es otro tema crítico que ha inspirado el desarrollo de enfoques alternativos y el refinamiento de las cadenas PS existentes. Junto con los esfuerzos por desarrollar métodos para mitigar los tres problemas antes mencionados, el trabajo presentado en esta tesis también aborda la presencia de una nueva señal en interferogramas multilooked que no puede explicarse por cambios de ruido, atmosféricos o topográficos de la superficie terrestre. Esta tesis describe un método para la estimación de la fase atmosférica utilizando datos meteorológicos adquiridos in-situ y tres aproximaciones diferentes basadas en datos Sentinel-1 para obtener mapas de deformación del terreno. Estos enfoques tienen como objetivo explotar los interferogramas altamente coherentes proporcionados por Sentinel-1 gracias a su corto tiempo de revisita. El primer método llamado división hace uso de filtros en el dominico frecuencial de los interferogramas para dividir las señales en alta y baja frecuencia, y siguiendo una cadena de procesamiento consecutiva independiente para cada clase. Este enfoque ha dado como resultado una mejora substancial de PS minimizando los errores debidos al desenrollado de fase. El segundo enfoque, llamado Integración Directa (DI), tiene como objetivo proporcionar un enfoque muy rápido y sencillo para examinar áreas amplias y detectar fácilmente áreas activas. Este enfoque aprovecha al máximo los interferogramas coherentes de las imágenes consecutivas proporcionadas por Sentinel-1, lo que da como resultado una densidad de muestreo muy alta. Sin embargo, carece de robustez y su usabilidad depende de la experiencia del operador. El tercer método, llamado PSIG (Persistent Scatterer Interferometry Geomatics) de línea de base temporal corta, proporciona una aplicación restringida de la cadena PSIG, el enfoque CTTC para el PSI. Utiliza interferogramas de línea base temporales cortos y no asume ningún modelo de deformación para la selección de puntos. Su uso es complementario al enfoque de integración directa proporcionando robustez en las zonas. Mejora el rendimiento del enfoque estándar de PSIG, aumentando la densidad de PS y proporcionando mediciones robustas. La efectividad de los enfoques se ilustra a través de análisis realizados en diferentes sitios de prueba.Postprint (published version

Insar measurements at high latitudes.

This thesis contributes towards understanding of Interferometric Synthetic Aperture Radar (InSAR) measurements at high latitudes. Luossavaara-Kiirunavaara Aktiebolag (LKAB) mining company, the sponsor of this Ph.D. research, intends to use the InSAR techniques for subsidence measurements around the Kiruna underground iron ore mine. The Kiruna underground iron ore mine is located in direct proximity to the city of Kiruna, with the active mining area currently about 1 km west of the city center (67°51'20" N, 20°13'30 E). At present LKAB is exploring the possibility of using InSAR measurements as an operational technique for subsidence measurements. High latitudes InSAR measurements are known to be particularly affected by long periods of ground snow cover that contributes to temporal de-correlation between subsequent radar images. The objectives of this Ph.D. research are (1) to quantify the seasonal effects in InSAR measurements and (2) to identify techniques to improve the high latitude InSAR measurements. In this research study, spatial coherence was used to quantify the seasonal effects in the Differential InSAR (DInSAR) measurements for Kiruna region. A comparison between Static Global Positioning System (Static-GPS) and Corner Reflector InSAR (CRInSAR) measurements were carried out to quantify the seasonal effects in CRInSAR measurements. Spatial ground deformation patterns were used to improve the DInSAR measurements. A theoretical analysis for compact active transponders (CAT) was carried out to improve the North-South InSAR measurements. DInSAR, CRInSAR and Coherent Target Monitoring (CTM) techniques were used to evaluate the applicability of InSAR techniques for high latitude mining induce deformation measurements. The results show seasonal variations in DInSAR, CTM and CRInSAR measurements. Furthermore, DInSAR measurements around the Kiruna iron ore mine can be improved up to sub-centimeter accuracies by using the spatial ground deformation patterns. Also, the compact active transponders identi ed as a possible candidate to improve the accuracy of the North-South InSAR measurements. It is concluded that, all InSAR techniques (DInSAR, CTM and CRInSAR) were affected by the winter snow condition, and only the summer (snow-free) months are suitable for ground deformation measurements. Moreover, the study shows that without the winter images still it is possible to achieve accurate CTM and CRInSAR time series estimations for Kiruna. The East-West vector is the least noisy deformation vector, and both East-West and vertical vectors can be used to determine the LKAB environmental criterion. At present, every year, LKAB is acquiring 45 Radarsat-2 images from three beam modes. However, this study shows that less than 35% of those data are useful for subsidence measurements

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

Surface deformation analysis in Northeast Italy by using PS-InSAR and GNSS data

In the present study, we exploited the potential of satellite-based geodetic data for detecting and measuring surface displacement in Northeast Italy. In this contest, we focused mainly on 1) the estimation of the interseismic deformation during the satellites’ observation period, 2) the detection and analysis of the main deformation patterns, and 3) the correlation of the signals to the active tectonic structures. Despite the low convergence rates (~ 1.5-3 mm/yr), Northeast Italy is an active tectonic area, as testified by the instrumental and historical seismicity. The Adria-Eurasia convergence is mainly accommodated by the thrusts and strike-slip faults of the Southeastern Alps and the External Dinarides, located in the northern and northeastern sectors of the study area. The Venetian-Friulian plain and the Adriatic coasts, affected by active subsidence, dominate the southern region. We used the Stanford Method for Persistent Scatterers (StaMPS) applied to Sentinel-1 SAR images acquired along the ascending and descending orbit tracks between 2015 and 2019. Based on a stack of single-master differential interferograms, we detected coherent and temporally stable pixels based on amplitude and phase noise analysis. After applying spatial-temporal filters and additional post-processing operations to refine the measurements, we used Adria-fixed GNSS velocities derived by permanent stations in the study area to calibrate the InSAR velocities. The outcome consists of Line-OF-Sight (LOS) mean ground velocity maps derived by displacement time series along the radar directions for each satellite track. The combination of the LOS datasets yields vertical and east-west velocity maps, which are mostly in agreement with GNSS data and previous geodetic studies. Based on our measurements, we observe a significant positive velocity gradient of 1 mm/yr across the westernmost sector of the Alpine system, suggesting an aseismic motion of the root of the Bassano-Valdobbiadene thrust. The positive vertical gradients (~1 and up to 2 mm/yr) across the Alpine-Dinaric systems in the central and eastern sectors and the eastward motion that increases northeastward (1-2 mm/yr) may be related to the active Alpine-Dinaric thrusts and strike-slip faults. We also suggest that the detected westward motion of the Friulian plain (around Udine) might be attributed to the presence of tectonic structures characterized by transcurrent-transpressive kinematics. Finally, we detect other signals, such as the significant subsidence (2-4 mm/yr) along the coasts and on the southern Venetian-Friulian plain, confirming the correlation between subsidence and the geological setting of the study area. In conclusion, our study confirms the potential of MT-InSAR and GNSS data for the estimation of the surface deformations in response to active tectonics, even in areas characterized by low deformation rates, such as Northeast Italy.In the present study, we exploited the potential of satellite-based geodetic data for detecting and measuring surface displacement in Northeast Italy. In this contest, we focused mainly on 1) the estimation of the interseismic deformation during the satellites’ observation period, 2) the detection and analysis of the main deformation patterns, and 3) the correlation of the signals to the active tectonic structures. Despite the low convergence rates (~ 1.5-3 mm/yr), Northeast Italy is an active tectonic area, as testified by the instrumental and historical seismicity. The Adria-Eurasia convergence is mainly accommodated by the thrusts and strike-slip faults of the Southeastern Alps and the External Dinarides, located in the northern and northeastern sectors of the study area. The Venetian-Friulian plain and the Adriatic coasts, affected by active subsidence, dominate the southern region. We used the Stanford Method for Persistent Scatterers (StaMPS) applied to Sentinel-1 SAR images acquired along the ascending and descending orbit tracks between 2015 and 2019. Based on a stack of single-master differential interferograms, we detected coherent and temporally stable pixels based on amplitude and phase noise analysis. After applying spatial-temporal filters and additional post-processing operations to refine the measurements, we used Adria-fixed GNSS velocities derived by permanent stations in the study area to calibrate the InSAR velocities. The outcome consists of Line-OF-Sight (LOS) mean ground velocity maps derived by displacement time series along the radar directions for each satellite track. The combination of the LOS datasets yields vertical and east-west velocity maps, which are mostly in agreement with GNSS data and previous geodetic studies. Based on our measurements, we observe a significant positive velocity gradient of 1 mm/yr across the westernmost sector of the Alpine system, suggesting an aseismic motion of the root of the Bassano-Valdobbiadene thrust. The positive vertical gradients (~1 and up to 2 mm/yr) across the Alpine-Dinaric systems in the central and eastern sectors and the eastward motion that increases northeastward (1-2 mm/yr) may be related to the active Alpine-Dinaric thrusts and strike-slip faults. We also suggest that the detected westward motion of the Friulian plain (around Udine) might be attributed to the presence of tectonic structures characterized by transcurrent-transpressive kinematics. Finally, we detect other signals, such as the significant subsidence (2-4 mm/yr) along the coasts and on the southern Venetian-Friulian plain, confirming the correlation between subsidence and the geological setting of the study area. In conclusion, our study confirms the potential of MT-InSAR and GNSS data for the estimation of the surface deformations in response to active tectonics, even in areas characterized by low deformation rates, such as Northeast Italy

VALIDATION OF FULL-RESOLUTION DINSAR-DERIVED VERTICAL DISPLACEMENT IN CULTURAL HERITAGE MONITORING: INTEGRATION WITH GEODETIC LEVELLING MEASUREMENTS

Towards revealing the potential of satellite Synthetic Aperture Radar (SAR) Interferometry (InSAR) for efficient detection and monitoring of Cultural Heritage (CH) encouraging resilient built CH, this study is devoted to the validation of InSAR-derived vertical displacements with a full-resolution perspective taking advantage of high-precision geodetic levelling measurements. Considering the Cathedral of Como, northern Italy, as the case study, two different Persistent Scatterer Interferometry (PSI) techniques have been applied to Cosmo-SkyMed high-resolution SAR images acquired in both ascending and descending orbit tacks within the time interval of 2010–2012. Besides using the simplified approach for obtaining the vertical displacement velocity from Line of Sight (LOS) velocity, a weighted, localized, multi-track Vertical Displacement Extraction (VDE) approach is proposed and evaluated, which uses the technical outcome of Differential InSAR (DInSAR) and spatial information. The results, using a proper PSI technique, showed that the accuracy level of extracted vertical displacement velocities in a full-resolution application is ca. 0.6 [mm/year] with a dense concentration of InSAR-Levelling absolute errors lower than 0.3 [mm/year] which are reliable and reasonable levels based on the employed validation framework in this study. Also, the weighted localized VDE can significantly decrease the InSAR-Levelling errors, adding to the reliability of the InSAR application for CH monitoring and condition assessment in practice

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.The different research areas included in this paper has been supported by the projects: CGL2005-05500-C02, CGL2008-06426-C01-01/BTE, AYA2 010-17448, IPT-2011-1234-310000, TEC-2008-06764, ACOMP/2010/082, AGL2009-08931/AGR, 2012GA-LC-036, 2003-03-4.3-I-014, CGL2006-05415, BEST-2011/225, CGL2010-16775, TEC2011-28201, 2012GA-LC-021 and the Banting Postdoctoral Fellowship to PJG