Unmatched Projector/Backprojector Pairs: Perturbation and Convergence Analysis

In tomographic reconstruction problems it is not uncommon that there are errors in the implementation of the forward projector and/or the backprojector, and hence we encounter a so-called unmatched projektor/backprojector pair. Consequently, the matrices that represent the two projectors are not each other's transpose. Surprisingly, the influence of such errors in algebraic iterative reconstruction methods has received little attention in the literature. The goal of this paper is to perform a rigorous first-order perturbation analysis of the minimization problems underlying the algebraic methods in order to understand the role played by the nonmatch of the matrices. We also study the convergence properties of linear stationary iterations based on unmatched matrix pairs, leading to insight into the behavior of some important row- and column-oriented algebraic iterative methods. We conclude with numerical examples that illustrate the perturbation and convergence results.

[1]  Jan Sijbers,et al.  Fast and flexible X-ray tomography using the ASTRA toolbox. , 2016, Optics express.

[2]  Per Christian Hansen,et al.  Multicore Performance of Block Algebraic Iterative Reconstruction Methods , 2014, SIAM J. Sci. Comput..

[3]  Per Christian Hansen,et al.  Semi-convergence properties of Kaczmarz’s method , 2014 .

[4]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[5]  Stefan Kindermann,et al.  CONVERGENCE RATES FOR KACZMARZ-TYPE REGULARIZATION METHODS , 2014 .

[6]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[7]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[8]  Anne Greenbaum,et al.  Iterative methods for solving linear systems , 1997, Frontiers in applied mathematics.

[9]  Per Christian Hansen,et al.  AIR Tools - A MATLAB package of algebraic iterative reconstruction methods , 2012, J. Comput. Appl. Math..

[10]  A. Lent,et al.  Iterative algorithms for large partitioned linear systems, with applications to image reconstruction , 1981 .

[11]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[12]  P. Joseph An Improved Algorithm for Reprojecting Rays through Pixel Images , 1983, IEEE Transactions on Medical Imaging.

[13]  Klaus Mueller,et al.  A comparative study of popular interpolation and integration methods for use in computed tomography , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[14]  Heinz H. Bauschke,et al.  Optimal Rates of Linear Convergence of Relaxed Alternating Projections and Generalized Douglas-Rachford Methods for Two Subspaces , 2015, Numerical Algorithms.

[15]  Gumersindo Verdú,et al.  System Matrix Analysis for Computed Tomography Imaging , 2015, PloS one.

[16]  Lea Fleischer,et al.  Regularization of Inverse Problems , 1996 .

[17]  T. Elfving,et al.  Properties of a class of block-iterative methods , 2009 .

[18]  Per Christian Hansen,et al.  Convergence analysis for column-action methods in image reconstruction , 2016, Numerical Algorithms.

[19]  Wen Zhang,et al.  Convergence of General Nonstationary Iterative Methods for Solving Singular Linear Equations , 2011, SIAM J. Matrix Anal. Appl..

[20]  Yair Censor,et al.  On Diagonally Relaxed Orthogonal Projection Methods , 2007, SIAM J. Sci. Comput..

[21]  K J Batenburg,et al.  Performance improvements for iterative electron tomography reconstruction using graphics processing units (GPUs). , 2011, Journal of structural biology.

[22]  Yair Censor,et al.  BICAV: An Inherently Parallel Algorithm for Sparse Systems with Pixel-Dependent Weighting , 2001, IEEE Trans. Medical Imaging.

[23]  J. Navarro-Pedreño Numerical Methods for Least Squares Problems , 1996 .

[24]  Tommy Elfving,et al.  A Class of Iterative Methods: Semi-Convergence, Stopping Rules, Inconsistency, and Constraining , 2010 .