Non-coplanar beam direction optimization for intensity-modulated radiotherapy.

An algorithm for the optimization of the direction of intensity-modulated beams is presented. Although the global optimum dose distribution cannot be predicted, usually a large number of equivalent beam configurations exists. This degeneracy facilitates beam direction optimization (BDO) through a number of possible approximations and because the target set of good beam configurations is very large. Usually, the target volume is accessible through a finite number of paths of little resistance, which are defined by the properties of the objective function and the global optimum dose distribution. Since these paths can be occupied by a finite number of beams, it is reasonable to assume that a minimum number of beams for a configuration that is degenerate to the global optimum exists. Efficiency of the BDO will be characterized by detecting this degeneracy threshold. Beam configurations are altered by adding and deleting beams. A fast exhaustive (up to 3500 non-coplanar orientations) search finds beam directions that improve a configuration. Redundant beams of a configuration can be identified by a fast criterion based on second-order derivative information of the objective function. This offers a fast means of iteratively substituting redundant beams from a configuration. Inferior stationary states can be evaded by adding more beams than the desired number to the current configuration, followed by the subsequent cancellation of superfluous beams. The significance of BDO is examined in a coplanar and a non-coplanar test case. The existence of a threshold number for the minimum configuration and its dependence on the complexity of the problem are shown. BDO outperforms manual configurations and equispaced coplanar beam arrangements in both example cases.

On the degeneracy of the IMRT optimization problem.

One approach to the computation of photon IMRT treatment plans is the formulation of an optimization problem with an objective function that derives from an objective density. An investigation of the second-order properties of such an objective function in a neighborhood of the minimizer opens an intuitive access to many traits of this approach. A general finding is that only a small subset of the parameter space has nonzero curvature, while the objective function is entirely flat in a neighborhood of the minimizer in most directions. The dimension of the subspace of vanishing curvature serves as a measure for the degeneracy of the solution. This finding is important both for algorithm design and evaluation of the mathematical model of clinical intuition, expressed by the objective function. The structure of the subspace of great curvature is found to be imposed on the problem by conflicts between objectives of target and critical structures. These conflicts and their corresponding modes of resolution form a common trait between all reasonable treatment plans of a given case. The high degree of degeneracy makes the use of a conjugate gradient optimization algorithm particularly favorable since the number of iterations to convergence is equivalent to the number of different eigenvalues of the curvature tensor and is hence independent from the number of optimization parameters. A high level of degeneracy of the fluence profiles implies that it should be possible to stipulate further delivery-related conditions without causing severe deterioration of the dose distribution.