No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
The full text of this article is not currently available.
Efficient scatter distribution estimation and correction in CBCT using concurrent Monte Carlo fitting
1.J. M. Boone and J. A. Seibert, “Monte Carlo simulation of the scattered radiation distribution in diagnostic radiology,” Med. Phys. 15, 713–720 (1988).
2.J. H. Siewerdsen and D. A. Jaffray, “Cone-beam computed tomography with a flat-panel imager: Magnitude and effects of x-ray scatter,” Med. Phys. 28, 220–231 (2001).
3.G. Jarry, S. A. Graham, D. J. Moseley, D. A. Jaffray, J. H. Siewerdsen, and F. Verhaegen, “Characterization of scattered radiation in kV CBCT images using Monte Carlo simulations,” Med. Phys. 33, 4320–4329 (2006).
4.G. J. Bootsma, F. Verhaegen, and D. A. Jaffray, “Spatial and angular frequency spectrum of the scatter distribution in CBCT,” Med. Phys. 40, 111901 (15pp.) (2013).
7.J. O. Deasy, M. V. Wickerhauser, and M. Picard, “Accelerating Monte Carlo simulations of radiation therapy dose distributions using wavelet threshold de-noising,” Med. Phys. 29, 2366–2373 (2002).
8.I. El Naqa, I. Kawrakow, M. Fippel, J. V. Siebers, P. E. Lindsay, M. V. Wickerhauser, M. Vicic, K. Zakarian, N. Kauffmann, and J. O. Deasy, “A comparison of Monte Carlo dose calculation denoising techniques,” Physics in Medicine and Biology 50, 909–922 (2005).
10.W. Zbijewski and F. J. Beekman, “Efficient Monte Carlo based scatter artifact reduction in cone-beam micro-CT,” IEEE Trans. Med. Imaging 25, 817–827 (2006).
13.G. J. Bootsma, F. Verhaegen, and D. A. Jaffray, “The effects of compensator and imaging geometry on the distribution of x-ray scatter in CBCT,” Med. Phys. 38, 897–914 (2011).
14.G. J. Bootsma, F. Vehaegen, and D. A. Jaffray, “The effects of compensator design on scatter distribution and magnitude: A Monte Carlo study,” Proc. SPIE 7961, 796109-1–796109-13 (2011).
15.S. A. Graham, D. J. Moseley, J. H. Siewerdsen, and D. A. Jaffray, “Compensator for dose and scatter management in cone-beam computed tomography,” Med. Phys. 34, 2691–2703 (2007).
16.N. Mail, D. J. Moseley, J. H. Siewerdsen, and D. A. Jaffray, “The influence of bowtie filtration on cone-beam CT image quality,” Med. Phys. 36, 22–32 (2009).
17.N. Freud, J. M. Letang, and D. Babot, “A hybrid approach to simulate multiple photon scattering in x-ray imaging,” Nucl. Instrum. Methods Phys. Res. B 227, 551–558 (2005).
18.G. Pludniowski, P. M. Evans, V. N. Hansen, and S. Webb, “An efficient Monte Carlo-based algorithm for scatter correction in keV cone-beam CT,” Phys. Med. Biol. 54, 3847–3864 (2009).
20.A. Badal and A. Badano, “Accelerating Monte Carlo simulations of photon transport in a voxelized geometry using a massively parallel graphics processing unit,” Med. Phys. 36, 4878–4880 (2009).
21.I. Kawrakow and D. W. O. Rogers, “The EGSnrc code system: Monte Carlo simulation of electron and photon transport,” National Research Council of Canada Technical Report No. PIRS-701 (National Research Council of Canada, Ottawa, 2000).
24.C. B. Chiarot, J. H. Siewerdsen, T. Haycocks, D. J. Moseley, and D. A. Jaffray, “An innovative phantom for quantitative and qualitative investigation of advanced x-ray imaging technologies,” Phys. Med. Biol. 50, N287–N297 (2005).
25.D. W. O. Rogers, B. Faddegon, G. X. Ding, C. M. Ma, J. Wei, and T. R. Mackie, “Beam: A Monte Carlo code to simulate radiotherapy treatment units,” Med. Phys. 22, 503–524 (1995).
26.H. Gao, R. Fahrig, N. R. Bennett, M. Sun, J. Star-Lack, and L. Zhu, “Scatter correction method for x-ray CT using primary modulation: Phantom studies,” Med. Phys. 37, 934–946 (2010).
27.L. Zhu, N. R. Bennett, and R. Fahrig, “Scatter correction method for x-ray CT using primary modulation: Theory and preliminary results,” IEEE Trans. Med. Imaging 25, 1573–1587 (2006).
29.J. H. Siewerdsen and D. A. Jaffray, “A ghost story: Spatio-temporal response characteristics of an indirect-detection flat-panel imager,” Med. Phys. 33, 1624–1641 (1999).
30.N. Mail, D. J. Moseley, J. H. Siewerdsen, and D. A. Jaffray, “An empirical method for lag correction in cone-beam CT,” Med. Phys. 35, 5187–5196 (2008).
32.B. Haas, T. Coradi, M. Scholz, P. Kunz, M. Huber, U. Oppitz, L. André, V. Lengkeek, D. Huyskens, A. van esch, and R. Reddick, “Automatic segmentation of thoracic and pelvic CT images for radiotherapy planning using implicit anatomic knowledge and organ-specific segmentation strategies,” Phys. Med. Biol. 53, 1751–1771 (2008).
34.M. Bazalova, L. Beaulieu, S. Palefsky, and F. Verhaegen, “Correction of CT artifacts and its influence on Monte Carlo dose calculations,” Med. Phys. 34, 2120–2132 (2007).
37.A. Sisniega, M. Abella, E. Lage, M. Desco, and J. J. Vaquero, “Automatic Monte-Carlo based scatter correction for x-ray cone-beam CT using general purpose graphic processing units (GP-GPU): A feasibility study,” IEEE Nuclear Science Symposium and Medical Imaging Conference (Valencia, 2011), pp. 3705–3709.
40.W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes: The Art of Scientific Computing, 3rd ed. (Cambridge University Press, Cambridge, England, 2007).
41.R. S. Thing and E. Mainegra-Hing, “Optimizing cone beam CT scatter estimation in egs_cbct for a clinical and virtual chest phantom,” Med. Phys. 41(7), 071902 (7pp) (2014).
Article metrics loading...
X-ray scatter is a significant impediment to image quality improvements in cone-beam CT
(CBCT). The authors present and demonstrate a novel scatter correction algorithm using a scatter estimation method that simultaneously combines multiple Monte Carlo
CBCT simulations through the use of a concurrently evaluated fitting function, referred to as concurrent MC fitting (CMCF).
The CMCF method uses concurrently run MC
scatter projection simulations that are a subset of the projection angles used in the projection set, P, to be corrected. The scattered
photons reaching the detector in each MC simulation are simultaneously aggregated by an algorithm which computes the scatter detector response, S
is fit to a function, SF
, and if the fit of SF
is within a specified goodness of fit (GOF), the simulations are terminated. The fit, SF
, is then used to interpolate the scatter distribution over all pixel locations for every projection angle in the set P. The CMCF algorithm was tested using a frequency limited sum of sines and cosines as the fitting function on both simulated and measured data. The simulated data consisted of an anthropomorphic head and a pelvis phantom created from CT data, simulated with and without the use of a compensator. The measured data were a pelvis scan of a phantom and patient taken on an Elekta Synergy platform. The simulated data were used to evaluate various GOF metrics as well as determine a suitable fitness value. The simulated data were also used to quantitatively evaluate the image quality improvements provided by the CMCF method. A qualitative analysis was performed on the measured data by comparing the CMCF scatter corrected reconstruction to the original uncorrected and corrected by a constant scatter correction reconstruction, as well as a reconstruction created using a set of projections taken with a small cone angle.
Pearson’s correlation, r, proved to be a suitable GOF metric with strong correlation with the actual error of the scatter fit, SF
. Fitting the scatter distribution to a limited sum of sine and cosine functions using a low-pass filtered fast Fourier transform provided a computationally efficient and accurate fit. The CMCF algorithm reduces the number of photon histories required by over four orders of magnitude. The simulated experiments showed that using a compensator reduced the computational time by a factor between 1.5 and 1.75. The scatter estimates for the simulated and measured data were computed between 35–93 s and 114–122 s, respectively, using 16 Intel Xeon cores (3.0 GHz). The CMCF scatter correction improved the contrast-to-noise ratio by 10%–50% and reduced the reconstruction error to under 3% for the simulated phantoms.
The novel CMCF algorithm significantly reduces the computation time required to estimate the scatter distribution by reducing the statistical noise in the MC
scatter estimate and limiting the number of projection angles that must be simulated. Using the scatter estimate provided by the CMCF algorithm to correct both simulated and real projection data showed improved reconstruction image quality.
Full text loading...
Most read this month