A data-driven approach to optimal beam/arc angle selection for liver stereotactic body radiation therapy treatment planning

BACKGROUND

Stereotactic body radiation therapy (SBRT) for liver cancer has shown promising therapeutic effects. Effective treatment relies not only on the precise delivery provided by image-guided radiation therapy (IGRT) but also high dose gradient formed around the treatment volume to spare functional liver tissue, which is highly dependent on the beam/arc angle selection. In this study, we aim to develop a decision support model to learn human planner's beam navigation approach for beam angle/arc angle selection for liver SBRT.

METHODS

A total of 27 liver SBRT/HIGRT patients (10 IMRT, 17 VMAT/DCA) were included in this study. A dosimetric budget index was defined for each beam angle/control point considering dose penetration through the patient body and liver tissue. Optimal beam angle setting (beam angles for IMRT and start/terminal angle for VMAT/DCA) was determined by minimizing the loss function defined as the sum of total dosimetric budget index and beam span penalty function. Leave-one-out validation was exercised on all 27 cases while weighting coefficients in the loss function was tuned in nested cross validation. To compare the efficacy of the model, a model plan was generated using automatically generated beam setting while retaining the original optimization constraints in the clinical plan. Model plan was normalized to the same planning target volume (PTV) V100% as the clinical plans. Dosimetric endpoints including PTV D98%, D2%, liver V20Gy and total MU were compared between two plan groups. Wilcoxon Signed-Rank test was performed with the null hypothesis being that no difference exists between two plan groups.

RESULTS

Beam setting prediction was instantaneous. Mean PTV D98% was 91.3% and 91.3% (P=0.566), while mean PTV D2% was 107.9% and 108.1% (P=0.164) for clinical plan and model plan respectively. Liver V20Gy showed no significant difference (P=0.590) with 23.3% for clinical plan and 23.4% for the model plan. Total MU is comparable (P=0.256) between the clinical plan (avg. 2,389.6) and model plan (avg. 2,319.6).

CONCLUSIONS

The evidence driven beam setting model yielded similar plan quality as hand-crafted clinical plan. It is capable of capturing human's knowledge in beam selection decision making. This model could facilitate decision making for beam angle selection while eliminating lengthy trial-and-error process of adjusting beam setting during liver SBRT treatment planning.

The Scatter Search Based Algorithm for Beam Angle Optimization in Intensity-Modulated Radiation Therapy

This article introduces a new framework for beam angle optimization (BAO) in intensity-modulated radiation therapy (IMRT) using the Scatter Search Based Algorithm. The potential benefits of plans employing the coplanar optimized beam sets are also examined. In the proposed beam angle selection algorithm, the problem is solved in two steps. Initially, the gantry angles are selected using the Scatter Search Based Algorithm, which is a global optimization method. Then, for each beam configuration, the intensity profile is calculated by the conjugate gradient method to score each beam angle set chosen. A simulated phantom case with obvious optimal beam angles was used to benchmark the validity of the presented algorithm. Two clinical cases (TG-119 phantom and prostate cases) were examined to prepare a dose volume histogram (DVH) and determine the dose distribution to evaluate efficiency of the algorithm. A clinical plan with the optimized beam configuration was compared with an equiangular plan to determine the efficiency of the proposed algorithm. The BAO plans yielded significant improvements in the DVHs and dose distributions compared to the equispaced coplanar beams for each case. The proposed algorithm showed its potential to effectively select the beam direction for IMRT inverse planning at different tumor sites.

Automatic learning-based beam angle selection for thoracic IMRT

PURPOSE

The treatment of thoracic cancer using external beam radiation requires an optimal selection of the radiation beam directions to ensure effective coverage of the target volume and to avoid unnecessary treatment of normal healthy tissues. Intensity modulated radiation therapy (IMRT) planning is a lengthy process, which requires the planner to iterate between choosing beam angles, specifying dose-volume objectives and executing IMRT optimization. In thorax treatment planning, where there are no class solutions for beam placement, beam angle selection is performed manually, based on the planner's clinical experience. The purpose of this work is to propose and study a computationally efficient framework that utilizes machine learning to automatically select treatment beam angles. Such a framework may be helpful for reducing the overall planning workload.

