GMRES methods for tomographic reconstruction with an unmatched back projector

DOI10.1016/J.CAM.2022.114352zbMATH Open1489.65060arXiv2110.01481OpenAlexW3203301659MaRDI QIDQ2146331FDOQ2146331

Per Christian Hansen, Keiichi Morikuni, Ken Hayami

Publication date: 16 June 2022

Published in: Journal of Computational and Applied Mathematics (Search for Journal in Brave)

Abstract: Unmatched pairs of forward and back projectors are common in X-ray CT computations for large-scale problems; they are caused by the need for fast algorithms that best utilize the computer hardware, and it is an interesting and challenging task to develop fast and easy-to-use algorithms for these cases. Our approach is to use preconditioned GMRES, in the form of the AB- and BA-GMRES algorithms, to handle the unmatched normal equations associated with an unmatched pair. These algorithms are simple to implement, they rely only on computations with the available forward and back projectors, and they do not require the tuning of any algorithm parameters. We show that these algorithms are equivalent to well-known LSQR and LSMR algorithms in the case of a matched projector. Our numerical experiments demonstrate that AB- and BA-GMRES exhibit a desired semi-convergence behavior that is comparable with LSQR/LSMR and that standard stopping rules work well. Hence, AB- and BA-GMRES are suited for large-scale CT reconstruction problems with noisy data and unmatched projector pairs.

Full work available at URL: https://arxiv.org/abs/2110.01481

zbMATH Keywords

GMRESsemi-convergenceregularizing iterationsstopping rulesCT-reconstructionunmatched backprojector

Mathematics Subject Classification ID

Ill-posedness and regularization problems in numerical linear algebra (65F22) Image processing (compression, reconstruction, etc.) in information and communication theory (94A08) Iterative numerical methods for linear systems (65F10)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Regularization tools version \(4.0\) for matlab \(7.3\)
• LSMR: An Iterative Algorithm for Sparse Least-Squares Problems
• Perturbation bounds in connection with singular value decomposition
• The mathematics of computerized tomography
• Improved image deblurring with anti-reflective boundary conditions and re-blurring
• Old and new parameter choice rules for discrete ill-posed problems
• GMRES Methods for Least Squares Problems
• On the regularizing properties of the GMRES method
• Discrete Inverse Problems
• Iterative regularization with minimum-residual methods
• The Moore-Penrose inverse of a partitioned matrix \(M=\begin{pmatrix} D&D\\ B&C\end{pmatrix}\)
• Arnoldi methods for image deblurring with anti-reflective boundary conditions
• Inheritance of the discrete Picard condition in Krylov subspace methods
• Convergence of Inner-Iteration GMRES Methods for Rank-Deficient Least Squares Problems
• Multicore Performance of Block Algebraic Iterative Reconstruction Methods
• Semi-convergence properties of Kaczmarz’s method
• A Twin Error Gauge for Kaczmarz's Iterations
• AIR tools II: algebraic iterative reconstruction methods, improved implementation
• Regularization properties of Krylov iterative solvers CGME and LSMR for linear discrete ill-posed problems with an application to truncated randomized SVDs
• Fixing Nonconvergence of Algebraic Iterative Reconstruction with an Unmatched Backprojector
• Unmatched Projector/Backprojector Pairs: Perturbation and Convergence Analysis
• Convergence of proximal gradient algorithm in the presence of adjoint mismatch ^*
• Computed Tomography: Algorithms, Insight, and Just Enough Theory
• Kaczmarz-Type Inner-Iteration Preconditioned Flexible GMRES Methods for Consistent Linear Systems

Cited In (2)

• Convergence Results for Primal-Dual Algorithms in the Presence of Adjoint Mismatch
• Solution of mismatched monotone+Lipschitz inclusion problems

Uses Software

• Regularization tools
• UTV
• AIR tools
• LSMR
• ASTRA

This page was built for publication: GMRES methods for tomographic reconstruction with an unmatched back projector

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2146331)