Multiobjective optimization with a modified simulated annealing algorithm for external beam radiotherapy treatment planning

Personalized treatment planning with a model of radiation therapy outcomes for use in multiobjective optimization of IMRT plans for prostate cancer

PURPOSE

To build a new treatment planning approach that extends beyond radiation transport and IMRT optimization by modeling the radiation therapy process and prognostic indicators for more outcome-focused decision making.

METHODS

An in-house treatment planning system was modified to include multiobjective inverse planning, a probabilistic outcome model, and a multi-attribute decision aid. A genetic algorithm generated a set of plans embodying trade-offs between the separate objectives. An influence diagram network modeled the radiation therapy process of prostate cancer using expert opinion, results of clinical trials, and published research. A Markov model calculated a quality adjusted life expectancy (QALE), which was the endpoint for ranking plans.

RESULTS

The Multiobjective Evolutionary Algorithm (MOEA) was designed to produce an approximation of the Pareto Front representing optimal tradeoffs for IMRT plans. Prognostic information from the dosimetrics of the plans, and from patient-specific clinical variables were combined by the influence diagram. QALEs were calculated for each plan for each set of patient characteristics. Sensitivity analyses were conducted to explore changes in outcomes for variations in patient characteristics and dosimetric variables. The model calculated life expectancies that were in agreement with an independent clinical study.

CONCLUSIONS

The radiation therapy model proposed has integrated a number of different physical, biological and clinical models into a more comprehensive model. It illustrates a number of the critical aspects of treatment planning that can be improved and represents a more detailed description of the therapy process. A Markov model was implemented to provide a stronger connection between dosimetric variables and clinical outcomes and could provide a practical, quantitative method for making difficult clinical decisions.

A self-adaptive case-based reasoning system for dose planning in prostate cancer radiotherapy

PURPOSE

Prostate cancer is the most common cancer in the male population. Radiotherapy is often used in the treatment for prostate cancer. In radiotherapy treatment, the oncologist makes a trade-off between the risk and benefit of the radiation, i.e., the task is to deliver a high dose to the prostate cancer cells and minimize side effects of the treatment. The aim of our research is to develop a software system that will assist the oncologist in planning new treatments.

METHODS

A nonlinear case-based reasoning system is developed to capture the expertise and experience of oncologists in treating previous patients. Importance (weights) of different clinical parameters in the dose planning is determined by the oncologist based on their past experience, and is highly subjective. The weights are usually fixed in the system. In this research, the weights are updated automatically each time after generating a treatment plan for a new patient using a group based simulated annealing approach.

RESULTS

The developed approach is analyzed on the real data set collected from the Nottingham University Hospitals NHS Trust, City Hospital Campus, UK. Extensive experiments show that the dose plan suggested by the proposed method is coherent with the dose plan prescribed by an experienced oncologist or even better.

CONCLUSIONS

The developed case-based reasoning system enables the use of knowledge and experience gained by the oncologist in treating new patients. This system may play a vital role to assist the oncologist in making a better decision in less computational time; it utilizes the success rate of the previously treated patients and it can also be used in teaching and training processes.

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.

Functional avoidance of lung in plan optimization with an aperture-based inverse planning system

PURPOSE

To implement SPECT-based optimization in an anatomy-based aperture inverse planning system for the functional avoidance of lung in thoracic irradiation.

MATERIAL AND METHODS

SPECT information has been introduced as a voxel-by-voxel modulation of lung importance factors proportionally to the local perfusion count. Fifteen cases of lung cancer have been retrospectively analyzed by generating angle-optimized non-coplanar plans, comparing a purely anatomical approach and our functional approach. Planning target volume coverage and lung sparing have been compared. Statistical significance was assessed by a Wilcoxon matched pairs test.

RESULTS

For similar target coverage, perfusion-weighted volume receiving 10 Gy was reduced by a median of 2.2% (p=0.022) and mean perfusion-weighted lung dose, by a median of 0.9 Gy (p=0.001). A separate analysis of patients with localized or non-uniform hypoperfusion could not show which would benefit more from SPECT-based treatment planning. Redirection of dose sometimes created overdosage regions in the target volume. Plans consisted of a similar number of segments and monitor units.

CONCLUSIONS

Angle optimization and SPECT-based modulation of importance factors allowed for functional avoidance of the lung while preserving target coverage. The technique could be also applied to implement PET-based modulation inside the target volume, leading to a safer dose escalation.

An approach to multiobjective optimization of rotational therapy. II. Pareto optimal surfaces and linear combinations of modulated blocked arcs for a prostate geometry

PURPOSE

The purpose of this work is twofold: To further develop an approach to multiobjective optimization of rotational therapy treatments recently introduced by the authors [J. Pardo-Montero and J. D. Fenwick, "An approach to multiobjective optimization of rotational therapy," Med. Phys. 36, 3292-3303 (2009)], especially regarding its application to realistic geometries, and to study the quality (Pareto optimality) of plans obtained using such an approach by comparing them with Pareto optimal plans obtained through inverse planning.

METHODS

In the previous work of the authors, a methodology is proposed for constructing a large number of plans, with different compromises between the objectives involved, from a small number of geometrically based arcs, each arc prioritizing different objectives. Here, this method has been further developed and studied. Two different techniques for constructing these arcs are investigated, one based on image-reconstruction algorithms and the other based on more common gradient-descent algorithms. The difficulty of dealing with organs abutting the target, briefly reported in previous work of the authors, has been investigated using partial OAR unblocking. Optimality of the solutions has been investigated by comparison with a Pareto front obtained from inverse planning. A relative Euclidean distance has been used to measure the distance of these plans to the Pareto front, and dose volume histogram comparisons have been used to gauge the clinical impact of these distances. A prostate geometry has been used for the study.

