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

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

TomoSAR Mapping of 3D Forest Structure: Contributions of L-Band Configurations

Synthetic Aperture Radar (SAR) measurements are unique for mapping forest 3D structure and its changes in time. Tomographic SAR (TomoSAR) configurations exploit this potential by reconstructing the 3D radar reflectivity. The frequency of the SAR measurements is one of the main parameters determining the information content of the reconstructed reflectivity in terms of penetration and sensitivity to the individual vegetation elements. This paper attempts to review and characterize the structural information content of L-band TomoSAR reflectivity reconstructions, and their potential to forest structure mapping. First, the challenges in the accurate TomoSAR reflectivity reconstruction of volume scatterers (which are expected to dominate at L-band) and to extract physical structure information from the reconstructed reflectivity is addressed. Then, the L-band penetration capability is directly evaluated by means of the estimation performance of the sub-canopy ground topography. The information content of the reconstructed reflectivity is then evaluated in terms of complementary structure indices. Finally, the dependency of the TomoSAR reconstruction and of its structural information to both the TomoSAR acquisition geometry and the temporal change of the reflectivity that may occur in the time between the TomoSAR measurements in repeat-pass or bistatic configurations is evaluated. The analysis is supported by experimental results obtained by processing airborne acquisitions performed over temperate forest sites close to the city of Traunstein in the south of Germany

Estimation of the Degree of Polarization in Polarimetric SAR Imagery : Principles and Applications

Les radars à synthèse d’ouverture (RSO) polarimétriques sont devenus incontournables dans le domaine de la télédétection, grâce à leur zone de couverture étendue, ainsi que leur capacité à acquérir des données dans n’importe quelles conditions atmosphériques de jour comme de nuit. Au cours des trois dernières décennies, plusieurs RSO polarimétriques ont été utilisés portant une variété de modes d’imagerie, tels que la polarisation unique, la polarisation double et également des modes dits pleinement polarimétriques. Grâce aux recherches récentes, d’autres modes alternatifs, tels que la polarisation hybride et compacte, ont été proposés pour les futures missions RSOs. Toutefois, un débat anime la communauté de la télédétection quant à l’utilité des modes alternatifs et quant au compromis entre la polarimétrie double et la polarimétrie totale. Cette thèse contribue à ce débat en analysant et comparant ces différents modes d’imagerie RSO dans une variété d’applications, avec un accent particulier sur la surveillance maritime (la détection des navires et de marées noires). Pour nos comparaisons, nous considérons un paramètre fondamental, appelé le degré de polarisation (DoP). Ce paramètre scalaire a été reconnu comme l’un des paramètres les plus pertinents pour caractériser les ondes électromagnétiques partiellement polarisées. A l’aide d’une analyse statistique détaillée sur les images polarimétriques RSO, nous proposons des estimateurs efficaces du DoP pour les systèmes d’imagerie cohérente et incohérente. Ainsi, nous étendons la notion de DoP aux différents modes d’imagerie polarimétrique hybride et compacte. Cette étude comparative réalisée dans différents contextes d’application dégage des propriétés permettant de guider le choix parmi les différents modes polarimétriques. Les expériences sont effectuées sur les données polarimétriques provenant du satellite Canadian RADARSAT-2 et le RSO aéroporté Américain AirSAR, couvrant divers types de terrains tels que l’urbain, la végétation et l’océan. Par ailleurs nous réalisons une étude détaillée sur les potentiels du DoP pour la détection et la reconnaissance des marées noires basée sur les acquisitions récentes d’UAVSAR, couvrant la catastrophe de Deepwater Horizon dans le golfe du Mexique. ABSTRACT : Polarimetric Synthetic Aperture Radar (SAR) systems have become highly fruitful thanks to their wide area coverage and day and night all-weather capabilities. Several polarimetric SARs have been flown over the last few decades with a variety of polarimetric SAR imaging modes; traditional ones are linear singleand dual-pol modes. More sophisticated ones are full-pol modes. Other alternative modes, such as hybrid and compact dual-pol, have also been recently proposed for future SAR missions. The discussion is vivid across the remote sensing society about both the utility of such alternative modes, and also the trade-off between dual and full polarimetry. This thesis contributes to that discussion by analyzing and comparing different polarimetric SAR modes in a variety of geoscience applications, with a particular focus on maritime monitoring and surveillance. For our comparisons, we make use of a fundamental, physically related discriminator called the Degree of Polarization (DoP). This scalar parameter has been recognized as one of the most important parameters characterizing a partially polarized electromagnetic wave. Based on a detailed statistical analysis of polarimetric SAR images, we propose efficient estimators of the DoP for both coherent and in-coherent SAR systems. We extend the DoP concept to different hybrid and compact SAR modes and compare the achieved performance with different full-pol methods. We perform a detailed study of vessel detection and oil-spill recognition, based on linear and hybrid/compact dual-pol DoP, using recent data from the Deepwater Horizon oil-spill, acquired by the National Aeronautics and Space Administration (NASA)/Jet Propulsion Laboratory (JPL) Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR). Extensive experiments are also performed over various terrain types, such as urban, vegetation, and ocean, using the data acquired by the Canadian RADARSAT-2 and the NASA/JPL Airborne SAR (AirSAR) system

