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

DOI10.1016/J.CAM.2014.12.044zbMATH Open1311.65162DBLPjournals/jcam/CoolsGASV15arXiv1310.0956OpenAlexW2117347707WikidataQ62663805 ScholiaQ62663805MaRDI QIDQ2018482FDOQ2018482

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

Publication date: 24 March 2015

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

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.

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)

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

Cited In (4)

• Preconditioning conjugate gradient with symmetric algebraic reconstruction technique (ART) in computerized tomography
• A Diagonal Sweeping Domain Decomposition Method with Source Transfer for the Helmholtz Equation
• A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
• Fast Approximation of Algebraic Reconstruction Methods for Tomography

Uses Software

• TVR-DART

Recommendations

• Title not available (Why is that?) 👍 👎
• Title not available (Why is that?) 👍 👎
• Title not available (Why is that?) 👍 👎
• Title not available (Why is that?) 👍 👎
• Preconditioned alternating projection algorithms for maximum a posteriori ECT reconstruction 👍 👎
• A computationally efficient version of the algebraic method for computer tomography 👍 👎
• Preconditioning conjugate gradient with symmetric algebraic reconstruction technique (ART) in computerized tomography 👍 👎
• Superiorization of Preconditioned Conjugate Gradient Algorithms for Tomographic Image Reconstruction 👍 👎
• Fast Approximation of Algebraic Reconstruction Methods for Tomography 👍 👎
• A multilevel domain decomposition algorithm for fast \(O(N^2\log N)\) reprojection of tomographic images 👍 👎

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)