Consistency in contouring of organs at risk by artificial intelligence vs oncologists in head and neck cancer patients.

BACKGROUND: In the Danish Head and Neck Cancer Group (DAHANCA) 35 trial, patients are selected for proton treatment based on simulated reductions of Normal Tissue Complication Probability (NTCP) for proton compared to photon treatment at the referring departments. After inclusion in the trial, immobilization, scanning, contouring and planning are repeated at the national proton centre. The new contours could result in reduced expected NTCP gain of the proton plan, resulting in a loss of validity in the selection process. The present study evaluates if contour consistency can be improved by having access to AI (Artificial Intelligence) based contours. MATERIALS AND METHODS: The 63 patients in the DAHANCA 35 pilot trial had a CT from the local DAHANCA centre and one from the proton centre. A nationally validated convolutional neural network, based on nnU-Net, was used to contour OARs on both scans for each patient. Using deformable image registration, local AI and oncologist contours were transferred to the proton centre scans for comparison. Consistency was calculated with the Dice Similarity Coefficient (DSC) and Mean Surface Distance (MSD), comparing contours from AI to AI and oncologist to oncologist, respectively. Two NTCP models were applied to calculate NTCP for xerostomia and dysphagia. RESULTS: The AI contours showed significantly better consistency than the contours by oncologists. The median and interquartile range of DSC was 0.85 [0.78 - 0.90] and 0.68 [0.51 - 0.80] for AI and oncologist contours, respectively. The median and interquartile range of MSD was 0.9 mm [0.7 - 1.1] mm and 1.9 mm [1.5 - 2.6] mm for AI and oncologist contours, respectively. There was no significant difference in ΔNTCP. CONCLUSIONS: The study showed that OAR contours made by the AI algorithm were more consistent than those made by oncologists. No significant impact on the ΔNTCP calculations could be discerned.

High-resolution three-dimensional dosimetry in clinically relevant volumes utilizing optically stimulated luminescence.

BACKGROUND: The continued development of new radiotherapy techniques requires dosimetry systems that satisfy increasingly rigorous requirements, such as high sensitivity, wide dose range, and high spatial resolution. An emerging requirement is the ability to read out doses in three dimensions (3D) with high precision and spatial resolution. A few dosimetry systems with 3D capabilities are available, but their application in a clinical workflow is limited for various reasons, primarily originating from their chemical nature. The search for a 3D dosimetry system with potential for clinical implementation is thus ongoing. PURPOSE: To demonstrate the capabilities of a novel optically-stimulated-luminescence (OSL)-based 3D dosimetry system capable of measuring radiation doses in clinically relevant volumes. METHODS: A laser-based readout system was used to measure dose distributions delivered by both photons and protons, utilizing the OSL from a 50 × 50 × 50 $50\times 50\times 50$  mm 3 $^3$ YSO:Ce crystal. A homogeneous treatment plan consisting of two opposing photon fields was used to establish an inhomogeneity correction map of the crystal response and demonstrated the accuracy and precision of the system. The crystal was additionally irradiated with a photon treatment plan consisting of three overlapping 10 × 10 $10\times 10$  mm 2 $^2$ fields delivered from different angles, and a proton treatment plan consisting of four pencil beams with energies 90 MeV ( × 2 $\times 2$ ), 115 MeV, and 140 MeV. The system abilities were quantified by comparing the 3D-resolved measurements to Monte Carlo simulations. RESULTS: The dose map reproducibility of the system was found to be within 2% including both statistical and systematic errors. The measurements yielded integrated doses from a volume of 50 × 50 × 40 $50\times 50\times 40$  mm 3 $^3$ with voxel volumes of just 0.28 × 0.28 × 0.50 $0.28\times 0.28\times 0.50$  mm 3 $^3$ . An excellent agreement between the 3D-resolved measurements and the simulations was found for both photon- and proton-irradiation. CONCLUSIONS: The capabilities of the devised system for measuring clinically relevant fields of photons and proton pencil beams within a clinically relevant volume were demonstrated. The system poses as a promising candidate for clinical applications, and enables future research in the field of OSL-based tissue-equivalent 3D dosimetry.

