Guillaume Landry

Sensitivity of low energy brachytherapy Monte Carlo dose calculations to uncertainties in human tissue composition.

PURPOSE The objective of this work is to assess the sensitivity of Monte Carlo (MC) dose calculations to uncertainties in human tissue composition for a range of low photon energy brachytherapy sources:I125, P103d, C131s, and an electronic brachytherapy source (EBS). The low energy photons emitted by these sources make the dosimetry sensitive to variations in tissue atomic number due to the dominance of the photoelectric effect. This work reports dose to a small mass of water in medium Dw,m as opposed to dose to a small mass of medium in medium Dm,m. METHODS Mean adipose, mammary gland, and breast tissues (as uniform mixture of the aforementioned tissues) are investigated as well as compositions corresponding to one standard deviation from the mean. Prostate mean compositions from three different literature sources are also investigated. Three sets of MC simulations are performed with theGEANT4 code: (1) Dose calculations for idealized TG-43-like spherical geometries using point sources. Radial dose profiles obtained in different media are compared to assess the influence of compositional uncertainties. (2) Dose calculations for four clinical prostate LDR brachytherapy permanent seed implants using I125 seeds (Model 2301, Best Medical, Springfield, VA). The effect of varying the prostate composition in the planning target volume (PTV) is investigated by comparing PTV D90 values. (3) Dose calculations for four clinical breast LDR brachytherapy permanent seed implants using P103d seeds (Model 2335, Best Medical). The effects of varying the adipose/gland ratio in the PTV and of varying the elemental composition of adipose and gland within one standard deviation of the assumed mean composition are investigated by comparing PTV D90 values. For (2) and (3), the influence of using the mass density from CT scans instead of unit mass density is also assessed. RESULTS Results from simulation (1) show that variations in the mean compositions of tissues affect low energy brachytherapy dosimetry. Dose differences between mean and one standard deviation of the mean composition increasing with distance from the source are observed. It is established that theI125 and C131s sources are the least sensitive to variations in elemental compositions while P103d is most sensitive. The EBS falls in between and exhibits complex behavior due to significant spectral hardening. Results from simulation (2) show that two prostate compositions are dosimetrically equivalent to water while the third shows D90 differences of up to 4%. Results from simulation (3) show that breast is more sensitive than prostate with dose variations of up to 30% from water for 70% adipose/30% gland breast. The variability of the breast composition adds a ±10% dose variation. CONCLUSIONS Low energy brachytherapy dose distributions in tissue differ from water and are influenced by density, mean tissue composition, and patient-to-patient composition variations. The results support the use of a dose calculation algorithm accounting for heterogeneities such as MC. Since this work shows that variations in mean tissue compositions affect MC dosimetry and result in increased dose uncertainties, the authors conclude that imaging tools providing more accurate estimates of elemental compositions such as dual energy CT would be beneficial.

Extracting atomic numbers and electron densities from a dual source dual energy CT scanner: experiments and a simulation model.

BACKGROUND AND PURPOSE Dual energy CT (DECT) imaging can provide both the electron density ρ(e) and effective atomic number Z(eff), thus facilitating tissue type identification. This paper investigates the accuracy of a dual source DECT scanner by means of measurements and simulations. Previous simulation work suggested improved Monte Carlo dose calculation accuracy when compared to single energy CT for low energy photon brachytherapy, but lacked validation. As such, we aim to validate our DECT simulation model in this work. MATERIALS AND METHODS A cylindrical phantom containing tissue mimicking inserts was scanned with a second generation dual source scanner (SOMATOM Definition FLASH) to obtain Z(eff) and ρ(e). A model of the scanner was designed in ImaSim, a CT simulation program, and was used to simulate the experiment. RESULTS Accuracy of measured Z(eff) (labelled Z) was found to vary from -10% to 10% from low to high Z tissue substitutes while the accuracy on ρ(e) from DECT was about 2.5%. Our simulation reproduced the experiments within ±5% for both Z and ρ(e). CONCLUSIONS A clinical DECT scanner was able to extract Z and ρ(e) of tissue substitutes. Our simulation tool replicates the experiments within a reasonable accuracy.

