Fast Inverse Planning of Beam Directions and Weights for Small Animal Radiotherapy

Current preclinical irradiation systems, such as the small animal radiation research platform, provide many capabilities of clinical systems, including on-board cone beam computed tomography (CBCT) imaging, which motivates the development of treatment planning systems. The time for preclinical treatment planning is often limited, however, due to the large number of subjects and the need for anesthesia for imaging and radiation delivery. This paper presents a 3-D inverse planning solution that optimizes beam directions and weights within a timeframe suitable for small animal radiation research. The system begins with a large number of beams and takes advantage of a graphics processing unit-accelerated implementation of the superposition-convolution method to quickly compute the dose for each beam. Optimization is performed by a linear programming solver, with a minimum bound on dose to the target and a maximum bound on dose to organs at risk (OAR). It is extended to dose shells by using hollow cylinders instead of solid beams. The method is demonstrated on mouse CBCT images with arbitrarily defined targets and OAR. Inverse planning is performed in about 5 min, which enables a preclinical workflow that includes CBCT acquisition, treatment planning, and radiation delivery in a single session.

[1]  I I Rosen,et al.  Treatment planning optimization using constrained simulated annealing. , 1991, Physics in medicine and biology.

[2]  Horst W. Hamacher,et al.  Mathematical optimization in intensity modulated radiation therapy , 2010, Ann. Oper. Res..

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

[4]  M Monz,et al.  Pareto navigation—algorithmic foundation of interactive multi-criteria IMRT planning , 2008, Physics in medicine and biology.

[5]  Russell H. Taylor,et al.  Towards real-time radiation therapy: GPU accelerated superposition/convolution , 2010, Comput. Methods Programs Biomed..

[6]  Peter Kazanzides,et al.  Image-guided small animal radiation research platform: calibration of treatment beam alignment , 2009, Physics in medicine and biology.

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

[8]  A T Redpath,et al.  A new technique for radiotherapy planning using quadratic programming. , 1976, Physics in medicine and biology.

[9]  I. Rosen,et al.  Treatment plan optimization using linear programming. , 1991, Medical physics.

[10]  L. Khachiyan.  Polynomial algorithms in linear programming , 1980 .

[11]  J. Dai,et al.  Beam orientation optimization for intensity-modulated radiation therapy using mixed integer programming , 2006, Physics in medicine and biology.

[12]  Andras Lasso,et al.  SlicerRT: radiation therapy research toolkit for 3D Slicer. , 2012, Medical physics.

[13]  Arvind Kumar,et al.  A New Linear Programming Approach to Radiation Therapy Treatment Planning Problems , 2006, Oper. Res..

[14]  Frank Verhaegen,et al.  A review of treatment planning for precision image-guided photon beam pre-clinical animal radiation studies. , 2014, Zeitschrift fur medizinische Physik.

[15]  Eva K. Lee,et al.  Optimization of radiosurgery treatment planning via mixed integer programming. , 2000, Medical physics.

[16]  G. Dantzig.  Programming of Interdependent Activities: II Mathematical Model , 1949 .

[17]  J. Kereiakes,et al.  The method of linear programming applied to radiation treatment planning. , 1968, Radiology.

[18]  David A Jaffray,et al.  Two-dimensional inverse planning and delivery with a preclinical image guided microirradiator. , 2013, Medical physics.

[19]  Frank Verhaegen,et al.  A framework for inverse planning of beam-on times for 3D small animal radiotherapy using interactive multi-objective optimisation , 2015, Physics in medicine and biology.

[20]  Russell H. Taylor,et al.  Real-time dose computation: GPU-accelerated source modeling and superposition/convolution. , 2010, Medical physics.

[21]  Wufan Chen,et al.  A particle swarm optimization algorithm for beam angle selection in intensity-modulated radiotherapy planning , 2005, Physics in medicine and biology.

[22]  Justin David Durfee,et al.  Comparison of open-source linear programming solvers. , 2013 .

[23]  P. Kazanzides,et al.  Validation of GPU-accelerated superposition-convolution dose computations for the Small Animal Radiation Research Platform. , 2018, Medical physics.

[24]  Michael C. Ferris,et al.  Radiosurgery Treatment Planning via Nonlinear Programming , 2003, Ann. Oper. Res..

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

[26]  G A Ezzell,et al.  Genetic and geometric optimization of three-dimensional radiation therapy treatment planning. , 1996, Medical physics.

[27]  William E. Lorensen,et al.  The NA-MIC Kit: ITK, VTK, pipelines, grids and 3D slicer as an open platform for the medical image computing community , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[28]  Benoît Ozell,et al.  A convolution-superposition dose calculation engine for GPUs. , 2010, Medical physics.

[29]  Peter Kazanzides,et al.  Robotic delivery of complex radiation volumes for small animal research , 2010, 2010 IEEE International Conference on Robotics and Automation.

[30]  Horst W. Hamacher,et al.  Inverse Radiation Therapy Planning: A Multiple Objective Optimisation Approach , 1999 .