Evaluation of decentralised model-based selection of head and neck cancer patients for a proton treatment study. DAHANCA 35.

INTRODUCTION: Proton treatment can potentially spare patients with H&N cancer for substantial treatment-related toxicities. The current study investigated the reproducibility of a decentralised model-based selection of patients for a proton treatment study when the selection plans were compared to the clinical treatment plans performed at the proton centre. METHODS: Sixty-three patients were selected for proton treatment in the six Danish Head and Neck Cancer (DAHANCA) centres. The patients were selected based on normal tissue complication probability (NTCP) estimated from local photon and proton treatment plans, which showed a ΔNTCP greater than 5%-point for either grade 2 + dysphagia or grade 2 + xerostomia at six months. The selection plans were compared to the clinical treatment plans performed at the proton centre. RESULTS: Of the 63 patients, 49 and 25 were selected based on an estimated benefit in risk of dysphagia and xerostomia, respectively. Eleven patients had a potential gain in both toxicities. The mean ΔNTCP changed from the local selection plan comparison to the clinical comparison from 6.9 to 5.3 %-points (p = 0.01) and 7.3 to 4.9 %-points (p = 0.03) for dysphagia and xerostomia, respectively. Volume differences in both CTV and OAR could add to the loss in ΔNTCP. 61 of the 63 clinical plans had a positive ΔNTCP, and 38 had a ΔNTCP of 5%-points for at least one of the two endpoints. CONCLUSION: A local treatment plan comparison can be used to select candidates for proton treatment. The local comparative proton plan overestimates the potential benefit of the clinical proton plan. Continuous quality assurance of the delineation procedures and planning is crucial in the subsequent randomised clinical trial setting.

Dosimetric verification of complex radiotherapy with a 3D optically based dosimetry system: dose painting and target tracking.

BACKGROUND: The increasing complexity of radiotherapy (RT) has motivated research into three-dimensional (3D) dosimetry. In this study we investigate the use of 3D dosimetry with polymerizing gels and optical computed tomography (optical CT) as a verification tool for complex RT: dose painting and target tracking. MATERIALS AND METHODS: For the dose painting studies, two dosimeters were irradiated with a seven-field intensity modulated radiotherapy (IMRT) plan with and without dose prescription based on a hypoxia image dataset of a head and neck patient. In the tracking experiments, two dosimeters were irradiated with a volumetric modulated arc therapy (VMAT) plan with and without clinically measured prostate motion and a third with both motion and target tracking. To assess the performance, 3D gamma analyses were performed between measured and calculated stationary dose distributions. RESULTS: Gamma pass-rates of 95.3% and 97.3% were achieved for the standard and dose-painted IMRT plans. Gamma pass-rates of 91.4% and 54.4% were obtained for the stationary and moving dosimeter, respectively, while tracking increased the pass-rate for the moving dosimeter to 90.4%. CONCLUSIONS: This study has shown that the 3D dosimetry system can reproduce and thus verify complex dose distributions, also when influenced by motion.

Technical note: Temporal and thermal stability of optical response for silicone-based 3D radiochromic dosimeters.

BACKGROUND: Radiochromic silicone-based dosimeters are flexible 3D dosimeters, which at appropriate concentration of leucomalachite green (LMG) and curing agent are dose-rate independent for clinical photon beams. However, their dose response is based on chemical processes that can be influenced by temporal and thermal conditions, impacting measurement stability. PURPOSE: The aim of this study was to investigate the temporal stability of the dose response of radiochromic dosimeters for different curing times and post-irradiation storage temperatures. METHODS: Six cylindrical dosimeters (5 cm diameter, 5 cm length) were produced in a single batch and separated into two groups that were irradiated 72 and 118 h after production. The same photon plan, consisting of two 10 × 1.6 cm2 opposing fields, was delivered to all dosimeters. After irradiation, the dosimeters were separated into three groups, stored at 5°C, 15°C, and 20°C, and read out for five consecutive days. RESULTS: Storage temperature influenced the measurement stability, and changes in the optical response with time differed between irradiated and non-irradiated parts of the dosimeters. The relative change between signal and background was greater than 10% for all measurements performed 24 h or more after irradiation, except for dosimeters stored at 5°C, which changed by 2%-5% after 24 h. The dosimeter temporal stability was not influenced by curing time. CONCLUSIONS: For room temperature storage (15°C and 20°C), readout should take place as soon as possible after irradiation since the background color increased rapidly for both curing times (72 and 118 h), whereas the dosimeters are stored at 5°C, readout can be performed up to 24 h after.

