Rapid scatter estimation for CBCT using the Boltzmann transport equation

Scatter in cone-beam computed tomography (CBCT) is a significant problem that degrades image contrast, uniformity and CT number accuracy. One means of estimating and correcting for detected scatter is through an iterative deconvolution process known as scatter kernel superposition (SKS). While the SKS approach is efficient, clinically significant errors on the order 2-4% (20-40 HU) still remain. We have previously shown that the kernel method can be improved by perturbing the kernel parameters based on reference data provided by limited Monte Carlo simulations of a first-pass reconstruction. In this work, we replace the Monte Carlo modeling with a deterministic Boltzmann solver (AcurosCTS) to generate the reference scatter data in a dramatically reduced time. In addition, the algorithm is improved so that instead of adjusting kernel parameters, we directly perturb the SKS scatter estimates. Studies were conducted on simulated data and on a large pelvis phantom scanned on a tabletop system. The new method reduced average reconstruction errors (relative to a reference scan) from 2.5% to 1.8%, and significantly improved visualization of low contrast objects. In total, 24 projections were simulated with an AcurosCTS execution time of 22 sec/projection using an 8-core computer. We have ported AcurosCTS to the GPU, and current run-times are approximately 4 sec/projection using two GPU’s running in parallel.

[1]  Mark Oldham,et al.  Cone-beam-CT guided radiation therapy: A model for on-line application. , 2005, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[2]  Rebecca Fahrig,et al.  A nonlinear lag correction algorithm for a-Si flat-panel x-ray detectors. , 2012, Medical physics.

[3]  K. Klingenbeck,et al.  A general framework and review of scatter correction methods in cone beam CT. Part 2: scatter estimation approaches. , 2011, Medical physics.

[4]  Firas Mourtada,et al.  Validation of a new grid-based Boltzmann equation solver for dose calculation in radiotherapy with photon beams , 2010, Physics in medicine and biology.

[5]  J H Siewerdsen,et al.  Intraoperative cone-beam CT for guidance of head and neck surgery: Assessment of dose and image quality using a C-arm prototype. , 2006, Medical physics.

[6]  Klaus Klingenbeck,et al.  A general framework and review of scatter correction methods in x-ray cone-beam computerized tomography. Part 1: Scatter compensation approaches. , 2011, Medical physics.

[7]  Edward W. Larsen,et al.  Convergence Rates of Spatial Difference Equations for the Discrete-Ordinates Neutron Transport Equations in Slab Geometry , 1980 .

[8]  Shawn D. Pautz,et al.  Discontinuous Finite Element SN Methods on Three-Dimensional Unstructured Grids , 2001 .

[9]  D. Jaffray,et al.  Advances in image-guided radiation therapy. , 2007, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[10]  Josh Star-Lack,et al.  User-friendly, ultra-fast simulation of detector DQE(f) , 2013, Medical Imaging.

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

[12]  K. Bush,et al.  Dosimetric validation of Acuros XB with Monte Carlo methods for photon dose calculations. , 2011, Medical physics.

[13]  Rebecca Fahrig,et al.  Investigation into the optimal linear time-invariant lag correction for radar artifact removal. , 2011, Medical physics.

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

[15]  Firas Mourtada,et al.  Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations , 2006, Physics in medicine and biology.

[16]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .