Variance reduction techniques for fast Monte Carlo CBCT scatter correction calculations

Several variance reduction techniques improving the efficiency of the Monte Carlo estimation of the scatter contribution to a cone beam computed tomography (CBCT) scan were implemented in egs_ctct, an EGSnrc-based application for CBCT-related calculations. The largest impact on the efficiency comes from the splitting + Russian Roulette techniques which are described in detail. The fixed splitting technique is outperformed by both the position-dependent importance splitting (PDIS) and the region-dependent importance splitting (RDIS). The superiority of PDIS over RDIS observed for a water phantom with bone inserts is not observed when applying these techniques to a more realistic human chest phantom. A maximum efficiency improvement of several orders of magnitude over an analog calculation is obtained. A scatter calculation combining the reported efficiency gain with a smoothing algorithm is already in the proximity of being of practical use if a medium size computer cluster is available.

[1]  J. H. Hubbell,et al.  XCOM: Photon cross sections on a personal computer , 1987 .

[2]  F. H. Attix,et al.  5 – Monte Carlo Techniques of Electron and Photon Transport for Radiation Dosimetry , 1990 .

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

[4]  D. Rogers,et al.  EGS4 code system , 1985 .

[5]  D. Jaffray,et al.  Cone-beam computed tomography with a flat-panel imager: magnitude and effects of x-ray scatter. , 2001, Medical physics.

[6]  I. Lux.  Monte Carlo Particle Transport Methods: Neutron and Photon Calculations , 1991 .

[7]  Michael Sandborg,et al.  CTmod - A toolkit for Monte Carlo simulation of projections including scatter in computed tomography , 2008, Comput. Methods Programs Biomed..

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

[9]  S. Webb,et al.  Rayleigh scatter in kilovoltage x-ray imaging: is the independent atom approximation good enough? , 2009, Physics in medicine and biology.

[10]  Freek J. Beekman,et al.  Efficient Monte Carlo based scatter artifact reduction in cone-beam micro-CT , 2006, IEEE Transactions on Medical Imaging.

[11]  Kenneth R. Kase,et al.  The Dosimetry of Ionizing Radiation , 1986 .

[12]  W. Kalender,et al.  Combining deterministic and Monte Carlo calculations for fast estimation of scatter intensities in CT , 2006, Physics in medicine and biology.

[13]  Frank Verhaegen,et al.  Characterization of scattered radiation in kV CBCT images using Monte Carlo simulations. , 2006, Medical physics.

[14]  J. W. Vieira,et al.  MAX06 and FAX06: update of two adult human phantoms for radiation protection dosimetry , 2006, Physics in medicine and biology.

[15]  R. Cloutier.  Tissue Substitutes in Radiation Dosimetry and Measurement. , 1989 .

[16]  Freek J. Beekman,et al.  Accelerated simulation of cone beam X-ray scatter projections , 2004, IEEE Transactions on Medical Imaging.

[17]  I Kawrakow,et al.  On the de-noising of Monte Carlo calculated dose distributions. , 2002, Physics in medicine and biology.

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

[19]  F. H. Attix,et al.  The Dosimetry of Ionizing Radiation, Vol. 2 , 1987 .