Dataset Information

Computed tomography dose index and dose length product for cone-beam CT: Monte Carlo simulations.

ABSTRACT: Dosimetry in kilovoltage cone beam computed tomography (CBCT) is a challenge due to the limitation of physical measurements. To address this, we used a Monte Carlo (MC) method to estimate the CT dose index (CTDI) and the dose length product (DLP) for a commercial CBCT system. As Dixon and Boone showed that CTDI concept can be applicable to both CBCT and conventional CT, we evaluated weighted CT dose index (CTDI(w)) and DLP for a commercial CBCT system. Two extended CT phantoms were created in our BEAMnrc/EGSnrc MC system. Before the simulations, the beam collimation of a Varian On-Board Imager (OBI) system was measured with radiochromic films (model: XR-QA). The MC model of the OBI X-ray tube, validated in a previous study, was used to acquire the phase space files of the full-fan and half-fan cone beams. Then, DOSXYZnrc user code simulated a total of 20 CBCT scans for the nominal beam widths from 1 cm to 10 cm. After the simulations, CBCT dose profiles at center and peripheral locations were extracted and integrated (dose profile integral, DPI) to calculate the CTDI per each beam width. The weighted cone-beam CTDI (CTDI(w,l)) was calculated from DPI values and mean CTDI(w,l) (CTDI(w,l)) and DLP were derived. We also evaluated the differences of CTDI(w) values between MC simulations and point dose measurements using standard CT phantoms. In results, it was found that CTDI(w,600) was 8.74 ± 0.01 cGy for head and CTDI(w,900) was 4.26 ± 0.01 cGy for body scan. The DLP was found to be proportional to the beam collimation. We also found that the point dose measurements with standard CT phantoms can estimate the CTDI within 3% difference compared to the full integrated CTDI from the MC method. This study showed the usability of CTDI as a dose index and DLP as a total dose descriptor in CBCT scans.

SUBMITTER: Kim S

PROVIDER: S-EPMC5718669 | biostudies-literature | 2011 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Computed tomography dose index and dose length product for cone-beam CT: Monte Carlo simulations.

Kim Sangroh S Song Haijun H Samei Ehsan E Yin Fang-Fang FF Yoshizumi Terry T TT

Journal of applied clinical medical physics 20110119 2

Dosimetry in kilovoltage cone beam computed tomography (CBCT) is a challenge due to the limitation of physical measurements. To address this, we used a Monte Carlo (MC) method to estimate the CT dose index (CTDI) and the dose length product (DLP) for a commercial CBCT system. As Dixon and Boone showed that CTDI concept can be applicable to both CBCT and conventional CT, we evaluated weighted CT dose index (CTDI(w)) and DLP for a commercial CBCT system. Two extended CT phantoms were created in ou ...[more]

PMID: 21587186

Similar Datasets

Project description:The variety of cone-beam computed tomography (CBCT) machines and their applications has rapidly increased in recent years, making the dose evaluation of individual devices an important issue. Patient doses from CBCT were assessed with two different methods: optically stimulated luminescence dosimeter (OSLD) measurements and Monte Carlo (MC) simulations, in four different examination modes. Based on an analysis of the measurement process and the obtained values, a recommendation is made regarding which method is more practical and efficient for acquiring the effective dose of CBCT. Twenty-two OSLDs were calibrated and equipped in human phantoms of head and neck organs. They were exposed to radiation from two CBCT units-CS9300 (Carestream Dental LLC, Atlanta, Georgia) and RAYSCAN α+ (Ray Co. Ltd, Hwaseong-si, Korea)-using two different examination modes. The dose recorded using the OSLDs was used to calculate the organ dose and the effective dose for each unit in each examination mode. These values were also calculated using MC software, PCXMC (STUK, Helsinki, Finland). The organ doses and effective doses obtained using both methods were compared for each examination mode of the individual units. The OSLD-measured effective dose value was higher than that obtained using the MC method for each examination mode, except the dual jaw mode of CS9300. The percent difference of the effective dose between the two methods ranged from 4.0% to 14.3%. The dose difference between the methods decreased as the field of view became smaller. The organ dose values varied according to the method, although the overall trend was similar for both methods. The organs showing high doses were mostly consistent for both methods. In this study, the effective dose obtained by OSLD measurements and MC simulations were compared, and both methods were described in detail. As a relatively efficient and easy-to-perform method, we cautiously suggest using MC simulations for dose evaluations in the future.

