FFT-based preconditioners for Toeplitz-block least squares problems
SIAM Journal on Numerical Analysis
Conjugate Gradient Methods for Toeplitz Systems
SIAM Review
Digital Image Restoration
Improved Methods of Maximum a Posteriori Restoration
IEEE Transactions on Computers
Non-negatively constrained image deblurring with an inexact interior point method
Journal of Computational and Applied Mathematics
Hi-index | 0.00 |
Nonlinear image restoration finds applications in a wide variety of research areas. In this paper, we consider nonlinear space-invariant imaging system with additive noise. The restored images can be found by solving weighted Toeplitz least squares problems. Since the normal equations matrices are non-Toeplitz in general, the fast Fourier transforms (FFTs) cannot be utilized in the evaluation of their inverses. We employ the preconditioned conjugate gradient method (PCG) with the FFT-based preconditioners to solve regularized linear systems arising from nonlinear image restoration problems. Thus we precondition these linear systems in the Fourier domain, while iterating in the spatial domain. Numerical examples are reported on a ground-based atmospheric imaging problem to demonstrate the fast convergence of the FFT-based PCG method.