METHODS

The authors introduce an automated beam selection method, based on learning the relationships between beam angles and anatomical features. Using a large set of clinically approved IMRT plans, a random forest regression algorithm is trained to map a multitude of anatomical features into an individual beam score. An optimization scheme is then built to select and adjust the beam angles, considering the learned interbeam dependencies. The validity and quality of the automatically selected beams evaluated using the manually selected beams from the corresponding clinical plans as the ground truth.

RESULTS

The analysis included 149 clinically approved thoracic IMRT plans. For a randomly selected test subset of 27 plans, IMRT plans were generated using automatically selected beams and compared to the clinical plans. The comparison of the predicted and the clinical beam angles demonstrated a good average correspondence between the two (angular distance 16.8° ± 10°, correlation 0.75 ± 0.2). The dose distributions of the semiautomatic and clinical plans were equivalent in terms of primary target volume coverage and organ at risk sparing and were superior over plans produced with fixed sets of common beam angles. The great majority of the automatic plans (93%) were approved as clinically acceptable by three radiation therapy specialists.

CONCLUSIONS

The results demonstrated the feasibility of utilizing a learning-based approach for automatic selection of beam angles in thoracic IMRT planning. The proposed method may assist in reducing the manual planning workload, while sustaining plan quality.

A surrogate-based metaheuristic global search method for beam angle selection in radiation treatment planning

An important element of radiation treatment planning for cancer therapy is the selection of beam angles (out of all possible coplanar and non-coplanar angles in relation to the patient) in order to maximize the delivery of radiation to the tumor site and minimize radiation damage to nearby organs-at-risk. This category of combinatorial optimization problem is particularly difficult because direct evaluation of the quality of treatment corresponding to any proposed selection of beams requires the solution of a large-scale dose optimization problem involving many thousands of variables that represent doses delivered to volume elements (voxels) in the patient. However, if the quality of angle sets can be accurately estimated without expensive computation, a large number of angle sets can be considered, increasing the likelihood of identifying a very high quality set. Using a computationally efficient surrogate beam set evaluation procedure based on single-beam data extracted from plans employing equallyspaced beams (eplans), we have developed a global search metaheuristic process based on the nested partitions framework for this combinatorial optimization problem. The surrogate scoring mechanism allows us to assess thousands of beam set samples within a clinically acceptable time frame. Tests on difficult clinical cases demonstrate that the beam sets obtained via our method are of superior quality.

Intensity modulated radiation therapy with field rotation--a time-varying fractionation study

This paper proposes a novel mathematical approach to the beam selection problem in intensity modulated radiation therapy (IMRT) planning. The approach allows more beams to be used over the course of therapy while limiting the number of beams required in any one session. In the proposed field rotation method, several sets of beams are interchanged throughout the treatment to allow a wider selection of beam angles than would be possible with fixed beam orientations. The choice of beamlet intensities and the number of identical fractions for each set are determined by a mixed integer linear program that controls jointly for the distribution per fraction and the cumulative dose distribution delivered to targets and critical structures. Trials showed the method allowed substantial increases in the dose objective and/or sparing of normal tissues while maintaining cumulative and fraction size limits. Trials for a head and neck site showed gains of 25%-35% in the objective (average tumor dose) and for a thoracic site gains were 7%-13%, depending on how strict the fraction size limits were set. The objective did not rise for a prostate site significantly, but the tolerance limits on normal tissues could be strengthened with the use of multiple beam sets.

PARETO: A novel evolutionary optimization approach to multiobjective IMRT planning

PURPOSE

