Experimental evaluation of a robust optimization method for IMRT of moving targets

Robustness of sweeping-window arc therapy treatment sequences against intrafractional tumor motion

PURPOSE

Due to the potentially periodic collimator dynamic in volumetric modulated arc therapy (VMAT) dose deliveries with the sweeping-window arc therapy (SWAT) technique, additional manifestations of dosimetric deviations in the presence of intrafractional motion may occur. With a fast multileaf collimator (MLC), and a flattening filter free dose delivery, treatment times close to 60 s per fraction are clinical reality. For these treatment sequences, the human breathing period can be close to the collimator sweeping period. Compared to a random arrangement of the segments, this will cause a further degradation of the dose homogeneity.

METHODS

Fifty VMAT sequences of potentially moving target volumes were delivered on a two dimensional ionization chamber array. In order to detect interplay effects along all three coordinate axes, time resolved measurements were performed twice--with the detector aligned in vertical (V) or horizontal (H) orientation. All dose matrices were then moved within a simulation software by a time-dependent motion vector. The minimum relative equivalent uniform dose EUDr,m for all breathing starting phases was determined for each amplitude and period. Furthermore, an estimation of periods with minimum EUD was performed. Additionally, LINAC logfiles were recorded during plan delivery. The MLC, jaw, gantry angle, and monitor unit settings were continuously saved and used to calculate the correlation coefficient between the target motion and the dose weighed collimator motion component for each direction (CC, LR, AP) separately.

RESULTS

The resulting EUDr,m were EUDr,m(CCV) = (98.3 ± 0.6)%, EUDr,m(CCH) = (98.6 ± 0.5)%, EUDr,m(APV) = (97.7 ± 0.9)%, and EUDr,m(LRH) = (97.8 ± 0.9)%. The overall minimum relative EUD observed for 360(∘) arc midventilation treatments was 94.6%. The treatment plan with the shortest period and a minimum relative EUD of less than 97% was found at T = 6.1 s. For a partial 120(∘) arc, an EUDr,m = 92.0% was found. In all cases, a correlation coefficient above 0.5 corresponded to a minimum in EUD.

CONCLUSIONS

With the advent of fast VMAT delivery techniques, nonrobust treatment sequences for human breathing patterns can be generated. These sequences are characterized by a large correlation coefficient between a target motion component and the corresponding collimator dynamic. By iteratively decreasing the maximum allowed dose rate, a low correlation coefficient and consequentially a robust treatment sequence are ensured.

Adaptive and robust radiation therapy in the presence of drift

Combining adaptive and robust optimization in radiation therapy has the potential to mitigate the negative effects of both intrafraction and interfraction uncertainty over a fractionated treatment course. A previously developed adaptive and robust radiation therapy (ARRT) method for lung cancer was demonstrated to be effective when the sequence of breathing patterns was well-behaved. In this paper, we examine the applicability of the ARRT method to less well-behaved breathing patterns. We develop a novel method to generate sequences of probability mass functions that represent different types of drift in the underlying breathing pattern. Computational results derived from applying the ARRT method to these sequences demonstrate that the ARRT method is effective for a much broader class of breathing patterns than previously demonstrated.

Quantitative metrics for assessing plan quality

Despite many studies over the last 3 decades that have attempted to explicitly quantify the decision-making process for radiotherapy treatment plan evaluation, judgments of an individual plan's degree of quality are still largely subjective and can show inter- and intra-practitioner variability even if the clinical treatment goals are the same. Several factors conspire to confound the full quantification of treatment plan quality, including uncertainties in dose response of cancerous and normal tissue, the rapid pace of new technology adoption, and the human component of treatment planning. However, new developments in clinical informatics and automation are lowering the bar for developing and implementing quantitative metrics into the treatment planning process. This review discusses general strategies for using quantitative metrics in the treatment planning process and presents a case study in intensity-modulated radiation therapy planning whereby control was established on a variable system via such techniques.

Bridging the gap between IMRT and VMAT: dense angularly sampled and sparse intensity modulated radiation therapy

PURPOSE

To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency.

METHODS

In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency.

RESULTS

Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.

CONCLUSION

DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

Intensity modulation under geometrical uncertainty: a deconvolution approach to robust fluence

A deconvolution algorithm has been developed to obtain robust fluence for external beam radiation treatment under geometrical uncertainties. Usually, the geometrical uncertainty is incorporated in the dose optimization process for inverse treatment planning to determine the additional intensity modulation of the beam to counter the geometrical uncertainty. Most of these approaches rely on dose convolution which is subject to the error caused by patient surface curvature and internal inhomogeneity. In this work, based on an 1D deconvolution algorithm developed by Ulmer and Kaissl, a fluence-deconvolution approach was developed to obtain robust fluence through the deconvolution of the nominal static one given by any treatment planning system. It incorporates the geometrical uncertainty outside the dose optimization procedure and therefore avoids the error of dose convolution. Robust fluences were calculated for a 4 x 4 cm flat field, a prostate IMRT and a head and neck IMRT plan in a commercial treatment planning system. The corresponding doses were simulated for 30 fractions with the random Gaussian distribution of the iso-centers showing good agreement with the nominal static doses. The feasibility of this deconvolution approach for clinical IMRT planning has been demonstrated. Because it is separated from the optimization procedure, this method is more flexible and easier to integrate into different existing treatment planning systems to obtain robust fluence.