Wavelet Operators and Multiplicative Observation Models - Application to Change-Enhanced Regularization of SAR Image Time Series

This paper first provides statistical properties of wavelet operators when the observation model can be seen as the product of a deterministic piecewise regular function (signal) and a stationary random field (noise). This multiplicative observation model is analyzed in two standard frameworks by considering either (1) a direct wavelet transform of the model or (2) a log-transform of the model prior to wavelet decomposition. The paper shows that, in Framework (1), wavelet coefficients of the time series are affected by intricate correlation structures which affect the signal singularities. Framework (2) is shown to be associated with a multiplicative (or geometric) wavelet transform and the multiplicative interactions between wavelets and the model highlight both sparsity of signal changes near singularities (dominant coefficients) and decorrelation of speckle wavelet coefficients. The paper then derives that, for time series of synthetic aperture radar data, geometric wavelets represent a more intuitive and relevant framework for the analysis of smooth earth fields observed in the presence of speckle. From this analysis, the paper proposes a fast-and-concise geometric wavelet based method for joint change detection and regularization of synthetic aperture radar image time series. In this method, geometric wavelet details are first computed with respect to the temporal axis in order to derive generalized-ratio change-images from the time series. The changes are then enhanced and speckle is attenuated by using spatial bloc sigmoid shrinkage. Finally, a regularized time series is reconstructed from the sigmoid shrunken change-images. An application of this method highlights the relevancy of the method for change detection and regularization of SENTINEL-1A dual-polarimetric image time series over Chamonix-Mont-Blanc test site

Wavelet Operators and Multiplicative Observation Models -Application to SAR Image Time Series Analysis

International audienceThis paper first provides statistical properties of wavelet operators when the observation model can be seen as the product of a deterministic piece-wise regular function (signal) and a stationary random field (noise). This multiplicative observation model is analyzed in two standard frameworks by considering either (1) a direct wavelet transform of the model or (2) a log-transform of the model prior to wavelet decomposition. The paper shows that, in Framework (1), wavelet coefficients of the time series are affected by intricate correlation structures which blur signal singularities. Framework (2) is shown to be associated with a multiplicative (or geometric) wavelet transform and the multiplicative interactions between wavelets and the model highlight both sparsity of signal changes near singularities (dominant coefficients) and decorre-lation of speckle wavelet coefficients. The paper then derives that, for time series of synthetic aperture radar data, geometric wavelets represent a more intuitive and relevant framework for the analysis of smooth earth fields observed in the presence of speckle. From this analysis, the paper proposes a fast-and-concise geometric wavelet based method for joint change detection and regularization of synthetic aperture radar image time series. In this method, geometric wavelet details are first computed with respect to the temporal axis in order to derive generalized-ratio change-images from the time series. The changes are then enhanced and speckle is attenuated by using spatial block sigmoid shrinkage. Finally, a regularized time series is reconstructed from the sigmoid shrunken change-images. Some applications highlight relevancy of the method for the analysis of SENTINEL-1A and TerraSAR-X image time series over Chamonix-Mont-Blanc

On the extension of multidimensional speckle noise model from single-look to multilook SAR imagery

Speckle noise represents one of the major problems when synthetic aperture radar (SAR) data are considered. Despite the fact that speckle is caused by the scattering process itself, it must be considered as a noise source due to the complexity of the scattering process. The presence of speckle makes data interpretation difficult, but it also affects the quantitative retrieval of physical parameters. In the case of one-dimensional SAR systems, speckle is completely determined by a multiplicative noise component. Nevertheless, for multidimensional SAR systems, speckle results from the combination of multiplicative and additive noise components. This model has been first developed for single-look data. The objective of this paper is to extend the single-look data model to define a multilook multidimensional speckle noise model. The asymptotic analysis of this extension, for a large number of averaged samples, is also considered to assess the model properties. Details and validation of the multilook multidimensional speckle noise model are provided both theoretically and by means of experimental SAR data acquired by the experimental synthetic aperture radar system, operated by the German Aerospace Center.Peer Reviewe

CFAR Ship Detection in Polarimetric Synthetic Aperture Radar Images Based on Whitening Filter

Polarimetric whitening filter (PWF) can be used to filter polarimetric synthetic aperture radar (PolSAR) images to improve the contrast between ships and sea clutter background. For this reason, the output of the filter can be used to detect ships. This paper deals with the setting of the threshold over PolSAR images filtered by the PWF. Two parameter-constant false alarm rate (2P-CFAR) is a common detection method used on whitened polarimetric images. It assumes that the probability density function (PDF) of the filtered image intensity is characterized by a log-normal distribution. However, this assumption does not always hold. In this paper, we propose a systemic analytical framework for CFAR algorithms based on PWF or multi-look PWF (MPWF). The framework covers the entire log-cumulants space in terms of the textural distributions in the product model, including the constant, gamma, inverse gamma, Fisher, beta, inverse beta, and generalized gamma distributions (GΓDs). We derive the analytical forms of the PDF for each of the textural distributions and the probability of false alarm (PFA). Finally, the threshold is derived by fixing the false alarm rate (FAR). Experimental results using both the simulated and real data demonstrate that the derived expressions and CFAR algorithms are valid and robust

Fuzzy superpixels for polarimetric SAR images classification

Superpixels technique has drawn much attention in computer vision applications. Each superpixels algorithm has its own advantages. Selecting a more appropriate superpixels algorithm for a specific application can improve the performance of the application. In the last few years, superpixels are widely used in polarimetric synthetic aperture radar (PolSAR) image classification. However, no superpixel algorithm is especially designed for image classification. It is believed that both mixed superpixels and pure superpixels exist in an image.Nevertheless, mixed superpixels have negative effects on classification accuracy. Thus, it is necessary to generate superpixels containing as few mixed superpixels as possible for image classification. In this paper, first, a novel superpixels concept, named fuzzy superpixels, is proposed for reducing the generation of mixed superpixels.In fuzzy superpixels ,not al lpixels are assigned to a corresponding superpixel. We would rather ignore the pixels than assigning them to improper superpixels. Second,a new algorithm, named FuzzyS(FS),is proposed to generate fuzzy superpixels for PolSAR image classification. Three PolSAR images are used to verify the effect of the proposed FS algorithm. Experimental results demonstrate the superiority of the proposed FS algorithm over several state-of-the-art superpixels algorithms