Fixing nonconvergence of algebraic iterative reconstruction with an unmatched backprojector

DOI10.1137/18M1206448zbMATH Open1420.65031arXiv1902.04282WikidataQ127753302 ScholiaQ127753302MaRDI QIDQ5230620FDOQ5230620

Authors: Yiqiu Dong, Per Christian Hansen, M. E. Hochstenbach, Nicolai André Brogaard Riis 

Publication date: 28 August 2019

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: We consider algebraic iterative reconstruction methods with applications in image reconstruction. In particular, we are concerned with methods based on an unmatched projector/backprojector pair; i.e., the backprojector is not the exact adjoint or transpose of the forward projector. Such situations are common in large-scale computed tomography, and we consider the common situation where the method does not converge due to the nonsymmetry of the iteration matrix. We propose a modified algorithm that incorporates a small shift parameter, and we give the conditions that guarantee convergence of this method to a fixed point of a slightly perturbed problem. We also give perturbation bounds for this fixed point. Moreover, we discuss how to use Krylov subspace methods to efficiently estimate the leftmost eigenvalue of a certain matrix to select a proper shift parameter. The modified algorithm is illustrated with test problems from computed tomography.

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

Recommendations

• Unmatched projector/backprojector pairs: perturbation and convergence analysis
• Projected non-stationary simultaneous iterative methods
• GMRES methods for tomographic reconstruction with an unmatched back projector
• Convergence of the simultaneous algebraic reconstruction technique (SART)
• Some results on simultaneous algebraic techniques in image reconstruction from projections

zbMATH Keywords

perturbation theorycomputed tomographyalgebraic iterative reconstructionleftmost eigenvalue estimationunmatched transpose

Mathematics Subject Classification ID

Eigenvalues, singular values, and eigenvectors (15A18) Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Ill-posedness and regularization problems in numerical linear algebra (65F22) Iterative numerical methods for linear systems (65F10) Norms of matrices, numerical range, applications of functional analysis to matrix theory (15A60)

Cites Work

• Regularization tools version \(4.0\) for matlab \(7.3\)
• Implicit Application of Polynomial Filters in a k-Step Arnoldi Method
• Title not available (Why is that?)
• Title not available (Why is that?)
• The mathematics of computerized tomography
• A Jacobi–Davidson Iteration Method for Linear Eigenvalue Problems
• The discrete Picard condition for discrete ill-posed problems
• Iterative solution of large sparse systems of equations
• A Krylov--Schur algorithm for large eigenproblems
• Matrix transformations for computing rightmost eigenvalues of large sparse non-symmetric eigenvalue problems
• Semiconvergence and Relaxation Parameters for Projected SIRT Algorithms
• AIR tools II: algebraic iterative reconstruction methods, improved implementation
• Oblique projections and standard-form transformations for discrete inverse problems
• The randomized Kaczmarz method with mismatched adjoint
• Test Matrices for Regularization Methods
• Using Krylov approximations to the matrix exponential operator in Davidson's method
• Unmatched projector/backprojector pairs: perturbation and convergence analysis

Cited In (10)

• Some transpose-free CG-like solvers for nonsymmetric ill-posed problems
• Convergence Results for Primal-Dual Algorithms in the Presence of Adjoint Mismatch
• Convergence analyses based on frequency decomposition for the randomized row iterative method
• Spherical Image Inpainting with Frame Transformation and Data-Driven Prior Deep Networks
• Unmatched projector/backprojector pairs: perturbation and convergence analysis
• Convergence of proximal gradient algorithm in the presence of adjoint mismatch
• Mathematical imaging and surface processing. Abstracts from the workshop held August 21--27, 2022
• GMRES methods for tomographic reconstruction with an unmatched back projector
• Solution of mismatched monotone+Lipschitz inclusion problems
• Chambolle-Pock's primal-dual method with mismatched adjoint

Uses Software

• Regularization tools
• eigs
• Spot
• IRAM
• AIR tools
• ASTRA
• TIGRE
• GitHub

This page was built for publication: Fixing nonconvergence of algebraic iterative reconstruction with an unmatched backprojector

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