Dose-volume constrained optimization in intensity-modulated radiation therapy treatment planning

We present a novel optimization method to handle dose-volume constraints (DVCs) directly in intensity-modulated radiation therapy (IMRT) treatment planning based on the idea of continuous dynamical methods. Most of the conventional methods are constructed for solving inconsistent inverse problems with, e.g., dose-volume based objective functions, and one expects to obtain a feasible solution that minimally violates the DVCs. We introduce the concept of ‘acceptable’, meaning that there exists a nonempty set of radiation beam weights satisfying the given DVCs, and we resolve the issue that the objective and evaluation are different in the conventional planning approach. We apply the initial-value problem of the proposed dynamical system to an acceptable and inconsistent inverse problem and prove that the convergence to an equilibrium in the acceptable set of solutions is theoretically guaranteed by using the Lyapunov theorem. Indeed, we confirmed that we can obtain acceptable beam weights through numerical experiments using phantom data simulating a clinical setup for an acceptable and inconsistent IMRT planning system.

[1]  Q. Hou,et al.  An optimization algorithm for intensity modulated radiotherapy--the simulated dynamics with dose-volume constraints. , 2002, Medical physics.

[2]  A. M. Li︠a︡punov Stability of Motion , 2016 .

[3]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[4]  Y. Censor,et al.  A unified approach for inversion problems in intensity-modulated radiation therapy , 2006, Physics in medicine and biology.

[5]  C. Byrne,et al.  Iterative oblique projection onto convex sets and the split feasibility problem , 2002 .

[6]  Johannes Schropp,et al.  Using dynamical systems methods to solve minimization problems , 1995 .

[7]  Rajat J Kudchadker,et al.  A comprehensive comparison of IMRT and VMAT plan quality for prostate cancer treatment. , 2012, International journal of radiation oncology, biology, physics.

[8]  N. Wouw,et al.  Stability and Convergence of Mechanical Systems with Unilateral Constraints , 2008 .

[9]  Bo Han,et al.  A dynamical system method for solving nonlinear ill-posed problems , 2008, Appl. Math. Comput..

[10]  George A. Anastassiou Applied mathematics reviews , 2000 .

[11]  Michael Lahanas,et al.  Multiobjective inverse planning for intensity modulated radiotherapy with constraint-free gradient-based optimization algorithms. , 2003, Physics in medicine and biology.

[12]  Tetsuya Yoshinaga,et al.  Continuous-time method and its discretization to inverse problem of intensity-modulated radiation therapy treatment planning , 2014, Commun. Nonlinear Sci. Numer. Simul..

[13]  M J Ackerman,et al.  The Visible Human Project , 1998, Proc. IEEE.

[14]  Omar Mohammad Abdel Latif Abou Al Ola Common Lyapunov function based on Kullback-Leibler divergence for a switched nonlinear system , 2011 .

[15]  Yoshiya Yamada,et al.  Ultra-high dose (86.4 Gy) IMRT for localized prostate cancer: toxicity and biochemical outcomes. , 2008, International journal of radiation oncology, biology, physics.

[16]  Tetsuya Yoshinaga,et al.  Continuous-time image reconstruction for binary tomography , 2013, Commun. Nonlinear Sci. Numer. Simul..

[17]  I. Kawrakow Improved modeling of multiple scattering in the Voxel Monte Carlo model. , 1997, Medical physics.

[18]  B. G. Fallone,et al.  Intensity-Modulated Radiation Therapy: The State of the Art , 2004 .

[19]  T. Bortfeld IMRT: a review and preview , 2006, Physics in medicine and biology.

[20]  D Baltas,et al.  A multiobjective gradient-based dose optimization algorithm for external beam conformal radiotherapy. , 2001, Physics in medicine and biology.

[21]  Alexander G. Ramm,et al.  Dynamical systems method for solving operator equations , 2004 .

[22]  I. Csiszár Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems , 1991 .

[23]  E. Yorke,et al.  Reduced-order constrained optimization in IMRT planning , 2008, Physics in medicine and biology.

[24]  Tetsuya Yoshinaga,et al.  Continuous-time image reconstruction using differential equations for computed tomography , 2010 .

[25]  Poom Kumam,et al.  Split feasibility and fixed-point problems for asymptotically quasi-nonexpansive mappings , 2013 .

[26]  R Mohan,et al.  Algorithms and functionality of an intensity modulated radiotherapy optimization system. , 2000, Medical physics.

[27]  A. Fuller,et al.  Stability of Motion , 1976, IEEE Transactions on Systems, Man, and Cybernetics.

[28]  Alexander G. Ramm,et al.  DYNAMICAL SYSTEMS AND DISCRETE METHODS FOR SOLVING NONLINEAR ILL-POSED PROBLEMS , 2000 .

[29]  Alexander G. Ramm,et al.  CONTINUOUS ANALOG OF THE GAUSS-NEWTON METHOD , 1999 .

[30]  Lei Dong,et al.  Speed and convergence properties of gradient algorithms for optimization of IMRT. , 2004, Medical physics.

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

[32]  C. Pelizzari,et al.  Late toxicity after intensity-modulated radiation therapy for localized prostate cancer: an exploration of dose-volume histogram parameters to limit genitourinary and gastrointestinal toxicity. , 2012, International journal of radiation oncology, biology, physics.

[33]  Ruban Airapetyan,et al.  Continuous newton method and its modification , 1999 .

[34]  Yair Censor,et al.  A multiprojection algorithm using Bregman projections in a product space , 1994, Numerical Algorithms.