Abstract:As an alternative to the use of traditional parallel hole collimators, SPECT imaging can be performed using rotating slat collimators. While maintaining the spatial resolution, a gain in image quality could be expected from the higher photon collection efficiency of this type of collimator. However, the use of iterative methods to do fully three-dimensional (3D) reconstruction is computationally much more expensive and furthermore involves slow convergence compared to a classical SPECT reconstruction. It has been proposed to do 3D reconstruction by splitting the system matrix into two separate matrices, forcing the reconstruction to first estimate the sinograms from the rotating slat SPECT data before estimating the image. While alleviating the computational load by one order of magnitude, this split matrix approach would result in fast computation of the projections in an iterative algorithm, but does not solve the problem of slow convergence. There is thus a need for an algorithm which speeds up convergence while maintaining image quality for rotating slat collimated SPECT cameras. Therefore, we developed a reconstruction algorithm based on the split matrix approach which allows both a fast calculation of the forward and backward projection and a fast convergence. In this work, an algorithm of the maximum likelihood expectation maximization (MLEM) type, obtained from a split system matrix MLEM reconstruction, is proposed as a reconstruction method for rotating slat collimated SPECT data. Here, we compare this new algorithm to the conventional split system matrix MLEM method and to a gold standard fully 3D MLEM reconstruction algorithm on the basis of computational load, convergence and contrast-to-noise. Furthermore, ordered subsets expectation maximization (OSEM) implementations of these three algorithms are compared. Calculation of computational load and convergence for the different algorithms shows a speedup for the new method of 38 and 426 compared to the split matrix MLEM approach and the fully 3D MLEM respectively and a speedup of 16 and 21 compared to the split matrix OSEM and the fully 3D OSEM respectively. A contrast-to-noise study based on simulated data shows that our new approach has comparable accuracy as the fully 3D reconstruction method. The algorithm developed in this study allows iterative image reconstruction of rotating slat collimated SPECT data with equal image quality in a comparable amount of computation time as a classical SPECT reconstruction.

Maximum likelihood reconstruction for pinhole SPECT with a displaced center-of-rotation

Monte Carlo-based quantitative pinhole SPECT reconstruction using a ray-tracing back-projector

Projection Process Modelling for Iterative Reconstruction of Pinhole SPECT

SPECT reconstruction of combined cone beam and parallel hole collimation with experimental data

Three-dimensional SPECT reconstruction of combined cone beam and parallel beam data

Ordered Subset Expectation Maximization Algorithm for Positron Emission Tomographic Image Reconstruction Using Belief Kernels

Joint Reconstruction of Dynamic PET Activity and Kinetic Parametric Images Using Total Variation Constrained Dictionary Sparse Coding

Fast 3D iterative image reconstruction for SPECT with rotating slat collimators

LMIRA: list-mode iterative reconstruction algorithm for SPECT

Uncertainty Penalized Weighted Least Squares Framework For Pet Reconstruction Under Uncertain System Models

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.

A Robust Least Squares Framework for PET Image Reconstruction

Superiorization of EM Algorithm and Its Application in Single-Photon Emission Computed Tomography(SPECT)

Robust Framework for PET Image Reconstruction Incorporating System and Measurement Uncertainties.

Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem

A generalization of the maximum likelihood expectation maximization (MLEM) method: Masked-MLEM

System Calibration and Statistical Image Reconstruction for Ultra-High Resolution Stationary Pinhole SPECT

Fast maximum entropy approximation in SPECT using the RBI-MAP algorithm

Maximum a Posteriori Reconstruction Algorithms in Emission CT

The Simulation of Centroid Guided Maximum Likelihood Positioning Algorithm