A Metropolis Monte Carlo Simulation Scheme for Fast Scattering Calculation in Cone Beam CT

Yusi Chen, Jianhui Ma, Bin Li, Zhen Tian, Linghong Zhou, Xun Jia, Yuan Xu


Published in:Fully3D 2017 Proceedings


Monte Carlo simulation, metropolis hasting sampling, path-by-path scheme, scattering calculation
Low computational efficiency is the largest obstacle hindering practical applications of Monte Carlo particle transport simulation. We develop a new MC scheme for photon transport that employs Metropolis-Hasting sampling algorithm to conduct the simulations via a path-by-path scheme in this GPU-based Metropolis MC (gMMC) package. By using the Metropolis algorithm, gMMC is able to sample an entire path of a photon each time, extending from the x-ray source to the detector. The sampled paths over a long run will follow a distribution governed by the particle transport physics. gMMC was benchmarked against an in-house developed GPU-based Monte Carol simulation tool gMCDRR that performs simulations in the conventional particle-by-particle scheme. An example problem of x-ray scatter calculation was studied. For pencil beam case, speed-up factors of 200 to 300 times for first order scatter and 12 to 18 times for second order scatter can be achieved, while the average differences comparing to the conventional approach are within 1%.
Yusi Chen
Southern Medical University
Jianhui Ma
Southern Medical University
Bin Li
Southern Medical University
Zhen Tian
UT Southern Medical Center
Linghong Zhou
Southern Medical University
Xun Jia
UT Southern Medical Center
Yuan Xu
Southern Medical University
  1. Y. Xu, T. Bai, H. Yan, L. Ouyang, A. Pompos, J. Wang, L. Zhou, S. B. Jiang, and X. Jia, "A practical cone-beam CT scatter correction method with optimized Monte Carlo simulations for image-guided radiation therapy," Physics in Medicine & Biology, vol. 60, pp. 3567-87, 2015.
  2. G. Pratx and L. Xing, "GPU computing in medical physics: a review," Medical Physics, vol. 38, pp. 2685-97, 2011.
  3. X. Jia, P. Ziegenhein, and S. B. Jiang, "GPU-based high-performance computing for radiation therapy,“ Physics in Medicine and Biology, vol. 59, p. R151, 2014.
  4. X. Jia, H. Yan, L. Cervino, M. Folkerts, and S. B. Jiang, "A GPU tool for efficient, accurate, and realistic simulation of cone beam CT projections," Med Phys, vol. 39, pp. 7368-78, Dec 2012.
  5. X. Jia, H. Yan, X. Gu, and S. B. Jiang, "Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation," Phys Med Biol, vol. 57, pp. 577-90, Feb 07 2012.
  6. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, "Equation of State Calculations by Fast Computing Machines," Journal of Biochemical & Biophysical Methods, vol. 21, pp. 1087-1092, 1952.
  7. W. K. Hastings, "Monte Carlo sampling methods using Markov chains and their applications," Biometrika, vol. 57, pp. 97-109, 1970.
  8. C. P. Robert and G. Casella, The Metropolis—Hastings Algorithm: Springer New York, 2004.