Inverse radiation treatment planning using the Dynamically Penalized Likelihood method.

In this paper we present a new method of solving the inverse radiation treatment planning problem. The method is based on a Maximum Likelihood Estimator with dynamically changing penalization terms. The resulting Dynamically Penalized Likelihood (DPL) algorithm achieves a dose distribution of excellent uniformity in a tumor volume and a much lower dose in regions containing sensitive volumes. A simple model of a patient and of energy deposition has been used for the initial results presented: a two-dimensional computer generated phantom and monochromatic x rays, without scattering. Three two-dimensional problems are solved with the DPL algorithm, corresponding to different size and spatial relationships between the tumor and sensitive tissue volumes. The results show that the DPL algorithm is robust and flexible; it only requires moderate computation times and leads to promising solutions, even in rather difficult problems. The results encourage the extension of the present work to more realistic therapy situations.

[1]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[2]  L. Lucy An iterative technique for the rectification of observed distributions , 1974 .

[3]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

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

[5]  M D Altschuler,et al.  Semi-automated radiotherapy treatment planning with a mathematical model to satisfy treatment goals. , 1989, International journal of radiation oncology, biology, physics.

[6]  T. Bortfeld,et al.  Optimization of beam orientations in radiation therapy: some theoretical considerations. , 1993, Physics in medicine and biology.

[7]  B. M. ter Haar Romeny,et al.  Representation of medical images by visual response functions , 1993, IEEE Engineering in Medicine and Biology Magazine.

[8]  Jorge Llacer,et al.  A GENERAL BAYESIAN IMAGE RECONSTRUCTION ALGORITHM WITH ENTROPY PRIOR. PRELIMINARY APPLICATION TO HST DATA , 1993 .

[9]  Jorge Núñez de Murga,et al.  Statistical analysis of maximum likelihood estimator images of human brain FDG PET studies , 1993, IEEE Trans. Medical Imaging.

[10]  A. Brahme,et al.  A generalized pencil beam algorithm for optimization of radiation therapy. , 1994, Medical physics.

[11]  T R Mackie,et al.  An iterative filtered backprojection inverse treatment planning algorithm for tomotherapy. , 1995, International journal of radiation oncology, biology, physics.

[12]  Tariq S. Durrani,et al.  Conformal therapy using maximum entropy optimization , 1995, Int. J. Imaging Syst. Technol..

[13]  Thomas Bortfeld,et al.  The exponential radon transform and projection filtering in radiotherapy planning , 1995, Int. J. Imaging Syst. Technol..

[14]  Steve Webb Optimizing radiation therapy inverse treatment planning using the simulated annealing technique , 1995, Int. J. Imaging Syst. Technol..