We present computational approaches for optimizing beam angles and fluence maps in Intensity Modulated Radiation Therapy (IMRT) planning. We assume that the number of angles to be used for the treatment is given by the treatment planner. A mixed integer programming (MIP) model and a linear programming (LP) model are used to find an optimal set of beam angles and their corresponding fluence maps. The MIP model is solved using the branch-and-bound method while the LP model is solved using the interior point method. In order to reduce the computational burden for solving the optimization models, we introduce iterative beam angle elimination algorithms in which an insignificant beam angle is eliminated in each iteration. Other techniques are also explored including feasible set reduction for LP and data reduction. Experiments are made to show the computational advantage of the iterative methods for optimizing angles using real patient cases.
This paper has been accepted for publication. OR Spectrum. 2007. IE Technical Report 10-01, October 2006, Department of Industrial Engineering, University of Houston