Paper
2 February 2006 Preconditioned conjugate gradient without linesearch: a comparison with the half-quadratic approach for edge-preserving image restoration
Christian Labat, Jérôme Idier
Author Affiliations +
Proceedings Volume 6065, Computational Imaging IV; 60650I (2006) https://doi.org/10.1117/12.641663
Event: Electronic Imaging 2006, 2006, San Jose, California, United States
Abstract
Our contribution deals with image restoration. The adopted approach consists in minimizing a penalized least squares (PLS) criterion. Here, we are interested in the search of efficient algorithms to carry out such a task. The minimization of PLS criteria can be addressed using a half-quadratic approach (HQ). However, the nontrivial inversion of a linear system is needed at each iteration. In practice, it is often proposed to approximate this inversion using a truncated preconditioned conjugate gradient (PCG) method. However, we point out that theoretical convergence is not proved for such approximate HQ algorithms, referred here as HQ+PCG. In the proposed contribution, we rely on a different scheme, also based on PCG and HQ ingredients and referred as PCG+HQ1D. General linesearch methods ensuring convergence of PCG type algorithms are difficult to code and to tune. Therefore, we propose to replace the linesearch step by a truncated scalar HQ algorithm. Convergence is established for any finite number of HQ1D sub-iterations. Compared to the HQ+PCG approach, we show that our scheme is preferable on both the theoretical and practical grounds.
© (2006) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Christian Labat and Jérôme Idier "Preconditioned conjugate gradient without linesearch: a comparison with the half-quadratic approach for edge-preserving image restoration", Proc. SPIE 6065, Computational Imaging IV, 60650I (2 February 2006); https://doi.org/10.1117/12.641663
Lens.org Logo
CITATIONS
Cited by 2 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Image restoration

Signal to noise ratio

Algorithms

Matrices

Denoising

Fourier transforms

MATLAB

Back to Top