A Variational Model for PolSAR Data Speckle Reduction Based on the Wishart Distribution

In this paper, we propose a variational model for polarimetric synthetic aperture radar (PolSAR) data speckle reduction, which is based on the complex Wishart distribution of the covariance or coherency matrix and multichannel total variation (TV) regularization defined for complex-valued matrices. By assuming the TV regularization to be a prior and taking the statistical distribution of the covariance matrix in each resolution element into account, the variational model for PolSAR covariance data speckle suppression, named WisTV-C, is derived from the maximum a posteriori estimate. A similar variational model for PolSAR coherency data speckle reduction, named WisTV-T, is also obtained. As far as we know, this is the first variational model for the whole PolSAR covariance or coherency matrix data despeckling. Since the model is nonconvex, a convex relaxation iterative algorithm is designed to solve the variational problem, based on the variable splitting and alternating minimization techniques. Experimental results on both simulated and real PolSAR data demonstrate that the proposed approach notably removes speckles in the extended uniform areas and, meanwhile, better preserves the spatial resolution, the details such as edges and point scatterers, and the polarimetric scattering characteristics, compared with other methods.

[1]  Wentao An,et al.  Nonlocal Filtering for Polarimetric SAR Data: A Pretest Approach , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[2]  Alexandre d'Aspremont,et al.  Model Selection Through Sparse Max Likelihood Estimation Model Selection Through Sparse Maximum Likelihood Estimation for Multivariate Gaussian or Binary Data , 2022 .

[3]  H. Zebker,et al.  Imaging radar polarization signatures: Theory and observation , 1987 .

[4]  Carlos López-Martínez,et al.  Model-Based Polarimetric SAR Speckle Filter , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[5]  Ke Chen,et al.  On High-Order Denoising Models and Fast Algorithms for Vector-Valued Images , 2010, IEEE Transactions on Image Processing.

[6]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[7]  Samuel Foucher,et al.  Polarimetric SAR image filtering with trace-based partial differential equations , 2007 .

[8]  Corina da Costa Freitas,et al.  Speckle reduction in polarimetric SAR imagery with stochastic distances and nonlocal means , 2013, Pattern Recognit..

[9]  T. Chan,et al.  Fast dual minimization of the vectorial total variation norm and applications to color image processing , 2008 .

[10]  Fawwaz T. Ulaby,et al.  SAR speckle reduction using wavelet denoising and Markov random field modeling , 2002, IEEE Trans. Geosci. Remote. Sens..

[11]  Gilles Aubert,et al.  A Variational Approach to Removing Multiplicative Noise , 2008, SIAM J. Appl. Math..

[12]  Christin Wirth,et al.  Geometric Level Set Methods In Imaging Vision And Graphics , 2016 .

[13]  Stanley Osher,et al.  Multiplicative Denoising and Deblurring: Theory and Algorithms , 2003 .

[14]  He Li,et al.  PolSAR image speckle reduction algorithm based on TV-PDE , 2011, International Symposium on Multispectral Image Processing and Pattern Recognition.

[15]  Florence Tupin,et al.  NL-SAR: A Unified Nonlocal Framework for Resolution-Preserving (Pol)(In)SAR Denoising , 2015, IEEE Transactions on Geoscience and Remote Sensing.

[16]  Daniel Cremers,et al.  The Natural Vectorial Total Variation Which Arises from Geometric Measure Theory , 2012, SIAM J. Imaging Sci..

[17]  N. R. Goodman.  Statistical analysis based on a certain multivariate complex Gaussian distribution , 1963 .

[18]  Jong-Sen Lee,et al.  Polarimetric SAR speckle filtering and its implication for classification , 1999, IEEE Trans. Geosci. Remote. Sens..

[19]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[20]  Maoguo Gong,et al.  SAR Image Despeckling Based on Local Homogeneous-Region Segmentation by Using Pixel-Relativity Measurement , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[21]  N. L. Johnson,et al.  Linear Statistical Inference and Its Applications , 1966 .

[22]  Samuel Foucher,et al.  Analysis, Evaluation, and Comparison of Polarimetric SAR Speckle Filtering Techniques , 2014, IEEE Transactions on Image Processing.

[23]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[24]  Carlos López-Martínez,et al.  Polarimetric SAR speckle noise model , 2003, IEEE Trans. Geosci. Remote. Sens..

[25]  Jian Yang,et al.  Speckle filtering in polarimetric SAR data based on the subspace decomposition , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[26]  Laurent Ferro-Famil,et al.  Scattering-model-based speckle filtering of polarimetric SAR data , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[27]  Gabriel Vasile,et al.  Intensity-driven adaptive-neighborhood technique for polarimetric and interferometric SAR parameters estimation , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[28]  José M. Bioucas-Dias,et al.  Multiplicative Noise Removal Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

[29]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[30]  Ridha Touzi,et al.  The principle of speckle filtering in polarimetric SAR imagery , 1994, IEEE Trans. Geosci. Remote. Sens..

[31]  Junfeng Yang,et al.  A Fast Algorithm for Edge-Preserving Variational Multichannel Image Restoration , 2009, SIAM J. Imaging Sci..

[32]  Jong-Sen Lee,et al.  Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery , 1994, IEEE Trans. Geosci. Remote. Sens..

[33]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[34]  Afzal Godil,et al.  Investigating the Bag-of-Words Method for 3D Shape Retrieval , 2010, EURASIP J. Adv. Signal Process..

[35]  Yihong Gong,et al.  Deep Learning with Kernel Regularization for Visual Recognition , 2008, NIPS.

[36]  Xiaoxia Huang,et al.  A Level Set Filter for Speckle Reduction in SAR Images , 2010, EURASIP J. Adv. Signal Process..

[37]  B. A. D. H. Brandwood.  A complex gradient operator and its applica-tion in adaptive array theory , 1983 .

[38]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[39]  E. Pottier,et al.  Polarimetric Radar Imaging: From Basics to Applications , 2009 .

[40]  Shiqian Ma,et al.  Sparse Inverse Covariance Selection via Alternating Linearization Methods , 2010, NIPS.

[41]  Tony F. Chan,et al.  Color TV: total variation methods for restoration of vector-valued images , 1998, IEEE Trans. Image Process..

[42]  Ron Kwok,et al.  Classification of multi-look polarimetric SAR imagery based on complex Wishart distribution , 1994 .

[43]  Jubo Zhu,et al.  Adaptive vectorial total variation models for multi-channel synthetic aperture radar images despeckling with fast algorithms , 2013 .

[44]  Calyampudi R. Rao,et al.  Linear Statistical Inference and Its Applications. , 1975 .

[45]  Nian Cai,et al.  Anisotropic diffusion for speckle filtering of polarimetric synthetic aperture radar imagery , 2013, J. Electronic Imaging.

[46]  T. Chan,et al.  Convergence of the alternating minimization algorithm for blind deconvolution , 2000 .

[47]  Hyenkyun Woo,et al.  Alternating Minimization Algorithm for Speckle Reduction With a Shifting Technique , 2012, IEEE Transactions on Image Processing.