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.

Quick reconstruction method of double-camera spectral imaging system based on GPU

Fast Image Super-Resolution for a Dual-Resolution Camera

Superpixel-Based Linear Reconstruction Method for Dual-Camera Compressed Hyperspectral Imaging System

Deeply Learned Broadband Encoding Stochastic Hyperspectral Imaging.

Residual image recovery method based on the dual-camera design of a compressive hyperspectral imaging system

Hyperspectral image recovery based on fusion of coded aperture snapshot spectral imaging and RGB images by guided filtering

Compressive Single-Pixel Hyperspectral Imaging Using RGB Sensors

Fast Multispectral Imaging by Spatial Pixel-Binning and Spectral Unmixing

Fast non-iterative blind restoration of hyperspectral images with PSFs

GPU-based Acceleration of the Hyperspectral Band Selection by SNR Estimation Using Wavelet Transform

Parallel Image Resample Algorithm Based on GPU for Land Remote Sensing Data Management

method for accelerating decoding of a camera video stream through a programmable GPU in a virtual-real fusion system

High-resolution Image Reconstruction Based on Independent Equations of Pixel Grey Level

A flexible and accurate total variation and cascaded denoisers-based image reconstruction algorithm for hyperspectrally compressed ultrafast photography

Accelerated CT Reconstruction Using GPU SIMD Parallel Computing with Bilinear Warping Method

An improved method based on twIst and dual-camera design for hyperspectral image reconstruction

Method for GPU-based spectral data cube reconstruction of integral field snapshot imaging spectrometers

Dual camera snapshot hyperspectral imaging system via physics informed learning

Near Real-Time Automatic Sub-Pixel Registration of Panchromatic and Multispectral Images for Pan-Sharpening

Development of GPU Based Image Reconstruction Method for Clinical SPECT

Fast high quality computational ghost imaging based on saliency variable sampling detection