Towards real-time EPID-based 3D in vivo dosimetry for IMRT with Deep Neural Networks: A feasibility study.

We investigate the potential of the Deep Dose Estimate (DDE) neural network to predict 3D dose distributions inside patients with Monte Carlo (MC) accuracy, based on transmitted EPID signals and patient CTs. The network was trained using as input patient CTs and first-order dose approximations (FOD). Accurate dose distributions (ADD) simulated with MC were given as training targets. 83 pelvic CTs were used to simulate ADDs and respective EPID signals for subfields of prostate IMRT plans (gantry at 0∘). FODs were produced as backprojections from the EPID signals. 581 ADD-FOD sets were produced and divided into training and test sets. An additional dataset simulated with gantry at 90∘ (lateral set) was used for evaluating the performance of the DDE at different beam directions. The quality of the FODs and DDE-predicted dose distributions (DDEP) with respect to ADDs, from the test and lateral sets, was evaluated with gamma analysis (3%,2 mm). The passing rates between FODs and ADDs were as low as 46%, while for DDEPs the passing rates were above 97% for the test set. Meaningful improvements were also observed for the lateral set. The high passing rates for DDEPs indicate that the DDE is able to convert FODs into ADDs. Moreover, the trained DDE predicts the dose inside a patient CT within 0.6 s/subfield (GPU), in contrast to 14 h needed for MC (CPU-cluster). 3D in vivo dose distributions due to clinical patient irradiation can be obtained within seconds, with MC-like accuracy, potentially paving the way towards real-time EPID-based in vivo dosimetry.

Measurement-based range evaluation for quality assurance of CBCT-based dose calculations in adaptive proton therapy.

PURPOSE: The implementation of volumetric in-room imaging for online adaptive radiotherapy makes extensive testing of this image data for treatment planning necessary. Especially for proton beams the higher sensitivity to stopping power properties of the tissue results in more stringent requirements. Current approaches mainly focus on recalculation of the plans on the new image data, lacking experimental verification, and ignoring the impact on the plan re-optimization process. The aim of this study was to use gel and film dosimetry coupled with a three-dimensional (3D) printed head phantom (based on the planning CT of the patient) for 3D range verification of intensity-corrected cone beam computed tomography (CBCT) image data for adaptive proton therapy. METHODS: Single field uniform dose pencil beam scanning proton plans were optimized for three different patients on the patients' planning CT (planCT) and the patients' intensity-corrected CBCT (scCBCT) for the same target volume using the same optimization constraints. The CBCTs were corrected on projection level using the planCT as a prior. The dose optimized on planCT and recalculated on scCBCT was compared in terms of proton range differences (80% distal fall-off, recalculation). Moreover, the dose distribution resulting from recalculation of the scCBCT-optimized plan on the planCT and the original planCT dose distribution were compared (simulation). Finally, the two plans of each patient were irradiated on the corresponding patient-specific 3D printed head phantom using gel dosimetry inserts for one patient and film dosimetry for all three patients. Range differences were extracted from the measured dose distributions. The measured and the simulated range differences were corrected for range differences originating from the initial plans and evaluated. RESULTS: The simulation approach showed high agreement with the standard recalculation approach. The median values of the range differences of these two methods agreed within 0.1 mm and the interquartile ranges (IQRs) within 0.3 mm for all three patients. The range differences of the film measurement were accurately matching with the simulation approach in the film plane. The median values of these range differences deviated less than 0.1 mm and the IQRs less than 0.4 mm. For the full 3D evaluation of the gel range differences, the median value and IQR matched those of the simulation approach within 0.7 and 0.5 mm, respectively. scCBCT- and planCT-based dose distributions were found to have a range agreement better than 3 mm (median and IQR) for all considered scenarios (recalculation, simulation, and measurement). CONCLUSIONS: The results of this initial study indicate that an online adaptive proton workflow based on scatter-corrected CBCT image data for head irradiations is feasible. The novel presented measurement- and simulation-based method was shown to be equivalent to the standard literature recalculation approach. Additionally, it has the capability to catch effects of image differences on the treatment plan optimization. This makes the measurement-based approach particularly interesting for quality assurance of CBCT-based online adaptive proton therapy. The observed uncertainties could be kept within those of the registration and positioning. The proposed validation could also be applied for other alternative in-room images, e.g. for MR-based pseudoCTs.

Real-time 4DMRI-based internal target volume definition for moving lung tumors.

