Optimal Multileaf Collimator Leaf Sequencing in IMRT Treatment Planning

We consider a problem dealing with the efficient delivery of intensity modulated radiation therapy (IMRT) to individual patients. IMRT treatment planning is usually performed in three phases. The first phase determines a set of beam angles through which radiation is delivered, followed by a second phase that determines an optimal radiation intensity profile (or fluence map). This intensity profile is selected to ensure that certain targets receive a required amount of dose while functional organs are spared. To deliver these intensity profiles to the patient, a third phase must decompose them into a collection of apertures and corresponding intensities. In this paper, we investigate this last problem. Formally, an intensity profile is represented as a nonnegative integer matrix; an aperture is represented as a binary matrix whose ones appear consecutively in each row. A feasible decomposition is one in which the original desired intensity profile is equal to the sum of a number of feasible binary matrices multiplied by corresponding intensity values. To most efficiently treat a patient, we wish to minimize a measure of total treatment time, which is given as a weighted sum of the number of apertures and the sum of the aperture intensities used in the decomposition. We develop the first exact algorithm capable of solving real-world problem instances to optimality within practicable computational limits, using a combination of integer programming decomposition and combinatorial search techniques. We demonstrate the efficacy of our approach on a set of 25 test instances derived from actual clinical data and on 100 randomly generated instances.

[1]  David S. Johnson,et al.  Computers and Intractability: A Guide to the Theory of NP-Completeness , 1978 .

[2]  T. Bortfeld,et al.  X-ray field compensation with multileaf collimators. , 1994, International journal of radiation oncology, biology, physics.

[3]  J. V. van Santvoort,et al.  Dynamic multileaf collimation without 'tongue-and-groove' underdosage effects. , 1996, Physics in medicine and biology.

[4]  P. Xia,et al.  Multileaf collimator leaf sequencing algorithm for intensity modulated beams with multiple static segments. , 1998, Medical physics.

[5]  W. Que.  Comparison of algorithms for multileaf collimator field segmentation. , 1999, Medical physics.

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

[7]  R. Siochi,et al.  Minimizing static intensity modulation delivery time using an intensity solid paradigm. , 1999, International journal of radiation oncology, biology, physics.

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

[9]  H. Hamacher,et al.  A mixed integer programming approach to the multileaf collimator problem , 2000 .

[10]  T. Pawlicki,et al.  The MLC tongue-and-groove effect on IMRT dose distributions , 2000, Proceedings of the 22nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (Cat. No.00CH37143).

[11]  Hanif D. Sherali,et al.  Improving Discrete Model Representations via Symmetry Considerations , 2001, Manag. Sci..

[12]  M. Langer,et al.  Improved leaf sequencing reduces segments or monitor units needed to deliver IMRT using multileaf collimators. , 2001, Medical physics.

[13]  J. Dai,et al.  Minimizing the number of segments in a delivery sequence for intensity-modulated radiation therapy with a multileaf collimator. , 2001, Medical physics.

[14]  D M Shepard,et al.  Direct aperture optimization: a turnkey solution for step-and-shoot IMRT. , 2002, Medical physics.

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

[16]  Sartaj Sahni,et al.  Leaf sequencing algorithms for segmented multileaf collimation. , 2003, Physics in medicine and biology.

[17]  C. Thieke,et al.  Intensity-Modulated Radiotherapy : A Large Scale Multi-Criteria Programming Problem Berichte , 2003 .

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

[19]  T. Solberg,et al.  Segmental and dynamic intensity-modulated radiotherapy delivery techniques for micro-multileaf collimator. , 2003, Medical physics.

[20]  Sartaj Sahni,et al.  Algorithms for optimal sequencing of dynamic multileaf collimators. , 2004, Physics in medicine and biology.

[21]  S. Sahni,et al.  A comparison of step-and-shoot leaf sequencing algorithms that eliminate tongue-and-groove effects. , 2004, Physics in medicine and biology.

[22]  Horst W. Hamacher,et al.  Minimizing beam‐on time in cancer radiation treatment using multileaf collimators , 2004, Networks.

[23]  S. Sahni,et al.  Optimal field splitting for large intensity-modulated fields. , 2004, Medical physics.

[24]  R. Siochi,et al.  Modifications to the IMFAST leaf sequencing optimization algorithm. , 2004, Medical physics.

[25]  Ronald L. Rardin,et al.  A coupled column generation, mixed integer approach to optimal planning of intensity modulated radiation therapy for cancer , 2004, Math. Program..

[26]  J Dai,et al.  'Tongue-and-groove' effect in intensity modulated radiotherapy with static multileaf collimator fields. , 2004, Physics in medicine and biology.

[27]  S. Sahni,et al.  Optimal leaf sequencing with elimination of tongue-and-groove underdosage. , 2004, Physics in medicine and biology.

[28]  Thomas Kalinowski,et al.  The algorithmic complexity of the minimization of the number of segments in multileaf collimator eld segmentation , 2004 .

[29]  Arvind Kumar,et al.  A Column Generation Approach to Radiation Therapy Treatment Planning Using Aperture Modulation , 2005, SIAM J. Optim..

[30]  Gerhard J. Woeginger,et al.  Decomposition of integer matrices and multileaf collimator sequencing , 2005, Discret. Appl. Math..

[31]  Thomas Kalinowski,et al.  A duality based algorithm for multileaf collimator field segmentation with interleaf collision constraint , 2005, Discret. Appl. Math..

[32]  T. Kalinowski.  Reducing the number of monitor units in multileaf collimator field segmentation , 2005, Physics in medicine and biology.

[33]  Davaatseren Baatar.  Matrix Decomposition with Times and Cardinality Objectives: Theory, Algorithms and Application to Multileaf Collimator Sequencing , 2005 .

[34]  Konrad Engel,et al.  A new algorithm for optimal multileaf collimator field segmentation , 2005, Discret. Appl. Math..

[35]  R. Ahuja,et al.  A network flow algorithm to minimize beam-on time for unconstrained multileaf collimator problems in cancer radiation therapy , 2005 .

[36]  Lauren Slater.  Tongue and Groove , 2006 .

[37]  Arvind Kumar,et al.  A New Linear Programming Approach to Radiation Therapy Treatment Planning Problems , 2006, Oper. Res..

[38]  Peter J. Stuckey,et al.  Minimum Cardinality Matrix Decomposition into Consecutive-Ones Matrices: CP and IP Approaches , 2007, CPAIOR.

[39]  R Alfredo C Siochi,et al.  Variable depth recursion algorithm for leaf sequencing. , 2007, Medical physics.

[40]  M. Ehrgott,et al.  Decomposition of matrices and static multileaf collimators: a survey , 2008 .

[41]  Andreas T. Ernst,et al.  An Exact Method for the Minimum Cardinality Problem in the Treatment Planning of Intensity-Modulated Radiotherapy , 2009, INFORMS J. Comput..