A multi-level preconditioned Krylov method for the efficient solution of algebraic tomographic reconstruction problems

DOI10.1016/J.CAM.2014.12.044MaRDI QIDQ2018482zbMATHOpenAlexDBLPWikidataFDO

Authors Siegfried Cools, Pieter Ghysels, Wim van Aarle, Jan Sijbers, Wim Vanroose

Publication date 24 March 2015

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

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

zbMATH Keywords

waveletsconvergencepreconditioningnumerical experimenttomographymultigridKrylov methodslinear iterative methodalgebraic reconstructionlinear inversion problem

Mathematics Subject Classification ID

Numerical methods for wavelets (65T60) Preconditioners for iterative methods (65F08) Radon transform (44A12) Numerical methods for integral transforms (65R10) Numerical methods for inverse problems for integral equations (65R32)

Abstract: Classical iterative methods for tomographic reconstruction include the class of Algebraic Reconstruction Techniques (ART). Convergence of these stationary linear iterative methods is however notably slow. In this paper we propose the use of Krylov solvers for tomographic linear inversion problems. These advanced iterative methods feature fast convergence at the expense of a higher computational cost per iteration, causing them to be generally uncompetitive without the inclusion of a suitable preconditioner. Combining elements from standard multigrid (MG) solvers and the theory of wavelets, a novel wavelet-based multi-level (WMG) preconditioner is introduced, which is shown to significantly speed-up Krylov convergence. The performance of the WMG-preconditioned Krylov method is analyzed through a spectral analysis, and the approach is compared to existing methods like the classical Simultaneous Iterative Reconstruction Technique (SIRT) and unpreconditioned Krylov methods on a 2D tomographic benchmark problem. Numerical experiments are promising, showing the method to be competitive with the classical Algebraic Reconstruction Techniques in terms of convergence speed and overall performance (CPU time) as well as precision of the reconstruction.

Recommendations

• Fast Approximation of Algebraic Reconstruction Methods for Tomography
• Superiorization of preconditioned conjugate gradient algorithms for tomographic image reconstruction
• A computationally efficient version of the algebraic method for computer tomography
• A multilevel domain decomposition algorithm for fast \(O(N^2\log N)\) reprojection of tomographic images
• scientific article; zbMATH DE number 418829
• Preconditioned alternating projection algorithms for maximum \textit{a posteriori} ECT reconstruction
• Preconditioning conjugate gradient with symmetric algebraic reconstruction technique (ART) in computerized tomography
• Numerical approach for a relaxed minimization problem arising in tomographic reconstruction
• scientific article; zbMATH DE number 869701
• scientific article; zbMATH DE number 624827

Cites work

• scientific article; zbMATH DE number 3938445 (Why is no real title available?)
• scientific article; zbMATH DE number 4001918 (Why is no real title available?)
• scientific article; zbMATH DE number 1561761 (Why is no real title available?)
• scientific article; zbMATH DE number 3797117 (Why is no real title available?)
• scientific article; zbMATH DE number 803211 (Why is no real title available?)
• scientific article; zbMATH DE number 936298 (Why is no real title available?)
• scientific article; zbMATH DE number 967332 (Why is no real title available?)
• A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
• A Multigrid Tutorial, Second Edition
• A general framework for nonlinear multigrid inversion
• A multigrid method enhanced by Krylov subspace iteration for discrete Helmholtz equations
• Algebraic multilevel Krylov methods
• Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer.
• Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
• DART: A Practical Reconstruction Algorithm for Discrete Tomography
• Discrete inverse problems. Insight and algorithms.
• Fourier Analysis of GMRES(m) Preconditioned by Multigrid
• Linear and nonlinear inverse problems with practical applications
• Local Fourier analysis of the complex shifted Laplacian preconditioner for Helmholtz problems.
• Multi-Level Adaptive Solutions to Boundary-Value Problems
• Multigrid solution of a linearized, regularized least-squares problem in electrical impedance tomography
• Multilevel Image Reconstruction with Natural Pixels
• Optical tomography in medical imaging
• Orthonormal bases of compactly supported wavelets
• Principles of computerized tomography imaging
• Recent computational developments in Krylov subspace methods for linear systems
• The hierarchical basis multigrid method
• The speed of convergence of one iterative process
• Tikhonov Regularization and Total Least Squares
• Two-level preconditioned Krylov subspace methods for the solution of three-dimensional heterogeneous Helmholtz problems in seismics
• Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction?

Cited in

(4)

• A diagonal sweeping domain decomposition method with source transfer for the Helmholtz equation
• Preconditioning conjugate gradient with symmetric algebraic reconstruction technique (ART) in computerized tomography
• A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
• Fast Approximation of Algebraic Reconstruction Methods for Tomography

Describes a project that uses

Uses Software

• TVR-DART

This page was built for publication: A multi-level preconditioned Krylov method for the efficient solution of algebraic tomographic reconstruction problems

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