Interactively exploring optimized treatment plans.

PURPOSE A new paradigm for treatment planning is proposed that embodies the concept of interactively exploring the space of optimized plans. In this approach, treatment planning ignores the details of individual plans and instead presents the physician with clinical summaries of sets of solutions to well-defined clinical goals in which every solution has been optimized in advance by computer algorithms. METHODS AND MATERIALS Before interactive planning, sets of optimized plans are created for a variety of treatment delivery options and critical structure dose-volume constraints. Then, the dose-volume parameters of the optimized plans are fit to linear functions. These linear functions are used to show in real time how the target dose-volume histogram (DVH) changes as the DVHs of the critical structures are changed interactively. A bitmap of the space of optimized plans is used to restrict the feasible solutions. The physician selects the critical structure dose-volume constraints that give the desired dose to the planning target volume (PTV) and then those constraints are used to create the corresponding optimized plan. RESULTS The method is demonstrated using prototype software, Treatment Plan Explorer (TPEx), and a clinical example of a patient with a tumor in the right lung. For this example, the delivery options included 4 open beams, 12 open beams, 4 wedged beams, and 12 wedged beams. Beam directions and relative weights were optimized for a range of critical structure dose-volume constraints for the lungs and esophagus. Cord dose was restricted to 45 Gy. Using the interactive interface, the physician explored how the tumor dose changed as critical structure dose-volume constraints were tightened or relaxed and selected the best compromise for each delivery option. The corresponding treatment plans were calculated and compared with the linear parameterization presented to the physician in TPEx. The linear fits were best for the maximum PTV dose and worst for the minimum PTV dose. Based on the root-mean-square error between the fit values and their corresponding data values, the linear fit appears to be adequate, although higher order polynomials could give better results. Some of the variance in fit is due to the stochastic nature of the simulated annealing optimization algorithm, which does not reproduce the exact same results in repetitions of the same calculation. Using a directed search algorithm for plan optimization should produce better parameter fits and, therefore, better predictions of plan characteristics by TPEx. CONCLUSIONS Using TPEx, the physician can easily select the optimum plan for a patient, with no imposed arbitrary definition of the "best" plan. More importantly, the physician can readily see what can be achieved for the patient with a given delivery technique. There is no more uncertainty about whether or not a better plan exists. By comparing the "best" plans for different delivery options (e.g., three-dimensional conformal radiotherapy versus intensity-modulated radiation therapy), the physician can gauge the clinical benefits of greater technical complexity. However, before the TPEx process can be clinical useful, faster computers and/or algorithms are needed and more studies are needed to better model the spaces of optimized solutions.

[1]  I I Rosen,et al.  Treatment planning optimization using constrained simulated annealing. , 1991, Physics in medicine and biology.

[2]  K S Lam,et al.  Comparison of simulated annealing algorithms for conformal therapy treatment planning. , 1995, International journal of radiation oncology, biology, physics.

[3]  L. Xing,et al.  Incorporating prior knowledge into beam orientation optimization in IMRT. , 2002, International journal of radiation oncology, biology, physics.

[4]  I. Rosen,et al.  Collimator scatter in modeling radiation beam profiles. , 1990, Medical physics.

[5]  D Baltas,et al.  A multiobjective gradient-based dose optimization algorithm for external beam conformal radiotherapy. , 2001, Physics in medicine and biology.

[6]  C. Cotrutz,et al.  Segment-based dose optimization using a genetic algorithm. , 2003, Physics in medicine and biology.

[7]  R. Lane,et al.  A generic genetic algorithm for generating beam weights. , 1996, Medical physics.

[8]  Radhe Mohan,et al.  Algorithm and performance of a clinical IMRT beam-angle optimization system. , 2003, Physics in medicine and biology.

[9]  S. Spirou,et al.  A gradient inverse planning algorithm with dose-volume constraints. , 1998, Medical physics.

[10]  Y. Yu,et al.  Decision theoretic steering and genetic algorithm optimization: application to stereotactic radiosurgery treatment planning. , 1997, Medical physics.

[11]  J. Deasy Multiple local minima in radiotherapy optimization problems with dose-volume constraints. , 1997, Medical physics.

[12]  Ying Xiao,et al.  The use of mixed-integer programming for inverse treatment planning with pre-defined field segments. , 2001, Physics in medicine and biology.

[13]  R. Lane,et al.  A comparison of mixed integer programming and fast simulated annealing for optimizing beam weights in radiation therapy. , 1996, Medical physics.

[14]  A L Boyer,et al.  Optimization of importance factors in inverse planning. , 1999, Physics in medicine and biology.

[15]  L. Xing,et al.  Multiobjective evolutionary optimization of the number of beams, their orientations and weights for intensity-modulated radiation therapy , 2004, Physics in Medicine and Biology.

[16]  Eva K. Lee,et al.  Operations research applied to radiotherapy, an NCI-NSF-sponsored workshop February 7-9, 2002. , 2003, International journal of radiation oncology, biology, physics.

[17]  G Starkschall,et al.  Evaluation and scoring of radiotherapy treatment plans using an artificial neural network. , 1996, International journal of radiation oncology, biology, physics.

[18]  Lei Dong,et al.  Speed and convergence properties of gradient algorithms for optimization of IMRT. , 2004, Medical physics.

[19]  Yair Censor,et al.  The dose–volume constraint satisfaction problem for inverse treatment planning with field segments , 2004, Physics in medicine and biology.

[20]  Jürgen Meyer,et al.  Application of influence diagrams to prostate intensity-modulated radiation therapy plan selection , 2004, Physics in medicine and biology.

[21]  C G Rowbottom,et al.  Simultaneous optimization of beam orientations and beam weights in conformal radiotherapy. , 2001, Medical physics.

[22]  A. Pollack,et al.  Treatment planning optimization based on response-surface modeling of cost function versus multiple constraints , 2000, Proceedings of the 22nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (Cat. No.00CH37143).

[23]  George Starkschall,et al.  Plan space: representation of treatment plans in multidimensional space. , 2002, International journal of radiation oncology, biology, physics.

[24]  T. Bortfeld,et al.  Number and orientations of beams in intensity-modulated radiation treatments. , 1997, Medical physics.

[25]  Radhe Mohan,et al.  Intensity-modulated radiotherapy optimization with gEUD-guided dose-volume objectives. , 2003, Physics in medicine and biology.

[26]  Steve Webb,et al.  Elimination of importance factors for clinically accurate selection of beam orientations, beam weights and wedge angles in conformal radiation therapy. , 2003, Medical physics.

[27]  G Starkschall,et al.  Treatment planning using a dose-volume feasibility search algorithm. , 2001, International journal of radiation oncology, biology, physics.

[28]  I. Rosen,et al.  Treatment plan optimization using linear programming. , 1991, Medical physics.

[29]  M. Alber,et al.  Tools for the analysis of dose optimization: II. Sensitivity analysis , 2002, Physics in medicine and biology.