Bayesian Analysis

Bayesian Multiscale Smoothing of Gaussian Noised Images

Meng Li and Subhashis Ghosal

Full-text: Open access


We propose a multiscale model for Gaussian noised images under a Bayesian framework for both 2-dimensional (2D) and 3-dimensional (3D) images. We use a Chinese restaurant process prior to randomly generate ties among intensity values at neighboring pixels in the image. The resulting Bayesian estimator enjoys some desirable asymptotic properties for identifying precise structures in the image. The proposed Bayesian denoising procedure is completely data-driven. A conditional conjugacy property allows analytical computation of the posterior distribution without involving Markov chain Monte Carlo (MCMC) methods, making the method computationally efficient. Simulations on Shepp-Logan phantom and Lena test images confirm that our smoothing method is comparable with the best available methods for light noise and outperforms them for heavier noise both visually and numerically. The proposed method is further extended for 3D images. A simulation study shows that the proposed method is numerically better than most existing denoising approaches for 3D images. A 3D Shepp-Logan phantom image is used to demonstrate the visual and numerical performance of the proposed method, along with the computational time. MATLAB toolboxes are made available online (both 2D and 3D) to implement the proposed method and reproduce the numerical results.

Article information

Bayesian Anal., Volume 9, Number 3 (2014), 733-758.

First available in Project Euclid: 5 September 2014

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Chinese Restaurant Process MCMC-free computation 3-dimensional image


Li, Meng; Ghosal, Subhashis. Bayesian Multiscale Smoothing of Gaussian Noised Images. Bayesian Anal. 9 (2014), no. 3, 733--758. doi:10.1214/14-BA871.

Export citation


  • Baddeley, A. (1992). “An error metric for binary images.” Proceedings of the IEEE Workshop on Robust Computer Vision, 59–78.
  • Castro, R., Willett, R., and Nowak, R. (2004). “Coarse-to-fine manifold learning.” In IEEE International Conference on Acoustics, Speech, and Signal Processing, 2004. (ICASSP’04).
  • Coifman, R. and Donoho, D. (1995). “Translation-Invariant De-Noising.” In Antoniadis, A. and Oppenheim, G. (eds.), Wavelets and Statistics, volume 103 of Lecture Notes in Statistics, 125–150. Springer New York.
  • Coupé, P., Yger, P., Prima, S., Hellier, P., Kervrann, C., and Barillot, C. (2008). “An optimized blockwise nonlocal means denoising filter for 3-D magnetic resonance images.” IEEE Transactions on Medical Imaging, 27(4): 425–441.
  • Donoho, D. L. (1999). “Wedgelets: Nearly Minimax Estimation of Edges.” The Annals of Statistics, 27(3): pp. 859–897.
  • Ferreira, M. A. and Lee, H. K. (2007). Multiscale modeling: a Bayesian perspective. Springer Verlag.
  • Hartigan, J. A. and Wong, M. A. (1979). “Algorithm AS 136: A k-means clustering algorithm.” Journal of the Royal Statistical Society. Series C (Applied Statistics), 28(1): 100–108.
  • Huttenlocher, D., Klanderman, G., and Rucklidge, W. (1993). “Comparing images using the Hausdorff distance.” IEEE Transactions on Pattern Analysis and Machine Intelligence, 15(9): 850–863.
  • Jain, A. (1989). Fundamentals of Digital Image Processing. Prentice-Hall Information and System Sciences Series. Prentice Hall.
  • Kolaczyk, E. D. (1999). “Bayesian Multiscale Models for Poisson Processes.” Journal of the American Statistical Association, 94: 920–933.
  • Kolaczyk, E. D. and Nowak, R. D. (2004). “Multiscale likelihood analysis and complexity penalized estimation.” Annals of Statistics, 32: 500–527.
  • Lagarias, J. C., Reeds, J. A., Wright, M. H., and Wright, P. E. (1998). “Convergence properties of the Nelder–Mead simplex method in low dimensions.” SIAM Journal on Optimization, 9(1): 112–147.
  • Lindeberg, T. (1993). Scale-Space Theory in Computer Vision. Kluwer International Series in Engineering and Computer Science: Robotics: Vision, manipulation and sensors. Springer.
  • McLachlan, G. and Peel, D. (2000). Finite Mixture Models. John Wiley & Sons.
  • Mukherjee, P. and Qiu, P. (2011). “3-D image denoising by local smoothing and nonparametric regression.” Technometrics, 53(2): 196–208.
  • Perona, P. and Malik, J. (1990). “Scale-space and edge detection using anisotropic diffusion.” IEEE Transactions on Pattern Analysis and Machine Intelligence, 12(7): 629–639.
  • Rudin, L., Osher, S., and Fatemi, E. (1992). “Nonlinear total variation based noise removal algorithms.” Physica D: Nonlinear Phenomena, 60(1): 259–268.
  • Sanyal, N. and Ferreira, M. A. (2012). “Bayesian hierarchical multi-subject multiscale analysis of functional MRI data.” NeuroImage, 63(3): 1519–1531.
  • Schabel, M. (2005). “3D Shepp-Logan phantom.” MATLAB Central File Exchange.
  • White, J. T. and Ghosal, S. (2011). “Bayesian smoothing of photon-limited images with applications in astronomy.” Journal of the Royal Statistical Society, Series B, 73: 579–599.
  • (2013). “Denoising three-Dimensional and colored images using a Bayesian multi-scale model for photon counts.” Signal Processing, 93: 2906–2914.
  • Willett, R. and Nowak, R. (2003). “Platelets: a multiscale approach for recovering edges and surfaces in photon-limited medical imaging.” IEEE Transactions on Medical Imaging, 22(3): 332–350.
  • (2004). “Fast multiresolution photon-limited image reconstruction.” In IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2004, 1192–1195.
  • Wilson, D., Baddeley, A., and Owens, R. (1997). “A new metric for grey-scale image comparison.” International Journal of Computer Vision, 24(1): 5–17.
  • Zhou, M., Chen, H., Paisley, J., Ren, L., Sapiro, G., and Carin, L. (2012). “Non-parametric Bayesian dictionary learning for sparse image representations.” IEEE Transactions on Image Processing, 21: 13—144.