Accurate Despeckling and Estimation of Polarimetric Features by Means of a Spatial Decorrelation of the Noise in Complex PolSAR Data

In this work, we extended a procedure for the spatial decorrelation of fully-developed speckle, originally developed for single-polarization SAR data, to fully-polarimetric SAR data. The spatial correlation of the noise depends on the tapering window in the Fourier domain used by the SAR processor to avoid defocusing of targets caused by Gibbs effects. Since each polarimetric channel is focused independently of the others, the noise-whitening procedure can be performed applying the decorrelation stage to each channel separately. Equivalently, the noise-whitening stage is applied to each element of the scattering matrix before any multilooking operation, either coherent or not, is performed. In order to evaluate the impact of a spatial decorrelation of the noise on the performance of polarimetric despeckling filters, we make use of simulated PolSAR data, having user-defined polarimetric features. We optionally introduce a spatial correlation of the noise in the simulated complex data by means of a 2D separable Hamming window in the Fourier domain. Then, we remove such a correlation by using the whitening procedure and compare the accuracy of both despeckling and polarimetric features estimation for the three following cases: uncorrelated, correlated, and decorrelated images. Simulation results showed a steady improvement of performance scores, most notably the equivalent number of looks (ENL), which increased after decorrelation and closely attained the value of the uncorrelated case. Besides ENL, the benefits of the noise decorrelation hold also for polarimetric features, whose estimation accuracy is diminished by the correlation. Also, the trends of simulations were confirmed by qualitative results of experiments carried out on a true Radarsat-2 image

Magma and fluid migration at Yellowstone Caldera in the last three decades inferred from InSAR, leveling and gravity measurements

We studied the Yellowstone caldera geological unrest between 1977 and 2010 by investigating temporal changes in differential Interferometric Synthetic Aperture Radar (InSAR), precise spirit leveling and gravity measurements. The analysis of the 1992–2010 displacement time series, retrieved by applying the SBAS InSAR technique, allowed the identification of three areas of deformation: (i) the Mallard Lake (ML) and Sour Creek (SC) resurgent domes, (ii) a region close to the Northern Caldera Rim (NCR), and (iii) the eastern Snake River Plain (SRP). While the eastern SRP shows a signal related to tectonic deformation, the other two regions are influenced by the caldera unrest. We removed the tectonic signal from the InSAR displacements, and we modeled the InSAR, leveling, and gravity measurements to retrieve the best fitting source parameters. Our findings confirmed the existence of different distinct sources, beneath the brittle-ductile transition zone, which have been intermittently active during the last three decades. Moreover, we interpreted our results in the light of existing seismic tomography studies. Concerning the SC dome, we highlighted the role of hydrothermal fluids as the driving force behind the 1977–1983 uplift; since 1983–1993 the deformation source transformed into a deeper one with a higher magmatic component. Furthermore, our results support the magmatic nature of the deformation source beneath ML dome for the overall investigated period. Finally, the uplift at NCR is interpreted as magma accumulation, while its subsidence could either be the result of fluids migration outside the caldera or the gravitational adjustment of the source from a spherical to a sill-like geometr

Forest Structure Characterization From SAR Tomography at L-Band

Synthetic aperture radar (SAR) remote sensing configurations are able to provide continuous measurements on global scales sensitive to the vertical structure of forests with a high spatial and temporal resolution. Furthermore, the development of tomographic SAR techniques allows the reconstruction of the three-dimensional (3-D) radar reflectivity opening the door for 3-D forest monitoring. However, the link between 3-D radar reflectivity and 3-D forest structure is not yet established. In this sense, this paper introduced a framework that allows a qualitative and quantitative interpretation of physical forest structure from tomographic SAR data at L-band. For this, forest structure is parameterized into a set of a horizontal and a vertical structure index. From inventory data, both indices can be derived from the spatial distribution and the dimensions of the trees. Similarly, two structure indices are derived from the 3-D spatial distribution of the local maxima of the reconstructed 3-D radar reflectivity profiles at L-band. The proposed methodology is tested by means of experimental tomographic L-band data acquired over the temperate forest site of Traunstein in Germany. The obtained horizontal and vertical structure indices are validated against the corresponding estimates obtained from inventory measurements and against the same indices derived from the vertical profiles of airborne Lidar data. The high correlation between the forest structure indices obtained from these three different data sources (expressed by correlation coefficients between 0.75 and 0.87) indicates the potential of the proposed framework

