This paper reports on the analysis of intensity modulated radiation treatment optimization problems in the presence of non-convex feasible parameter spaces caused by the specification of dose-volume constraints for the organs-at-risk (OARs). The main aim was to determine whether the presence of those non-convex spaces affects the optimization of clinical cases in any significant way. This was done in two phases: (1) Using a carefully designed two-dimensional mathematical phantom that exhibits two controllable minima and with randomly initialized beamlet weights, we developed a methodology for exploring the nature of the convergence characteristics of quadratic cost function optimizations (deterministic or stochastic). The methodology is based on observing the statistical behaviour of the residual cost at the end of optimizations in which the stopping criterion is progressively more demanding and carrying out those optimizations to very small error changes per iteration. (2) Seven clinical cases were then analysed with dose-volume constraints that are stronger than originally used in the clinic. The clinical cases are two prostate cases differently posed, a meningioma case, two head-and-neck cases, a spleen case and a spine case. Of the 14 different sets of optimizations (with and without the specification of maximum doses allowed for the OARs), 12 fail to show any effect due to the existence of non-convex feasible spaces. The remaining two sets of optimizations show evidence of multiple minima in the solutions, but those minima are very close to each other in cost and the resulting treatment plans are practically identical, as measured by the quality of the dose-volume histograms (DVHs). We discuss the differences between fluence maps resulting from those similar treatment plans. We provide a possible reason for the observed results and conclude that, although the study is necessarily limited, the annealing characteristics of a simulated annealing method may not be justified in clinical optimization in the presence of dose-volume constraints. The results of optimizations by the Newton gradient (NG) method with a quadratic cost function are reported in detail. An adaptive simulated annealing method, optimizing the same function, and the dynamically penalized likelihood method, optimizing a log likelihood function, have also been used in the study. The results of the latter two methods have only been discussed briefly, as they yielded the same conclusions as the NG method.

This paper presents a description of tests carried out to compare the behaviour of five algorithms in inverse radiation therapy planning: (1) The Dynamically Penalized Likelihood (DPL), an algorithm based on statistical estimation theory; (2) an accelerated version of the same algorithm: (3) a new fast adaptive simulated annealing (ASA) algorithm; (4) a conjugate gradient method; and (5) a Newton gradient method. A three-dimensional mathematical phantom and two clinical cases have been studied in detail. The phantom consisted of a U-shaped tumour with a partially enclosed 'spinal cord'. The clinical examples were a cavernous sinus meningioma and a prostate case. The algorithms have been tested in carefully selected and controlled conditions so as to ensure fairness in the assessment of results. It has been found that all five methods can yield relatively similar optimizations, except when a very demanding optimization is carried out. For the easier cases. the differences are principally in robustness, ease of use and optimization speed. In the more demanding case, there are significant differences in the resulting dose distributions. The accelerated DPL emerges as possibly the algorithm of choice for clinical practice. An appendix describes the differences in behaviour between the new ASA method and the one based on a patent by the Nomos Corporation.

This paper proposes a method for automatic selection of beam orientations in non-coplanar cranial IMRT. Methods of computer vision, beam's eye view techniques and neural networks are used to define a new geometry-based methodology that leads to treatment plans for cranial lesions that are comparable in quality to those generated by experienced radiation physicists. The automatic beam selection (ABS) process can be carried out in clinically useful computation times, in 1 min or less for most cases. In the process of describing the ABS process, it is shown that the cranial beam orientation optimization problem is mathematically ill posed, with the expectation that a large number of solutions will lead to similar results. Nevertheless, there are better and worse solutions and we show that the proposed ABS process, by its design, has to lead to one of the better ones. We have carried out extensive tests with 14 patients with beam selection tasks ranging from the rather simple to quite complex. The ABS process has always yielded optimizations with results that are considered good for clinic use. Seven-beam coplanar optimizations for some of the patients have also been investigated. Comparisons with non-coplanar optimizations indicate in which cases the simpler coplanar plans can be used to advantage. Parameters used in the comparisons are dose-volume histograms, minimum and maximum PTV doses, equivalent uniform doses for the PTV and OARs, and treatment volume, conformity and normal tissue indices. It is felt that the current ABS methodology is ready for extensive clinical tests.

This paper attempts to provide an answer to some questions that remain either poorly understood, or not well documented in the literature, on basic issues related to intensity modulated radiation therapy (IMRT). The questions examined are: the relationship between degeneracy and frequency response of optimizations, effects of initial beamlet fluence assignment and stopping point, what does filtering of an optimized beamlet map actually do and how could image analysis help to obtain better optimizations? Two target functions are studied, a quadratic cost function and the log likelihood function of the dynamically penalized likelihood (DPL) algorithm. The algorithms used are the conjugate gradient, the stochastic adaptive simulated annealing and the DPL. One simple phantom is used to show the development of the analysis tools used and two clinical cases of medium and large dose matrix size (a meningioma and a prostate) are studied in detail. The conclusions reached are that the high number of iterations that is needed to avoid degeneracy is not warranted in clinical practice, as the quality of the optimizations, as judged by the DVHs and dose distributions obtained, does not improve significantly after a certain point. It is also shown that the optimum initial beamlet fluence assignment for analytical iterative algorithms is a uniform distribution, but such an assignment does not help a stochastic method of optimization. Stopping points for the studied algorithms are discussed and the deterioration of DVH characteristics with filtering is shown to be partially recoverable by the use of space-variant filtering techniques.

No abstract