Title Fast image reconstruction algorithms combining half-quadraticregularization and preconditioning

In this paper, we focus on image deconvolution and image reconstruction problems where a sought image is recovered from degraded observed data. The solution is defined to be the minimizer of an objective function combining a data-fidelity term and a edge-preserving, convex regularization term. Our objective is to speed up the calculation of the solution in a wide range of situations. To this end, we propose a method applying pertinent preconditioning to an adapted half-quadratic equivalent form of the objective function. The optimal solution is then found using an alternating minimization (AM) scheme. We focus specifically on Huber regularization. We exhibit the possibility get very fast calculations while preserving the edges in the solution. Preliminary numerical results are reported to illustrate the effectiveness of our method.

[1]  Rachid Deriche,et al.  Image Sequence Analysis via Partial Differential Equations , 1999, Journal of Mathematical Imaging and Vision.

[2]  Laurent D. Cohen,et al.  Auxiliary variables and two-step iterative algorithms in computer vision problems , 2004, Journal of Mathematical Imaging and Vision.

[3]  Mila Nikolova,et al.  Local Strong Homogeneity of a Regularized Estimator , 2000, SIAM J. Appl. Math..

[4]  Raymond H. Chan,et al.  A Fast Algorithm for Deblurring Models with Neumann Boundary Conditions , 1999, SIAM J. Sci. Comput..

[5]  Yoram Bresler,et al.  Globally convergent edge-preserving regularized reconstruction: an application to limited-angle tomography , 1998, IEEE Trans. Image Process..

[6]  Nirmal K. Bose,et al.  High‐resolution image reconstruction with multisensors , 1998, Int. J. Imaging Syst. Technol..

[7]  Michel Barlaud,et al.  Deterministic edge-preserving regularization in computed imaging , 1997, IEEE Trans. Image Process..

[8]  Raymond H. Chan,et al.  Conjugate Gradient Methods for Toeplitz Systems , 1996, SIAM Rev..

[9]  Fadil Santosa,et al.  Recovery of Blocky Images from Noisy and Blurred Data , 1996, SIAM J. Appl. Math..

[10]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[11]  Michael J. Black,et al.  On the unification of line processes , 1996 .

[12]  Donald Geman,et al.  Nonlinear image recovery with half-quadratic regularization , 1995, IEEE Trans. Image Process..

[13]  Ken D. Sauer,et al.  A generalized Gaussian image model for edge-preserving MAP estimation , 1993, IEEE Trans. Image Process..

[14]  Donald Geman,et al.  Constrained Restoration and the Recovery of Discontinuities , 1992, IEEE Trans. Pattern Anal. Mach. Intell..