Evaluation of a novel triple-channel radiochromic film analysis procedure using EBT2

A novel approach to read out radiochromic film was introduced recently by the manufacturer of GafChromic film. In this study, the performance of this triple-channel film dosimetry method was compared against the conventional single-red-channel film dosimetry procedure, with and without inclusion of a pre-irradiation (pre-IR) film scan, using EBT2 film and kilo- and megavoltage photon beams up to 10 Gy. When considering regions of interest averaged doses, the triple-channel method and both single-channel methods produced equivalent results. Absolute dose discrepancies between the triple-channel method, both single-channel methods and the treatment planning system calculated dose values, were no larger than 5 cGy for dose levels up to 2.2 Gy. Signal to noise in triple-channel dose images was found to be similar to signal to noise in single-channel dose images. The accuracy of resulting dose images from the triple- and single-channel methods with inclusion of pre-IR film scan was found to be similar. Results of a comparison of EBT2 data from a kilovoltage depth dose experiment to corresponding Monte Carlo depth dose data produced dose discrepancies of 9.5 ± 12 cGy and 7.6 ± 6 cGy for the single-channel method with inclusion of a pre-IR film scan and the triple-channel method, respectively. EBT2 showed to be energy sensitive at low kilovoltage energies with response differences of 11.9% and 15.6% in the red channel at 2 Gy between 50-225 kVp and 80-225 kVp photon spectra, respectively. We observed that the triple-channel method resulted in non-uniformity corrections of ±1% and consistency values of 0-3 cGy for the batches and dose levels studied. Results of this study indicate that the triple-channel radiochromic film read-out method performs at least as well as the single-channel method with inclusion of a pre-IR film scan, reduces film non-uniformity and saves time with elimination of a pre-IR film scan.

Simulation study on potential accuracy gains from dual energy CT tissue segmentation for low-energy brachytherapy Monte Carlo dose calculations

This work compares Monte Carlo (MC) dose calculations for (125)I and (103)Pd low-dose rate (LDR) brachytherapy sources performed in virtual phantoms containing a series of human soft tissues of interest for brachytherapy. The geometries are segmented (tissue type and density assignment) based on simulated single energy computed tomography (SECT) and dual energy (DECT) images, as well as the all-water TG-43 approach. Accuracy is evaluated by comparison to a reference MC dose calculation performed in the same phantoms, where each voxels material properties are assigned with exactly known values. The objective is to assess potential dose calculation accuracy gains from DECT. A CT imaging simulation package, ImaSim, is used to generate CT images of calibration and dose calculation phantoms at 80, 120, and 140 kVp. From the high and low energy images electron density ρ(e) and atomic number Z are obtained using a DECT algorithm. Following a correction derived from scans of the calibration phantom, accuracy on Z and ρ(e) of ±1% is obtained for all soft tissues with atomic number Z ∊ [6,8] except lung. GEANT4 MC dose calculations based on DECT segmentation agreed with the reference within ±4% for (103)Pd, the most sensitive source to tissue misassignments. SECT segmentation with three tissue bins as well as the TG-43 approach showed inferior accuracy with errors of up to 20%. Using seven tissue bins in our SECT segmentation brought errors within ±10% for (103)Pd. In general (125)I dose calculations showed higher accuracy than (103)Pd. Simulated image noise was found to decrease DECT accuracy by 3-4%. Our findings suggest that DECT-based segmentation yields improved accuracy when compared to SECT segmentation with seven tissue bins in LDR brachytherapy dose calculation for the specific case of our non-anthropomorphic phantom. The validity of our conclusions for clinical geometry as well as the importance of image noise in the tissue segmentation procedure deserves further experimental investigation.

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation. METHODS Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT. RESULTS The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy]. CONCLUSIONS In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

