One of the main topics in synthetic aperture radar (SAR) tomography (TomoSAR) is the estimation of the vertical structures‚ location, which scatter the field back towards the sensor, constrained to a reduced number of passes. Moreover, the introduction of artifacts and the increase of the ambiguity levels due to irregular sampling, consequence of non-uniform acquisition constellations, complicate the accurate estimation of the source parameters. Pursuing the alleviation of such drawbacks, the use of statistical regularization approaches, based on the maximum-likelihood (ML) estimation theory, have been successfully demonstrated in the previous related studies. However, these techniques are constrained to the assumption that the probability density function (pdf) of the observed data is Gaussian. In this paper, in order to solve the ill-posed non-linear TomoSAR inverse problem, we relax this assumption and apply the weighted covariance fitting (WCF) criterion instead. The latter alleviates the previously mentioned drawbacks, and retrieves a power spectrum pattern (PSP) with an outline more similar to the expected one, i.e., recovered using matched spatial filtering (MSF) with a higher number of tracks. First, we present the mathematical background of the related regularization methods, adapted to solve the TomoSAR inverse problem, from which we derive our novel technique, named WISE (WCF-based Iterative Spectral Estimator). Then, the differences and similarities between the addressed regularization approaches are discussed, besides their main advantages and disadvantages. Finally, the implementation details of WISE are treated, along with simulated examples and experimental results gotten from a forested test site.
14 March 2025 in IEEE Geoscience and Remote Sensing Magazine
The main goal of the TerraSAR-X Add-On for Digital Elevation Measurements (TanDEM-X) mission is the generation of a global digital elevation model (DEM) of unprecedented accuracy and coverage. The global TanDEM-X DEM product became available in 2016, surpassed all expectations, and became a reference for a wide range of Earth science, commercial, and geospatial applications. In addition, new information products, such as DEM change maps (DCMs), have been developed and are available to the geoscience and remote sensing community. Beyond the operational products, new science applications have been demonstrated and are summarized in this article, along with experimental data acquisitions. This article also aims to provide an overview of science activities with TanDEM-X data and science data acquisitions planned for the coming years.
Among the objectives of the upcoming space missions Tandem-L and BIOMASS, is the 3-D representation of the global forest structure via synthetic aperture radar (SAR) tomography (TomoSAR). To achieve such a goal, modern approaches suggest solving the TomoSAR inverse problems by exploiting polarimetric diversity and structural model properties of the different scattering mechanisms. This way, the related tomographic imaging problems are treated in descriptive regularization settings, applying modern non-parametric spatial spectral analysis (SSA) techniques. Nonetheless, the achievable resolution of the commonly performed SSA-based estimators highly depends on the span of the tomographic aperture; furthermore, irregular sampling and non-uniform constellations sacrifice the attainable resolution, introduce artifacts and increase ambiguity. Overcoming these drawbacks, in this paper, we address a new multi-stage iterative technique for feature-enhanced TomoSAR imaging that aggregates the virtual adaptive beamforming (VAB)-based SSA approach, with the wavelet domain thresholding (WDT) regularization framework, which we refer to as WAVAB (WDT-refined VAB). First, high resolution imagery is recovered applying the descriptive experiment design regularization (DEDR)-inspired reconstructive processing. Next, the additional resolution enhancement with suppression of artifacts is performed, via the WDT-based sparsity promoting refinement in the wavelet transform (WT) domain. Additionally, incorporation of the sum of Kronecker products (SKP) decomposition technique at the pre-processing stage, improves ground and canopy separation and allows for the utilization of different better adapted TomoSAR imaging techniques, on the ground and canopy structural components, separately. The feature enhancing capabilities of the novel robust WAVAB TomoSAR imaging technique are corroborated through the processing of airborne data of the German Aerospace Center (DLR), providing detailed volume height profiles reconstruction, as an alternative to the competing non-parametric SSA-based methods.
6 October 2021 in IEEE Transactions on Geoscience and Remote Sensing
Synthetic aperture radar (SAR) polarimetry can improve change detection in terms of detection capabilities. In this work, we are proposing to extend the idea of target decomposition to changes affecting partial targets. This will allow the separation of polarimetric-dependent changes, providing extra information that can be used to better understand the processes affecting the targets. Three models for changes are proposed and compared. The methodologies are based on Lagrangian optimizations of distinct operators built using quadratic forms for a power ratio and a power difference. The optimizations can be accomplished by diagonalizations of specific matrices derived from polarimetric covariance matrices. These are, therefore, spectral decompositions of an appropriate matrix which we define as change matrix. The theoretical validity of the models is assessed using Monte Carlo simulations. Additionally, we perform real data validation exploiting L-band quad-polarimetric data from the E-SAR (DLR) SARTOM 2006 campaign and ALOS PALSAR (JAXA) acquisitions in Morecombe Bay (U.K.). We observed that the two algorithms based on power difference allow to decompose the change into the minimal set of scattering mechanisms (SMs) that have been added or removed from the scene. The two algorithms differ on the initial assumption on the change. On the other hand, the ratio operator provides a better detection performance although the eigenvalues do not correspond to meaningful SMs. A combination of the three methodologies can, therefore, improve detection and classification of changes.
Synthetic aperture radar (SAR) tomography (TomoSAR) is a multibaseline interferometric technique that estimates the power spectrum pattern (PSP) along the perpendicular to the line-of-sight (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.
The aim of this work is to exploit the MUSIC algorithm performance in order to enhance target separability in range and azimuth, i.e. achieve point targets separation inside a resolution cell. Simulations have been done in order to plan and check the feasibility of a super-resolution experiment that took place in September 2006 on the test site of Oberpfaffenhofen (Germany). The data set has been acquired with the E-SAR system of the DLR in X-band. The targets to be separated were seven small corner reflectors that have been placed in a way that their response falls in one or, at maximum, two resolution cells of the standard Fourier SAR image. A post-processing implementation of the MUSIC algorithm has been proposed allowing, in the already focused SAR image, to retrieve the targets geometry. Conditions and analysis of the results have been carried out.