Project description:The purpose of this study is to adapt an equivalent source model originally developed for conventional CT Monte Carlo dose quantification to the radiation oncology context and validate its application for evaluating concomitant dose incurred by a kilovoltage (kV) cone-beam CT (CBCT) system integrated into a linear accelerator.In order to properly characterize beams from the integrated kV CBCT system, the authors have adapted a previously developed equivalent source model consisting of an equivalent spectrum module that takes into account intrinsic filtration and an equivalent filter module characterizing the added bowtie filtration. An equivalent spectrum was generated for an 80, 100, and 125 kVp beam with beam energy characterized by half-value layer measurements. An equivalent filter description was generated from bowtie profile measurements for both the full- and half-bowtie. Equivalent source models for each combination of equivalent spectrum and filter were incorporated into the Monte Carlo software package MCNPX. Monte Carlo simulations were then validated against in-phantom measurements for both the radiographic and CBCT mode of operation of the kV CBCT system. Radiographic and CBCT imaging dose was measured for a variety of protocols at various locations within a body (32 cm in diameter) and head (16 cm in diameter) CTDI phantom. The in-phantom radiographic and CBCT dose was simulated at all measurement locations and converted to absolute dose using normalization factors calculated from air scan measurements and corresponding simulations. The simulated results were compared with the physical measurements and their discrepancies were assessed quantitatively.Strong agreement was observed between in-phantom simulations and measurements. For the radiographic protocols, simulations uniformly underestimated measurements by 0.54%-5.14% (mean difference = -3.07%, SD = 1.60%). For the CBCT protocols, simulations uniformly underestimated measurements by 1.35%-5.31% (mean difference = -3.42%, SD = 1.09%).This work demonstrates the feasibility of using a measurement-based kV CBCT source model to facilitate dose calculations with Monte Carlo methods for both the radiographic and CBCT mode of operation. While this initial work validates simulations against measurements for simple geometries, future work will involve utilizing the source model to investigate kV CBCT dosimetry with more complex anthropomorphic phantoms and patient specific models.

Project description:The proliferation of cone-beam CT (CBCT) has created interest in performance optimization, with x-ray scatter identified among the main limitations to image quality. CBCT often contends with elevated scatter, but the wide variety of imaging geometry in different CBCT configurations suggests that not all configurations are affected to the same extent. Graphics processing unit (GPU) accelerated Monte Carlo (MC) simulations are employed over a range of imaging geometries to elucidate the factors governing scatter characteristics, efficacy of antiscatter grids, guide system design, and augment development of scatter correction.A MC x-ray simulator implemented on GPU was accelerated by inclusion of variance reduction techniques (interaction splitting, forced scattering, and forced detection) and extended to include x-ray spectra and analytical models of antiscatter grids and flat-panel detectors. The simulator was applied to small animal (SA), musculoskeletal (MSK) extremity, otolaryngology (Head), breast, interventional C-arm, and on-board (kilovoltage) linear accelerator (Linac) imaging, with an axis-to-detector distance (ADD) of 5, 12, 22, 32, 60, and 50 cm, respectively. Each configuration was modeled with and without an antiscatter grid and with (i) an elliptical cylinder varying 70-280 mm in major axis; and (ii) digital murine and anthropomorphic models. The effects of scatter were evaluated in terms of the angular distribution of scatter incident upon the detector, scatter-to-primary ratio (SPR), artifact magnitude, contrast, contrast-to-noise ratio (CNR), and visual assessment.Variance reduction yielded improvements in MC simulation efficiency ranging from ∼17-fold (for SA CBCT) to ∼35-fold (for Head and C-arm), with the most significant acceleration due to interaction splitting (∼6 to ∼10-fold increase in efficiency). The benefit of a more extended geometry was evident by virtue of a larger air gap-e.g., for a 16 cm diameter object, the SPR reduced from 1.5 for ADD = 12 cm (MSK geometry) to 1.1 for ADD = 22 cm (Head) and to 0.5 for ADD = 60 cm (C-arm). Grid efficiency was higher for configurations with shorter air gap due to a broader angular distribution of scattered photons-e.g., scatter rejection factor ∼0.8 for MSK geometry versus ∼0.65 for C-arm. Grids reduced cupping for all configurations but had limited improvement on scatter-induced streaks and resulted in a loss of CNR for the SA, Breast, and C-arm. Relative contribution of forward-directed scatter increased with a grid (e.g., Rayleigh scatter fraction increasing from ∼0.15 without a grid to ∼0.25 with a grid for the MSK configuration), resulting in scatter distributions with greater spatial variation (the form of which depended on grid orientation).A fast MC simulator combining GPU acceleration with variance reduction provided a systematic examination of a range of CBCT configurations in relation to scatter, highlighting the magnitude and spatial uniformity of individual scatter components, illustrating tradeoffs in CNR and artifacts and identifying the system geometries for which grids are more beneficial (e.g., MSK) from those in which an extended geometry is the better defense (e.g., C-arm head imaging). Compact geometries with an antiscatter grid challenge assumptions of slowly varying scatter distributions due to increased contribution of Rayleigh scatter.