PURPOSE: In photon radiotherapy, respiratory-induced target motion can be accounted for by internal target volumes (ITV) or mid-ventilation target volumes (midV) defined on the basis of four-dimensional computed tomography (4D-CT). Intrinsic limitations of these approaches can result in target volumes that are not representative for the gross tumor volume (GTV) motion over the course of treatment. To address these limitations, we propose a novel patient-specific ITV definition method based on real-time 4D magnetic resonance imaging (rt-4DMRI). METHODS: Three lung cancer patients underwent weekly rt-4DMRI scans. A total of 24 datasets were included in this retrospective study. The GTV was contoured on breath-hold MR images and propagated to all rt-4DMRI images by deformable image registration. Different targets were created for the first (reference) imaging sessions: ITVs encompassing all GTV positions over the complete (ITV 80 s ) or partial acquisition time ( ITV 10 s ), ITVs including only voxels with a GTV probability-of-presence (POP) of at least 5% ( ITV 5 % ) or 10% ( ITV 10 % ), and the mid-ventilation GTV position. Reference planning target volumes ( PTV r ) were created by adding margins around the ITVs and midV target volumes. The geometrical overlap of the PTV r with ITV n 5 % from the six to eight subsequent imaging sessions on days n was quantified in terms of the Dice similarity coefficient (DSC), sensitivity [SE: ( PTV r ∩ ITV n 5 % )/ ITV n 5 % ] and precision [PRE: ( PTV r ∩ ITV n 5 % )/ PTV r ] as surrogates for target coverage and normal tissue sparing. RESULTS: Patient-specific analysis yielded a high variance of the overlap values of PTV r 10 s , when different periods within the reference imaging session were sampled. The mid-ventilation-based PTVs were smaller than the ITV-based PTVs. While the SE was high for patients with small breathing pattern variations, changes of the median breathing amplitudes in different imaging sessions led to inferior SE values for the mid-ventilation PTV for one patient. In contrast, PTV r 5 % and PTV r 10 % showed higher SE values with a higher robustness against interfractional changes, at the cost of larger target volumes. CONCLUSIONS: The results indicate that rt-4DMRI could be valuable for the definition of target volumes based on the GTV POP to achieve a higher robustness against interfractional changes than feasible with today's 4D-CT-based target definition concepts.

Optimization of Phase Space files from clinical linear accelerators.

This work proposes a methodology to produce an optimized phase-space (PhSp) for the Elekta Synergy linac by tuning the energy and direction of particles inside the 6-MV Elekta Precise PhSp, provided by the International Atomic Energy Agency (IAEA), for Monte Carlo (MC) simulations. First, the energies of the particles emerging from the original PhSp were increased by different factors, producing new PhSps. Percentage depth dose (PDD) profiles were simulated and compared to measured data from a Synergy linac for 6-MV photon beam. This process was repeated until a minimum difference was reached. Particles' directions were then manipulated following identified correlations to lateral profiles, resulting in two distinct perturbation factors based on inline and crossline profiles. Both factors were merged into one single optimal factor. For energy optimization, an increase of 0.32 MeV applied to all particles inside the original PhSp, but to 0.511 MeV annihilation photons, provided the best results. The direction optimization factor was the combination of the individual factors for inline (0.605%) and crossline (0.051%). The agreement between measured and simulated profiles, when using the optimized PhSp, improved considerably in comparison to simulations performed with the original IAEA PhSp. For all fields and depths analyzed, the discrepancies for PDD, inline and crossline profiles dropped from 11.2%, 15.7% and 27.5% to under 1.4%, 4.7% and 13.2%, respectively. The optimized PhSp should not replace the full linac modelling, however it offers an alternative for MC dose calculations when neither geometric details nor validated IAEA PhSp are available to the user.

Evaluation of proton and photon dose distributions recalculated on 2D and 3D Unet-generated pseudoCTs from T1-weighted MR head scans.

Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.

Feasibility of 4DCBCT-based proton dose calculation: An ex vivo porcine lung phantom study.

Inter-fractional variations of breathing pattern and patient anatomy introduce dose uncertainties in proton therapy. One approach to monitor these variations is to utilize the cone-beam computed tomography (CT, CBCT) scans routinely taken for patient positioning, reconstruct them as 4DCBCTs, and generate 'virtual CTs' (vCTs), combining the accurate CT numbers of the diagnostic 4DCT and the geometry of the daily 4DCBCT by using deformable image registration (DIR). In this study different algorithms for 4DCBCT reconstruction and DIR were evaluated. For this purpose, CBCT scans of a moving ex vivo porcine lung phantom with 663 and 2350 projections respectively were acquired, accompanied by an additional 4DCT as reference. The CBCT projections were sorted in 10 phase bins with the Amsterdam-shroud method and reconstructed phase-by-phase using first a FDK reconstruction from the Reconstruction Toolkit (RTK) and again an iterative reconstruction algorithm implemented in the Gadgetron Toolkit. The resulting 4DCBCTs were corrected by DIR of the corresponding 4DCT phases, using both a morphons algorithm from REGGUI and a b-spline deformation from Plastimatch. The resulting 4DvCTs were compared to the 4DCT by visual inspection and by calculating water equivalent thickness (WET) maps from the phantom's surface to the distal edge of a target from various angles. The optimized procedure was successfully repeated with mismatched input phases and on a clinical patient dataset. Proton treatment plans were simulated on the 4DvCTs and the dose distributions compared to the reference based on the 4DCT via gamma pass rate analysis. A combination of iterative reconstruction and morphons DIR yielded the most accurate 4DvCTs, with median WET differences under 2mm and 3%/3mm gamma pass rates per phase between 89% and 99%. These results suggest that image correction of iteratively reconstructed 4DCBCTs with a morphons DIR of the planning CT may yield sufficiently accurate 4DvCTs for daily time resolved proton dose calculations.