In radiation therapy treatment planning, the clinical objectives of uniform high dose to the planning target volume (PTV) and low dose to the organs-at-risk (OARs) are invariably in conflict, often requiring compromises to be made between them when selecting the best treatment plan for a particular patient. In this work, the authors introduce Pareto-Aware Radiotherapy Evolutionary Treatment Optimization (pareto), a multiobjective optimization tool to solve for beam angles and fluence patterns in intensity-modulated radiation therapy (IMRT) treatment planning.

METHODS

pareto is built around a powerful multiobjective genetic algorithm (GA), which allows us to treat the problem of IMRT treatment plan optimization as a combined monolithic problem, where all beam fluence and angle parameters are treated equally during the optimization. We have employed a simple parameterized beam fluence representation with a realistic dose calculation approach, incorporating patient scatter effects, to demonstrate feasibility of the proposed approach on two phantoms. The first phantom is a simple cylindrical phantom containing a target surrounded by three OARs, while the second phantom is more complex and represents a paraspinal patient.

RESULTS

pareto results in a large database of Pareto nondominated solutions that represent the necessary trade-offs between objectives. The solution quality was examined for several PTV and OAR fitness functions. The combination of a conformity-based PTV fitness function and a dose-volume histogram (DVH) or equivalent uniform dose (EUD) -based fitness function for the OAR produced relatively uniform and conformal PTV doses, with well-spaced beams. A penalty function added to the fitness functions eliminates hotspots. Comparison of resulting DVHs to those from treatment plans developed with a single-objective fluence optimizer (from a commercial treatment planning system) showed good correlation. Results also indicated that pareto shows promise in optimizing the number of beams.

CONCLUSIONS

This initial evaluation of the evolutionary optimization software tool pareto for IMRT treatment planning demonstrates feasibility and provides motivation for continued development. Advantages of this approach over current commercial methods for treatment planning are many, including: (1) fully automated optimization that avoids human controlled iterative optimization and potentially improves overall process efficiency, (2) formulation of the problem as a true multiobjective one, which provides an optimized set of Pareto nondominated solutions refined over hundreds of generations and compiled from thousands of parameter sets explored during the run, and (3) rapid exploration of the final nondominated set accomplished by a graphical interface used to select the best treatment option for the patient.

An algorithm for fast beam angle selection in intensity modulated radiotherapy

PURPOSE

This article aims to introduce a novel algorithm for fast beam angle selection in intensity modulated radiotherapy (IMRT).

METHODS

The algorithm models the optimization problem as a beam angle ranking problem and chooses suitable beam angles according to their rank. A new parameter called "beam intensity profile perturbation score (BIPPS)" is used for ranking the beam angles. The BIPPS-based beam angle ranking implicitly accounts for the dose-volume effects of the involved structures. A simulated phantom case with obvious optimal beam angles is used to verify the validity of the presented technique. In addition, the efficiency of the algorithm was examined in three clinical cases (prostate, pancreas, and head and neck) in terms of DVH and dose distribution. In all cases, the judgment of the algorithm's efficiency was based on the comparison between plans with equidistant beams (equal-angle-plan) and plans with beams obtained using the algorithm (suitable-angle-plan).

RESULTS

It is observed from the study that the beam angle ranking function over BIPPS instantly picks up a suitable set of beam angles for a specific case. It takes only about 15 min for choosing the suitable beam angles even for the most complicated cases. The DVHs and dose distributions confirm that the proposed algorithm can efficiently reduce the mean or maximum dose to OARs, while guaranteeing the target coverage and dose uniformity. On the average, about 17% reduction in the mean dose to critical organs, such as rectum, bladder, kidneys and parotids, is observed. Also, about 12% (averaged) reduction in the maximum dose to critical organs (spinal cord) is observed in the clinical cases presented in this study.

CONCLUSIONS

This study demonstrates that the algorithm can be effectively applied to IMRT scenarios to get fast and case specific beam angle configurations.