Comparison of planned dose on different CT image sets to four-dimensional Monte Carlo dose recalculation using the patient's actual breathing trace for lung stereotactic body radiation therapy

Optimized scoring of end-to-end dosimetry audits for passive motion management - A simulation study using the IROC thorax phantom

Dosimetry audits for passive motion management require dynamically-acquired measurements in a moving phantom to be compared to statically calculated planned doses. This study aimed to characterise the relationship between planning and delivery errors, and the measured dose in the Imaging and Radiation Oncology Core (IROC) thorax phantom, to assess different audit scoring approaches. Treatment plans were created using a 4DCT scan of the IROC phantom, equipped with film and thermoluminescent dosimeters (TLDs). Plans were created on the average intensity projection from all bins. Three levels of aperture complexity were explored: dynamic conformal arcs (DCAT), low-, and high-complexity volumetric modulated arcs (VMATLo, VMATHi). Simulated-measured doses were generated by modelling motion using isocenter shifts. Various errors were introduced including incorrect setup position and target delineation. Simulated-measured film doses were scored using gamma analysis and compared within specific regions of interest (ROIs) as well as the entire film plane. Positional offsets were estimated based on isodoses on the film planes, and point doses within TLD contours were compared. Motion-induced differences between planned and simulated-measured doses were evident even without introduced errors Gamma passing rates within target-centred ROIs correlated well with error-induced dose differences, while whole film passing rates did not. Isodose-based setup position measurements demonstrated high sensitivity to errors. Simulated point doses at TLD locations yielded erratic responses to introduced errors. ROI gamma analysis demonstrated enhanced sensitivity to simulated errors compared to whole film analysis. Gamma results may be further contextualized by other metrics such as setup position or maximum gamma.

Quantifying the Effect of 4-Dimensional Computed Tomography-Based Deformable Dose Accumulation on Representing Radiation Damage for Patients with Locally Advanced Non-Small Cell Lung Cancer Treated with Standard-Fractionated Intensity-Modulated Radiation Therapy

PURPOSE

The aim of this study was to investigate the dosimetric and clinical effects of 4-dimensional computed tomography (4DCT)-based longitudinal dose accumulation in patients with locally advanced non-small cell lung cancer treated with standard-fractionated intensity-modulated radiation therapy (IMRT).

METHODS AND MATERIALS

Sixty-seven patients were retrospectively selected from a randomized clinical trial. Their original IMRT plan, planning and verification 4DCTs, and ∼4-month posttreatment follow-up CTs were imported into a commercial treatment planning system. Two deformable image registration algorithms were implemented for dose accumulation, and their accuracies were assessed. The planned and accumulated doses computed using average-intensity images or phase images were compared. At the organ level, mean lung dose and normal-tissue complication probability (NTCP) for grade ≥2 radiation pneumonitis were compared. At the region level, mean dose in lung subsections and the volumetric overlap between isodose intervals were compared. At the voxel level, the accuracy in estimating the delivered dose was compared by evaluating the fit of a dose versus radiographic image density change (IDC) model. The dose-IDC model fit was also compared for subcohorts based on the magnitude of NTCP difference (|ΔNTCP|) between planned and accumulated doses.

RESULTS

Deformable image registration accuracy was quantified, and the uncertainty was considered for the voxel-level analysis. Compared with planned doses, accumulated doses on average resulted in <1-Gy lung dose increase and <2% NTCP increase (up to 8.2 Gy and 18.8% for a patient, respectively). Volumetric overlap of isodose intervals between the planned and accumulated dose distributions ranged from 0.01 to 0.93. Voxel-level dose-IDC models demonstrated a fit improvement from planned dose to accumulated dose (pseudo-R2 increased 0.0023) and a further improvement for patients with ≥2% |ΔNTCP| versus for patients with <2% |ΔNTCP|.

CONCLUSIONS

With a relatively large cohort, robust image registrations, multilevel metric comparisons, and radiographic image-based evidence, we demonstrated that dose accumulation more accurately represents the delivered dose and can be especially beneficial for patients with greater longitudinal response.

4D Monte Carlo dose reconstructions using surface motion measurements

PURPOSE

To investigate the feasibility of a 4D Monte Carlo based dose reconstruction method to study the dosimetric impact of respiratory motion using surface motion measurements for patients undergoing VMAT treatments for Non-Small Cell Lung Cancer.