A novel approach to EPID-based 3D volumetric dosimetry for IMRT and VMAT QA.

Intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) are relatively complex treatment delivery techniques and require quality assurance (QA) procedures. Pre-treatment dosimetric verification represents a fundamental QA procedure in daily clinical routine in radiation therapy. The purpose of this study is to develop an EPID-based approach to reconstruct a 3D dose distribution as imparted to a virtual cylindrical water phantom to be used for plan-specific pre-treatment dosimetric verification for IMRT and VMAT plans. For each depth, the planar 2D dose distributions acquired in air were back-projected and convolved by depth-specific scatter and attenuation kernels. The kernels were obtained by making use of scatter and attenuation models to iteratively estimate the parameters from a set of reference measurements. The derived parameters served as a look-up table for reconstruction of arbitrary measurements. The summation of the reconstructed 3D dose distributions resulted in the integrated 3D dose distribution of the treatment delivery. The accuracy of the proposed approach was validated in clinical IMRT and VMAT plans by means of gamma evaluation, comparing the reconstructed 3D dose distributions with Octavius measurement. The comparison was carried out using (3%, 3 mm) criteria scoring 99% and 96% passing rates for IMRT and VMAT, respectively. An accuracy comparable to the one of the commercial device for 3D volumetric dosimetry was demonstrated. In addition, five IMRT and five VMAT were validated against the 3D dose calculation performed by the TPS in a water phantom using the same passing rate criteria. The median passing rates within the ten treatment plans was 97.3%, whereas the lowest was 95%. Besides, the reconstructed 3D distribution is obtained without predictions relying on forward dose calculation and without external phantom or dosimetric devices. Thus, the approach provides a fully automated, fast and easy QA procedure for plan-specific pre-treatment dosimetric verification.

Investigating deformable image registration and scatter correction for CBCT-based dose calculation in adaptive IMPT.

PURPOSE: This work aims at investigating intensity corrected cone-beam x-ray computed tomography (CBCT) images for accurate dose calculation in adaptive intensity modulated proton therapy (IMPT) for prostate and head and neck (H&N) cancer. A deformable image registration (DIR)-based method and a scatter correction approach using the image data obtained from DIR as prior are characterized and compared on the basis of the same clinical patient cohort for the first time. METHODS: Planning CT (pCT) and daily CBCT data (reconstructed images and measured projections) of four H&N and four prostate cancer patients have been considered in this study. A previously validated Morphons algorithm was used for DIR of the planning CT to the current CBCT image, yielding a so-called virtual CT (vCT). For the first time, this approach was translated from H&N to prostate cancer cases in the scope of proton therapy. The warped pCT images were also used as prior for scatter correction of the CBCT projections for both tumor sites. Single field uniform dose and IMPT (only for H&N cases) treatment plans have been generated with a research version of a commercial planning system. Dose calculations on vCT and scatter corrected CBCT (CBCTcor) were compared by means of the proton range and a gamma-index analysis. For the H&N cases, an additional diagnostic replanning CT (rpCT) acquired within three days of the CBCT served as additional reference. For the prostate patients, a comprehensive contour comparison of CBCT and vCT, using a trained physician's delineation, was performed. RESULTS: A high agreement of vCT and CBCTcor was found in terms of the proton range and gamma-index analysis. For all patients and indications between 95% and 100% of the proton dose profiles in beam's eye view showed a range agreement of better than 3 mm. The pass rate in a (2%,2 mm) gamma-comparison was between 96% and 100%. For H&N patients, an equivalent agreement of vCT and CBCTcor to the reference rpCT was observed. However, for the prostate cases, an insufficient accuracy of the vCT contours retrieved from DIR was found, while the CBCTcor contours showed very high agreement to the contours delineated on the raw CBCT. CONCLUSIONS: For H&N patients, no considerable differences of vCT and CBCTcor were found. For prostate cases, despite the high dosimetric agreement, the DIR yields incorrect contours, probably due to the more pronounced anatomical changes in the abdomen and the reduced soft-tissue contrast in the CBCT. Using the vCT as prior, these inaccuracies can be overcome and images suitable for accurate delineation and dose calculation in CBCT-based adaptive IMPT can be retrieved from scatter correction of the CBCT projections.