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

The worsening impacts of land reclamation assessed with Sentinel-1: The Rize (Turkey) test case

Massive amounts of land are being reclaimed to build airports, new cities, ports, and highways. Hundreds of kilometers are added each year, as coastlines are extended further out to the sea. In this paper, this urbanization approach is monitored by Persistent Scatterer Interferometry (PSI) technique with Sentinel-1 SAR data. The study aims to explore this technology in order to support local authorities to detect and evaluate subtle terrain displacements. For this purpose, a large 3-years Sentinel-1 stack composed by 92 images acquired between 07/01/2015 to 27/01/2018 is employed and stacking techniques are chosen to assess ground motion. The test site of this study, Rize, Turkey, has been declared at high risk of collapse and radical solutions such as the relocation of the entire city in another area are been taken into consideration. A media fact-checking approach, i.e. evaluating national and international press releases on the test site, is considered for the paper and this work presents many findings in different areas of the city. For instance, alerts are confirmed by inspecting several buildings reported by the press. Critical infrastructures are monitored as well. Portions of the harbor show high displacement rates, up to 1 cm/year, proving reported warnings. Rural villages belonging to the same municipality are also investigated and a mountainous village affected by landslide is considered in the study. Sentinel-1 is demonstrated to be a suitable system to detect and monitor small changes or buildings and infrastructures for these scenarios. These changes may be highly indicative of imminent damage which can lead to the loss of the structural integrity and subsequent failure of the structure in the long-term. In Rize, only a few known motion-critical structures are monitored daily with in-situ technologies. SAR interferometry can assist to save expensive inspection and monitoring services, especially in highly critical cases such as the one studied in this paper

Analysis of permanent scatterers properties by land cover type in Ida-Virumaa

http://tartu.ester.ee/record=b2693107~S1*es

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

Simulating SAR geometric distortions and predicting Persistent Scatterer densities for ERS-1/2 and ENVISAT C-band SAR and InSAR applications: nationwide feasibility assessment to monitor the landmass of Great Britain with SAR imagery

We assess the feasibility of monitoring the landmass of Great Britain with satellite Synthetic Aperture Radar (SAR) imagery, by analysing ERS-1/2 SAR and ENVISAT IS2 Advanced SAR (ASAR) archive data availability, geometric distortions and land cover control on the success of (non-)interferometric analyses. Our assessment both addresses the scientific and operational question of whether a nationwide SAR-based monitoring of ground motion would succeed in Great Britain, and helps to understand controlling factors and possible solutions to overcome the limitations of undertaking SAR-based imaging of the landmass. This is the first time such a nationwide assessment is performed in preparation for acquisition and processing of SAR data in the United Kingdom, and any other country in the world. Analysis of the ERS-1/2 and ENVISAT archives reveals potential for multi-interferogram SAR Interferometry (InSAR) for the entirety of Britain using ERS-1/2 in descending mode, with 100% standard image frames showing at least 20 archive scenes available. ERS-1/2 ascending and both ENVISAT modes show potential for non-interferometric and single-pair InSAR for the vast majority of Britain, and multi-interferogram only for 13% to 38% of the available standard frames. Based on NEXTMap® Britain Digital Terrain Model (DTM) we simulate SAR layover, foreshortening and shadow to the ERS-1/2 and ENVISAT Lines-Of-Sight (LOS), and quantify changes of SAR distortions with variations in mode, LOS incidence angles and ground track angles, local terrain orientation, and the effect of scale due to the input DTM resolution. The simulation is extended to the ~ 230,000 km2 landmass, and shows limited control of local topography on the radar terrain visibility. According to the 50 m to 5 m DTM-based simulations, ~ 1.0–1.4% of Great Britain could potentially be affected by shadow and layover in each mode. Only ~ 0.02–0.04% overlapping between ascending and descending mode distortions is found, this indicating the negligible proportion of the landmass that cannot be monitored using either imaging mode. We calibrate the CORINE Land Cover 2006 (CLC2006) using Persistent Scatterer (PS) datasets available for London, Stoke-On-Trent, Newcastle and Bristol, to quantify land cover control on the PS distribution and characterise the CLC2006 classes in terms of the potential PS density they could provide. Despite predominance of rural land cover types, we predict potential for over 12.8 M monitoring targets for each acquisition mode using a set of image frames covering the entire landmass. We validate our assessment by processing with the Interferometric Point Target Analysis (IPTA) 55 ERS-1/2 SAR scenes depicting South Wales between 1992 and 1999. Although absolute differences between predicted and observed target density are revealed, relative densities and rankings among the various CLC2006 classes are found constant across the calibration and validation datasets. Rescaled predictions for Britain show potential for a total of 2.5 M monitoring targets across the landmass. We examine the use of the topographic and land cover feasibility maps for landslide studies in relation to the British Geological Survey's National Landslide Database and DiGMapGB mass movement layer. Building upon recent literature, we finally discuss future perspectives relating to the replication of our feasibility assessment to account for higher resolution SAR imagery, new Earth explorers (e.g., Sentinel-1) and improved processing techniques, showing potential to generate invaluable sources of information on land motions and geohazards in Great Britai

Imaging multi-age construction settlement behaviour by advanced SAR interferometry

This paper focuses on the application of Advanced Satellite Synthetic Aperture Radar Interferometry (A-DInSAR) to subsidence-related issues, with particular reference to ground settlements due to external loads. Beyond the stratigraphic setting and the geotechnical properties of the subsoil, other relevant boundary conditions strongly influence the reliability of remotely sensed data for quantitative analyses and risk mitigation purposes. Because most of the Persistent Scatterer Interferometry (PSI) measurement points (Persistent Scatterers, PSs) lie on structures and infrastructures, the foundation type and the age of a construction are key factors for a proper interpretation of the time series of ground displacements. To exemplify a methodological approach to evaluate these issues, this paper refers to an analysis carried out in the coastal/deltaic plain west of Rome (Rome and Fiumicino municipalities) affected by subsidence and related damages to structures. This region is characterized by a complex geological setting (alternation of recent deposits with low and high compressibilities) and has been subjected to different urbanisation phases starting in the late 1800s, with a strong acceleration in the last few decades. The results of A-DInSAR analyses conducted from 1992 to 2015 have been interpreted in light of high-resolution geological/geotechnical models, the age of the construction, and the types of foundations of the buildings on which the PSs are located. Collection, interpretation, and processing of geo-thematic data were fundamental to obtain high-resolution models; change detection analyses of the land cover allowed us to classify structures/infrastructures in terms of the construction period. Additional information was collected to define the types of foundations, i.e., shallow versus deep foundations. As a result, we found that only by filtering and partitioning the A-DInSAR datasets on the basis of the above-mentioned boundary conditions can the related time series be considered a proxy of the consolidation process governing the subsidence related to external loads as confirmed by a comparison with results from a physically based back analysis based on Terzaghi's theory. Therefore, if properly managed, the A-DInSAR data represents a powerful tool for capturing the evolutionary stage of the process for a single building and has potential for forecasting the behaviour of the terrain-foundation-structure combination

Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work