Dataset Information

Monte Carlo study of the effects of system geometry and antiscatter grids on cone-beam CT scatter distributions.

ABSTRACT: 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.

SUBMITTER: Sisniega A

PROVIDER: S-EPMC3651212 | biostudies-other | 2013 May

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

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.

Project description:X-ray scatter is a significant barrier to image quality improvements in cone-beam computed tomography (CBCT). A moving blocker-based strategy was previously proposed to simultaneously estimate scatter and reconstruct the complete volume within the field of view (FOV) from a single CBCT scan. A blocker consisting of lead stripes is inserted between the X-ray source and the imaging object, and moves back and forth along the rotation axis during gantry rotation. While promising results were obtained in our previous studies, the geometric design and moving speed of the blocker were set empirically. The goal of this work is to optimize the geometry and speed of the moving block system.Performance of the blocker was examined through Monte Carlo (MC) simulation and experimental studies with various geometry designs and moving speeds. All hypothetical designs employed an anthropomorphic pelvic phantom. The scatter estimation accuracy was quantified by using lead stripes ranging from 5 to 100 pixels on the detector plane. An iterative reconstruction based on total variation minimization was used to reconstruct CBCT images from unblocked projection data after scatter correction. The reconstructed image was evaluated under various combinations of lead strip width and interspace (ranging from 10 to 60 pixels) and different moving speed (ranging from 1 to 30 pixels per projection).MC simulation showed that the scatter estimation error varied from 0.8% to 5.8%. Phantom experiment showed that CT number error in the reconstructed CBCT images varied from 13 to 35. Highest reconstruction accuracy was achieved when the strip width was 20 pixels and interspace was 60 pixels and the moving speed was 15 pixels per projection.Scatter estimation can be achieved in a large range of lead strip width and interspace combinations. The moving speed does not have a very strong effect on reconstruction result if it is above 5 pixels per projection. Geometry design of the blocker affected image reconstruction accuracy more. The optimal geometry of the blocker has a strip width of 20 pixels and an interspace three times the strip width, which means 25% detector is covered by the blocker, while the optimal moving speed is 15 pixels per projection.

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 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:X-ray scatter is a major detriment to image quality in cone-beam CT (CBCT). Existing geometries exhibit strong differences in scatter susceptibility with more compact geometries, e.g., dental or musculoskeletal, benefiting from antiscatter grids, whereas in more extended geometries, e.g., IGRT, grid use carries tradeoffs in image quality per unit dose. This work assesses the tradeoffs in dose and image quality for grids applied in the context of low-dose CBCT on a mobile C-arm for image-guided surgery.Studies were performed on a mobile C-arm equipped with a flat-panel detector for high-quality CBCT. Antiscatter grids of grid ratio (GR) 6:1-12:1, 40 lp∕cm, were tested in "body" surgery, i.e., spine, using protocols for bone and soft-tissue visibility in the thoracic and abdominal spine. Studies focused on grid orientation, CT number accuracy, image noise, and contrast-to-noise ratio (CNR) in quantitative phantoms at constant dose.There was no effect of grid orientation on possible gridline artifacts, given accurate angle-dependent gain calibration. Incorrect calibration was found to result in gridline shadows in the projection data that imparted high-frequency artifacts in 3D reconstructions. Increasing GR reduced errors in CT number from 31%, thorax, and 37%, abdomen, for gridless operation to 2% and 10%, respectively, with a 12:1 grid, while image noise increased by up to 70%. The CNR of high-contrast objects was largely unaffected by grids, but low-contrast soft-tissues suffered reduction in CNR, 2%-65%, across the investigated GR at constant dose.While grids improved CT number accuracy, soft-tissue CNR was reduced due to attenuation of primary radiation. CNR could be restored by increasing dose by factors of ~1.6-2.5 depending on GR, e.g., increase from 4.6 mGy for the thorax and 12.5 mGy for the abdomen without antiscatter grids to approximately 12 mGy and 30 mGy, respectively, with a high-GR grid. However, increasing the dose poses a significant impediment to repeat intraoperative CBCT and can cause the cumulative intraoperative dose to exceed that of a single diagnostic CT scan. This places the mobile C-arm in the category of extended CBCT geometries with sufficient air gap for which the tradeoffs between CNR and dose typically do not favor incorporation of an antiscatter grid.

Dataset Information

Monte Carlo study of the effects of system geometry and antiscatter grids on cone-beam CT scatter distributions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets