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

SAR Tomography via Nonlinear Blind Scatterer Separation

Layover separation has been fundamental to many synthetic aperture radar applications, such as building reconstruction and biomass estimation. Retrieving the scattering profile along the mixed dimension (elevation) is typically solved by inversion of the SAR imaging model, a process known as SAR tomography. This paper proposes a nonlinear blind scatterer separation method to retrieve the phase centers of the layovered scatterers, avoiding the computationally expensive tomographic inversion. We demonstrate that conventional linear separation methods, e.g., principle component analysis (PCA), can only partially separate the scatterers under good conditions. These methods produce systematic phase bias in the retrieved scatterers due to the nonorthogonality of the scatterers' steering vectors, especially when the intensities of the sources are similar or the number of images is low. The proposed method artificially increases the dimensionality of the data using kernel PCA, hence mitigating the aforementioned limitations. In the processing, the proposed method sequentially deflates the covariance matrix using the estimate of the brightest scatterer from kernel PCA. Simulations demonstrate the superior performance of the proposed method over conventional PCA-based methods in various respects. Experiments using TerraSAR-X data show an improvement in height reconstruction accuracy by a factor of one to three, depending on the used number of looks.Comment: This work has been accepted by IEEE TGRS for publicatio

Multiresolution Detection of Persistent Scatterers: A Performance Comparison Between Multilook GLRT and CAESAR

Persistent scatterers (PS) interferometry tools are extensively used for the monitoring of slow, long-term ground deformation. High spatial resolution is typically required in urban areas to cope with the variability of the signal, whereas in rural regions, multilook shall be implemented to improve the coverage of monitored areas. Along this line, SqueeSAR and later Component extrAction and sElection SAR (CAESAR) were introduced for the monitoring of both persistent and (decorrelating) distributed scatterers (DS). Multilook generalized likelihood ratio test (MGLRT) is a detector derived in the context of tomographic SAR processing that has been investigated for a fixed multilook degree. In this work, we address MGLRT and CAESAR in the multiresolution context characterized by a spatially variable multilook degree. We compare the two schemes for the multiresolution selection of PS and DS, highlighting the pros and cons of each scheme, particularly the peculiarities of CAESAR that have important implications at the implementation stage. A performance analysis of both detectors in case of model mismatch is also addressed. Experiments carried out with data acquired by the COSMO-SkyMed constellation support both the theoretical argumentation and the results achieved by resorting to Monte Carlo simulations

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Very High Resolution Tomographic SAR Inversion for Urban Infrastructure Monitoring — A Sparse and Nonlinear Tour

The topic of this thesis is very high resolution (VHR) tomographic SAR inversion for urban infrastructure monitoring. To this end, SAR tomography and differential SAR tomography are demonstrated using TerraSAR-X spotlight data for providing 3-D and 4-D (spatial-temporal) maps of an entire high rise city area including layover separation and estimation of deformation of the buildings. A compressive sensing based estimator (SL1MMER) tailored to VHR SAR data is developed for tomographic SAR inversion by exploiting the sparsity of the signal. A systematic performance assessment of the algorithm is performed regarding elevation estimation accuracy, super-resolution and robustness. A generalized time warp method is proposed which enables differential SAR tomography to estimate multi-component nonlinear motion. All developed methods are validated with both simulated and extensive processing of large volumes of real data from TerraSAR-X

Polarimetric Synthetic Aperture Radar

This open access book focuses on the practical application of electromagnetic polarimetry principles in Earth remote sensing with an educational purpose. In the last decade, the operations from fully polarimetric synthetic aperture radar such as the Japanese ALOS/PalSAR, the Canadian Radarsat-2 and the German TerraSAR-X and their easy data access for scientific use have developed further the research and data applications at L,C and X band. As a consequence, the wider distribution of polarimetric data sets across the remote sensing community boosted activity and development in polarimetric SAR applications, also in view of future missions. Numerous experiments with real data from spaceborne platforms are shown, with the aim of giving an up-to-date and complete treatment of the unique benefits of fully polarimetric synthetic aperture radar data in five different domains: forest, agriculture, cryosphere, urban and oceans

