A Tutorial on Speckle Reduction in Synthetic Aperture Radar Images

Speckle is a granular disturbance, usually modeled as a multiplicative noise, that affects synthetic aperture radar (SAR) images, as well as all coherent images. Over the last three decades, several methods have been proposed for the reduction of speckle, or despeckling, in SAR images. Goal of this paper is making a comprehensive review of despeckling methods since their birth, over thirty years ago, highlighting trends and changing approaches over years. The concept of fully developed speckle is explained. Drawbacks of homomorphic filtering are pointed out. Assets of multiresolution despeckling, as opposite to spatial-domain despeckling, are highlighted. Also advantages of undecimated, or stationary, wavelet transforms over decimated ones are discussed. Bayesian estimators and probability density function (pdf) models in both spatial and multiresolution domains are reviewed. Scale-space varying pdf models, as opposite to scale varying models, are promoted. Promising methods following non-Bayesian approaches, like nonlocal (NL) filtering and total variation (TV) regularization, are reviewed and compared to spatial- and wavelet-domain Bayesian filters. Both established and new trends for assessment of despeckling are presented. A few experiments on simulated data and real COSMO-SkyMed SAR images highlight, on one side the costperformance tradeoff of the different methods, on the other side the effectiveness of solutions purposely designed for SAR heterogeneity and not fully developed speckle. Eventually, upcoming methods based on new concepts of signal processing, like compressive sensing, are foreseen as a new generation of despeckling, after spatial-domain and multiresolution-domain method

An interferometric phase noise reduction method based on modified denoising convolutional neural network

Traditional interferometric synthetic aperture radar (InSAR) denoising methods normally try to estimate the phase fringes directly from the noisy interferogram. Since the statistics of phase noise are more stable than the phase corresponding to complex terrain, it could be easier to estimate the phase noise. In this paper, phase noises rather than phase fringes are estimated first, and then they are subtracted from the noisy interferometric phase for denoising. The denoising convolutional neural network (DnCNN) is introduced to estimate phase noise and then a modified network called IPDnCNN is constructed for the problem. Based on the IPDnCNN, a novel interferometric phase noise reduction algorithm is proposed, which can reduce phase noise while protecting fringe edges and avoid the use of filter windows. Experimental results using simulated and real data are provided to demonstrate the effectiveness of the proposed method

Range Spectral Filtering in SAR Interferometry: Methods and Limitations

A geometrical decorrelation constitutes one of the sources of noise present in Synthetic Aperture Radar (SAR) interferograms. It comes from the different incidence angles of the two images used to form the interferograms, which cause a spectral (frequency) shift between them. A geometrical decorrelation must be compensated by a specific filtering technique known as range filtering, the goal of which is to estimate this spectral displacement and retain only the common parts of the images’ spectra, reducing the noise and improving the quality of the interferograms. Multiple range filters have been proposed in the literature. The most widely used methods are an adaptive filter approach, which estimates the spectral shift directly from the data; a method based on orbital information, which assumes a constant-slope (or flat) terrain; and slope-adaptive algorithms, which consider both orbital information and auxiliary topographic data. Their advantages and limitations are analyzed in this manuscript and, additionally, a new, more refined approach is proposed. Its goal is to enhance the filtering process by automatically adapting the filter to all types of surface variations using a multi-scale strategy. A pair of RADARSAT-2 images that mapped the mountainous area around the Etna volcano (Italy) are used for the study. The results show that filtering accuracy is improved with the new method including the steepest areas and vegetation-covered regions in which the performance of the original methods is limited.This work was supported by the Spanish Ministry of Science and Innovation (State Agency of Research, AEI) and the European Funds for Regional Development (ERFD) under Projects PID2020-117303GB-C21 and PID2020-117303-C22

Improved Goldstein Interferogram Filter Based on Local Fringe Frequency Estimation

The quality of an interferogram, which is limited by various phase noise, will greatly affect the further processes of InSAR, such as phase unwrapping. Interferometric SAR (InSAR) geophysical measurements’, such as height or displacement, phase filtering is therefore an essential step. In this work, an improved Goldstein interferogram filter is proposed to suppress the phase noise while preserving the fringe edges. First, the proposed adaptive filter step, performed before frequency estimation, is employed to improve the estimation accuracy. Subsequently, to preserve the fringe characteristics, the estimated fringe frequency in each fixed filtering patch is removed from the original noisy phase. Then, the residual phase is smoothed based on the modified Goldstein filter with its parameter alpha dependent on both the coherence map and the residual phase frequency. Finally, the filtered residual phase and the removed fringe frequency are combined to generate the filtered interferogram, with the loss of signal minimized while reducing the noise level. The effectiveness of the proposed method is verified by experimental results based on both simulated and real data

Multiresolution Detection of Persistent Scatterers: A Performance Comparison Between Multilook GLRT and CAESAR

Persistent scatterers (PS) interferometry tools are extensively used for the monitoring of slow, long-term ground deformation. High spatial resolution is typically required in urban areas to cope with the variability of the signal, whereas in rural regions, multilook shall be implemented to improve the coverage of monitored areas. Along this line, SqueeSAR and later Component extrAction and sElection SAR (CAESAR) were introduced for the monitoring of both persistent and (decorrelating) distributed scatterers (DS). Multilook generalized likelihood ratio test (MGLRT) is a detector derived in the context of tomographic SAR processing that has been investigated for a fixed multilook degree. In this work, we address MGLRT and CAESAR in the multiresolution context characterized by a spatially variable multilook degree. We compare the two schemes for the multiresolution selection of PS and DS, highlighting the pros and cons of each scheme, particularly the peculiarities of CAESAR that have important implications at the implementation stage. A performance analysis of both detectors in case of model mismatch is also addressed. Experiments carried out with data acquired by the COSMO-SkyMed constellation support both the theoretical argumentation and the results achieved by resorting to Monte Carlo simulations

A Sparsity-Based InSAR Phase Denoising Algorithm Using Nonlocal Wavelet Shrinkage

An interferometric synthetic aperture radar (InSAR) phase denoising algorithm using the local sparsity of wavelet coefficients and nonlocal similarity of grouped blocks was developed. From the Bayesian perspective, the double-l1 norm regularization model that enforces the local and nonlocal sparsity constraints was used. Taking advantages of coefficients of the nonlocal similarity between group blocks for the wavelet shrinkage, the proposed algorithm effectively filtered the phase noise. Applying the method to simulated and acquired InSAR data, we obtained satisfactory results. In comparison, the algorithm outperformed several widely-used InSAR phase denoising approaches in terms of the number of residues, root-mean-square errors and other edge preservation indexes