RESULTS

For geometries where a blocked OAR abuts the target, moderate OAR unblocking can substantially improve target dose distribution and minimize hot spots while not overly compromising dose sparing of the organ. Image-reconstruction type and gradient-descent blocked-arc computations generate similar results. The Pareto front for the prostate geometry, reconstructed using a large number of inverse plans, presents a hockey-stick shape comprising two regions: One where the dose to the target is close to prescription and trade-offs can be made between doses to the organs at risk and (small) changes in target dose, and one where very substantial rectal sparing is achieved at the cost of large target underdosage. Plans computed following the approach using a conformal arc and four blocked arcs generally lie close to the Pareto front, although distances of some plans from high gradient regions of the Pareto front can be greater. Only around 12% of plans lie a relative Euclidean distance of 0.15 or greater from the Pareto front. Using the alternative distance measure of Craft ["Calculating and controlling the error of discrete representations of Pareto surfaces in convex multi-criteria optimization," Phys. Medica (to be published)], around 2/5 of plans lie more than 0.05 from the front. Computation of blocked arcs is quite fast, the algorithms requiring 35%-80% of the running time per iteration needed for conventional inverse plan computation.

CONCLUSIONS

The geometry-based arc approach to multicriteria optimization of rotational therapy allows solutions to be obtained that lie close to the Pareto front. Both the image-reconstruction type and gradient-descent algorithms produce similar modulated arcs, the latter one perhaps being preferred because it is more easily implementable in standard treatment planning systems. Moderate unblocking provides a good way of dealing with OARs which abut the PTV. Optimization of geometry-based arcs is faster than usual inverse optimization of treatment plans, making this approach more rapid than an inverse-based Pareto front reconstruction.

The minimum knowledge base for predicting organ-at-risk dose-volume levels and plan-related complications in IMRT planning

IMRT treatment planning requires consideration of two competing objectives: achieving the required amount of radiation for the planning target volume and minimizing the amount of radiation delivered to all other tissues. It is important for planners to understand the tradeoff between competing factors so that the time-consuming human interaction loop (plan-evaluate-modify) can be eliminated. Treatment-plan-surface models have been proposed as a decision support tool to aid treatment planners and clinicians in choosing between rival treatment plans in a multi-plan environment. In this paper, an empirical approach is introduced to determine the minimum number of treatment plans (minimum knowledge base) required to build accurate representations of the IMRT plan surface in order to predict organ-at-risk (OAR) dose-volume (DV) levels and complications as a function of input DV constraint settings corresponding to all involved OARs in the plan. We have tested our approach on five head and neck patients and five whole pelvis/prostate patients. Our results suggest that approximately 30 plans were sufficient to predict DV levels with less than 3% relative error in both head and neck and whole pelvis/prostate cases. In addition, approximately 30-60 plans were sufficient to predict saliva flow rate with less than 2% relative error and to classify rectal bleeding with an accuracy of 90%.

Optimization of photon beam energy in aperture-based inverse planning

Optimal choice of beam energy in radiation therapy is easy in many well‐documented cases, but less obvious in some others. Low‐energy beams may provide better conformity around the target than their high‐energy counterparts due to reduced lateral scatter, but they also contribute to overdosage of peripheral normal tissue. Beam energy was added as an optimization parameter in an automatic aperture‐based inverse planning system. We have investigated a total of six cases for two sites (prostate and lung), representative of deep‐seated and moderately deep‐seated tumors. For one case for each site, different numbers of beam incidences were considered. The other cases for each site were optimized using a fixed number of incidences. Four types of plans were optimized: 6 MV, 23 MV, and mixed energy plans with one or two energies per incidence. Each plan was scored with a dose‐volume cost function. Cost function values, number of segments, monitor units, dose‐volume parameters, and isodose distributions were compared. For the prostate and lung cases, energy mixing improved plans in terms of cost function values, with a more important reduction for a small number of beam incidences. Use of high energy allowed better peripheral tissue sparing, while keeping similar target coverage and sensitive structures avoidance. Low energy contribution to monitor units usually increased with the number of beam incidences. Thus, for deep‐seated and moderately deep‐seated tumors, energy optimization can produce interesting plans with less peripheral dose and monitor units than for low energy alone.

PACS numbers: 87.55.de, 87.55.dk, 87.56.N‐

An approach for practical multiobjective IMRT treatment planning

PURPOSE

To introduce and demonstrate a practical multiobjective treatment planning procedure for intensity-modulated radiation therapy (IMRT) planning.

METHODS AND MATERIALS

The creation of a database of Pareto optimal treatment plans proceeds in two steps. The first step solves an optimization problem that finds a single treatment plan which is close to a set of clinical aspirations. This plan provides an example of what is feasible, and is then used to determine mutually satisfiable hard constraints for the subsequent generation of the plan database. All optimizations are done using linear programming.

RESULTS

The two-step procedure is applied to a brain, a prostate, and a lung case. The plan databases created allow for the selection of a final treatment plan based on the observed tradeoffs between the various organs involved.

CONCLUSIONS

The proposed method reduces the human iteration time common in IMRT treatment planning. Additionally, the database of plans, when properly viewed, allows the decision maker to make an informed final plan selection.