METHODS

The 4Ddefdosxyznrc/EGSnrc algorithm was used to reconstruct VMAT doses delivered to the patients using machine log files and respiratory traces measured with the RADPOS 4D dosimetry system. The RADPOS sensor was adhered to the patient's abdomen prior to each treatment fraction and its position was used as a surrogate for tumour motion. Treatment log files were synchronized with the patient respiratory traces. Patient specific respiratory models were generated from deformable registration of the inhale and exhale 4DCT images and the respiratory traces. The reconstructed doses were compared to planned doses calculated with DOSXYZnrc/EGSnrc on the average-intensity and the exhale phase CT images.

RESULTS

Respiratory motion measurements and log files were acquired for 2 patients over 5 treatment fractions each. The motion was predominantly along the anterior/posterior direction (A/P). The average respiratory amplitudes were 8.7 ± 2.7 mm and 10.0 ± 1.2 mm for Patient 1 and 2, respectively. Both patients displayed inter- and intra-fractional variations in the baseline position. Small inter-fractional differences were observed in the reconstructed doses for each patient. Differences between the reconstructed and planned doses were attributed to differences in organ volumes.

CONCLUSION

The 4D reconstruction method was successfully implemented for the two patients studied. Small differences between the planned and reconstructed doses were observed due to the small tumour motion of these patients.

Adoption of respiratory motion management in radiation therapy

Background and Purpose

A survey on the patterns of practice of respiratory motion management (MM) was distributed to 111 radiation therapy facilities to inform the development of an end-to-end dosimetry audit including respiratory motion.

Materials and methods

The survey (distributed via REDCap) asked facilities to provide information specific to the combinations of MM techniques (breath-hold gating – BHG, internal target volume – ITV, free-breathing gating – FBG, mid-ventilation – MidV, tumour tracking – TT), sites treated (thorax, upper abdomen, lower abdomen), and fractionation regimes (conventional, stereotactic ablative body radiation therapy – SABR) used in their clinic.

Results

The survey was completed by 78% of facilities, with 98% of respondents indicating that they used at least one form of MM. The ITV approach was common to all MM-users, used for thoracic treatments by 89% of respondents, and upper and lower abdominal treatments by 38%. BHG was the next most prevalent (41% of MM users), with applications in upper abdominal and thoracic treatment sites (28% vs 25% respectively), but minimal use in the lower abdomen (9%). FBG and TT were utilised sparingly (17%, 7% respectively), and MidV was not selected at all.

Conclusions

Two distinct treatment workflows (including use of motion limitation, imaging used for motion assessment, dose calculation, and image guidance procedures) were identified for the ITV and BHG MM techniques, to form the basis of the initial audit. Thoracic SABR with the ITV approach was common to nearly all respondents, while upper abdominal SABR using BHG stood out as more technically challenging. Other MM techniques were sparsely used, but may be considered for future audit development.

Statistical breathing curve sampling to quantify interplay effects of moving lung tumors in a 4D Monte Carlo dose calculation framework

PURPOSE

The interplay between respiratory tumor motion and dose application by intensity modulated radiotherapy (IMRT) techniques can potentially lead to undesirable and non-intuitive deviations from the planned dose distribution. We developed a 4D Monte Carlo (MC) dose recalculation framework featuring statistical breathing curve sampling, to precisely simulate the dose distribution for moving target volumes aiming at a comprehensive assessment of interplay effects.

METHODS

We implemented a dose accumulation tool that enables dose recalculations of arbitrary breathing curves including the actual breathing curve of the patient. This MC dose recalculation framework is based on linac log-files, facilitating a high temporal resolution up to 0.1 s. By statistical analysis of 128 different breathing curves, interplay susceptibility of different treatment parameters was evaluated for an exemplary patient case. To facilitate prospective clinical application in the treatment planning stage, in which patient breathing curves or linac log-files are not available, we derived a log-file free version with breathing curves generated by a random walk approach. Interplay was quantified by standard deviations σ in D_5%, D_50% and D_95%.

RESULTS

Interplay induced dose deviations for single fractions were observed and evaluated for IMRT and volumetric arc therapy (σ_D95% up to 1.3 %) showing a decrease with higher fraction doses and an increase with higher MU rates. Interplay effects for conformal treatment techniques were negligible (σ<0.1%). The log-file free version and the random walk generated breathing curves yielded similar results (deviations in σ< 0.1 %) and can be used as substitutes for interplay assessment.

