Trace Coherence: A New Operator for Polarimetric and Interferometric SAR images

Quadratic forms play an important role in the development of several Polarimetric and Interferometric Synthetic Aperture Radar (Pol-InSAR) methodologies, which are very powerful tools for Earth Observation. This work investigates integrals of Pol-InSAR operators based on quadratic forms, with special interest for the Pol-InSAR coherence. A new operator is introduced, namely Trace Coherence, that provides an approximation for the center of mass of the Coherence Region (CoRe). The latter is the locus of points on the polar plot containing all the possible coherence values. Such center of mass can be calculated as the integral of Pol-InSAR coherences over the scattering mechanisms. The Trace Coherence provides a synthetic information regarding the partial target as one single entity. Therefore, it provides a representation, which is not dependent on the selection of one specific polarization channel. It may find application in change detection (e.g. Coherent Change Detection and differential DEM), classification (e.g. building structure parameters) and modeling (e.g. for the retrieval of forest height). In calculating the integral of the Pol-InSAR coherences, an approximate Trace Coherence expression is derived and shown to improve the calculation speed by several orders of magnitude. The Trace Coherence approximation is investigated using Monte Carlo simulations and validated ESA (DLR) L-band quad-polarimetric data acquired during the AGRISAR 2006 campaign. The result of the analysis using simulated and real data is that the average error in approximating the integral of the Coherence Region is 0.025 in magnitude and 3° in phase (in scenarios with sufficiently high coherence)

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

Pol-InSAR-Island - A benchmark dataset for multi-frequency Pol-InSAR data land cover classification

This paper presents Pol-InSAR-Island, the first publicly available multi-frequency Polarimetric Interferometric Synthetic Aperture Radar (Pol-InSAR) dataset labeled with detailed land cover classes, which serves as a challenging benchmark dataset for land cover classification. In recent years, machine learning has become a powerful tool for remote sensing image analysis. While there are numerous large-scale benchmark datasets for training and evaluating machine learning models for the analysis of optical data, the availability of labeled SAR or, more specifically, Pol-InSAR data is very limited. The lack of labeled data for training, as well as for testing and comparing different approaches, hinders the rapid development of machine learning algorithms for Pol-InSAR image analysis. The Pol-InSAR-Island benchmark dataset presented in this paper aims to fill this gap. The dataset consists of Pol-InSAR data acquired in S- and L-band by DLR\u27s airborne F-SAR system over the East Frisian island Baltrum. The interferometric image pairs are the result of a repeat-pass measurement with a time offset of several minutes. The image data are given as 6 × 6 coherency matrices in ground range on a 1 m × 1m grid. Pixel-accurate class labels, consisting of 12 different land cover classes, are generated in a semi-automatic process based on an existing biotope type map and visual interpretation of SAR and optical images. Fixed training and test subsets are defined to ensure the comparability of different approaches trained and tested prospectively on the Pol-InSAR-Island dataset. In addition to the dataset, results of supervised Wishart and Random Forest classifiers that achieve mean Intersection-over-Union scores between 24% and 67% are provided to serve as a baseline for future work. The dataset is provided via KITopenData: https://doi.org/10.35097/170

Monitoring Water Level Change and Seasonal Vegetation Change in the Coastal Wetlands of Louisiana Using L-Band Time-Series