Project description:Current CT dosimetry employing CTDI methodology has come under fire in recent years, partially in response to the increasing width of collimated x-ray fields in modern CT scanners. This study was conducted to provide a better understanding of the radiation dose distributions in CT.Monte Carlo simulations were used to evaluate radiation dose distributions along the z axis arising from CT imaging in cylindrical phantoms. Mathematical cylinders were simulated with compositions of water, polymethyl methacrylate (PMMA), and polyethylene. Cylinder diameters from 10 to 50 cm were studied. X-ray spectra typical of several CT manufacturers (80, 100, 120, and 140 kVp) were used. In addition to no bow tie filter, the head and body bow tie filters from modern General Electric and Siemens CT scanners were evaluated. Each cylinder was divided into three concentric regions of equal volume such that the energy deposited is proportional to dose for each region. Two additional dose assessment regions, central and edge locations 10 mm in diameter, were included for comparisons to CTDI100 measurements. Dose spread functions (DSFs) were computed for a wide number of imaging parameters.DSFs generally exhibit a biexponential falloff from the z=0 position. For a very narrow primary beam input (<< 1 mm), DSFs demonstrated significant low amplitude long range scatter dose tails. For body imaging conditions (30 cm diameter in water), the DSF at the center showed 160 mm at full width at tenth maximum (FWTM), while at the edge the FWTM was approximately 80 mm. Polyethylene phantoms exhibited wider DSFs than PMMA or water, as did higher tube voltages in any material. The FWTM were 80, 180, and 250 mm for 10, 30, and 50 cm phantom diameters, respectively, at the center in water at 120 kVp with a typical body bow tie filter. Scatter to primary dose ratios (SPRs) increased with phantom diameter from 4 at the center (1 cm diameter) for a 16 cm diameter cylinder to approximately 12.5 for a 32 cm diameter cylinder. The SPRs increased dramatically at the center of the phantom compared to the edge. For the three equal area regions, the edge to center SPRs for a 32 cm diameter phantom were approximately 1.8, 3.5, and 6.3, respectively.DSFs demonstrate low amplitude long ranging tails which reach considerable distances in cylindrical phantoms. The buildup that results from these long-ranged tails increases at the center of the field (at z=0) with increasing scan length. The DSF distributions lend a better understanding of the trends in CT dose deposition over a range of relevant imaging parameters. The DSFs as well as other related data are available to interested parties using EPAPS at http://www.aip.org/pubservs/epaps.html.

Dataset Information

Computed tomography dose index and dose length product for cone-beam CT: Monte Carlo simulations.

Publications

Computed tomography dose index and dose length product for cone-beam CT: Monte Carlo simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets