Application of constrained least-squares techniques to IMRT treatment planning.

PURPOSE The purpose of this work was to apply the method of constrained least-squares to inverse treatment planning and to explore its potential for providing a fast interactive planning environment for intensity-modulated radiation therapy (IMRT). METHODS AND MATERIALS The description of the dose inside a patient is a linear matrix transformation of beamlet weights. The constrained least-squares method adds additional matrix operators and produces beamlet weights by a direct linear transformation. These matrix operators contain a priori knowledge about the radiation distribution. The constrained least-squares technique was applied to obtain IMRT plans for prostate and paraspinal cancer patients and compared with the corresponding plans optimized using the CORVUS inverse planning system. RESULTS It was demonstrated that a constrained least-squares technique is suitable for IMRT plan optimization with significantly increased computing speed. For the two cases we have tested, the constrained least-squares method was an order of magnitude faster than conventional iterative techniques because of the avoidance of the iterative calculations. We also found that the constrained least-squares method is capable of generating clinically acceptable treatment plans with less trial-and-error adjustments of system variables, and with improved target volume coverage as well as sensitive structure sparing in comparison with that obtained using CORVUS. CONCLUSIONS The constrained least-squares method has the advantage that it does not require iterative calculation and thus significantly speeds up the therapeutic plan optimization process. Besides shedding important insight into the inverse planning problem, the technique has strong potential to provide a fast and interactive environment for IMRT treatment planning.

[1]  Steven J. Leon Linear Algebra With Applications , 1980 .

[2]  Gerhard Winkler,et al.  Image analysis, random fields and dynamic Monte Carlo methods: a mathematical introduction , 1995, Applications of mathematics.

[3]  T. Holmes,et al.  A filtered backprojection dose calculation method for inverse treatment planning. , 1994, Medical physics.

[4]  B G Fallone,et al.  An active set algorithm for treatment planning optimization. , 1997, Medical physics.

[5]  D M Shepard,et al.  Application of constrained optimization to radiotherapy planning. , 1999, Medical physics.

[6]  A. Boyer,et al.  Beam orientation optimization in intensity-modulated radiation treatment planning. , 2000, Medical physics.

[7]  K S Lam,et al.  Comparison of simulated annealing algorithms for conformal therapy treatment planning. , 1995, International journal of radiation oncology, biology, physics.

[8]  G Starkschall,et al.  A constrained least-squares optimization method for external beam radiation therapy treatment planning. , 1984, Medical physics.

[9]  A L Boyer,et al.  Optimization of importance factors in inverse planning. , 1999, Physics in medicine and biology.

[10]  T. Bortfeld,et al.  Methods of image reconstruction from projections applied to conformation radiotherapy. , 1990, Physics in medicine and biology.

[11]  J. Roeske,et al.  Initial clinical experience with intensity-modulated whole-pelvis radiation therapy in women with gynecologic malignancies. , 2001, Gynecologic oncology.

[12]  F. R. Gantmakher The Theory of Matrices , 1984 .

[13]  L. Mirsky,et al.  The Theory of Matrices , 1961, The Mathematical Gazette.

[14]  T. Mackie,et al.  Iterative approaches to dose optimization in tomotherapy. , 2000, Physics in medicine and biology.

[15]  G T Chen,et al.  Fast iterative algorithms for three-dimensional inverse treatment planning. , 1998, Medical physics.

[16]  J. G. Li,et al.  Estimation theory and model parameter selection for therapeutic treatment plan optimization. , 1999, Medical physics.

[17]  I. Rosen,et al.  Very fast simulated reannealing in radiation therapy treatment plan optimization. , 1995, International journal of radiation oncology, biology, physics.

[18]  L. Xing,et al.  Iterative methods for inverse treatment planning. , 1996, Physics in medicine and biology.

[19]  Steven Kay,et al.  Fundamentals Of Statistical Signal Processing , 2001 .

[20]  G Starkschall,et al.  An interactive beam-weight optimization tool for three-dimensional radiotherapy treatment planning. , 1992, Medical physics.

[21]  R Mohan,et al.  Conformal radiation treatment of prostate cancer using inversely-planned intensity-modulated photon beams produced with dynamic multileaf collimation. , 1996, International journal of radiation oncology, biology, physics.

[22]  I I Rosen,et al.  Dose-volume considerations with linear programming optimization. , 1991, Medical physics.

[23]  S. Webb Optimisation of conformal radiotherapy dose distributions by simulated annealing. , 1989, Physics in medicine and biology.

[24]  J. Llacer Inverse radiation treatment planning using the Dynamically Penalized Likelihood method. , 1997, Medical physics.

[25]  Michael C. Ferris,et al.  Optimizing the Delivery of Radiation Therapy to Cancer Patients , 1999, SIAM Rev..

[26]  P Xia,et al.  Comparison of treatment plans involving intensity-modulated radiotherapy for nasopharyngeal carcinoma. , 2000, International journal of radiation oncology, biology, physics.

[27]  R. Viswanathan,et al.  An introduction to statistical signal processing with applications , 1979 .

[28]  S. Kay Fundamentals of statistical signal processing: estimation theory , 1993 .