Gabriel Vasile

Intensity-driven adaptive-neighborhood technique for polarimetric and interferometric SAR parameters estimation

In this paper, a new method to filter coherency matrices of polarimetric or interferometric data is presented. For each pixel, an adaptive neighborhood (AN) is determined by a region growing technique driven exclusively by the intensity image information. All the available intensity images of the polarimetric and interferometric terms are fused in the region growing process to ensure the validity of the stationarity assumption. Afterward, all the pixels within the obtained AN are used to yield the filtered values of the polarimetric and interferometric coherency matrices, which can be derived either by direct complex multilooking or from the locally linear minimum mean-squared error (LLMMSE) estimator. The entropy/alpha/anisotropy decomposition is then applied to the estimated polarimetric coherency matrices, and coherence optimization is performed on the estimated polarimetric and interferometric coherency matrices. Using this decomposition, unsupervised classification for land applications by an iterative algorithm based on a complex Wishart density function is also applied. The method has been tested on airborne high-resolution polarimetric interferometric synthetic aperture radar (POL-InSAR) images (Oberpfaffenhofen area-German Space Agency). For comparison purposes, the two estimation techniques (complex multilooking and LLMMSE) were tested using three different spatial supports: a fix-sized symmetric neighborhood (boxcar filter), directional nonsymmetric windows, and the proposed AN. Subjective and objective performance analysis, including coherence edge detection, receiver operating characteristics plots, and bias reduction tables, recommends the proposed algorithm as an effective POL-InSAR postprocessing technique.

Hierarchical Segmentation of Polarimetric SAR Images Using Heterogeneous Clutter Models

In this paper, heterogeneous clutter models are used to describe polarimetric synthetic aperture radar (PolSAR) data. The KummerU distribution is introduced to model the PolSAR clutter. Then, a detailed analysis is carried out to evaluate the potential of this new multivariate distribution. It is implemented in a hierarchical maximum likelihood segmentation algorithm. The segmentation results are shown on both synthetic and high-resolution PolSAR data at the X- and L-bands. Finally, some methods are examined to determine automatically the “optimal” number of segments in the final partition.

Combining Airborne Photographs and Spaceborne SAR Data to Monitor Temperate Glaciers: Potentials and Limits

Monitoring temperate glacier activity has become more and more necessary for economical and security reasons and as an indicator of the local effects of global climate change. Remote sensing data provide useful information on such complex geophysical objects, but they require specific processing techniques to cope with the difficult context of moving and changing features in high-relief areas. This paper presents the first results of a project involving four laboratories developing and combining specific methods to extract information from optical and synthetic aperture radar (SAR) data. Two different information sources are processed, namely: 1) airborne photography and 2) spaceborne C-band SAR interferometry. The difficulties and limitations of their processing in the context of Alpine glaciers are discussed and illustrated on two glaciers located in the Mont-Blanc area. The results obtained by aerial triangulation techniques provide digital terrain models with an accuracy that is better than 30 cm, which is compatible with the computation of volume balance and useful for precise georeferencing and slope measurement updating. The results obtained by SAR differential interferometry using European Remote Sensing Satellite images show that it is possible to measure temperate glacier surface velocity fields from October to April in one-day interferograms with approximately 20-m ground sampling. This allows to derive ice surface strain rate fields required to model the glacier flow. These different measurements are complementary to results obtained during the summer from satellite optical data and ground measurements that are available only in few accessible points

Monitoring Temperate Glacier Displacement by Multi-Temporal TerraSAR-X Images and Continuous GPS Measurements

A new generation of space-borne SAR sensors were launched in 2006-2007 with ALOS, TerraSAR-X, COSMO-Sky-Med and RadarSat-2 satellites. The data available in different bands (L, C and X bands), with High Resolution (HR) or multi-polarization modes offer new possibilities to monitor glacier displacement and surface evolution by SAR remote sensing. In this paper, the first results obtained with TerraSAR-X HR SAR image time series acquired over the temperate glaciers of the Chamonix Mont-Blanc test site are presented. This area involves well-known temperate glaciers which have been monitored and instrumented i.e. stakes for annual displacement/ablation, GPS for surface displacement and cavitometer for basal displacement, for more than 50 years. The potential of 11-day repeated X-band HR SAR data for Alpine glacier monitoring is investigated by a combined use of in situ measurements and multi-temporal images. Interpretations of HR images, analysis of interferometric pairs and performance assessments of target/texture tracking methods for glacier motion estimation are presented. The results obtained with four time series covering the Chamonix Mont-Blanc glaciers over one year show that the phase information is rarely preserved after 11 days on such glaciers, whereas the high resolution intensity information allows the main glacier features to be observed and displacement fields on the textured areas to be derived.

Statistical Classification for Heterogeneous Polarimetric SAR Images

This paper presents a general approach for high- resolution polarimetric SAR data classification in heterogeneous clutter, based on a statistical test of equality of covariance matrices. The Spherically Invariant Random Vector (SIRV) model is used to describe the clutter. Several distance measures, including classical ones used in standard classification methods, can be derived from the general test. The new approach provide a threshold over which pixels are rejected from the image, meaning they are not sufficiently “close” from any existing class. A distance measure using this general approach is derived and tested on a high-resolution polarimetric data set acquired by the ONERA RAMSES system. It is compared to the results of the classical H-α decomposition and Wishart classifier under Gaussian and SIRV assumption. Results show that the new approach rejects all pixels from heterogeneous parts of the scene and classifies its Gaussian parts.