Temperature and temporal dependence of the optical response for a radiochromic dosimeter.

PURPOSE: Both temporal and thermal dependencies of the dose response have been observed in radiochromic dosimeters. As these dependencies may be influenced by the dose level, the present study investigates the temperature dependence during irradiation and the temporal change of the optical response following irradiation of radiochromic dosimeters at a range of doses. METHODS: Cuvette samples of the PRESAGE™ radiochromic dosimeter were irradiated within a dose range of 0-10 Gy at irradiation temperatures within 5-35 °C and postirradiation storage within 6-30 °C. The optical response due to irradiation was measured using a standard spectrophotometer and the data were analyzed in terms of thermal and temporal change. RESULTS: The initial dose response was linear over the applied dose range independent of irradiation temperature. However, the optical response to a specific dose increased exponentially with irradiation temperature corresponding to an activation energy of 0.114 ± 0.007 eV. The temporal change in dose response after irradiation consisted of an offset, an auto-oxidation rate with activation energy 0.84 ± 0.03 eV, and an initial exponential increase in optical response (1.6 ± 0.2 eV) followed by an exponential decrease in optical response (0.98 ± 0.08 eV). These contributions depended on both storage temperature and the dose given, leading to a nonlinear dose response with time at low storage temperatures and a high auto-oxidation rate at high storage temperatures. CONCLUSIONS: Thermal equilibration is important to the radiochromic dosimeter investigated due to an exponential change in dose response with irradiation temperature and a considerable postirradiation temporal change in response. For the dosimeter version investigated in this study, a compromise in storage temperature has to be made between increasing the nonlinearity of the dose response with time and inducing a high auto-oxidation rate.

Towards range-guidance in proton therapy to detect organ motion-induced dose degradations.

Introduction.Internal organ motion and deformations may cause dose degradations in proton therapy (PT) that are challenging to resolve using conventional image-guidance strategies. This study aimed to investigate the potential ofrange guidanceusing water-equivalent path length (WEPL) calculations to detect dose degradations occurring in PT.Materials and methods. Proton ranges were estimated using WEPL calculations. Field-specific isodose surfaces in the planning CT (pCT), from robustly optimised five-field proton plans (opposing lateral and three posterior/posterior oblique beams) for locally advanced prostate cancer patients, were used as starting points. WEPLs to each point on the field-specific isodoses in the pCT were calculated. The corresponding range for each point was found in the repeat CTs (rCTs). The spatial agreement between the resulting surfaces in the rCTs (hereafter referred to as iso-WEPLs) and the isodoses re-calculated in rCTs was evaluated for different dose levels and Hausdorff thresholds (2-5 mm). Finally, the sensitivity and specificity of detecting target dose degradation (V95% < 95%) using spatial agreement measures between the iso-WEPLs and isodoses in the pCT was evaluated.Results. The spatial agreement between the iso-WEPLs and isodoses in the rCTs depended on the Hausdorff threshold. The agreement was 65%-88% for a 2 mm threshold, 83%-96% for 3 mm, 90%-99% for 4 mm, and 94%-99% for 5 mm, across all fields and isodose levels. Minor differences were observed between the different isodose levels investigated. Target dose degradations were detected with 82%-100% sensitivity and 75%-80% specificity using a 2 mm Hausdorff threshold for the lateral fields.Conclusion. Iso-WEPLs were comparable to isodoses re-calculated in the rCTs. The proposed strategy could detect target dose degradations occurring in the rCTs and could be an alternative to a fully-fledged dose re-calculation to detect anatomical variations severely influencing the proton range.