Coastal wetlands are productive ecosystems driven by highly dynamic hydrological processes such as tides and river discharge, which operate at daily to seasonal timescales, respectively. The scientific community has been calling for landscape-scale measurements of hydrological variables that could help understand the flow of water and transport of sediment across coastal wetlands. While in situ water level gauge data have enabled significant advances, they are limited in coverage and largely unavailable in many parts of the world. In preparation for the NISAR mission, we investigate the use of spaceborne Interferometric Synthetic Aperture Radar (InSAR) observations of phase and coherence at L-band for landscape-scale monitoring of water level change and vegetation cover in coastal wetlands across seasons. We use L-band SAR images acquired by ALOS/PALSAR from 2007 to 2011 to study the impact of seasonal changes in vegetation cover on InSAR sensitivity to water level change in the wetlands of the Atchafalaya basin located in coastal Louisiana, USA. Seasonal variations are observed in the interferometric coherence (γ) time-series over wetlands, with higher coherence during the winter and lower coherence during the summer. We show with InSAR time-series that coherence is inversely correlated with Normalized Difference Vegetation Index (NDVI). Our analysis of polarimetric scattering mechanisms demonstrates that double-bounce is the dominant mechanism in swamps while its weakness in marshes hinders estimation of water level changes. In swamps, water level change maps derived from InSAR are highly correlated (r2 = 0.83) with in situ data from the Coastwide Reference Monitoring System (CRMS). From October to December, we observed that the water level may be below wetland elevation and thus not inundating wetlands significantly. Our analysis shows that water level can only be retrieved when both images used for InSAR are acquired when wetlands are inundated. The L-band derived-maps of water level change show large scale gradients originating from the Gulf Intracoastal Waterway rather than the main delta trunk channel, confirming its significant role as a source of hydrologic connectivity across these coastal wetlands. These results indicate that NISAR, with its InSAR observations every 12 days, will provide the measurements necessary to reveal large scale hydrodynamic processes that occur in swamps across seasons

Monitoring Water Level Change and Seasonal Vegetation Change in the Coastal Wetlands of Louisiana Using L-Band Time-Series

Coastal wetlands are productive ecosystems driven by highly dynamic hydrological processes such as tides and river discharge, which operate at daily to seasonal timescales, respectively. The scientific community has been calling for landscape-scale measurements of hydrological variables that could help understand the flow of water and transport of sediment across coastal wetlands. While in situ water level gauge data have enabled significant advances, they are limited in coverage and largely unavailable in many parts of the world. In preparation for the NISAR mission, we investigate the use of spaceborne Interferometric Synthetic Aperture Radar (InSAR) observations of phase and coherence at L-band for landscape-scale monitoring of water level change and vegetation cover in coastal wetlands across seasons. We use L-band SAR images acquired by ALOS/PALSAR from 2007 to 2011 to study the impact of seasonal changes in vegetation cover on InSAR sensitivity to water level change in the wetlands of the Atchafalaya basin located in coastal Louisiana, USA. Seasonal variations are observed in the interferometric coherence (γ) time-series over wetlands, with higher coherence during the winter and lower coherence during the summer. We show with InSAR time-series that coherence is inversely correlated with Normalized Difference Vegetation Index (NDVI). Our analysis of polarimetric scattering mechanisms demonstrates that double-bounce is the dominant mechanism in swamps while its weakness in marshes hinders estimation of water level changes. In swamps, water level change maps derived from InSAR are highly correlated (r2 = 0.83) with in situ data from the Coastwide Reference Monitoring System (CRMS). From October to December, we observed that the water level may be below wetland elevation and thus not inundating wetlands significantly. Our analysis shows that water level can only be retrieved when both images used for InSAR are acquired when wetlands are inundated. The L-band derived-maps of water level change show large scale gradients originating from the Gulf Intracoastal Waterway rather than the main delta trunk channel, confirming its significant role as a source of hydrologic connectivity across these coastal wetlands. These results indicate that NISAR, with its InSAR observations every 12 days, will provide the measurements necessary to reveal large scale hydrodynamic processes that occur in swamps across seasons

Method for landslides detection with semi-automatic procedures: The case in the zone center-east of Cauca department, Colombia

