Statistical Modeling of Low SNR Magnetic Resonance Images in Wavelet Domain Using Laplacian Prior and Two-Sided Rayleigh Noise for Visual Quality Improvement

Statistical Modeling of Low SNR Magnetic Resonance Images in Wavelet Domain Using Laplacian Prior and Two-Sided Rayleigh Noise for Visual Quality Improvement In this paper we introduce a new wavelet-based image denoising algorithm using maximum a posteriori (MAP) criterion. For this reason we propose Laplace distribution with local variance for clean image and two-sided Rayleigh model for noise in wavelet domain. The local Laplace probability density function (pdf) is able to simultaneously model the heavy-tailed nature of marginal distribution and intrascale dependency between spatial adjacent coefficients. Using local Laplace prior and two-sided Rayleigh noise, we derive a new shrinkage function for image denoising in the wavelet domain. We propose our new spatially adaptive wavelet-based image denoising algorithm for several low signal-to-noise ratio (SNR) magnetic resonance (MR) images and compare the results with other methods. The simulation results show that this algorithm is able to truly improve the visual quality of noisy MR images with very low computational cost. In case the input MR image is blurred, a blind deconvolution (BD) algorithm is necessary for visual quality improvement. Since BD techniques are usually sensitive to noise, in this paper we also apply a BD algorithm to an appropriate subband in the wavelet domain to eliminate the effect of noise in the BD procedure and to further improve visual quality.

[1]  Haluk Derin,et al.  Modeling and Segmentation of Noisy and Textured Images Using Gibbs Random Fields , 1987, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[2]  Richard Baraniuk,et al.  The Dual-tree Complex Wavelet Transform , 2007 .

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

[4]  Hayit Greenspan,et al.  MRI inter-slice reconstruction using super-resolution , 2002 .

[5]  Hossein Rabbani,et al.  Statistical modeling of low SNR magnetic resonance images in wavelet domain using Laplacian prior and two-sided Rayleigh noise for visual quality improvement , 2008 .

[6]  Jan Flusser,et al.  Multichannel blind deconvolution of spatially misaligned images , 2005, IEEE Transactions on Image Processing.

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

[8]  Saeed Gazor,et al.  Image Denoising Based on a Mixture of Laplace Distributions with Local Parameters in Complex Wavelet Domain , 2006, 2006 International Conference on Image Processing.

[9]  Aggelos K. Katsaggelos,et al.  Parameter Estimation in TV Image Restoration Using Variational Distribution Approximation , 2008, IEEE Transactions on Image Processing.

[10]  T. M. Cannon,et al.  Blind deconvolution through digital signal processing , 1975, Proceedings of the IEEE.

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

[12]  D. A. Fish,et al.  Blind deconvolution by means of the Richardson-Lucy algorithm. , 1995 .

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

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

[15]  M. Cannon.  Blind deconvolution of spatially invariant image blurs with phase , 1976 .

[16]  Reza Nezafat,et al.  Wavelet-Domain Medical Image Denoising Using Bivariate Laplacian Mixture Model , 2009, IEEE Transactions on Biomedical Engineering.

[17]  Nick G. Kingsbury,et al.  A dual-tree complex wavelet transform with improved orthogonality and symmetry properties , 2000, Proceedings 2000 International Conference on Image Processing (Cat. No.00CH37101).

[18]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[19]  Robert D. Nowak,et al.  Wavelet-based Rician noise removal for magnetic resonance imaging , 1999, IEEE Trans. Image Process..

[20]  Pierre Moulin,et al.  Analysis of Multiresolution Image Denoising Schemes Using Generalized Gaussian and Complexity Priors , 1999, IEEE Trans. Inf. Theory.

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

[22]  Hossein Rabbani,et al.  Image/video denoising based on a mixture of Laplace distributions with local parameters in multidimensional complex wavelet domain , 2008, Signal Process..

[23]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[24]  Aggelos K. Katsaggelos,et al.  A Bayesian approach to blind deconvolution based on Dirichlet distributions , 1997, 1997 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[25]  Frédo Durand,et al.  Understanding and evaluating blind deconvolution algorithms , 2009, CVPR.

[26]  L. Lucy.  An iterative technique for the rectification of observed distributions , 1974 .