CONCLUSION

It is feasible to combine statistically sampled breathing curves with MC dose calculations. The universality of the presented framework allows comprehensive assessment of interplay effects in retrospective and prospective clinically relevant scenarios.

Biologically consistent dose accumulation using daily patient imaging

BACKGROUND

This work addresses a basic inconsistency in the way dose is accumulated in radiotherapy when predicting the biological effect based on the linear quadratic model (LQM). To overcome this inconsistency, we introduce and evaluate the concept of the total biological dose, bEQD_d.

METHODS

Daily computed tomography imaging of nine patients treated for prostate carcinoma with intensity-modulated radiotherapy was used to compute the delivered deformed dose on the basis of deformable image registration (DIR). We compared conventional dose accumulation (DA) with the newly introduced bEQD_d, a new method of accumulating biological dose that considers each fraction dose and tissue radiobiology. We investigated the impact of the applied fractionation scheme (conventional/hypofractionated), uncertainties induced by the DIR and by the assigned α/β-value.

RESULTS

bEQD_d was systematically higher than the conventionally accumulated dose with difference hot spots of 3.3-4.9 Gy detected in six out of nine patients in regions of high dose gradient in the bladder and rectum. For hypofractionation, differences are up to 8.4 Gy. The difference amplitude was found to be in a similar range to worst-case uncertainties induced by DIR and was higher than that induced by α/β.

CONCLUSION

Using bEQD_d for dose accumulation overcomes a potential systematic inaccuracy in biological effect prediction based on accumulated dose. Highest impact is found for serial-type late responding organs at risk in dose gradient regions and for hypofractionation. Although hot spot differences are in the order of several Gray, in dose-volume parameters there is little difference compared with using conventional or biological DA. However, when local dose information is used, e.g. dose surface maps, difference hot spots can potentially change outcomes of dose-response modelling and adaptive treatment strategies.

Integrating CVH and LVH metrics into an optimization strategy for the selection of Iris collimator for Cyberknife Xsight lung tracking treatment

PURPOSE

We conducted this study to construct a target coverage-volume histogram (CVH) and leakage-volume histogram (LVH) metrics and optimization strategy for the selection of the Iris collimator in Cyberknife Xsight lung tracking treatment through a retrospective analysis of target structures and clinical data.

METHODS AND MATERIALS

CVH and LVH metrics were retrospectively analyzed for 37 lung cancer patients. CVH and LVH were the same as dose-volume histogram (DVH), but with a coverage and leakage replacing dose. For each patient, Iris collimator was optimized and selected based on CVH and LVH metrics. The CVH and LVH metrics were then compared to ascertain differences in 95% (C95) or 90% (C90) of the target coverage thresholds. The planning target volume (PTV) C95 and C90 coverage, absolute mean leakage value, leakage/coverage ratio, selected collimator diameter (Φ), Φ/length of the long axis of PTV (A_max ), and Φ/length of the short axis (A_min ) of PTV were compared. The correlation of the absolute mean leakage value, leakage/coverage ratio, Φ/A_min and Φ/A_max were evaluated.

RESULTS

For each patient, the PTV C95 coverage (70.45 vs 63.19) and C90 coverage (77.25 vs 69.96) were higher in the C95 coverage threshold group compared to the C90 coverage threshold group. The leakage/coverage ratio (0.56 vs 0.69) and absolute mean leakage value (0.56 vs 0.61) were lower in C90 coverage threshold group than in C95 coverage threshold group. The Spearmen correlation test showed the Φ/A_min were significantly correlated with leakage/coverage ratio and absolute mean leakage value. Upon analysis of the selected collimator diameters, the mean value of Φ/A_min of the optimized collimator diameters was found to be 1.10.

CONCLUSION

The CVH and LVH analysis is able to quantitatively evaluate the tradeoff between target coverage and normal tissue sparing.

Monte Carlo evaluation of target dose coverage in lung stereotactic body radiation therapy with flattening filter-free beams

Aim:

Previous studies showed that replacing conventional flattened beams (FF) with flattening filter-free (FFF) beams improves the therapeutic ratio in lung stereotactic body radiation therapy (SBRT), but these findings could have been impacted by dose calculation uncertainties caused by the heterogeneity of the thoracic anatomy and by respiratory motion, which were particularly high for target coverage. In this study, we minimised such uncertainties by calculating doses using high-spatial-resolution Monte Carlo and four-dimensional computed tomography (4DCT) images. We aimed to evaluate more reliably the benefits of using FFF beams for lung SBRT.

