Optimization Tools for Radiation Treatment Planning in Matlab

by Gino (Jinho) Lim , Michael C. Ferris and David M. Shepard

We provide a suite of optimization tools for radiation treatment planning within the Matlab programming environment. The data included with these tools was computed for real patient cases using a Monte Carlo dose engine. The formulation of a series of optimization models is described that utilizes this data within a modeling system. Furthermore, visualization techniques are provided that assist in validating the quality of each solution. The versatility and utility of the tools are shown using a sequence of optimization techniques designed to generate a practical solution.

Before downloading any files, you must have:

A licensed Matlab (release 13 or higher) installed on your computer.
A fullly licensed GAMS (20.7 or higher preferred) with appropriate solvers for LP, NLP, and MIP.
Matlab/GAMS interface. All necessary files and documents for the setup can be found and downloaded from here.

MATLAB utility files (Matlab M- & MEX-files ):

Model Setup: inputs.m
Model Data Generation: gendata.m, gendata.mexglx (Mex file), and genData (an executable file for gendata)
Generation of the rind of an organ structure: genrind.m and (For linux): genrind.mexglx (For pc): genrind.dll
Read 3D structure information from a GDX file: readstruct.m and readstruct.mexglx and readstruct.dll
Write 3D structure information into a GDX file: writestruct.m and writestruct.mexglx and writestruct.dll
Treatment Plan Generation: rungms.m
DVH plot: dvh.m
Dose distribution plot: doseplot.m and vol2surf.m
Solution refinement utility: extract.m
Data sampling utility: gensamp.m
Creating arbitrary organ shapes:neworgans.m

All of the above files can be downloaded here (toolutils.zip, 501K)

GAMS Optimization Models

Beam Weight Optimization using linear programming: lp.gms
Beam Weight Optimization using quadratic programming: qp.gms
Beam Weight & Angle Optimization using mixed integer programming: optangle.gms
Beam Weight & Angle & Wedge Orientation Optimization using a mixed integer programming: optwedge.gms

All GAMS optimization files can be downloaded here (models.zip, 5K)

Data

ferris@cs.wisc.edu

jin-ho@cs.wisc.edu

The three cases are:

Pelvis
- Dosimetrist Plan: (6 beam angles: 45, 90, 125, 235, 270, 315)
- Prescription:
  - PTV: prescription = 1.0 min = 0.97 max = 1.05
  - Rectum: min = 0.1 max = 1.0; DVH = let no more that 30% of volume exceed a dose of 0.6
  - Bladder: min = 0.1 max = 1.0; DVH = let no more that 35% of volume exceed a dose of 0.6
Breast
- Dosimetrist Plan: (2 beam angles: 317, 137 - wedge for 137 only)
- Prescription:
  - PTV: (Left breast) prescription = 1.0 min = 0.97 max = 1.05
  - Heart: min = 0.1 max = 1.0; DVH = let no more that 15% of volume exceed a dose of 0.3
  - Left Lung: min = 0.1 max = 1.0; DVH = let no more that 20% of volume exceed a dose of 0.2
Pancreas
- Dosimetrist Plan: (4 beam angles: 345, 260, 90, 170 - wedges used on 3 of 4 beams)
- Prescription:
  - PTV: prescription = 1.0 min = 0.99 max = 1.05
  - Spinal Cord: min = 0.0 max = 0.4; nothing above 0.4
  - Left Kidney: min = 0.0 max = 0.55; DVH = let no more that 20% of volume exceed a dose of 0.15
  - Rt. Kidney: min = 0.0 max = 0.55; DVH = let no more that 20% of volume exceed a dose of 0.20
  - Liver: min = 0.0 max = 0.8; DVH = let no more than 25% of volume exceed a dose of 0.15

User Guide and related reference

Michael C. Ferris, Jinho Lim, and David Shepard, "Optimization tools for radiation treatment planning in Matlab", , Optimization Technical Report 03-01, Computer Sciences Department, University of Wisconsin, Madison, Wisconsin, 2003.
J.-H. Lim, M. C. Ferris, S. J. Wright, D. M. Shepard, and M. A. Earl., An optimization framework for conformal radiation treatment planning, , Optimization Technical Report 02-10, Computer Sciences Department, University of Wisconsin, Madison, Wisconsin, 2002.
Jinho Lim, "Optimization in radiation treatment planning", , PhD Thesis, Optimization Technical Report 02-11,Computer Sciences Department, University of Wisconsin, Madison, Wisconsin, December 2002.