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

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

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

• A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
• Title not available (Why is that?)
• 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
• Title not available (Why is that?)
• 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
• Title not available (Why is that?)
• DART: A Practical Reconstruction Algorithm for Discrete Tomography
• Recent computational developments in Krylov subspace methods for linear systems
• Title not available (Why is that?)
• A Multigrid Tutorial, Second Edition
• Title not available (Why is that?)
• Discrete inverse problems. Insight and algorithms.
• Title not available (Why is that?)
• Title not available (Why is that?)
• 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

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)