Robust and Flexible Persistent Scatterer Interferometry for Long-Term and Large-Scale Displacement Monitoring

Die Persistent Scatterer Interferometrie (PSI) ist eine Methode zur Überwachung von Verschiebungen der Erdoberfläche aus dem Weltraum. Sie basiert auf der Identifizierung und Analyse von stabilen Punktstreuern (sog. Persistent Scatterer, PS) durch die Anwendung von Ansätzen der Zeitreihenanalyse auf Stapel von SAR-Interferogrammen. PS Punkte dominieren die Rückstreuung der Auflösungszellen, in denen sie sich befinden, und werden durch geringfügige Dekorrelation charakterisiert. Verschiebungen solcher PS Punkte können mit einer potenziellen Submillimetergenauigkeit überwacht werden, wenn Störquellen effektiv minimiert werden. Im Laufe der Zeit hat sich die PSI in bestimmten Anwendungen zu einer operationellen Technologie entwickelt. Es gibt jedoch immer noch herausfordernde Anwendungen für die Methode. Physische Veränderungen der Landoberfläche und Änderungen in der Aufnahmegeometrie können dazu führen, dass PS Punkte im Laufe der Zeit erscheinen oder verschwinden. Die Anzahl der kontinuierlich kohärenten PS Punkte nimmt mit zunehmender Länge der Zeitreihen ab, während die Anzahl der TPS Punkte zunimmt, die nur während eines oder mehrerer getrennter Segmente der analysierten Zeitreihe kohärent sind. Daher ist es wünschenswert, die Analyse solcher TPS Punkte in die PSI zu integrieren, um ein flexibles PSI-System zu entwickeln, das in der Lage ist mit dynamischen Veränderungen der Landoberfläche umzugehen und somit ein kontinuierliches Verschiebungsmonitoring ermöglicht. Eine weitere Herausforderung der PSI besteht darin, großflächiges Monitoring in Regionen mit komplexen atmosphärischen Bedingungen durchzuführen. Letztere führen zu hoher Unsicherheit in den Verschiebungszeitreihen bei großen Abständen zur räumlichen Referenz. Diese Arbeit befasst sich mit Modifikationen und Erweiterungen, die auf der Grund lage eines bestehenden PSI-Algorithmus realisiert wurden, um einen robusten und flexiblen PSI-Ansatz zu entwickeln, der mit den oben genannten Herausforderungen umgehen kann. Als erster Hauptbeitrag wird eine Methode präsentiert, die TPS Punkte vollständig in die PSI integriert. In Evaluierungsstudien mit echten SAR Daten wird gezeigt, dass die Integration von TPS Punkten tatsächlich die Bewältigung dynamischer Veränderungen der Landoberfläche ermöglicht und mit zunehmender Zeitreihenlänge zunehmende Relevanz für PSI-basierte Beobachtungsnetzwerke hat. Der zweite Hauptbeitrag ist die Vorstellung einer Methode zur kovarianzbasierten Referenzintegration in großflächige PSI-Anwendungen zur Schätzung von räumlich korreliertem Rauschen. Die Methode basiert auf der Abtastung des Rauschens an Referenzpixeln mit bekannten Verschiebungszeitreihen und anschließender Interpolation auf die restlichen PS Pixel unter Berücksichtigung der räumlichen Statistik des Rauschens. Es wird in einer Simulationsstudie sowie einer Studie mit realen Daten gezeigt, dass die Methode überlegene Leistung im Vergleich zu alternativen Methoden zur Reduktion von räumlich korreliertem Rauschen in Interferogrammen mittels Referenzintegration zeigt. Die entwickelte PSI-Methode wird schließlich zur Untersuchung von Landsenkung im Vietnamesischen Teil des Mekong Deltas eingesetzt, das seit einigen Jahrzehnten von Landsenkung und verschiedenen anderen Umweltproblemen betroffen ist. Die geschätzten Landsenkungsraten zeigen eine hohe Variabilität auf kurzen sowie großen räumlichen Skalen. Die höchsten Senkungsraten von bis zu 6 cm pro Jahr treten hauptsächlich in städtischen Gebieten auf. Es kann gezeigt werden, dass der größte Teil der Landsenkung ihren Ursprung im oberflächennahen Untergrund hat. Die präsentierte Methode zur Reduzierung von räumlich korreliertem Rauschen verbessert die Ergebnisse signifikant, wenn eine angemessene räumliche Verteilung von Referenzgebieten verfügbar ist. In diesem Fall wird das Rauschen effektiv reduziert und unabhängige Ergebnisse von zwei Interferogrammstapeln, die aus unterschiedlichen Orbits aufgenommen wurden, zeigen große Übereinstimmung. Die Integration von TPS Punkten führt für die analysierte Zeitreihe von sechs Jahren zu einer deutlich größeren Anzahl an identifizierten TPS als PS Punkten im gesamten Untersuchungsgebiet und verbessert damit das Beobachtungsnetzwerk erheblich. Ein spezieller Anwendungsfall der TPS Integration wird vorgestellt, der auf der Clusterung von TPS Punkten basiert, die innerhalb der analysierten Zeitreihe erschienen, um neue Konstruktionen systematisch zu identifizieren und ihre anfängliche Bewegungszeitreihen zu analysieren

