A GPU tool for efficient, accurate, and realistic simulation of cone beam CT projections.

PURPOSE Simulation of x-ray projection images plays an important role in cone beam CT (CBCT) related research projects, such as the design of reconstruction algorithms or scanners. A projection image contains primary signal, scatter signal, and noise. It is computationally demanding to perform accurate and realistic computations for all of these components. In this work, the authors develop a package on graphics processing unit (GPU), called gDRR, for the accurate and efficient computations of x-ray projection images in CBCT under clinically realistic conditions. METHODS The primary signal is computed by a trilinear ray-tracing algorithm. A Monte Carlo (MC) simulation is then performed, yielding the primary signal and the scatter signal, both with noise. A denoising process specifically designed for Poisson noise removal is applied to obtain a smooth scatter signal. The noise component is then obtained by combining the difference between the MC primary and the ray-tracing primary signals, and the difference between the MC simulated scatter and the denoised scatter signals. Finally, a calibration step converts the calculated noise signal into a realistic one by scaling its amplitude according to a specified mAs level. The computations of gDRR include a number of realistic features, e.g., a bowtie filter, a polyenergetic spectrum, and detector response. The implementation is fine-tuned for a GPU platform to yield high computational efficiency. RESULTS For a typical CBCT projection with a polyenergetic spectrum, the calculation time for the primary signal using the ray-tracing algorithms is 1.2-2.3 s, while the MC simulations take 28.1-95.3 s, depending on the voxel size. Computation time for all other steps is negligible. The ray-tracing primary signal matches well with the primary part of the MC simulation result. The MC simulated scatter signal using gDRR is in agreement with EGSnrc results with a relative difference of 3.8%. A noise calibration process is conducted to calibrate gDRR against a real CBCT scanner. The calculated projections are accurate and realistic, such that beam-hardening artifacts and scatter artifacts can be reproduced using the simulated projections. The noise amplitudes in the CBCT images reconstructed from the simulated projections also agree with those in the measured images at corresponding mAs levels. CONCLUSIONS A GPU computational tool, gDRR, has been developed for the accurate and efficient simulations of x-ray projections of CBCT with realistic configurations.

[1]  J. Star-Lack,et al.  Improved scatter correction using adaptive scatter kernel superposition , 2010, Physics in medicine and biology.

[2]  Gene H. Golub,et al.  Matrix computations , 1983 .

[3]  Russell H. Taylor,et al.  Towards real-time radiation therapy: GPU accelerated superposition/convolution , 2010, Comput. Methods Programs Biomed..

[4]  Benoît Ozell,et al.  Fast convolution-superposition dose calculation on graphics hardware. , 2009, Medical physics.

[5]  Steve B Jiang,et al.  GPU-based ultra-fast dose calculation using a finite size pencil beam model. , 2009, Physics in medicine and biology.

[6]  Steve B. Jiang,et al.  Ultrafast treatment plan optimization for volumetric modulated arc therapy (VMAT) , 2010, 1005.4396.

[7]  J. Boone,et al.  An accurate method for computer-generating tungsten anode x-ray spectra from 30 to 140 kV. , 1997, Medical physics.

[8]  S Webb,et al.  An efficient Monte Carlo-based algorithm for scatter correction in keV cone-beam CT , 2009, Physics in medicine and biology.

[9]  J H Siewerdsen,et al.  Cone-beam computed tomography with a flat-panel imager: initial performance characterization. , 2000, Medical physics.

[10]  J. Baró,et al.  PENELOPE: An algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter , 1995 .

[11]  Iwan Kawrakow,et al.  Variance reduction techniques for fast Monte Carlo CBCT scatter correction calculations , 2010, Physics in medicine and biology.

[12]  Qiong Xu,et al.  Projection correlation based view interpolation for cone beam CT: primary fluence restoration in scatter measurement with a moving beam stop array , 2010, Physics in medicine and biology.

[13]  Zhengrong Liang,et al.  Noise properties of low-dose CT projections and noise treatment by scale transformations , 2001, 2001 IEEE Nuclear Science Symposium Conference Record (Cat. No.01CH37310).

[14]  Xuanqin Mou,et al.  X-ray projection simulation based on physical imaging model , 2006 .

[15]  Steve B. Jiang,et al.  Development of a GPU-based Monte Carlo dose calculation code for coupled electron–photon transport , 2009, Physics in medicine and biology.

[16]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation , 2010 .

[17]  Bjorn De Sutter,et al.  A Fast Algorithm to Calculate the Exact Radiological Path through a Pixel or Voxel Space , 1998 .

[18]  Thomas J. Asaki,et al.  A Variational Approach to Reconstructing Images Corrupted by Poisson Noise , 2007, Journal of Mathematical Imaging and Vision.

[19]  Junyi Xia,et al.  High performance computing for deformable image registration: Towards a new paradigm in adaptive radiotherapy. , 2008, Medical physics.

[20]  J. Wong,et al.  Flat-panel cone-beam computed tomography for image-guided radiation therapy. , 2002, International journal of radiation oncology, biology, physics.

[21]  I. Kawrakow.  Accurate condensed history Monte Carlo simulation of electron transport. I. EGSnrc, the new EGS4 version. , 2000, Medical physics.

[22]  Aldo Badano,et al.  Monte Carlo simulation of X-ray imaging using a graphics processing unit , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[23]  Steve B Jiang,et al.  Implementation and evaluation of various demons deformable image registration algorithms on a GPU. , 2010, Physics in medicine and biology.

[24]  Josh Star-Lack,et al.  Efficient scatter correction using asymmetric kernels , 2009, Medical Imaging.

[25]  Zhengrong Liang,et al.  An experimental study on the noise properties of x-ray CT sinogram data in Radon space , 2008, Physics in medicine and biology.

[26]  Aldo Badano,et al.  Accelerating Monte Carlo simulations of photon transport in a voxelized geometry using a massively parallel graphics processing unit. , 2009, Medical physics.

[27]  Weichung Wang,et al.  A fast forward projection using multithreads for multirays on GPUs in medical image reconstruction. , 2011, Medical physics.

[28]  Zhengrong Liang,et al.  Noise properties of low-dose x-ray CT sinogram data in radon space , 2008, SPIE Medical Imaging.

[29]  Xun Jia,et al.  GPU-based fast Monte Carlo simulation for radiotherapy dose calculation. , 2011, Physics in medicine and biology.

[30]  Benoît Ozell,et al.  GPUMCD: A new GPU-oriented Monte Carlo dose calculation platform. , 2011, Medical physics.

[31]  Xun Jia,et al.  MO‐FF‐A4‐05: Implementation and Evaluation of Various DRR Algorithms on GPU , 2010 .

[32]  Steve B. Jiang,et al.  Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation , 2011, Physics in Medicine and Biology.

[33]  Lei Zhu,et al.  Scatter correction for cone-beam CT in radiation therapy. , 2009 .

[34]  Xun Jia,et al.  A GPU-based finite-size pencil beam algorithm with 3D-density correction for radiotherapy dose calculation. , 2011, Physics in medicine and biology.

[35]  Iwan Kawrakow,et al.  Fast Monte Carlo calculation of scatter corrections for CBCT images , 2008 .

[36]  Steve B. Jiang,et al.  GPU-based ultrafast IMRT plan optimization , 2009, Physics in medicine and biology.

[37]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[38]  Steve B. Jiang,et al.  GPU-based iterative cone-beam CT reconstruction using tight frame regularization , 2010, Physics in medicine and biology.