Stochastic relaxation for MAP restoration of gray level images with multiplicative noise

This paper is concerned with developing alagorithms for maximum a posteriori (MAP) restoration of gray level images degraded by multiplicative noise. The MAP algorithm requires the probability density function of the original undegraded image which is rarely available and the probability density function of the corrupting noise. By assuming that the original image is represented by a 2-D noncausal Gaussian Markov random field (GMRF) model, the MAP algorithm is written in terms of GMRF model parameters. The computer implementation of the MAP estimator equations is realized by a stochastic relaxation (SR) algorithm. The SR algorithm generates a sequence of images which converges in probability to the global MAP estimate. Several examples of restoration of the gray level image degraded by multiplicative noise are included.

[1]  Leonard M. Silverman,et al.  Nonlinear Restoration of Noisy Images , 1982, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[2]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[3]  Bobby R. Hunt,et al.  Bayesian Methods in Nonlinear Digital Image Restoration , 1977, IEEE Transactions on Computers.

[4]  J. Besag Errors‐In‐Variables Estimation for Gaussian Lattice Schemes , 1977 .

[5]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[6]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[7]  Rama Chellappa,et al.  A nonrecursive filter for edge preserving image restoration , 1985, ICASSP '85. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[8]  John W. Woods,et al.  Two-dimensional discrete Markovian fields , 1972, IEEE Trans. Inf. Theory.

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