Single-Look SAR Tomography of Urban Areas

Synthetic aperture radar (SAR) tomography (TomoSAR) is a multibaseline interferometric technique that estimates the power spectrum pattern (PSP) along the perpendicular to the line-ofsight (PLOS) direction. TomoSAR achieves the separation of individual scatterers in layover areas, allowing for the 3D representation of urban zones. These scenes are typically characterized by buildings of different heights, with layover between the facades of the higher structures, the rooftop of the smaller edifices and the ground surface. Multilooking, as required by most spectral estimation techniques, reduces the azimuth-range spatial resolution, since it is accomplished through the averaging of adjacent values, e.g., via Boxcar filtering. Consequently, with the aim of avoiding the spatial mixture of sources due to multilooking, this article proposes a novel methodology to perform single-look TomoSAR over urban areas. First, a robust version of Capon is applied to focus the TomoSAR data, being robust against the rank-deficiencies of the data covariance matrices. Afterward, the recovered PSP is refined using statistical regularization, attaining resolution enhancement, suppression of artifacts and reduction of the ambiguity levels. The capabilities of the proposed methodology are demonstrated by means of strip-map airborne data of the Jet Propulsion Laboratory (JPL) and the National Aeronautics and Space Administration (NASA), acquired by the uninhabited aerial vehicle SAR (UAVSAR) system over the urban area of Munich, Germany in 2015. Making use of multipolarization data [horizontal/horizontal (HH), horizontal/vertical (HV) and vertical/vertical (VV)], a comparative analysis against popular focusing techniques for urban monitoring (i.e., matched filtering, Capon and compressive sensing (CS)) is addressed

High-resolution deformation measurement using "Persistent Scatterer Interferometry"

