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

An accurate method to correct atmospheric phase delay for InSAR with the ERA5 global atmospheric model

Differential SAR Interferometry (DInSAR) has proven its unprecedented ability and merits of monitoring ground deformation on a large scale with centimeter to millimeter accuracy. However, atmospheric artifacts due to spatial and temporal variations of the atmospheric state often affect the reliability and accuracy of its results. The commonly-known Atmospheric Phase Screen (APS) appears in the interferograms as ghost fringes not related to either topography or deformation. Atmospheric artifact mitigation remains one of the biggest challenges to be addressed within the DInSAR community. State-of-the-art research works have revealed that atmospheric artifacts can be partially compensated with empirical models, point-wise GPS zenith path delay, and numerical weather prediction models. In this study, we implement an accurate and realistic computing strategy using atmospheric reanalysis ERA5 data to estimate atmospheric artifacts. With this approach, the Line-of-Sight (LOS) path along the satellite trajectory and the monitored points is considered, rather than estimating it from the zenith path delay. Compared with the zenith delay-based method, the key advantage is that it can avoid errors caused by any anisotropic atmospheric phenomena. The accurate method is validated with Sentinel-1 data in three different test sites: Tenerife island (Spain), Almería (Spain), and Crete island (Greece). The effectiveness and performance of the method to remove APS from interferograms is evaluated in the three test sites showing a great improvement with respect to the zenith-based approach.Peer ReviewedPostprint (published version

A temporal phase coherence estimation algorithm and its application on DInSAR pixel selection

© 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.Pixel selection is a crucial step of all advanced Differential Interferometric Synthetic Aperture Radar (DInSAR) techniques that have a direct impact on the quality of the final DInSAR products. In this paper, a full-resolution phase quality estimator, i.e., the temporal phase coherence (TPC), is proposed for DInSAR pixel selection. The method is able to work with both distributed scatterers (DSs) and permanent scatterers (PSs). The influence of different neighboring window sizes and types of interferograms combinations [both the single-master (SM) and the multi-master (MM)] on TPC has been studied. The relationship between TPC and phase standard deviation (STD) of the selected pixels has also been derived. Together with the classical coherence and amplitude dispersion methods, the TPC pixel selection algorithm has been tested on 37 VV polarization Radarsat-2 images of Barcelona Airport. Results show the feasibility and effectiveness of TPC pixel selection algorithm. Besides obvious improvements in the number of selected pixels, the new method shows some other advantages comparing with the other classical two. The proposed pixel selection algorithm, which presents an affordable computational cost, is easy to be implemented and incorporated into any advanced DInSAR processing chain for high-quality pixels' identification.Peer ReviewedPostprint (author's final draft

Automatic refocus and feature extraction of single-look complex SAR signatures of vessels

In recent years, spaceborne synthetic aperture radar ( SAR) technology has been considered as a complement to cooperative vessel surveillance systems thanks to its imaging capabilities. In this paper, a processing chain is presented to explore the potential of using basic stripmap single-look complex ( SLC) SAR images of vessels for the automatic extraction of their dimensions and heading. Local autofocus is applied to the vessels' SAR signatures to compensate blurring artefacts in the azimuth direction, improving both their image quality and their estimated dimensions. For the heading, the orientation ambiguities of the vessels' SAR signatures are solved using the direction of their ground-range velocity from the analysis of their Doppler spectra. Preliminary results are provided using five images of vessels from SLC RADARSAT-2 stripmap images. These results have shown good agreement with their respective ground-truth data from Automatic Identification System ( AIS) records at the time of the acquisitions.Postprint (published version

An improvement of offset tracking for cross hair (CH) and patch like (PL) elimination and reliability estimation for large deformation monitoring with SAR data

SAR based offset tracking (OT) is an efficient tool for ground deformation observation, and signal to noise ratio (SNR) is its common error indicator. However, ground feature variations often weaken the accuracy of OT. In addition, SNR shows the signal reliability instead of result accuracy. Based on amplitude selection, cubic spline interpolation and double offset detecting, an improved OT method has been proposed in this paper. Subsidence caused by mining and GPS data have been employed to evaluate the performance of this scheme with TerraSAR-X data. The results indicate that patch like (PL) and cross hair(CH) are reduced efficiently, and the error estimated by the proposed method has a higher correlation with real error than SNR in the mountainous area.This research work has been supported by the China Scholarship Council (CSC NO. 201806420035), by the Spanish Ministry of Economy, Industry and Competitiveness (MINECO), the State Research Agency (AEI) and the European Funds for Regional Development (EFRD) under project TEC2017- 85244-C2-1-P. CommSensLab, which is Unidad de Excelencia Maria de Maeztu MDM-2016-0600 financed by the Agencia Estatal de Investigacion, ´ Spain. TerraSAR-X images have been provided by the German Aerospace Center (DLR) in the framework of Project GEO0389 of the TerraSAR-X scientific program.Peer ReviewedPostprint (author's final draft

Métodos numéricos para aplicaciones biomédicas: problemas directo e inverso electromagnéticos

El objetivo básico de la tesis ha sido el desarrollo de algoritmos de reconstrucción que permiten obtener imágenes de cuerpos biológicos de alto contraste. Los datos de entrada de los distintos algoritmos son los campos dispersados por el cuerpo para distintas ondas incidentes a frecuencias de microondas. Además, se han desarrollado algoritmos para el calculo de los campos dispersados por modelos numéricos, o resolución del problema directo, necesarios para el estudio del problema inverso, o reconstrucción. Se incluye una extensa recopilación de resultados, tanto simulados como obtenidos a partir de medidas reales, que permiten conocer el comportamiento de los distintos algoritmos en función de los parámetros de diseño de un tomógrafo.Postprint (published version

SMF-POLOPT: an adaptive multitemporal pol(DIn)SAR filtering and phase optimization algorithm for PSI applications

© 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.Speckle noise and decorrelation can hamper the application and interpretation of PolSAR images. In this paper, a new adaptive multitemporal Pol(DIn)SAR filtering and phase optimization algorithm is proposed to address these limitations. This algorithm first categorizes and adaptively filters permanent scatterer (PS) and distributed scatterer (DS) pixels according to their polarimetric scattering mechanisms [i.e., the scattering-mechanism-based filtering (SMF)]. Then, two different polarimetric DInSAR (POLDInSAR) phase OPTimization methods are applied separately on the filtered PS and DS pixels (i.e., POLOPT). Finally, an inclusive pixel selection approach is used to identify high-quality pixels for ground deformation estimation. Thirty-one full-polarization Radarsat-2 SAR images over Barcelona (Spain) and 31 dual-polarization TerraSAR-X images over Murcia (Spain) have been used to evaluate the performance of the proposed algorithm. The PolSAR filtering results show that the speckle of PolSAR images has been well reduced with the preservation of details by the proposed SMF. The obtained ground deformation monitoring results have shown significant improvements, about ×7.2 (the full-polarization case) and ×3.8 (the dual-polarization case) with respect to the classical full-resolution single-pol amplitude dispersion method, on the valid pixels' densities. The excellent PolSAR filtering and ground deformation monitoring results achieved by the adaptive Pol(DIn)SAR filtering and phase optimization algorithm (i.e., the SMF-POLOPT) have validated the effectiveness of this proposed scheme.Peer ReviewedPostprint (author's final draft

Sensitivity equations and calibration requirements on airborne interferometry

The usage of an airborne interferometric SAR to mass map production is strongly dependent on the capability of performing the data processing in a near automatic mode. One of the main limitations is related to both the accurate calibration of the system parameters, its stability from flight to flight and the correct in-flight recording of aircraft position and attitude. The unstable movements of an airborne SAR platform, if recorded in an accurate manner, can be corrected during the processing step. Measurement errors and time drifts lead to location errors in the final DEM. The system requirements in both accuracy and stability to fulfil a quality-mapping requirement can be resolved from the sensitivity equations. The same equations can be also used to calibrate the system parameters from the location errors on the imaging of well-known targets, usually corner reflectors.Peer ReviewedPostprint (published version

Estimation of Azimuth Phase Undulations with Multisquint Processing in Airborne Interferometric SAR Images

This letter presents a technique to detect and correct phase errors appearing in interferometric airborne synthetic aperture radar (SAR) systems due to the lack of precision in the navigation system. The technique is based on a multisquint processing approach, i.e., by processing the same image pairs with different squint angles we can combine the information of different interferograms to obtain the desired phase correction. Airborne single-pass interferometric data from the Deutsches Zentrum für Luft- und Raumfahrt (DLR) Experimental airborne SAR is used to validate the metho

Study of sea clutter influence in ship classification algorithms based on Polarimetric SAR Inteferometry

This paper is focused to evaluate the influence of sea clutter in the performance of ship classification algorithms based on single-pass Polarimetric SAR Interferometry (PolInSAR). For such purpose, series of numerical simulations have been carried out with GRECOSAR, the SAR simulator of complex targets developed by UPC. There, different types of vessels have been considered for a TerraSAR-X like sensor and a sea surface following the two-scale wave approach. The quality of ship discrimination has been quantitatively evaluated with a novel identification method that exploits the particular scattering properties of ships. The results show that the presence of clutter does not notably drop identification performance, despite negative matches can be observed in some particular situations. But the requirement of single-pass interferometric capabilities is not achieved by any of the existing orbital system. This drawback can difficult the validation of what has been observed in simulation environments and can be one of the most limiting factors for the practical implementation of these techniques. Ideas and possible solutions to relax the system requirements are preliminary discussed.Postprint (published version