matRad - a multi-modality open source 3D treatment planning toolkit

We present matRad, an open source software for three-dimensional radiation treatment planning of intensitymodulated photon, proton, and carbon ion therapy. matRad is developed for educational and research purposes; it is entirely written in MATLAB. A first beta release is available for download. The toolkit features a highly modular design with a set of individual functions modeling the entire treatment planning workflow based on a segmented patient CT. All algorithms, e.g. for ray tracing, photon/proton/carbon dose calculation, fluence optimization, and multileaf collimator sequencing, follow well-established approaches and operate on clinically adequate voxel and bixel resolution. Patient data as well as base data for all required computations is included in matRad. We achieve computation times of 60-100s (60-400s) for realistic patient cases including photon (particle) dose calculation and fluence optimization. Memory consumption ranges between 0.2GB and 2.2GB. Dose distributions of a treatment planning study for a phantom and prostate patient case considering multiple radiation modalities are shown. Both the computational and dosimetric results encourage a future use of matRad in an educational and scientific setting.

[1]  P. Xia,et al.  Multileaf collimator leaf sequencing algorithm for intensity modulated beams with multiple static segments. , 1998, Medical physics.

[2]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[3]  Y. X. Wang,et al.  Nuclear Instruments and Methods in Physics Research Section B : Beam Interactions with Materials and Atoms , 2018 .

[4]  Dávid Papp,et al.  Shared data for intensity modulated radiation therapy (IMRT) optimization research: the CORT dataset , 2014, GigaScience.

[5]  Robert J. Schneider,et al.  Multiple Coulomb scattering of 160 MeV protons , 1993 .

[6]  T. Bortfeld,et al.  Decomposition of pencil beam kernels for fast dose calculations in three-dimensional treatment planning. , 1993, Medical physics.

[7]  T Bortfeld,et al.  An analytical approximation of the Bragg curve for therapeutic proton beams. , 1997, Medical physics.

[8]  D Sánchez-Parcerisa,et al.  FoCa: a modular treatment planning system for proton radiotherapy with research and educational purposes , 2014, Physics in medicine and biology.

[9]  T. Bortfeld,et al.  Inverse planning for photon and proton beams. , 2001, Medical dosimetry : official journal of the American Association of Medical Dosimetrists.

[10]  Uwe Oelfke,et al.  A 3D isodose manipulation tool for interactive dose shaping , 2014 .

[11]  Simeon Nill.  Development and application of a multi-modality inverse treatment planning system , 2001 .

[12]  M Goitein,et al.  A pencil beam algorithm for proton dose calculations. , 1996, Physics in medicine and biology.

[13]  Philipp Hennig,et al.  Analytical probabilistic modeling for radiation therapy treatment planning , 2013, Physics in medicine and biology.

[14]  Joseph O Deasy,et al.  CERR: a computational environment for radiotherapy research. , 2003, Medical physics.

[15]  O Jäkel,et al.  Treatment planning for heavy-ion radiotherapy: physical beam model and dose optimization. , 2000, Physics in medicine and biology.

[16]  J. Mechalakos,et al.  IMRT commissioning: multiple institution planning and dosimetry comparisons, a report from AAPM Task Group 119. , 2009, Medical physics.