A fast and accurate basis pursuit denoising algorithm with application to super-resolving tomographic SAR

$L_1$ regularization is used for finding sparse solutions to an underdetermined linear system. As sparse signals are widely expected in remote sensing, this type of regularization scheme and its extensions have been widely employed in many remote sensing problems, such as image fusion, target detection, image super-resolution, and others and have led to promising results. However, solving such sparse reconstruction problems is computationally expensive and has limitations in its practical use. In this paper, we proposed a novel efficient algorithm for solving the complex-valued $L_1$ regularized least squares problem. Taking the high-dimensional tomographic synthetic aperture radar (TomoSAR) as a practical example, we carried out extensive experiments, both with simulation data and real data, to demonstrate that the proposed approach can retain the accuracy of second order methods while dramatically speeding up the processing by one or two orders. Although we have chosen TomoSAR as the example, the proposed method can be generally applied to any spectral estimation problems.Comment: 11 pages, IEEE Transactions on Geoscience and Remote Sensin

InSAR Deformation Analysis with Distributed Scatterers: A Review Complemented by New Advances

Interferometric Synthetic Aperture Radar (InSAR) is a powerful remote sensing technique able to measure deformation of the earth’s surface over large areas. InSAR deformation analysis uses two main categories of backscatter: Persistent Scatterers (PS) and Distributed Scatterers (DS). While PS are characterized by a high signal-to-noise ratio and predominantly occur as single pixels, DS possess a medium or low signal-to-noise ratio and can only be exploited if they form homogeneous groups of pixels that are large enough to allow for statistical analysis. Although DS have been used by InSAR since its beginnings for different purposes, new methods developed during the last decade have advanced the field significantly. Preprocessing of DS with spatio-temporal filtering allows today the use of DS in PS algorithms as if they were PS, thereby enlarging spatial coverage and stabilizing algorithms. This review explores the relations between different lines of research and discusses open questions regarding DS preprocessing for deformation analysis. The review is complemented with an experiment that demonstrates that significantly improved results can be achieved for preprocessed DS during parameter estimation if their statistical properties are used

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

BUILDING RECONSTRUCTION BASED ON A SMALL NUMBER OF TRACKS USING NONPARAMETRIC SAR TOMOGRAPHIC METHODS

Nowadays, the synthetic aperture radar (SAR) tomography (TomoSAR) technique plays a notable role in the 3D reconstruction of urban buildings through several SAR acquisitions with slightly different positions. Nonparametric-based TomoSAR spectral estimation algorithms usually work well when a large number of SAR observations. In this study, with a limited number of SAR images, we have assessed the efficiency of the nonparametric spectral estimation methods, including maximum entropy (ME), singular value decomposition (SVD), linear prediction (LP), Capon, minimum norm (MN), and beamforming (BF) in the reconstruction of the third dimension of urban buildings. The experiments are conducted on both simulated and TerraSAR-X stripmap images to indicate the effectiveness of the LP proposed estimation algorithm. The analysis of the results proves that by minimizing the average output signal power over the antenna array elements, the LP spectral estimation achieves the discrimination of distinct scatterers inside an image pixel. In addition, this low computational estimator improves the sidelobe suppression and the height estimates of the scatterers in the complex multiple-scattering urban environment. Compared to SVD, maximum entropy, Capon, minimum norm, and beamforming, the height of the Eskan tower in Tehran, Iran, obtained with the LP technique, is considerably near to field-based measurement

Non-local methods for InSAR parameters estimation

In the thesis work the nonlocal paradigm has been investigated in the framework of Multitemporal SAR Interferometry, e.g. Differential Interferometry, Tomography, etc., and single InSAR pair, e.g. DEM generation. In the former, Adaptive Multi-Looking methods have been developed for the generation of interferometric data-stacks. Following the nonlocal approach, the proposed methods rely only on similar pixels according to a suitable similarity measure that exploits the stack's temporal information. An hybrid approach that jointly uses the nonlocal paradigm and transform domain filtering has been investigated for InSAR pair phase estimation. On the track of the BM3D and SARBM3D algorithms, different approaches to the filtering in the transform domain are investigated. Furthermore, a novel approach to the similarity computation and filtering, based on a relative-topography content of the interferometric phase rather than its absolute value, is proposed

A novel statistical model for differential synthetic aperture radar tomography

A deterministic differential tomographic synthetic aperture radar (D-TomoSAR) model, based on geometrical derivations and the assumption of accurate phase calibration, has been widely employed for spatially locating and temporally monitoring the point-like scatterers. In this work, we model the phase miscalibration effects of the extended scatters caused by partial correlation, i.e. the decorrelation effects from temporal and spatial changes as well as the residual atmospheric and deformation effects after preprocessing. Starting from the origin of four-dimensional SAR focusing, correlation of the target is analysed and a statistical D-TomoSAR model accounting for partial correlation effects is proposed. Based on the proposed model, a D-TomoSAR stack simulator is designed using Cholesky decomposition. Moreover, a linear minimum mean square error estimator based on the proposed model is developed for estimation of the height and deformation velocity of extended scatterers. Reconstruction results with both simulated data and real data acquired by TerraSAR-X/Tandem-X sensors are provided to demonstrate the effectiveness of the proposed model

Approches tomographiques structurelles pour l'analyse du milieu urbain par tomographie SAR THR : TomoSAR

SAR tomography consists in exploiting multiple images from the same area acquired from a slightly different angle to retrieve the 3-D distribution of the complex reflectivity on the ground. As the transmitted waves are coherent, the desired spatial information (along with the vertical axis) is coded in the phase of the pixels. Many methods have been proposed to retrieve this information in the past years. However, the natural redundancies of the scene are generally not exploited to improve the tomographic estimation step. This Ph.D. presents new approaches to regularize the estimated reflectivity density obtained through SAR tomography by exploiting the urban geometrical structures.La tomographie SAR exploite plusieurs acquisitions d'une même zone acquises d'un point de vue légerement différent pour reconstruire la densité complexe de réflectivité au sol. Cette technique d'imagerie s'appuyant sur l'émission et la réception d'ondes électromagnétiques cohérentes, les données analysées sont complexes et l'information spatiale manquante (selon la verticale) est codée dans la phase. De nombreuse méthodes ont pu être proposées pour retrouver cette information. L'utilisation des redondances naturelles à certains milieux n'est toutefois généralement pas exploitée pour améliorer l'estimation tomographique. Cette thèse propose d'utiliser l'information structurelle propre aux structures urbaines pour régulariser les densités de réflecteurs obtenues par cette technique