Temperature dependence of the dose response for a solid-state radiochromic dosimeter during irradiation and storage.

PURPOSE: The dose response of radiochromic dosimeters is based on radiation-induced chemical reactions and is thus likely to be thermally influenced. In this study we have therefore investigated the temperature dependence of the dose response for such dosimeters, regarding both irradiation and storage conditions. METHODS: Dosimeter samples in cuvettes were irradiated to 5 Gy. The temperature for the different cuvettes during irradiation and post-irradiation storage was varied in the range of 3-30 degrees C in order to quantify the temperature dependence of the dosimeter response. The optical properties of the dosimeter samples were measured using a spectrophotometer before irradiation as well as at several times after irradiation to quantify the temporal variation of dose response (expressed as the optical density change induced by irradiation) as a function of storage temperature. RESULTS: The measurements show considerable temperature dependencies of dose response both during irradiation and storage. Fit to an Arrhenius equation revealed an activation energy of 1.4 +/- 0.2 eV for the variation in irradiation temperature, indicating a contribution from a thermally activated process. Variation in dose response at different storage temperatures showed an exponential increase with time followed by a decrease in optical density. Exponential Arrhenius fits to rate constants gave activation energies of 1.7 +/- 0.2 eV for the increase in dose response and 2.3 +/- 0.5 eV for the subsequent decrease, in this case dominated by thermally activated processes. CONCLUSIONS: Due to the exponential dependencies, stabilization of the dosimeter during irradiation at low temperatures (e.g., 5 degrees C) is preferable in clinical use to optimize the accuracy of the dose response. In addition, a low storage temperature is recommended in order to minimize the post-irradiation temporal change in dose response and thereby increase the post-irradiation stability of the dosimeter. The measurements in this study show that if the observed temperature and temporal dependencies are not considered, this could potentially deteriorate the accuracy of the dosimeter.

Characterization of the optical properties and stability of Presage™ following irradiation with photons and carbon ions.

BACKGROUND: The on-going development of both intensity-modulated radiotherapy (IMRT), including the more recent intensity-modulated arc therapy, as well as particle beam therapy, has created a clear need for accurate verification of dose distributions in three dimensions (3D). Presage™ is a new 3D dosimetry material that exhibits a radiochromic response when exposed to ionizing radiation. In this study we have 1) developed an improved optical set-up for measurements of changes in OD of Presage™ point dosimeters, 2) investigated the dose response of Presage™ for photons and carbon ions in the therapy range, 3) investigated the dose response of Presage™ for photons in the kGy range and 4) investigated the fading (i.e. bleaching) of Presage™ postirradiation. MATERIALS AND METHODS: Presage™ was examined in 1 × 1 × 4.5 cm(3) optical cuvettes; a cuvette holder assured accurate repositioning, and the optical setup included a reference detector to take into account laser intensity fluctuations. The cuvettes were measured pre- and postirradiation for a two week period. RESULTS: A linear response was observed between dose and optical response between 0 Gy and 100 Gy for γ-radiation from Co-60 and for carbon ions (both plateau and SOBP) from 0 to 20 Gy. The dosimeter was found to have a saturation dose of approximately 100 Gy for photons. A linear energy transfer (LET) effect was not observed in the dose response of different LET radiation. The postirradiation change in optical fading was found to be 0.5% ΔOD/day. CONCLUSIONS: Our study shows that Presage™ remains a dosimeter of interest for radiation therapy with other particles as well as photons in the therapy dose range.

Empirical quenching correction in radiochromic silicone-based three-dimensional dosimetry of spot-scanning proton therapy.

