Hyperparameter estimation for satellite image restoration using a MCMC maximum-likelihood method. (English) Zbl 0993.68136

Summary: The satellite image deconvolution problem is ill-posed and must be regularized. Herein, we use an edge-preserving regularization model using a \(\varphi\) function, involving two hyperparameters. Our goal is to estimate the optimal parameters in order to automatically reconstruct images. We propose to use the maximum-likelihood estimator, applied to the observed image. We need sampling from prior and posterior distributions. Since the convolution prevents use of standard samplers, we have developed a modified Geman-Yang algorithm, using an auxiliary variable and a cosine transform. We present a Markov chain Monte Carlo maximum-likelihood technique which is able to simultaneously achieve the estimation and the reconstruction.


68U10 Computing methodologies for image processing
68T10 Pattern recognition, speech recognition
Full Text: DOI


[1] Hadamard, J., Lectures on Cauchy’s Problem in Linear Partial Differential Equations (1923), Yale University Press: Yale University Press New Haven · JFM 49.0725.04
[2] Chalmond, B., Image restoration using an estimated Markov model, Signal Process., 15, 115-129 (1988)
[3] Demoment, G., Image reconstruction and restoration: overview of common estimation structures and problems, IEEE Trans. ASSP, 37, 12, 2024-2036 (1989)
[4] Geman, S.; Geman, D., Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images, IEEE Trans. PAMI, 6, 6, 721-741 (1984) · Zbl 0573.62030
[5] Tikhonov, A. N., Regularization of incorrectly posed problems, Sov. Math. Dokl., 4, 1624-1627 (1963) · Zbl 0183.11601
[6] Charbonnier, P.; Blanc-Féraud, L.; Aubert, G.; Barlaud, M., Deterministic edge-preserving regularization in computed imaging, IEEE Trans. IP, 6, 2, 298-311 (1997)
[7] Geman, D.; Reynolds, G., Constrained restoration and recovery of discontinuities, IEEE Trans. PAMI, 16, 3, 367-383 (1992)
[8] Aubert, G.; Vese, L., A variational method in image recovery, SIAM J. Numer. Anal., 34, 5, 1948-1979 (1997) · Zbl 0890.35033
[9] van Laarhoven, P. J.M.; Aarts, E. H.L., Simulated Annealing: Theory and Applications (1987), D. Reidel: D. Reidel Dordrecht · Zbl 0643.65028
[10] Higdon, D. M.; Bowsher, J. E.; Johnson, V. E.; Turkington, T. G.; Gilland, D. R.; Jaszczak, R. J., Fully Bayesian estimation of Gibbs hyperparameters for emission computed tomography data, IEEE Trans. MI, 16, 5, 516-526 (1997)
[11] Descombes, X.; Morris, R.; Zerubia, J.; Berthod, M., Estimation of Markov random field prior parameters using Markov chain Monte Carlo maximum-likelihood, IEEE Trans. IP, 8, 7, 954-963 (1999) · Zbl 0959.62082
[13] Lakshmanan, S.; Derin, H., Simultaneous parameter estimation and segmentation of Gibbs random fields using simulated annealing, IEEE Trans. PAMI, 11, 8, 799-813 (1989)
[14] Champagnat, F.; Goussard, Y.; Idier, J., Unsupervised deconvolution of sparse spike train using stochastic approximation, IEEE Trans. IP, 44, 12, 2988-2998 (1996)
[16] Younes, L., Estimation and annealing for Gibbsian fields, Ann. Inst. Poincaré, 24, 2, 269-294 (1988) · Zbl 0651.62091
[17] Younes, L., Parametric inference for imperfectly observed Gibbsian fields, Probab. Theory Related Fields, 82, 625-645 (1989) · Zbl 0659.62115
[18] Zhou, Z.; Leahy, R. M.; Qi, J., Approximate maximum-likelihood estimation for Gibbs priors, IEEE Trans. IP, 6, 6, 844-861 (1997)
[19] Saquib, S. S.; Bouman, C. A.; Sauer, K., ML parameter estimation for Markov random fields, with applications to Bayesian tomography, IEEE Trans. IP, 7, 7, 1029-1044 (1998)
[21] Metropolis, N.; Rosenbluth, A.; Rosenbluth, M.; Teller, A.; Teller, E., Equation of state calculations by fast computing machines, J. Chem. Phys., 21, 1087-1092 (1953) · Zbl 1431.65006
[22] Hastings, W., Monte Carlo sampling methods using Markov chains and their applications, Biometrika, 57, 97-109 (1970) · Zbl 0219.65008
[23] Geman, D.; Yang, C., Nonlinear image recovery with half-quadratic regularization and FFTs, IEEE Trans. IP, 4, 7, 932-946 (1995)
[24] Jain, A. K., A sinusoidal family of unitary transforms, IEEE Trans. PAMI, 1, 4, 356-365 (1979) · Zbl 0437.65110
[25] O Ruanaidh, J. J.K.; Fitzgerald, W. J., Numerical Bayesian methods applied to signal processing. Numerical Bayesian methods applied to signal processing, Statistics and Computing (1996), Springer: Springer Berlin · Zbl 0871.62025
[27] Tanner, M. A., Tools for Statistical Inference. Tools for Statistical Inference, Springer Series in Statistics (1996), Springer: Springer Berlin · Zbl 0846.62001
[28] Helstrom, C. W., Image restoration by the method of least squares, J. Opt. Soc. Amer., 57, 3, 297-303 (1967)
[30] Besag, J., Spatial interaction and statistical analysis of lattice systems, A. Roy. Stat. Soc. Ser. B, 36, 721-741 (1974)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.