An exact approach to direct aperture optimization in IMRT treatment planning

We consider the problem of intensity-modulated radiation therapy (IMRT) treatment planning using direct aperture optimization. While this problem has been relatively well studied in recent years, most approaches employ a heuristic approach to the generation of apertures. In contrast, we use an exact approach that explicitly formulates the fluence map optimization (FMO) problem as a convex optimization problem in terms of all multileaf collimator (MLC) deliverable apertures and their associated intensities. However, the number of deliverable apertures, and therefore the number of decision variables and constraints in the new problem formulation, is typically enormous. To overcome this, we use an iterative approach that employs a subproblem whose optimal solution either provides a suitable aperture to add to a given pool of allowable apertures or concludes that the current solution is optimal. We are able to handle standard consecutiveness, interdigitation and connectedness constraints that may be imposed by the particular MLC system used, as well as jaws-only delivery. Our approach has the additional advantage that it can explicitly account for transmission of dose through the part of an aperture that is blocked by the MLC system, yielding a more precise assessment of the treatment plan than what is possible using a traditional beamlet-based FMO problem. Finally, we develop and test two stopping rules that can be used to identify treatment plans of high clinical quality that are deliverable very efficiently. Tests on clinical head-and-neck cancer cases showed the efficacy of our approach, yielding treatment plans comparable in quality to plans obtained by the traditional method with a reduction of more than 75% in the number of apertures and a reduction of more than 50% in beam-on time, with only a modest increase in computational effort. The results also show that delivery efficiency is very insensitive to the addition of traditional MLC constraints; however, jaws-only treatment requires about a doubling in beam-on time and number of apertures used. Finally, we showed the importance of accounting for transmission effects when assessing or, preferably, optimizing treatment plan quality.