Persistent Scatterer Interferometry (PSI) is a group of advanced differential interferometric SAR techniques that are used to measure and monitor terrain deformation. Different PSI techniques have been proposed in the last two decades. In this thesis, the two PSI chains implemented and used at the Geomatics division of CTTC are described: the local area PSI and the PSIG chains. The first part of the thesis is devoted to the local area PSI chain, used to analyse the deformations over small areas. The chain includes a linear deformation model to directly deal with interferometric wrapped phases. Moreover, it does not directly involve the estimation of the APS, thus simplifying the procedure and its computational cost. The chain has been tested using different types of SAR data. The availability of high resolution X-band SAR data has led to an improvement of the PSI results with respect to C-band data. The higher image resolution and phase quality implies an increase of the PS density, an improvement in the estimation precision of the residual topographic error and a higher sensibility to very small deformations, including the displacements caused by thermal dilation. An extension of the classical PSI linear deformation model has been proposed, to account for the thermal dilation effects. This allows obtaining a new PSI outcome, the thermal dilation parameter, which opens new interesting applications since it provides information on the physical properties of single objects, i.e. the coefficient of thermal expansion, and the static structures of the same objects. The second part of the thesis describes the PSIG chain, whose aim was to extend the interferometric processing to wider areas. The ability to cover wide areas is essential to obtain a unique and consistent deformation monitoring for the available SAR image full scenes, i.e. typically 30 by 50 km for TerraSAR-X, 40 by 40 km for CosmoSkyMed and 100 by 100 km for ASAR ENVISAT and ERS. This is particularly important for the forthcoming C-band Sentinel SAR data that will cover 250 by 250 km with a single image scene. The key steps of the PSIG procedure include a new selection of candidate PSs based on a phase similitude criteria and a 2+1D phase unwrapping algorithm. The procedure offers different tools to control the quality of the processing steps. It has been successfully tested over urban, rural and vegetated areas using X-band PSI data. The performance of the PSIG chain is illustrated and discussed in detail, analysing the procedure step by step.Persistent Scatterer Interferometry (PSI) és un grup de tècniques avançades d'interferometria diferencial SAR que s'utilitzen per mesurar i monitoritzar deformacions del terreny. Durant les últimes dues dècades s’han proposat diverses tècniques PSI. En aquesta tesi es descriuen les dues cadenes PSI implementades i utilitzades en la divisió de Geomàtica del CTTC: la cadena PSI d’àrea local i la cadena PSIG. La primera part de la tesi està dedicada a la cadena PSI d’àrea local, que s'utilitza per analitzar deformacions en zones d’extensió limitada. La cadena inclou un model de deformació lineal per tractar directament amb les fases interferomètriques wrapped. En canvi, no estima directament la component atmosfèrica, cosa que simplifica el procediment i el seu cost computacional. La cadena s’ha provat sobre diferents tipus de dades SAR. La disponibilitat de dades SAR d’alta resolució en banda X ha donat lloc a una millora dels resultats del PSI respecte a les dades en banda C. La resolució més gran de la imatge i la qualitat de la fase impliquen un augment de la densitat de PS, una millora en la precisió de l'estimació de l'error topogràfic residual i una sensibilitat més alta a deformacions subtils, incloent-hi els desplaçaments causats per la dilatació tèrmica. Per tenir en compte els efectes de la dilatació tèrmica, s'ha proposat una extensió del model PSI clàssic que ens permet obtenir un nou producte PSI: el paràmetre de dilatació tèrmica. Aquest paràmetre obre noves aplicacions interessants: proporciona informació relacionada amb les propietats físiques dels objectes mesurats –com el coeficient d'expansió tèrmica– i amb la seva pròpia estructura estàtica. La segona part de la tesi descriu la cadena PSIG, l'objectiu de la qual és estendre el processament interferomètric a àrees més extenses. La capacitat de cobrir àrees grans és fonamental per obtenir un únic mapa global de deformacions que sigui consistent i cobreixi l’extensió sencera de les imatges SAR disponibles, de 30 km per 50 km per TerraSAR-X, de 40 km per 40 km per CosmoSkyMed i de 100 km per 100 km per ASAR-ENVISAT i ERS. Això és particularment important tenint en compte la propera disponibilitat de les dades del satèl•lit Sentinel, que opera en banda C i cobrirà 250 km per 250 km amb una sola imatge. Els passos clau del procediment PSIG són una nova selecció de PS candidats en base a un criteri de similitud de fase i un algoritme de 2+1D phase unwrapping. El procediment ofereix diferents eines per controlar la qualitat dels diferents passos del processament. La cadena PSIG s’ha utilitzat amb èxit en àrees urbanes, rurals i amb vegetació utilitzant dades PSI en banda X. El funcionament de la cadena PSIG s'il•lustra i es descriu en detall, analitzant el procediment pas a pas