Technical note: A GPU-based shared Monte Carlo method for fast photon transport in multi-energy x-ray exposures

Abstract

Background

The Monte Carlo (MC) method is an accurate technique for particle transport calculation due to the precise modeling of physical interactions. Nevertheless, the MC method still suffers from the problem of expensive computational cost, even with graphics processing unit (GPU) acceleration. Our previous works have investigated the acceleration strategies of photon transport simulation for single-energy CT. But for multi-energy CT, conventional individual simulation leads to unnecessary redundant calculation, consuming more time.

Purpose

This work proposes a novel GPU-based shared MC scheme (gSMC) to reduce unnecessary repeated simulations of similar photons between different spectra, thereby enhancing the efficiency of scatter estimation in multi-energy x-ray exposures.

Methods

The shared MC method selects shared photons between different spectra using two strategies. Specifically, we introduce spectral region classification strategy to select photons with the same initial energy from different spectra, thus generating energy-shared photon groups. Subsequently, the multi-directional sampling strategy is utilized to select energy-and-direction-shared photons, which have the same initial direction, from energy-shared photon groups. Energy-and-direction-shared photons perform shared simulations, while others are simulated individually. Finally, all results are integrated to obtain scatter distribution estimations for different spectral cases.

Results

The efficiency and accuracy of the proposed gSMC are evaluated on the digital phantom and clinical case. The experimental results demonstrate that gSMC can speed up the simulation in the digital case by ∼37.8% and the one in the clinical case by ∼20.6%, while keeping the differences in total scatter results within 0.09%, compared to the conventional MC package, which performs an individual simulation.

Conclusions

The proposed GPU-based shared MC simulation method can achieve fast photon transport calculation for multi-energy x-ray exposures.