Abstract:Purpose: Performance of the preconditioned alternating projection algorithm (PAPA) using relaxed ordered subsets (ROS) with a non-smooth penalty function was investigated in positron emission tomography (PET). A higher order total variation (HOTV) regularizer was applied and a method for unsupervised selection of penalty weights based on the measured data is introduced. Methods: A ROS version of PAPA with HOTV penalty (ROS-HOTV-PAPA) for PET image reconstruction was developed and implemented. Two-dimensional PET data were simulated using two synthetic phantoms (geometric and brain) in geometry similar to GE D690/710 PET/CT with uniform attenuation, and realistic scatter (25%) and randoms (25%). Three count levels (high/medium/low) corresponding to mean information densities (ID¯s) of 125, 25, and 5 noise equivalent counts (NEC) per support voxel were reconstructed using ROS-HOTV-PAPA. The patients' brain and whole body PET data were acquired at similar ID¯s on GE D690 PET/CT with time-of-fight and were reconstructed using ROS-HOTV-PAPA and available clinical ordered-subset expectation-maximization (OSEM) algorithms. A power-law model of the penalty weights' dependence on ID¯ was semi-empirically derived. Its parameters were elucidated from the data and used for unsupervised selection of the penalty weights within a reduced search space. The resulting image quality was evaluated qualitatively, including reduction of staircase artifacts, image noise, spatial resolution and contrast, and quantitatively using root mean squared error (RMSE) as a global metric. The convergence rates were also investigated. Results: ROS-HOTV-PAPA converged rapidly, in comparison to non-ROS-HOTV-PAPA, with no evidence of limit cycle behavior. The reconstructed image quality was superior to optimally post-filtered OSEM reconstruction in terms of noise, spatial resolution, and contrast. Staircase artifacts were not observed. Images of the measured phantom reconstructed using ROS-HOTV-PAPA showed reductions in RMSE of 5%-44% as compared with optimized OSEM. The greatest improvement occurred in the lowest count images. Further, ROS-HOTV-PAPA reconstructions produced images with RMSE similar to images reconstructed using optimally post-filtered OSEM but at one-quarter the NEC. Conclusion: Acceleration of HOTV-PAPA was achieved using ROS. This was accompanied by an improved RMSE metric and perceptual image quality that were both superior to that obtained with either clinical or optimized OSEM. This may allow up to a four-fold reduction of the radiation dose to the patients in a PET study, as compared with current clinical practice. The proposed unsupervised parameter selection method provided useful estimates of the penalty weights for the selected phantoms' and patients' PET studies. In sum, the outcomes of this research indicate that ROS-HOTV-PAPA is an appropriate candidate for clinical applications and warrants further research.

A relaxed OSEM algorithm applied in 60Co container CT image reconstruction

OSEM Reconstruction with Overrelaxation

A novel reconstruction approach combining OSEM and split Bregman method for low dose CT

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

Evaluation of the 3D IMF-OSEM Algorithm by Using Data from a High Resolution PET Scanner

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

Research on Some Key Problems in the OSEM Algorithm

Reconstruction Algorithm for Regions of Interest in γ-Photon Images

A Low-dose CT Reconstruction Network Based on TV-regularized OSEM Algorithm

Exploration of Optimal Reconstruction Parameters for 18F-FDG Total-body PET/CT with Ultra-low Activity Injection

Millisecond SPECT Image Reconstruction Acceleration Using OSEM Algorithm Based on CPU / GPU Hybrid Parallel Hardware Platform

Optimization of the Algorithm for the Implementation of Point Spread Function in the 3D-OSEM Reconstruction Algorithm Based on the List-Mode Micro PET Data

THE SIMULATED STUDY OF VARIABLE CONSTRAINT ORDERED SUBSETS EXPECTATION MAXIMIZATION (VCOSEM)

Communication Optimization and Auto Load Balancing in Parallel OSEM Algorithm for Fully 3-D SPECT Reconstruction

Pde Regularization for Bayesian Reconstruction of Emission Tomography

An Iterative Image Reconstruction Algorithm for Spect

Image Quality of Solid Tumors of Liver: A Comparison of Automatic Exposure Control and Organ Effective Modulation in Abdominopelvic CT

Relaxed ordered subset preconditioned alternating projection algorithm for PET reconstruction with automated penalty weight selection

A Robust Least Squares Framework for PET Image Reconstruction

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

Performance of the Iterative OSEM and HYPER Algorithm for Total-body PET at SUVmax with a Low 18F-FDG Activity, a Short Acquisition Time and Small Lesions