Short-Range Wideband FMCW Radar for Millimetric Displacement Measurements

The frequency-modulated continuous-wave (FMCW) radar is an alternative to the pulse radar when the distance to the target is short. Typical FMCW radar implementations have a homodyne architecture transceiver which limits the performances for short-range applications: The beat frequency can be relatively small and placed in the frequency range affected by the specific homodyne issues (dc offset, self-mixing, and 1/f noise). In addition, one classical problem of an FMCW radar is that the voltage-controlled oscillator adds a certain degree of nonlinearity which can cause a dramatic resolution degradation for wideband sweeps. This paper proposes a short-range X-band FMCW radar platform which solves these two problems by using a heterodyne transceiver and a wideband nonlinearity correction algorithm based on high-order ambiguity functions and time resampling. The platforms displacement measurement capability was tested on range profiles and synthetic aperture radar images acquired for various targets. The displacements were computed from the interferometric phase, and the measurement errors were situated below 0.1 mm for metal bar targets placed at a few meters from the radar.

General adaptive-neighborhood technique for improving synthetic aperture radar interferometric coherence estimation

A new method for filtering the coherence map issued from synthetic aperture radar (SAR) interferometric data is presented. For each pixel of the interferogram, an adaptive neighborhood is determined by a region-growing technique driven by the information provided by the amplitude images. Then pixels in the derived adaptive neighborhood are complex averaged to yield the filtered value of the coherence, after a phase-compensation step is performed. An extension of the algorithm is proposed for polarimetric interferometric SAR images. The proposed method has been applied to both European Remote Sensing (ERS) satellite SAR images and airborne high-resolution polarimetric interferometric SAR images. Both subjective and objective performance analysis, including coherence edge detection, shows that the proposed method provides better results than the standard phase-compensated fixed multilook filter and the Lee adaptive coherence filter.

Hierarchical segmentation of Polarimetric SAR images using heterogeneous clutter models

In this paper, heterogeneous clutter models are introduced to describe Polarimetric Synthetic Aperture Radar (PolSAR) data. Based on the Spherically Invariant Random Vectors (SIRV) estimation scheme, the scalar texture parameter and the normalized covariance matrix are extracted. If the texture parameter is modeled by a Fisher PDF, the observed target scattering vector follows a KummerU PDF. Then, this PDF is implemented in a hierarchical segmentation algorithm. Segmentation results are shown on high resolution PolSAR data at L and}band.

Polarimetric Incoherent Target Decomposition by Means of Independent Component Analysis

This paper presents an alternative approach for polarimetric incoherent target decomposition (ICTD) dedicated to the analysis of very high-resolution polarimetric synthetic aperture radar (POLSAR) images. Given the non-Gaussian nature of the heterogeneous POLSAR clutter due to the increase in spatial resolution, the conventional methods based on the eigenvector target decomposition can ensure uncorrelation of the derived backscattering components at most. By introducing the independent component analysis (ICA) in lieu of the eigenvector decomposition, our method is rather deriving statistically independent components. The adopted algorithm, i.e., FastICA, uses the non-Gaussianity of the components as the criterion for their independence. Considering the eigenvector decomposition as being analogs to the principal component analysis (PCA), we propose the generalization of the ICTD methods to the level of the blind source separation (BSS) techniques (comprising both PCA and ICA). The proposed method preserves the invariance properties of the conventional ones, appearing to be robust both with respect to the rotation around the line of sight and to the change of the polarization basis. The efficiency of the method is demonstrated comparatively using POLSAR RAMSES X-band and ALOS L-band data sets. The main differences with respect to the conventional methods are mostly found in the behavior of the second most dominant component, which is not necessarily orthogonal to the first one. The potential of retrieving nonorthogonal mechanisms is moreover demonstrated using synthetic data. On the expense of a negligible entropy increase, the proposed method is capable of retrieving the edge diffraction of an elementary trihedral by recognizing dipole as the second component.

High-Resolution SAR Interferometry: Estimation of Local Frequencies in the Context of Alpine Glaciers

Synthetic aperture radar (SAR) interferometric data offer the opportunity to measure temperate glacier surface topography and displacement. The increase of the resolution provided by the most recent SAR systems has some critical implications. For instance, a reliable estimate of the phase gradient can only be achieved by using interferogram local frequencies. In this paper, an original two-step method for estimating local frequencies is proposed. The 2-D phase signal is considered to have two deterministic components corresponding to low-resolution (LR) fringes and high-resolution (HR) patterns due to the local microrelief, respectively. The first step of the proposed algorithm consists in the LR phase flattening. In the second step, the local HR frequencies are estimated from the phase 2-D autocorrelation function computed on adaptive neighborhoods. This neighborhood is the set of connected pixels belonging to the same HR spatial feature and respecting the ldquolocal stationarityrdquo hypothesis. Results with both simulated TerraSAR-X interferograms and real airborne E-SAR images are presented to illustrate the potential of the proposed method.