SAR speckle reduction using wavelet denoising and Markov random field modeling

The granular appearance of speckle noise in synthetic aperture radar (SAR) imagery makes it very difficult to visually and automatically interpret SAR data. Therefore, speckle reduction is a prerequisite for many SAR image processing tasks. In this paper, we develop a speckle reduction algorithm by fusing the wavelet Bayesian denoising technique with Markov-random-field-based image regularization. Wavelet coefficients are modeled independently and identically by a two-state Gaussian mixture model, while their spatial dependence is characterized by a Markov random field imposed on the hidden state of Gaussian mixtures. The Expectation-Maximization algorithm is used to estimate hyperparameters and specify the mixture model, and the iterated-conditional-modes method is implemented to optimize the state configuration. The noise-free wavelet coefficients are finally estimated by a shrinkage function based on local weighted averaging of the Bayesian estimator. Experimental results show that the proposed method outperforms standard wavelet denoising techniques in terms of the signal-to-noise ratio and the equivalent-number-of-looks measures in most cases. It also achieves better performance than the refined Lee filter.

[1]  M. Sties,et al.  Efficient speckle filtering of SAR images , 1994, Proceedings of IGARSS '94 - 1994 IEEE International Geoscience and Remote Sensing Symposium.

[2]  Fawwaz T. Ulaby,et al.  Statistical properties of logarithmically transformed speckle , 2002, IEEE Trans. Geosci. Remote. Sens..

[3]  F. Ulaby,et al.  Handbook of radar scattering statistics for terrain , 1989 .

[4]  E. Nezry,et al.  Adaptive speckle filters and scene heterogeneity , 1990 .

[5]  Dirk H. Hoekman Speckle ensemble statistics of logarithmically scaled data [radar] , 1991, IEEE Trans. Geosci. Remote. Sens..

[6]  G. Beylkin On the representation of operators in bases of compactly supported wavelets , 1992 .

[7]  H. Cramér Mathematical methods of statistics , 1947 .

[8]  D. Donoho,et al.  Translation-Invariant De-Noising , 1995 .

[9]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[10]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[11]  Corina da Costa Freitas,et al.  A model for extremely heterogeneous clutter , 1997, IEEE Trans. Geosci. Remote. Sens..

[12]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[13]  C. McDiarmid SIMULATED ANNEALING AND BOLTZMANN MACHINES A Stochastic Approach to Combinatorial Optimization and Neural Computing , 1991 .

[14]  Ramesh A. Gopinath,et al.  Wavelet based speckle reduction with application to SAR based ATD/R , 1994, Proceedings of 1st International Conference on Image Processing.

[15]  I. Daubechies Ten Lectures on Wavelets , 1992 .

[16]  Alexandre Jouan,et al.  Speckle filtering of SAR images: a comparative study between complex-wavelet-based and standard filters , 1997, Optics & Photonics.

[17]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .

[18]  Anil K. Jain,et al.  A Markov random field model for classification of multisource satellite imagery , 1996, IEEE Trans. Geosci. Remote. Sens..

[19]  Martin Vetterli,et al.  Spatially adaptive wavelet thresholding with context modeling for image denoising , 1998, Proceedings 1998 International Conference on Image Processing. ICIP98 (Cat. No.98CB36269).

[20]  Victor S. Frost,et al.  A Model for Radar Images and Its Application to Adaptive Digital Filtering of Multiplicative Noise , 1982, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[21]  J. Goodman Statistical Properties of Laser Speckle Patterns , 1963 .

[22]  Stan Z. Li,et al.  Markov Random Field Modeling in Computer Vision , 1995, Computer Science Workbench.

[23]  Henri Maitre,et al.  Smoothing speckled synthetic aperture radar images by using maximum homgeneous region filters , 1992 .

[24]  Kohei Arai,et al.  A method of speckle noise reduction for SAR data , 1996 .

[25]  Pierre Moulin,et al.  Information-theoretic analysis of interscale and intrascale dependencies between image wavelet coefficients , 2001, IEEE Trans. Image Process..

[26]  Martin Vetterli,et al.  Spatial adaptive wavelet thresholding for image denoising , 1997, Proceedings of International Conference on Image Processing.

[27]  Alexander A. Sawchuk,et al.  Adaptive Noise Smoothing Filter for Images with Signal-Dependent Noise , 1985, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[28]  H. Arsenault,et al.  Properties of speckle integrated with a finite aperture and logarithmically transformed , 1976 .

[29]  Edward H. Adelson,et al.  Noise removal via Bayesian wavelet coring , 1996, Proceedings of 3rd IEEE International Conference on Image Processing.

[30]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[31]  Jong-Sen Lee,et al.  Refined filtering of image noise using local statistics , 1981 .

[32]  Eric L. Miller,et al.  Wavelet domain image restoration with adaptive edge-preserving regularization , 2000, IEEE Trans. Image Process..

[33]  Jong-Sen Lee Speckle suppression and analysis for synthetic aperture radar images , 1986 .

[34]  H. Chipman,et al.  Adaptive Bayesian Wavelet Shrinkage , 1997 .

[35]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[36]  Jia Jie Bayesian denoising of visual images in the wavelet domain , 2003 .

[37]  B. Silverman,et al.  The Stationary Wavelet Transform and some Statistical Applications , 1995 .

[38]  Kannan Ramchandran,et al.  Low-complexity image denoising based on statistical modeling of wavelet coefficients , 1999, IEEE Signal Processing Letters.

[39]  R. Redner,et al.  Mixture densities, maximum likelihood, and the EM algorithm , 1984 .

[40]  S. Quegan,et al.  Understanding Synthetic Aperture Radar Images , 1998 .

[41]  G. Franceschetti,et al.  Iterative homomorphic technique for speckle reduction in synthetic-aperture radar imaging , 1995 .

[42]  Martin Vetterli,et al.  Adaptive wavelet thresholding for image denoising and compression , 2000, IEEE Trans. Image Process..

[43]  Anil K. Jain,et al.  Random field models in image analysis , 1989 .

[44]  D. L. Donoho,et al.  Ideal spacial adaptation via wavelet shrinkage , 1994 .

[45]  D. Kaplan,et al.  On the statistical characteristics of log-compressed Rayleigh signals: theoretical formulation and experimental results , 1993 .

[46]  Dirk Roose,et al.  Wavelet-based image denoising using a Markov random field a priori model , 1997, IEEE Trans. Image Process..

[47]  Andrew F. Laine,et al.  Speckle reduction and contrast enhancement of echocardiograms via multiscale nonlinear processing , 1998, IEEE Transactions on Medical Imaging.

[48]  Robert D. Nowak,et al.  Wavelet-based statistical signal processing using hidden Markov models , 1998, IEEE Trans. Signal Process..

[49]  J. Benediktsson,et al.  Speckle reduction of SAR images using wavelet-domain hidden Markov models , 2000, IGARSS 2000. IEEE 2000 International Geoscience and Remote Sensing Symposium. Taking the Pulse of the Planet: The Role of Remote Sensing in Managing the Environment. Proceedings (Cat. No.00CH37120).