Algebraic Iterative Reconstruction Methods Research Articles

Tomographic task solution using a dichotomous discretization scheme in polar coordinates and partial system matrices invariant to rotations

Objectives. The purpose of this work was to create an effective iterative algorithm for the tomographic reconstruction of objects with large volumes of initial data. Unlike the convolutional projection algorithm, widely used in commercial industrial and medical tomographic devices, algebraic iterative reconstruction methods use significant amounts of memory and typically involve long reconstruction times. At the same time, iterative methods enable a wider range of diagnostic tasks to be resolved where greater accuracy of reconstruction is required, as well as in cases where a limited amount of data is used for sparse-view angle shooting or shooting with a limited angular range.Methods. A feature of the algorithm thus created is the use of a polar coordinate system in which the projection system matrices are invariant with respect to the rotation of the object. This enables a signification reduction of the amount of memory required for system matrices storage and the use of graphics processors for reconstruction. Unlike the simple polar coordinate system used earlier, we used a coordinate system with a dichotomous division of the reconstruction field enabling us to ensure invariance to rotations and at the same time a fairly uniform distribution of spatial resolution over the reconstruction field.Results. A reconstruction algorithm was developed on the basis of the use of partial system matrices corresponding to the dichotomous division of the image field into partial annular reconstruction regions. A 2D and 3D digital phantom was used to show the features of the proposed reconstruction algorithm and its applicability to solving tomographic problems.Conclusions. The proposed algorithm allows algebraic image reconstruction to be implemented using standard libraries for working with sparse matrices based on desktop computers with graphics processors.

Fixing Nonconvergence of Algebraic Iterative Reconstruction with an Unmatched Backprojector

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.

Unmatched Projector/Backprojector Pairs: Perturbation and Convergence Analysis

In tomographic reconstruction problems it is not uncommon that there are errors in the implementation of the forward projector and/or the backprojector, and hence we encounter a so-called unmatched projektor/backprojector pair. Consequently, the matrices that represent the two projectors are not each other's transpose. Surprisingly, the influence of such errors in algebraic iterative reconstruction methods has received little attention in the literature. The goal of this paper is to perform a rigorous first-order perturbation analysis of the minimization problems underlying the algebraic methods in order to understand the role played by the nonmatch of the matrices. We also study the convergence properties of linear stationary iterations based on unmatched matrix pairs, leading to insight into the behavior of some important row- and column-oriented algebraic iterative methods. We conclude with numerical examples that illustrate the perturbation and convergence results.

AIR Tools II: algebraic iterative reconstruction methods, improved implementation

We present a MATLAB software package with efficient, robust, and flexible implementations of algebraic iterative reconstruction (AIR) methods for computing regularized solutions to discretizations of inverse problems. These methods are of particular interest in computed tomography and similar problems where they easily adapt to the particular geometry of the problem. All our methods are equipped with stopping rules as well as heuristics for computing a good relaxation parameter, and we also provide several test problems from tomography. The package is intended for users who want to experiment with algebraic iterative methods and their convergence properties. The present software is a much expanded and improved version of the package AIR Tools from 2012, based on a new modular design. In addition to improved performance and memory use, we provide more flexible iterative methods, a column-action method, new test problems, new demo functions, and—perhaps most important—the ability to use function handles instead of (sparse) matrices, allowing larger problems to be handled.

Proton Computed Tomography: iterative image reconstruction and dose evaluation

Proton Computed Tomography (pCT) is a medical imaging method with a potential for increasing accuracy of treatment planning and patient positioning in hadron therapy. A pCT system based on a Silicon microstrip tracker and a YAG:Ce crystal calorimeter has been developed within the INFN Prima-RDH collaboration. The prototype has been tested with a 175 MeV proton beam at The Svedberg Laboratory (Uppsala, Sweden) with the aim to reconstruct and characterize a tomographic image. Algebraic iterative reconstruction methods (ART), together with the most likely path formalism, have been used to obtain tomographies of an inhomogeneous phantom to eventually extract density and spatial resolutions. These results will be presented and discussed together with an estimation of the average dose delivered to the phantom and the dependence of the image quality on the dose. Due to the heavy computation load required by the algebraic algorithms the reconstruction programs have been implemented to fully exploit the high calculation parallelism of Graphics Processing Units. An extended field of view pCT system is in an advanced construction stage. This apparatus will be able to reconstruct objects of the size of a human head making possible to characterize this pCT approach in a pre-clinical environment.

Multicore Performance of Block Algebraic Iterative Reconstruction Methods

Algebraic iterative methods are routinely used for solving the ill-posed sparse linear systems arising in tomographic image reconstruction. Here we consider the algebraic reconstruction technique (ART) and the simultaneous iterative reconstruction techniques (SIRT), both of which rely on semiconvergence. Block versions of these methods, based on a partitioning of the linear system, are able to combine the fast semiconvergence of ART with the better multicore properties of SIRT. These block methods separate into two classes: those that, in each iteration, access the blocks in a sequential manner, and those that compute a result for each block in parallel and then combine these results before the next iteration. The goal of this work is to demonstrate which block methods are best suited for implementation on modern multicore computers. To compare the performance of the different block methods, we use a fixed relaxation parameter in each method, namely, the one that leads to the fastest semiconvergence. Comp...

AIR Tools — A MATLAB package of algebraic iterative reconstruction methods

We present a MATLAB package with implementations of several algebraic iterative reconstruction methods for discretizations of inverse problems. These so-called row action methods rely on semi-convergence for achieving the necessary regularization of the problem. Two classes of methods are implemented: Algebraic Reconstruction Techniques (ART) and Simultaneous Iterative Reconstruction Techniques (SIRT). In addition we provide a few simplified test problems from medical and seismic tomography. For each iterative method, a number of strategies are available for choosing the relaxation parameter and the stopping rule. The relaxation parameter can be fixed, or chosen adaptively in each iteration; in the former case we provide a new “training” algorithm that finds the optimal parameter for a given test problem. The stopping rules provided are the discrepancy principle, the monotone error rule, and the NCP criterion; for the first two methods “training” can be used to find the optimal discrepancy parameter.