Dual energy CT in radiotherapy: Current applications and future outlook

Dual energy CT (DECT) scanners are nowadays available in many radiology departments. For radiotherapy purposes, new strategies using DECT imaging are investigated to optimize radiation treatment for multiple steps in the radiotherapy chain. This review describes how DECT based methods can be used for electron density estimation, effective atomic number decomposition and contrast material quantification. Clinical radiotherapy related applications for improved dose calculation accuracy of brachytherapy and proton therapy, metal artifact reduction techniques and normal tissue characterization are also summarized together with future perspectives on the use of DECT for radiotherapy purposes.

ALGEBRA: ALgorithm for the heterogeneous dosimetry based on GEANT4 for BRAchytherapy

Task group 43 (TG43)-based dosimetry algorithms are efficient for brachytherapy dose calculation in water. However, human tissues have chemical compositions and densities different than water. Moreover, the mutual shielding effect of seeds on each other (interseed attenuation) is neglected in the TG43-based dosimetry platforms. The scientific community has expressed the need for an accurate dosimetry platform in brachytherapy. The purpose of this paper is to present ALGEBRA, a Monte Carlo platform for dosimetry in brachytherapy which is sufficiently fast and accurate for clinical and research purposes. ALGEBRA is based on the GEANT4 Monte Carlo code and is capable of handling the DICOM RT standard to recreate a virtual model of the treated site. Here, the performance of ALGEBRA is presented for the special case of LDR brachytherapy in permanent prostate and breast seed implants. However, the algorithm is also capable of handling other treatments such as HDR brachytherapy.

A combined dose calculation and verification method for a small animal precision irradiator based on onboard imaging

PURPOSE Novel small animal precision microirradiators (micro-IR) are becoming available for preclinical use and are often equipped with onboard imaging (OBI) devices. We investigated the use of OBI as a means to infer the accuracy of the delivered treatment plan. METHODS Monte Carlo modeling of the micro-IR including an elliptical Gaussian electron beam incident on the x-ray tube was used to score dose and to continue photon transport to the plane of the OBI device. A model of the OBI detector response was used to generate simulated onboard images. Experimental OBI was performed at 225 kVp, gain∕offset and scatter-glare were corrected. Simulated and experimentally obtained onboard images of phantoms and a mouse specimen were compared for a range of photon beam sizes from 2.5 cm down to 0.1 cm. RESULTS Simulated OBI can be used in small animal radiotherapy to determine if a treatment plan was delivered according to the prescription within an uncertainty of 5% for beams as small as 4 mm in diameter. For collimated beams smaller than 4 mm, beam profile differences remain primarily in the penumbra region of the smallest beams, which may be tolerable for specific preclinical micro-IR investigations. CONCLUSIONS Comparing simulated to acquired OBI during small animal treatment radiotherapy represents a useful treatment delivery tool.PURPOSE Novel small animal precision microirradiators (micro-IR) are becoming available for preclinical use and are often equipped with onboard imaging (OBI) devices. We investigated the use of OBI as a means to infer the accuracy of the delivered treatment plan. METHODS Monte Carlo modeling of the micro-IR including an elliptical Gaussian electron beam incident on the x-ray tube was used to score dose and to continue photon transport to the plane of the OBI device. A model of the OBI detector response was used to generate simulated onboard images. Experimental OBI was performed at 225 kVp, gain/offset and scatter-glare were corrected. Simulated and experimentally obtained onboard images of phantoms and a mouse specimen were compared for a range of photon beam sizes from 2.5 cm down to 0.1 cm. RESULTS Simulated OBI can be used in small animal radiotherapy to determine if a treatment plan was delivered according to the prescription within an uncertainty of 5% for beams as small as 4 mm in diameter. For collimated beams smaller than 4 mm, beam profile differences remain primarily in the penumbra region of the smallest beams, which may be tolerable for specific preclinical micro-IR investigations. CONCLUSIONS Comparing simulated to acquired OBI during small animal treatment radiotherapy represents a useful treatment delivery tool.