BACKGROUND AND PURPOSE: Three-dimensional dosimetry of proton therapy (PT) with chemical dosimeters is challenged by signal quenching, which is a lower dose-response in regions with high ionization density due to high linear-energy-transfer (LET) and dose-rate. This study aimed to assess the viability of an empirical correction model for 3D radiochromic silicone-based dosimeters irradiated with spot-scanning PT, by parametrizing its LET and dose-rate dependency. MATERIALS AND METHODS: Ten cylindrical radiochromic dosimeters (Ø50 and Ø75 mm) were produced in-house, and irradiated with different spot-scanning proton beam configurations and machine-set dose rates ranging from 56 to 145 Gy/min. Beams with incident energies of 75, 95 and 120 MeV, a spread-out Bragg peak and a plan optimized to an irregular target volume were included. Five of the dosimeters, irradiated with 120 MeV beams, were used to estimate the quenching correction factors. Monte Carlo simulations were used to obtain dose and dose-averaged-LET (LETd) maps. Additionally, a local dose-rate map was estimated, using the simulated dose maps and the machine-set dose-rate information retrieved from the irradiation log-files. Finally, the correction factor was estimated as a function of LETd and local dose-rate and tested on the different fields. RESULTS: Gamma-pass-rates of the corrected measurements were >94% using a 3%-3 mm gamma analysis and >88% using 2%-2 mm, with a dose deviation of <5.6 ± 1.8%. Larger dosimeters showed a 20% systematic increase in dose-response, but the same quenching in signal when compared to the smaller dosimeters. CONCLUSION: The quenching correction model was valid for different dosimeter sizes to obtain relative dosimetric maps of complex dose distributions in PT.

Chemically tuned linear energy transfer dependent quenching in a deformable, radiochromic 3D dosimeter.

Most solid-state detectors, including 3D dosimeters, show lower signal in the Bragg peak than expected, a process termed quenching. The purpose of this study was to investigate how variation in chemical composition of a recently developed radiochromic, silicone-based 3D dosimeter influences the observed quenching in proton beams. The dependency of dose response on linear energy transfer, as calculated through Monte Carlo simulations of the dosimeter, was investigated in 60 MeV proton beams. We found that the amount of quenching varied with the chemical composition: peak-to-plateau ratios (1 cm into the plateau) ranged from 2.2 to 3.4, compared to 4.3 using an ionization chamber. The dose response, and thereby the quenching, was predominantly influenced by the curing agent concentration, which determined the dosimeter's deformation properties. The dose response was found to be linear at all depths. All chemical compositions of the dosimeter showed dose-rate dependency; however this was not dependent on the linear energy transfer. Track-structure theory was used to explain the observed quenching effects. In conclusion, this study shows that the silicone-based dosimeter has potential for use in measuring 3D-dose-distributions from proton beams.

Technical Note: Improving proton stopping power ratio determination for a deformable silicone-based 3D dosimeter using dual energy CT.

PURPOSE: The aim of this study was to investigate whether the stopping power ratio (SPR) of a deformable, silicone-based 3D dosimeter could be determined more accurately using dual energy (DE) CT compared to using conventional methods based on single energy (SE) CT. The use of SECT combined with the stoichiometric calibration method was therefore compared to DECT-based determination. METHODS: The SPR of the dosimeter was estimated based on its Hounsfield units (HUs) in both a SECT image and a DECT image set. The stoichiometric calibration method was used for converting the HU in the SECT image to a SPR value for the dosimeter while two published SPR calibration methods for dual energy were applied on the DECT images. Finally, the SPR of the dosimeter was measured in a 60 MeV proton by quantifying the range difference with and without the dosimeter in the beam path. RESULTS: The SPR determined from SECT and the stoichiometric method was 1.10, compared to 1.01 with both DECT calibration methods. The measured SPR for the dosimeter material was 0.97. CONCLUSIONS: The SPR of the dosimeter was overestimated by 13% using the stoichiometric method and by 3% when using DECT. If the stoichiometric method should be applied for the dosimeter, the HU of the dosimeter must be manually changed in the treatment planning system in order to give a correct SPR estimate. Using a wrong SPR value will cause differences between the calculated and the delivered treatment plans.