Dataset Information

Comparison of Monte Carlo methods for fluorescence molecular tomography-computational efficiency.

ABSTRACT: The Monte Carlo method is an accurate model for time-resolved quantitative fluorescence tomography. However, this method suffers from low computational efficiency due to the large number of photons required for reliable statistics. This paper presents a comparison study on the computational efficiency of three Monte Carlo-based methods for time-domain fluorescence molecular tomography.The methods investigated to generate time-gated Jacobians were the perturbation Monte Carlo (pMC) method, the adjoint Monte Carlo (aMC) method and the mid-way Monte Carlo (mMC) method. The effects of the different parameters that affect the computation time and statistics reliability were evaluated. Also, the methods were applied to a set of experimental data for tomographic application.In silico results establish that, the investigated parameters affect the computational time for the three methods differently (linearly, quadratically, or not significantly). Moreover, the noise level of the Jacobian varies when these parameters change. The experimental results in preclinical settings demonstrates the feasibility of using both aMC and pMC methods for time-resolved whole body studies in small animals within a few hours.Among the three Monte Carlo methods, the mMC method is a computationally prohibitive technique that is not well suited for time-domain fluorescence tomography applications. The pMC method is advantageous over the aMC method when the early gates are employed and large number of detectors is present. Alternatively, the aMC method is the method of choice when a small number of source-detector pairs are used.

SUBMITTER: Chen J

PROVIDER: S-EPMC3298560 | biostudies-other | 2011 Oct

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation.The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique.A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 ± 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 ± 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria.The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

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.

Dataset Information

Comparison of Monte Carlo methods for fluorescence molecular tomography-computational efficiency.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets