Monte Carlo evaluation of scatter mitigation strategies in cone-beam CT.
Ontology highlight
ABSTRACT:
SUBMITTER: Lazos D
PROVIDER: S-EPMC2962662 | biostudies-other | 2010 Oct
REPOSITORIES: biostudies-other
Similar Datasets
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.
| S-EPMC3651212 | biostudies-other
Project description: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.
| S-EPMC5718669 | biostudies-literature
Project description:Background and purposeThe scatter induced image quality degradation of cone-beam computed tomography (CBCT) prevents more advanced applications in radiotherapy. We evaluated the dose calculation accuracy on CBCT of various disease sites using different scatter mitigation strategies.Materials and methodsCBCT scans of two patient cohorts (C1, C2) were reconstructed using a uniform (USC) and an iterative scatter correction (ISC) method, combined with an anti-scatter grid (ASG). Head and neck (H&N), lung, pelvic region, and prostate patients were included. To achieve a high accuracy Hounsfield unit and physical density calibrations were performed. The dose distributions of the original treatment plans were analyzed with the γ evaluation method using criteria of 1%/2 mm using the planning CT as the reference. The investigated parameters were the mean γ (γmean), the points in agreement (Pγ≤1) and the 99th percentile (γ1%).ResultsSignificant differences between USC and ISC in C1 were found for the lung and prostate, where the latter using the ISC produced the best results with medians of 0.38, 98%, and 1.1 for γmean, Pγ≤1 and γ1%, respectively. For C2 the ISC with ASG showed an improvement for all imaging sites. The lung demonstrated the largest relative increase in accuracy with improvements between 48% and 54% for the medians of γmean, Pγ≤1 and γ1%.ConclusionsThe introduced method demonstrated high dosimetric accuracy for H&N, prostate and pelvic region if an ASG is applied. A significantly lower accuracy was seen for lung. The ISC yielded a higher robustness against scatter variations than the USC.
| S-EPMC7807872 | biostudies-literature
Project description:PURPOSE:In 90 Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative 90 Y imaging by either SPECT or PET is challenging. In 90 Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate 90 Y bremsstrahlung SPECT/CT imaging. METHODS:Since a fully Monte Carlo (MC) approach to 90 Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration. RESULTS:There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively. CONCLUSIONS:Quantitative accuracy of 90 Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.
| S-EPMC5734647 | biostudies-literature
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.
| S-EPMC3808443 | biostudies-other
Project description:The electron Monte Carlo (eMC) dose calculation algorithm of the Eclipse treatment planning system is based heavily upon Monte Carlo simulation of the linac head and modeling of the linac beam characteristics with minimal measurement of beam data. Commissioning of the eMC algorithm on multiple identical linacs provided a unique opportunity to systematically evaluate the algorithm with actual measurements of clinically relevant beam and dose parameters. In this study, measured and eMC calculated dose distributions were compared both along and perpendicular to electron beam direction for electron energy/applicator/depth combination using measurement data from four Varian 21EX CLINAC linear accelerator (Varian Medical System, Palo Alto, CA). Cutout factors for sizes down to 3 x 3 cm were also compared. Comparisons between the measurement and the eMC calculated values show that the R90, R80, R50, and R10 values mostly agree within 3 mm. Measure and Calculated bremsstrahlung dose Dx correlates well statistically although eMC calculated Dx values are consistently smaller than the measured, with maximum discrepancy of 1% for the 20 MeV electron beams. Surface dose agrees mostly within 2%. Field width and penumbra agree mostly within 3mm. Calculation grid size is found to have a significant effect on the dose calculation. A grid size of 5 mm can produce erroneous dose distributions. Using a grid size of 2.5 mm and a 3% accuracy specified for the eMC to stop calculation iteration, the absolute output agrees with measurements within 3% for field sizes of 5 x 5 cm or larger. For cutout of 3 x 3 cm, however, the output disagreement can reach 8%. Our result indicate that eMC algorithm in Eclipse provides acceptable agreement with measurement data for most clinical situations. Calculation grid size of 2.5 mm or smaller is recommended.
| S-EPMC5722292 | biostudies-literature
Project description:Conventional mammography can suffer from poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter slot scan imaging is an imaging technique which provides additional information and is compatible with conventional mammography. A Monte Carlo simulation of coherent scatter slot scan imaging was performed to assess its performance and provide system optimization. Coherent scatter could be exploited using a system similar to conventional slot scan mammography system with antiscatter grids tilted at the characteristic angle of cancerous tissues. System optimization was performed across several parameters, including source voltage, tilt angle, grid distances, grid ratio, and shielding geometry. The simulated carcinomas were detectable for tumors as small as 5 mm in diameter, so coherent scatter analysis using a wide-slot setup could be promising as an enhancement for screening mammography. Employing coherent scatter information simultaneously with conventional mammography could yield a conventional high spatial resolution image with additional coherent scatter information.
| S-EPMC4999523 | biostudies-other
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.
| S-EPMC7108699 | biostudies-literature
Project description:OBJECTIVES: The aim of this study was to analyse the position and relationship of juxta-apical radiolucency (JAR) to the mandibular canal and buccal and/or lingual cortical plates using cone beam CT (CBCT). METHODS: A retrospective study was carried out to analyse the JAR on CBCT for 27 patients. These findings were compared with 27 age- and sex-matched patients without the presence of JAR, which acted as the control group. The CBCT images were analysed according to a checklist, to evaluate the position of the JAR and its relationship to the mandibular canal. Then, any thinning or perforation of either the buccal or lingual cortical plate due to JAR was noted, and a classification to quantify the thinning of cortical plates was proposed. The findings in the two groups were analysed using a paired comparison by McNemar test. RESULTS: A statistical increased thinning of cortical plates was seen in the JAR group compared with the control group, and most of the cases were in the J3 group. None of the patients in either the JAR or the control group showed perforation of the buccal and/or lingual cortical plate on CBCT images. CONCLUSIONS: A classification to quantify the thinning of cortical plates was proposed, which may be used for objective evaluation of the thinning of the cortical plates in future studies. The present study gives an insight into the relationship of the juxta-apical area with the mandibular canal and cortical plates in the mandible using CBCT.
| S-EPMC4082266 | biostudies-other
Project description:Background:Proton pencil beam (PB) dose calculation algorithms have limited accuracy within heterogeneous tissues of lung cancer patients, which may be addressed by modern commercial Monte Carlo (MC) algorithms. We investigated clinical pencil beam scanning (PBS) dose differences between PB and MC-based treatment planning for lung cancer patients. Methods:With IRB approval, a comparative dosimetric analysis between RayStation MC and PB dose engines was performed on ten patient plans. PBS gantry plans were generated using single-field optimization technique to maintain target coverage under range and setup uncertainties. Dose differences between PB-optimized (PBopt), MC-recalculated (MCrecalc), and MC-optimized (MCopt) plans were recorded for the following region-of-interest metrics: clinical target volume (CTV) V95, CTV homogeneity index (HI), total lung V20, total lung VRX (relative lung volume receiving prescribed dose or higher), and global maximum dose. The impact of PB-based and MC-based planning on robustness to systematic perturbation of range (±3% density) and setup (±3 mm isotropic) was assessed. Pairwise differences in dose parameters were evaluated through non-parametric Friedman and Wilcoxon sign-rank testing. Results:In this ten-patient sample, CTV V95 decreased significantly from 99-100% for PBopt to 77-94% for MCrecalc and recovered to 99-100% for MCopt (P<10-5). The median CTV HI (D95/D5) decreased from 0.98 for PBopt to 0.91 for MCrecalc and increased to 0.95 for MCopt (P<10-3). CTV D95 robustness to range and setup errors improved under MCopt (?D95 =-1%) compared to MCrecalc (?D95 =-6%, P=0.006). No changes in lung dosimetry were observed for large volumes receiving low to intermediate doses (e.g., V20), while differences between PB-based and MC-based planning were noted for small volumes receiving high doses (e.g., VRX). Global maximum patient dose increased from 106% for PBopt to 109% for MCrecalc and 112% for MCopt (P<10-3). Conclusions:MC dosimetry revealed a reduction in target dose coverage under PB-based planning that was regained under MC-based planning along with improved plan robustness. MC-based optimization and dose calculation should be integrated into clinical planning workflows of lung cancer patients receiving actively scanned proton therapy.
| S-EPMC5960654 | biostudies-literature