Abstract:The aim of this work is to develop, validate and preliminarily evaluate two graphic processing units (GPU) implementations of the OS-EM iterative image reconstruction algorithms for a clinical SPECT system. The Monte Carlo software SIMSET was used to simulate the data acquisition of a typical SPECT camera with a parallel-hole collimator for the Shepp-Logan phantom. CUDA was adopted for the GPU programming of the OS-EM algorithm with high complexity. For the first implementation, the system matrix used for reconstruction was pre-computed by Monte Carlo simulation. A major task in the GPU programming in this case was to efficiently access every element of the pre-stored system matrix and compute the projection and back-projection for each image update. In the second implementation, the system matrix was calculated on-the-fly. A rotation-based projector and backprojector was employed, which reduced the complexity for the GPU programming and can achieve higher efficiency theoretically than other methods of computing projector and back-projectors. We modeled the depth-dependent collimator detector response in the matrix computation. The results showed that, for both GPU implementations, there are seldom differences of images reconstructed by the GPU and CPU implementation of the OS-EM algorithm. In the first matrix computation strategy, an acceleration ratio of 31.9 was achieved by using the GPU implementation compared with the CPU implementation of the same algorithm. In the second strategy, the acceleration ratio was 33.9. The images reconstructed with the first strategy demonstrated higher quantitative accuracy due to the more accurate modeling of the system matrix. But the corresponding reconstructed time was more than the second strategy . We also applied the second GPU implementation of the OS-EM method to clinical SPECT patient cardiac data and achieved the reconstruction of a 64×64×64 image within 8 seconds. Compared with the FBP reconstruction image, the OS-EM reconstruction image had a better performance. In conclusion, GPU implementation of the OS-EM algorithm has a great potential in SPECT image reconstruction.

Toward Parallel Optimal Computation of Ultrasound Computed Tomography Using GPU

Sound Speed Image Reconstruction of Ultrasound Computed Tomography Based on GPU Acceleration

Parallel processing of SAFT ultrasonic imaging technique based on GPU

GPU-Based Parallel Computing Methods for Medical Photoacoustic Image Reconstruction

Fast GPU-Based Automatic Time Gain Compensation for Ultrasound Imaging

A Graphics Processing Unit Implementation and Optimization for Parallel Double-Difference Seismic Tomography

Parallel Multi-threaded Gridrec Algorithm for Computer Tomography on GPU for Edge Computing

High Speed Parallel Processing of Biomedical Optics Data with PC Graphic Hardware

Development of GPU Based Image Reconstruction Method for Clinical SPECT

Accelerating Image Reconstruction in Three-Dimensional Optoacoustic Tomography on Graphics Processing Units

Research Progress of GPU Applications in Medical Ultrasound Imaging

Graphic Processing Unit Based Phase Retrieval and CT Reconstruction for Differential X-ray Phase Contrast Imaging

Fast Cone-Beam CT Image Reconstruction Using GPU Hardware

Accelerate Helical Cone-Beam CT with Graphics Hardware

Parallel Computing For Quantitative Blood Flow Imaging In Photoacoustic Microscopy

GPU-based ultra fast dose calculation using differential convolution/superposition algorithm

Cs-Based Fast Ultrasound Imaging With Improved Fista Algorithm

Feasibility of GPU-assisted Iterative Image Reconstruction for Mobile C-arm CT

GPU-based 3D cone-beam CT image reconstruction: application to micro CT

GPU-Based Real-Time Software Coincidence Processing for Digital PET System

Cone-beam computed tomography reconstruction accelerated with CUDA