A deep learning solution for height estimation on a forested area based on Pol-TomoSAR data

Forest height and underlying terrain reconstruction is one of the main aims in dealing with forested areas. Theoretically, synthetic aperture radar tomography (TomoSAR) offers the possibility to solve the layover problem, making it possible to estimate the elevation of scatters located in the same resolution cell. This article describes a deep learning approach, named tomographic SAR neural network (TSNN), which aims at reconstructing forest and ground height using multipolarimetric multibaseline (MPMB) SAR data and light detection and ranging (LiDAR)-based data. The reconstruction of the forest and ground height is formulated as a classification problem, in which TSNN, a feedforward network, is trained using covariance matrix elements as input vectors and quantized LiDAR-based data as the reference. In our work, TSNN is trained and tested with P-band MPMB data acquired by ONERA over Paracou region of French Guiana in the frame of the European Space Agency's campaign TROPISAR and LiDAR-based data provided by the French Agricultural Research Center. The novelty of the proposed TSNN is related to its ability to estimate the height with a high agreement with LiDAR-based measurement and actual height with no requirement for phase calibration. Experimental results of different covariance window sizes are included to demonstrate that TSNN conducts height measurement with high spatial resolution and vertical accuracy outperforming the other two TomoSAR methods. Moreover, the conducted experiments on the effects of phase errors in different ranges show that TSNN has a good tolerance for small errors and is still able to precisely reconstruct forest heights

Estado da Arte do Sensoriamento Remoto de Radar: Fundamentos, Sensores, Processamento de Imagens e Aplicações.

Este artigo aborda o estado da arte do sensoriamento remoto por radar e foi elaborado para fazer parte da edição especial de comemoração dos 50 anos desta revista. Neste estudo, é apresentada uma breve introdução sobre os fundamentos do sensoriamento remoto por radar, com destaque para os parâmetros mais importantes de imageamento e da superfície terrestre envolvidos no processo de obtenção de imagens de radar. Ênfase é dada para o comprimento de onda, polarização das ondas eletromagnéticas e geometria de obtenção de imagens (parâmetros de imageamento) e para a umidade de solos e da vegetação, rugosidade do terreno e estrutura da vegetação (parâmetros da superfície terrestre). Em seguida, são apresentados os principais sensores orbitais de radar de abertura sintética que estão atualmente em operação e os principais processamentos digitais de imagens de radar, destacando-se a conversão dos valores digitais para coeficientes de retroespalhamento, os filtros espaciais para redução do ruído speckle, as técnicas de decomposição de imagens e o processamento InSAR. Finalmente, é apresentada uma breve discussão sobre algumas aplicações potenciais, com especial atenção para o monitoramento de derrame de óleo em plataformas continentais, estimativa de biomassa aérea, monitoramento de desmatamento em coberturas florestais tropicais, detecção de áreas de plantio de arroz irrigado e estimativa de umidade de solos

Indoor experiments on polarimetric SAR interferometry

A coherence optimization method, which makes use of polarimetry to enhance the quality of SAR interferograms, has been experimentally tested under laboratory conditions in an anechoic chamber. By carefully selecting the polarization in both images, the resulting interferogram exhibits an improved coherence above the standard HH or VV channel. This higher coherence produces a lower phase variance, thus estimating the underlying topography more accurately. The potential improvement that this technique provides in the generation of digital elevation models (DEM) of non-vegetated natural surfaces has been observed for the first time on some artificial surfaces created with gravel. An experiment on a true outdoor DEM has not been accomplished yet, but the first laboratory results show that the height error for an almost planar surface can be drastically reduced within a wide range of baselines by using the optimization algorithm. This algorithm leads to three possible interferograms associated with statistically independent scattering mechanisms. The phase difference between those interferograms has been employed for extracting the height of vegetation samples. This retrieval technique has been tested on three different samples: maize, rice, and young fir trees. The inverted heights are compared with ground truth for different frequency bands. The estimates are quite variable with frequency, but their complete physical justification is still in progress. Finally, an alternative simplified scheme for the optimization is proposed. The new approach (called polarization subspace method) yields suboptimum results but is more intuitive and has been used for illustrating the working principle of the original optimization algorithm.Peer Reviewe