Landslides are a common natural hazard that causes human casualties, but also infrastructure damage and land-use degradation. Therefore, a quantitative assessment of their presence is required by means of detecting and recognizing the potentially unstable areas. This research aims to develop a method supported on semiautomatic methods to detect potential mass movements at a regional scale. Five techniques were studied: Morphometry, SAR interferometry (InSAR), Persistent Scatterer InSAR (PS-InSAR), SAR polarimetry (PolSAR) and NDVI composites of Landsat 5, Landsat 7, and Landsat 8. The case study was chosen within the mid-eastern area of the Cauca state, which is characterised by its mountainous terrain and the presence of slope instabilities, officially registered in the CGS-SIMMA landslide inventory. This inventory revealed that the type `slide' occurred with 77.4% from the entire registries, `fall' with 16.5%, followed by `creeps' with 3%, flows with 2.6%, and `lateral spread' with 0.43%. As a result, we obtained the morphometric variables: slope, CONVI, TWI, landform, which were highly associated with landslides. The effect of a DEM in the processing flow of the InSAR method was similar for the InSAR coherence variable using the DEMs ASTER, PALSAR RTC, Topo-map, and SRTM. Then, a multiInSAR analysis gave displacement velocities in the LOS direction between -10 and 10 mm/year. With the dual-PolSAR analysis (Sentinel-1), VH and VV C-band polarised radar energy emitted median values of backscatters, for landslides, about of -14.5 dB for VH polarisation and -8.5 dB for VV polarisation. Also, L-band fully polarimetric NASA-UAVSAR data allowed to nd the mechanism of dispersion of CGS landslide inventory: 39% for surface scattering, 46.4% for volume dispersion, and 14.6% for double-bounce scattering. The optical remote sensing provided NDVI composites derived from Landsat series between 2012 and 2016, showing that NDVI values between 0.40 and 0.70 had a high correlation to landslides. In summary, we found the highest categories related to landslides by Weight of Evidence method (WofE) for each spaceborne technique applied. Finally, these results were merged to generate the landslide detection model by using the supervised machine learning method of Random Forest. By taking training and test samples, the precision of the detection model was of about 70% for the rotational and translational types.Los deslizamientos son una amenaza natural que causa pérdidas humanas, daños a la infraestructura y degradación del suelo. Una evaluación cuantitativa de su presencia se requiere mediante la detección y el reconocimiento de potenciales áreas inestables. Esta investigación tuvo como alcance desarrollar un método soportado en métodos semi-automáticos para detectar potenciales movimientos en masa a escala regional. Cinco técnicas fueron estudiadas: Morfometría, Interferometría radar, Interferometría con Persistent Scatterers, Polarimetría radar y composiciones del NDVI con los satélites Landsat 5, Landsat 7 y Landsat 8. El caso de estudio se seleccionó dentro de la región intermedia al este del departamento del Cauca, la cual se caracteriza por terreno montañoso y la presencia de inestabilidades de la pendiente oficialmente registrados en el servicio SIMMA del Servicio Geológico Colombiano. Este inventario reveló que el tipo de movimiento deslizamiento ocurrió con una frecuencia relativa de 77.4%, caidos con el 16.5% de los casos y reptaciones con 3%, flujos con 2.6% y propagación lateral con 0.43%. Como resultado, se obtuvo las variables morfométricas: pendiente, convergencia, índice topográfico de humedad y forma del terreno altamente asociados con los deslizamientos. El efecto de un DEM en el procesamiento del método InSAR fue similar para la variable coherencia usando los DEMs: ASTER, PAlSAR RTC, Topo-map y SRTM. Un análisis Multi-InSAR estimó velocidades de desplazamiento en dirección de vista del radar entre -10 y 10 mm/año. El análisis de polarimetría dual del Sentinel-1 arrojó valores de retrodispersión promedio de -14.5 dB en la banda VH y -8.5dB en la banda VV. Las cuatro polarimetrías del sensor aéreo UAVSAR permitió caracterizar el mecanismo de dispersión del Inventario de Deslizamiento así: 39% en el mecanismo de superficie, 46.4% en el mecanismo de volumen y 14.6% en el mecanismo de doble rebote. La información generada en el rango óptico permitió obtener composiciones de NDVI derivados de la plataforma Landsat entre los años 2012 y 2016, mostrando que el rango entre 0.4 y 0.7 tuvieron una alta asociación con los deslizamientos. En esta investigación se determinaron las categorías de las variables de Teledetección más altamente relacionadas con los movimientos en masa mediante el método de Pesos de Evidencias (WofE). Finalmente, estos resultados se fusionaron para generar el modelo de detección de deslizamientos usando el método supervisado de aprendizaje de máquina Random Forest. Tomando muestras aleatorias para entrenar y validar el modelo en una proporción 70:30, el modelo de detección, especialmente los movimientos de tipo rotacional y traslacional fueron clasificados con una tasa general de éxito del 70%.Ministerio de CienciasConvocatoria 647 de 2014Research line: Geotechnics and Geoenvironmental HazardDoctorad

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

Application of Differential and Polarimetric Synthetic Aperture Radar (SAR) Interferometry for Studying Natural Hazards

In the following work, I address the problem of coherence loss in standard Differential Interferometric SAR (DInSAR) processing, which can result in incomplete or poor quality deformation measurements in some areas. I incorporate polarimetric information with DInSAR in a technique called Polarimetric SAR Interferometry (PolInSAR) in order to acquire more accurate and detailed maps of surface deformation. In Chapter 2, I present a standard DInSAR study of the Ahar double earthquakes (Mw=6.4 and 6.2) which occurred in northwest Iran, August 11, 2012. The DInSAR coseismic deformation map was affected by decorrelation noise. Despite this, I employed an advanced inversion technique, in combination with a Coulomb stress analysis, to find the geometry and the slip distribution on the ruptured fault plane. The analysis shows that the two earthquakes most likely occurred on a single fault, not on conjugate fault planes. This further implies that the minor strike-slip faults play more significant role in accommodating convergence stress accumulation in the northwest part of Iran. Chapter 3 presents results from the application of PolInSAR coherence optimization on quad-pol RADARSAT-2 images. The optimized solution results in the identification of a larger number of reliable measurement points, which otherwise are not recognized by the standard DInSAR technique. I further assess the quality of the optimized interferometric phase, which demonstrates an increased phase quality with respect to those phases recovered by applying standard DInSAR alone. Chapter 4 discusses results from the application of PolInSAR coherence optimization from different geometries to the study of creep on the Hayward fault and landslide motions near Berkeley, CA. The results show that the deformation rates resolved by PolInSAR are in agreement with those of standard DInSAR. I also infer that there is potential motion on a secondary fault, northeast and parallel to the Hayward fault, which may be creeping with a lower velocity

Investigation of polarimetric coherence optimization in persistent scatterer interferometry

"December 2010.""A Thesis presented to the Faculty of the Graduate School University of Missouri--Columbia In Partial Fulfillment of the Requirements for the Degree Master of Science."Thesis supervisor: Dr. Justin J. Legarsky.Interferometric Synthetic Aperture Radar (InSAR) provides a means of imaging small deformations of the Earth's surface (millimeter scale). Furthermore, interferometric point target analysis (IPTA) and time-series analyses can be enhanced using polarimetric and interferometric synthetic aperture radar (PolInSAR) processing techniques to improve detection of point targets or coherent scatterers (CS) that are characterized by a point-like scattering behavior. A mixture of single or multiple scattering mechanisms with a high range of backscattering amplitudes characterize many SAR images over a number of areas. However, the fact that a number of scatterers are characterized by a strongly polarized behavior and are located at different heights even within the same resolution cell make the combination of polarimetric and interferometric information a promising one. For the study site of Socorro, New Mexico, full polarization techniques provided a significant increase (i.e. more than doubled) in CSs compared to a non full polarization. Scattering mechanism and temporal behavior of the CSs were investigated. The CS scattering mechanism behavior was found consistent amongst the CS lists. Full polarization processing was found to increase the CS quantity significantly compared to single polarizations while maintaining high temporal coherence. Thus, this significant increase in full polarization CSs has the potential to increase the overall InSAR processing quality. For sites that challenge (e.g. vegetated areas and other low coherence environments) the single polarization InSAR processing, the full polarization increase in CSs may be essential for processing.Includes bibliographical references (pages 55-56)

Cryosphere Applications

Synthetic aperture radar (SAR) provides large coverage and high resolution, and it has been proven to be sensitive to both surface and near-surface features related to accumulation, ablation, and metamorphism of snow and firn. Exploiting this sensitivity, SAR polarimetry and polarimetric interferometry found application to land ice for instance for the estimation of wave extinction (which relates to sub surface ice volume structure) and for the estimation of snow water equivalent (which relates to snow density and depth). After presenting these applications, the Chapter proceeds by reviewing applications of SAR polarimetry to sea ice for the classification of different ice types, the estimation of thickness, and the characterisation of its surface. Finally, an application to the characterisation of permafrost regions is considered. For each application, the used (model-based) decomposition and polarimetric parameters are critically described, and real data results from relevant airborne campaigns and space borne acquisitions are reported