Materials and methods:

For a cohort of 15 patients with early-stage lung cancer that we investigated in a previous treatment planning study, we recalculated dose distributions with Monte Carlo using 4DCT images. This included 15 FF and 15 FFF treatment plans.

Results:

Compared to Monte Carlo, the treatment planning system (TPS) over-predicted doses in low-dose regions of the planning target volume (PTV). For most patients, replacing FF beams with FFF beams improved target coverage, tumour control, and uncomplicated tumour control probabilities.

Conclusions:

Monte Carlo tends to reveal deficiencies in target coverage compared to coverage predicted by the TPS. Our data support previously reported benefits of using FFF beams for lung SBRT.

Technological quality requirements for stereotactic radiotherapy : Expert review group consensus from the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy

This review details and discusses the technological quality requirements to ensure the desired quality for stereotactic radiotherapy using photon external beam radiotherapy as defined by the DEGRO Working Group Radiosurgery and Stereotactic Radiotherapy and the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy. The covered aspects of this review are 1) imaging for target volume definition, 2) patient positioning and target volume localization, 3) motion management, 4) collimation of the irradiation and beam directions, 5) dose calculation, 6) treatment unit accuracy, and 7) dedicated quality assurance measures. For each part, an expert review for current state-of-the-art techniques and their particular technological quality requirement to reach the necessary accuracy for stereotactic radiotherapy divided into intracranial stereotactic radiosurgery in one single fraction (SRS), intracranial fractionated stereotactic radiotherapy (FSRT), and extracranial stereotactic body radiotherapy (SBRT) is presented. All recommendations and suggestions for all mentioned aspects of stereotactic radiotherapy are formulated and related uncertainties and potential sources of error discussed. Additionally, further research and development needs in terms of insufficient data and unsolved problems for stereotactic radiotherapy are identified, which will serve as a basis for the future assignments of the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy. The review was group peer-reviewed, and consensus was obtained through multiple working group meetings.

4D Monte Carlo dose calculations for pre-treatment quality assurance of VMAT SBRT: a phantom-based feasibility study

Volumetric arc therapy (VMAT) for lung stereotactic body radiotherapy (SBRT) is challenging due to both breathing-induced motion and the dynamic components of the linear accelerator. In this study, a 4D Monte Carlo (4DMC) dose calculation method for VMAT SBRT is proposed and the feasibility of the method is evaluated. A rigidly-moving lung phantom was imaged using four dimensional computed tomography (4DCT). VMAT SBRT plans were generated on the average intensity projection dataset using the internal target volume (ITV) strategy (ITV-plan) and a single phase to simulate a dynamic treatment-couch tracking technique (TRACKING-plan). 4DMC simulations were performed and compared to 3D Monte Carlo (3DMC) and 3D- and 4D- calculations in the treatment planning system using the adaptive convolution (AC) algorithm. Dose metrics calculated for the ITV-plan showed an overestimation with 3D adaptive convolution (3DAC) for D[Formula: see text] (GTV) by 3.5% and by 2.0% for 3DMC, both compared to 4DMC. The TRACKING-plan D[Formula: see text] (GTV) calculated with the 3DAC method overestimated by 2.0% compared with 4DMC. Deviations between the calculation methods for D _mean (Lung) and D[Formula: see text] (PTV) were minimal. For both plans, measurements were taken with EBT3 film inside the phantom tumour. EBT3 film profiles showed good agreement with 4DMC for the TRACKING-plan giving a gamma pass rate of 97.2% for 3%/3 mm global and for 3DAC compared with measured, 95.8%. Whereas for the ITV-plan, the 3D profiles varied from film in the ITV periphery region with a pass rates of 50% and 48.6% for 3DAC and 3DMC, respectively. 4DMC agreed more closely to measurements for this plan with a pass rate of 95.8%. We have proposed an accurate method to perform 4D dose calculations for pre-treatment quality assurance of VMAT SBRT. The method was compared to experimental measurements and for both plans, 4DMC dose agreed with measurements more closely than other evaluated dose calculation methods. This study has demonstrated the feasibility of this 4DMC method.