A Multilevel Iteration Method for Solving a Coupled Integral Equation Model in Image Restoration

The problem of out-of-focus image restoration can be modeled as an ill-posed integral equation, which can be regularized as a second kind of equation using the Tikhonov method. The multiscale collocation method with the compression strategy has already been developed to discretize this well-posed equation. However, the integral computation and solution of the large multiscale collocation integral equation are two time-consuming processes. To overcome these difficulties, we propose a fully discrete multiscale collocation method using an integral approximation strategy to compute the integral, which efficiently converts the integral operation to the matrix operation and reduces costs. In addition, we also propose a multilevel iteration method (MIM) to solve the fully discrete integral equation obtained from the integral approximation strategy. Herein, the stopping criterion and the computation complexity that correspond to the MIM are shown. Furthermore, a posteriori parameter choice strategy is developed for this method, and the final convergence order is evaluated. We present three numerical experiments to display the performance and computation efficiency of our proposed methods.

[1]  Charles W. Groetsch,et al.  Convergence analysis of a regularized degenerate kernel method for Fredholm integral equations of the first kind , 1990 .

[2]  Raymond H. Chan,et al.  Wavelet Algorithms for High-Resolution Image Reconstruction , 2002, SIAM J. Sci. Comput..

[3]  Weifu Fang,et al.  A fast collocation method for an inverse boundary value problem , 2004 .

[4]  Charles A. Micchelli,et al.  Wavelet Galerkin methods for second-kind integral equations , 1997 .

[5]  Yuesheng Xu,et al.  Fast Collocation Methods for Second Kind Integral Equations , 2002, SIAM J. Numer. Anal..

[6]  L. Shampine,et al.  A collocation method for boundary value problems , 1972 .

[7]  C. Groetsch Uniform convergence of regularization methods for Fredholm equations of the first kind , 1985, Journal of the Australian Mathematical Society. Series A. Pure Mathematics and Statistics.

[8]  Yuesheng Xu,et al.  Reconstruction and Decomposition Algorithms for Biorthogonal Multiwavelets , 1997, Multidimens. Syst. Signal Process..

[9]  Yao Lu,et al.  Integral equation models for image restoration: high accuracy methods and fast algorithms , 2010 .

[10]  Sergei V. Pereverzyev,et al.  On the Adaptive Selection of the Parameter in Regularization of Ill-Posed Problems , 2005, SIAM J. Numer. Anal..

[11]  Yuesheng Xu,et al.  A construction of interpolating wavelets on invariant sets , 1999, Math. Comput..

[12]  Fanchun Li,et al.  Fast multi-level iteration methods with compression technique for solving ill-posed integral equations , 2014, J. Comput. Appl. Math..

[13]  Yuesheng Xu,et al.  A multilevel augmentation method for solving ill-posed operator equations , 2006 .

[14]  M. P. Rajan,et al.  Convergence analysis of a regularized approximation for solving Fredholm integral equations of the first kind , 2003 .

[15]  Yuesheng Xu,et al.  Multiscale collocation methods for ill-posed integral equations via a coupled system , 2012 .

[16]  Yuesheng Xu,et al.  A Higher-Order Polynomial Method for SPECT Reconstruction , 2019, IEEE Transactions on Medical Imaging.

[17]  Zhongying Chen,et al.  Fast collocation methods for solving ill-posed integral equations of the first kind , 2008 .

[18]  Q. Jin,et al.  Analysis of the iteratively regularized Gauss–Newton method under a heuristic rule , 2018 .

[19]  F. Bauer,et al.  The quasi-optimality criterion for classical inverse problems , 2008 .

[20]  Yuesheng Xu,et al.  Computing Integrals Involved the Gaussian Function with a Small Standard Deviation , 2018, J. Sci. Comput..

[21]  Toomas Raus,et al.  About the Balancing Principle for Choice of the Regularization Parameter , 2009 .