Estimation of regularization parameters in multiple-image deblurring

¹ Chip Computers Consulting s.r.l., Viale Don L. Sturzo 82, S. Liberale di Marcon, 30020 Venice, Italy ESA-VILSPA, Apartado 50727, 28080 Madrid, Spain e-mail: robertovio@tin.it
² Department of Statistics, Harvard University, Cambridge, MA 02138, USA e-mail: pingma@stat.harvard.edu
³ Department of Statistics, Purdue University, West Lafayette, IN 47907, USA e-mail: wenxuan@stat.purdue.edu
⁴ Department of Mathematics and Computer Science, Emory University, Atlanta, GA 30322, USA e-mail: nagy@mathcs.emory.edu
⁵ Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden CO 80401, USA e-mail: ltenorio@Mines.EDU
⁶ ESA-VILSPA, Apartado 50727, 28080 Madrid, Spain e-mail: willem.wamsteker@esa.int

Received: 21 January 2004
Accepted: 25 May 2004

Abstract

We consider the estimation of the regularization parameter for the simultaneous deblurring of multiple noisy images via Tikhonov regularization. We approach the problem in three ways. We first reduce the problem to a single-image deblurring for which the regularization parameter can be estimated through a classic generalized cross-validation () method. A modification of this function is used for correcting the undersmoothing typical of the original technique. With a second method, we minimize an average least-squares fit to the images and define a new  function. In the last approach, we use the classical  on a single higher-dimensional image obtained by concatenating all the images into a single vector. With a reliable estimator of the regularization parameter, one can fully exploit the excellent computational characteristics typical of direct deblurring methods, which, especially for large images, makes them competitive with the more flexible but much slower iterative algorithms. The performance of the techniques is analyzed through numerical experiments. We find that under the independent homoscedastic and Gaussian assumptions made on the noise, the three approaches provide almost identical results with the first single image providing the practical advantage that no new software is required and the same image can be used with other deblurring algorithms.