Comparison of proton therapy treatment planning for head tumors with a pencil beam algorithm on dual and single energy CT images

PURPOSE Dual energy CT (DECT) has recently been proposed as an improvement over single energy CT (SECT) for stopping power ratio (SPR) estimation for proton therapy treatment planning (TP), thereby potentially reducing range uncertainties. Published literature investigated phantoms. This study aims at performing proton therapy TP on SECT and DECT head images of the same patients and at evaluating whether the reported improved DECT SPR accuracy translates into clinically relevant range shifts in clinical head treatment scenarios. METHODS Two phantoms were scanned at a last generation dual source DECT scanner at 90 and 150 kVp with Sn filtration. The first phantom (Gammex phantom) was used to calibrate the scanner in terms of SPR while the second served as evaluation (CIRS phantom). DECT images of five head trauma patients were used as surrogate cancer patient images for TP of proton therapy. Pencil beam algorithm based TP was performed on SECT and DECT images and the dose distributions corresponding to the optimized proton plans were calculated using a Monte Carlo (MC) simulation platform using the same patient geometry for both plans obtained from conversion of the 150 kVp images. Range shifts between the MC dose distributions from SECT and DECT plans were assessed using 2D range maps. RESULTS SPR root mean square errors (RMSEs) for the inserts of the Gammex phantom were 1.9%, 1.8%, and 1.2% for SECT phantom calibration (SECTphantom), SECT stoichiometric calibration (SECTstoichiometric), and DECT calibration, respectively. For the CIRS phantom, these were 3.6%, 1.6%, and 1.0%. When investigating patient anatomy, group median range differences of up to -1.4% were observed for head cases when comparing SECTstoichiometric with DECT. For this calibration the 25th and 75th percentiles varied from -2% to 0% across the five patients. The group median was found to be limited to 0.5% when using SECTphantom and the 25th and 75th percentiles varied from -1% to 2%. CONCLUSIONS Proton therapy TP using a pencil beam algorithm and DECT images was performed for the first time. Given that the DECT accuracy as evaluated by two phantoms was 1.2% and 1.0% RMSE, it is questionable whether the range differences reported here are significant.

Dual-energy CT quantitative imaging: a comparison study between twin-beam and dual-source CT scanners

Purpose: To assess image quality and to quantify the accuracy of relative electron densities (ρe) and effective atomic numbers (Zeff) for three dual‐energy computed tomography (DECT) scanners: a novel single‐source split‐filter (i.e., twin‐beam) and two dual‐source scanners. Methods: Measurements were made with a second generation dual‐source scanner at 80/140Sn kVp, a third‐generation twin‐beam single‐source scanner at 120 kVp with gold (Au) and tin (Sn) filters, and a third‐generation dual‐source scanner at 90/150Sn kVp. Three phantoms with tissue inserts were scanned and used for calibration and validation of parameterized methods to extract ρe and Zeff, whereas iodine and calcium inserts were used to quantify Contrast‐to‐Noise‐Ratio (CNR). Spatial resolution in tomographic images was also tested. Results: The third‐generation scanners have an image resolution of 6.2, ˜0.5 lp/cm higher than the second generation scanner. The twin‐beam scanner has low imaging contrast for iodine materials due to its limited spectral separation. The parameterization methods resulted in calibrations with low fit residuals for the dual‐source scanners, yielding values of ρe and Zeff close to the reference values (errors within 1.2% for ρe and 6.2% for Zeff for a dose of 20 mGy, excluding lung substitute tissues). The twin‐beam scanner presented overall higher errors (within 3.2% for ρe and 28% for Zeff, also excluding lung inserts) and also larger variations for uniform inserts. Conclusions: Spatial resolution is similar for the three scanners. The twin‐beam is able to derive ρe and Zeff, but with inferior accuracy compared to both dual‐source scanners.