Degeneracy, frequency response and filtering in IMRT optimization.

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.

[1]  Robert Jeraj,et al.  The method of intercepts in parameter space for the analysis of local minima caused by dose-volume constraints. , 2003, Physics in medicine and biology.

[2]  M Gambaccini,et al.  Dual-energy tissue cancellation in mammography with quasi-monochromatic x-rays. , 2002, Physics in medicine and biology.

[3]  F Nüsslin,et al.  Intensity modulated photon beams subject to a minimal surface smoothing constraint. , 2000, Physics in medicine and biology.

[4]  T D Solberg,et al.  Comparative behaviour of the dynamically penalized likelihood algorithm in inverse radiation therapy planning. , 2001, Physics in medicine and biology.

[5]  Jorge Núñez de Murga,et al.  Statistically based image reconstruction for emission tomography , 1989, Int. J. Imaging Syst. Technol..

[6]  J. Llacer,et al.  Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints. , 2003, Physics in medicine and biology.

[7]  J. Deasy.  Multiple local minima in radiotherapy optimization problems with dose-volume constraints. , 1997, Medical physics.

[8]  C S Chui,et al.  Inverse planning algorithms for external beam radiation therapy. , 2001, Medical dosimetry : official journal of the American Association of Medical Dosimetrists.

[9]  Steve Webb,et al.  Configuration space analysis of common cost functions in radiotherapy beam-weight optimization algorithms. , 2002, Physics in medicine and biology.

[10]  R. Mohan,et al.  Multiple local minima in IMRT optimization based on dose-volume criteria. , 2002, Medical physics.

[11]  Bart M. ter Haar Romeny,et al.  Front-End Vision and Multi-Scale Image Analysis , 2003, Computational Imaging and Vision.

[12]  J. Llacer.  Tomographic Image Reconstruction by Eigenvector Decomposition: Its Limitations and Areas of Applicability , 1982 .

[13]  M. Alber,et al.  On the degeneracy of the IMRT optimization problem. , 2002, Medical physics.

[14]  S. Spirou,et al.  A gradient inverse planning algorithm with dose-volume constraints. , 1998, Medical physics.

[15]  R. Jeraj,et al.  Optimizer convergence and local minima errors and their clinical importance. , 2003, Physics in medicine and biology.