A New Linear Programming Approach to Radiation Therapy Treatment Planning Problems

We consider the problem of radiation therapy treatment planning for cancer patients. During radiation therapy, beams of radiation pass through a patient, killing both cancerous and normal cells. Thus, the radiation therapy must be carefully planned so that a clinically prescribed dose is delivered to targets containing cancerous cells, while nearby organs and tissues are spared. Currently, a technique called intensity-modulated radiation therapy (IMRT) is considered to be the most effective radiation therapy for many forms of cancer. In IMRT, the patient is irradiated from several beams, each of which is decomposed into hundreds of small beamlets, the intensities of which can be controlled individually. In this paper, we consider the problem of designing a treatment plan for IMRT when the orientations of the beams are given. We propose a new model that has the potential to achieve most of the goals with respect to the quality of a treatment plan that have been considered to date. However, in contrast with established mixed-integer and nonlinear programming formulations, we do so while retaining linearity of the optimization problem, which substantially improves the tractability of the optimization problem. Furthermore, we discuss how several additional quality and practical aspects of the problem that have been ignored to date can be incorporated into our linear model. We demonstrate the effectiveness of our approach on clinical data.

[1]  Horst W. Hamacher,et al.  Inverse Radiation Therapy Planning: A Multiple Objective Optimisation Approach , 1999 .

[2]  L. Xing,et al.  Iterative methods for inverse treatment planning. , 1996, Physics in medicine and biology.

[3]  J. Llacer,et al.  Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints. , 2003, Physics in medicine and biology.

[4]  R. Mohan,et al.  Multiple local minima in IMRT optimization based on dose-volume criteria. , 2002, Medical physics.

[5]  Cedric X. Yu,et al.  Guidance document on delivery, treatment planning, and clinical implementation of IMRT: report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee. , 2003, Medical physics.

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

[7]  P. Maddock Intensity modulated radiation therapy. , 2006, Medicine and health, Rhode Island.

[8]  J. Shapiro,et al.  Large scale optimization of beam weights under dose-volume restrictions. , 1990, International journal of radiation oncology, biology, physics.

[9]  Stan Uryasev,et al.  Conditional value-at-risk: optimization algorithms and applications , 2000, Proceedings of the IEEE/IAFE/INFORMS 2000 Conference on Computational Intelligence for Financial Engineering (CIFEr) (Cat. No.00TH8520).

[10]  Wlodzimierz Ogryczak,et al.  Conditional Median: A Parametric Solution Concept for Location Problems , 2002, Ann. Oper. Res..

[11]  S. Sutlief,et al.  Optimization of intensity modulated beams with volume constraints using two methods: cost function minimization and projections onto convex sets. , 1998, Medical physics.

[12]  Eva K. Lee,et al.  Optimization of radiosurgery treatment planning via mixed integer programming. , 2000, Medical physics.

[13]  S. Webb Optimisation of conformal radiotherapy dose distributions by simulated annealing. , 1989, Physics in medicine and biology.

[14]  Peter J. Slater,et al.  Centers to centroids in graphs , 1978, J. Graph Theory.

[15]  R. Rockafellar,et al.  Optimization of conditional value-at risk , 2000 .

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

[17]  G. Andreatta,et al.  k-eccentricity and absolute k-centrum of a probabilistic tree , 1985 .

[18]  T. Mackie,et al.  Iterative approaches to dose optimization in tomotherapy. , 2000, Physics in medicine and biology.

[19]  S Webb,et al.  There is no IMRT? , 2001, Physics in medicine and biology.

[20]  Robert J. Marks,et al.  Comparison of Algorithms for Intensity Modulated Beam Optimization : Projections Onto Convex Sets and Simulated Annealing , 2007 .

[21]  Arie Tamir,et al.  The k-centrum multi-facility location problem , 2001, Discret. Appl. Math..

[22]  G. Steel,et al.  Basic clinical radiobiology for radiation oncologists , 1993 .

[23]  T. Bortfeld,et al.  Methods of image reconstruction from projections applied to conformation radiotherapy. , 1990, Physics in medicine and biology.

[24]  I J Das,et al.  Treatment plan evaluation using dose-volume histogram (DVH) and spatial dose-volume histogram (zDVH). , 1999, International journal of radiation oncology, biology, physics.

[25]  T. Holmes,et al.  A filtered backprojection dose calculation method for inverse treatment planning. , 1994, Medical physics.

[26]  Eva K. Lee,et al.  Integer Programming Applied to Intensity-Modulated Radiation Therapy Treatment Planning , 2003, Ann. Oper. Res..

[27]  John Yarnold,et al.  Principles and practice of radiotherapy. , 2001 .

[28]  Fernando Alonso,et al.  Intensity-modulated radiotherapy – a large scale multi-criteria programming problem , 2003, OR Spectr..

[29]  S. Lippman,et al.  Head and Neck Cancer , 1993, Cancer treatment and research.

[30]  Michael C. Ferris,et al.  Optimizing the Delivery of Radiation Therapy to Cancer Patients , 1999, SIAM Rev..

[31]  S. Webb Optimization by simulated annealing of three-dimensional, conformal treatment planning for radiation fields defined by a multileaf collimator: II. Inclusion of two-dimensional modulation of the x-ray intensity. , 1992, Physics in medicine and biology.

[32]  Avraham Eisbruch,et al.  Intensity-modulated radiation therapy (IMRT) for locally advanced paranasal sinus tumors: incorporating clinical decisions in the optimization process. , 2003, International journal of radiation oncology, biology, physics.

[33]  R Mohan,et al.  Application of fast simulated annealing to optimization of conformal radiation treatments. , 1993, Medical physics.

[34]  Joseph O. Deasy,et al.  Treatment plan review tools incorporating spatial dose-volume information , 2001 .

[35]  Wlodzimierz Ogryczak,et al.  Minimizing the sum of the k largest functions in linear time , 2003, Inf. Process. Lett..

[36]  Stanislav Uryasev,et al.  Conditional Value-at-Risk for General Loss Distributions , 2002 .

[37]  C. Fiorino,et al.  Optimizing the movement of a single absorber for 1D non-uniform dose delivery by (fast) simulated annealing. , 1997, Physics in medicine and biology.

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

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

[40]  Robert Jeraj,et al.  The method of intercepts in parameter space for the analysis of local minima caused by dose-volume constraints. , 2003, Physics in medicine and biology.

[41]  J. Dempsey,et al.  Validation of a precision radiochromic film dosimetry system for quantitative two-dimensional imaging of acute exposure dose distributions. , 2000, Medical physics.

[42]  R Mohan,et al.  Algorithms and functionality of an intensity modulated radiotherapy optimization system. , 2000, Medical physics.

[43]  T. Bortfeld,et al.  Inverse planning for photon and proton beams. , 2001, Medical dosimetry : official journal of the American Association of Medical Dosimetrists.