Abstract:BackgroundFull Monte Carlo (MC)-based SPECT reconstructions have a strong potential for correcting for image degrading factors, but the reconstruction times are long. The objective of this study was to develop a highly parallel Monte Carlo code for fast, ordered subset expectation maximum (OSEM) reconstructions of SPECT/CT images. The MC code was written in the Compute Unified Device Architecture language for a computer with four graphics processing units (GPUs) (GeForce GTX Titan X, Nvidia, USA). This enabled simulations of parallel photon emissions from the voxels matrix (1283 or 2563). Each computed tomography (CT) number was converted to attenuation coefficients for photo absorption, coherent scattering, and incoherent scattering. For photon scattering, the deflection angle was determined by the differential scattering cross sections. An angular response function was developed and used to model the accepted angles for photon interaction with the crystal, and a detector scattering kernel was used for modeling the photon scattering in the detector. Predefined energy and spatial resolution kernels for the crystal were used. The MC code was implemented in the OSEM reconstruction of clinical and phantom 177Lu SPECT/CT images. The Jaszczak image quality phantom was used to evaluate the performance of the MC reconstruction in comparison with attenuated corrected (AC) OSEM reconstructions and attenuated corrected OSEM reconstructions with resolution recovery corrections (RRC).ResultThe performance of the MC code was 3200 million photons/s. The required number of photons emitted per voxel to obtain a sufficiently low noise level in the simulated image was 200 for a 1283 voxel matrix. With this number of emitted photons/voxel, the MC-based OSEM reconstruction with ten subsets was performed within 20 s/iteration. The images converged after around six iterations. Therefore, the reconstruction time was around 3 min. The activity recovery for the spheres in the Jaszczak phantom was clearly improved with MC-based OSEM reconstruction, e.g., the activity recovery was 88% for the largest sphere, while it was 66% for AC-OSEM and 79% for RRC-OSEM.ConclusionThe GPU-based MC code generated an MC-based SPECT/CT reconstruction within a few minutes, and reconstructed patient images of 177Lu-DOTATATE treatments revealed clearly improved resolution and contrast.

Dithiothreitol-induced alterations of blood pressure, vascular reactivity and aortic microsomal calcium uptake in spontaneously hypertensive rats.

Toward real-time Monte Carlo simulation using a commercial cloud computing infrastructure

Technical Note: Radiotherapy dose calculations using GEANT4 and the Amazon Elastic Compute Cloud

A Serverless Architecture for Efficient and Scalable Monte Carlo Markov Chain Computation

CloudMC: a cloud computing application for Monte Carlo simulation

Radiation therapy calculations using an on-demand virtual cluster via cloud computing

Monte Carlo Performance Study of Virtual High Performance Computing Cluster over Cloud

Cloud-based Geophysical Inversion Modeling Using GNU Octave and MatlabMPI on Amazon EC2.

Cross-metathesis reactions as an efficient tool in the synthesis of fluorinated cyclic beta-amino acids.

A fast GPU-based Monte Carlo simulation of proton transport with detailed modeling of non-elastic interactions

Solid-liquid interface free energy in binary systems: theory and atomistic calculations for the (110) Cu-Ag interface.

Monte Carlo simulation of two Siemens Biograph PET/CT system using GATE: Image quality performance

A GPU-accelerated Monte Carlo code, RT2for coupled transport of photon, electron/positron, and neutron

Serverless Approach to Sensitivity Analysis of Computational Models

Are Clouds Ready to Accelerate Ad hoc Financial Simulations?

Molecular Dynamics Simulations on Cloud Computing and Machine Learning Platforms

UMC-PET: a fast and flexible Monte Carlo PET simulator

Fast GPU-based Monte Carlo code for SPECT/CT reconstructions generates improved 177Lu images

GPU-based fast Monte Carlo simulation for radiotherapy dose calculation

GROMACS in the Cloud: A Global Supercomputer to Speed Up Alchemical Drug Design

Effect of streptomycin on the response of Escherichia coli ribosomes to the dissociation factor.