Development of the open-source dose calculation and optimization toolkit matRad.

Abstract:

PURPOSE:We report on the development of the open-source cross-platform radiation treatment planning toolkit matRad and its comparison against validated treatment planning systems. The toolkit enables three-dimensional intensity-modulated radiation therapy treatment planning for photons, scanned protons and scanned carbon ions. METHODS:matRad is entirely written in Matlab and is freely available online. It re-implements well-established algorithms employing a modular and sequential software design to model the entire treatment planning workflow. It comprises core functionalities to import DICOM data, to calculate and optimize dose as well as a graphical user interface for visualization. matRad dose calculation algorithms (for carbon ions this also includes the computation of the relative biological effect) are compared against dose calculation results originating from clinically approved treatment planning systems. RESULTS:We observe three-dimensional γ-analysis pass rates ≥ 99.67% for all three radiation modalities utilizing a distance to agreement of 2 mm and a dose difference criterion of 2%. The computational efficiency of matRad is evaluated in a treatment planning study considering three different treatment scenarios for every radiation modality. For photons, we measure total run times of 145 s-1260 s for dose calculation and fluence optimization combined considering 4-72 beam orientations and 2608-13597 beamlets. For charged particles, we measure total run times of 63 s-993 s for dose calculation and fluence optimization combined considering 9963-45574 pencil beams. Using a CT and dose grid resolution of 0.3 cm3 requires a memory consumption of 1.59 GB-9.07 GB and 0.29 GB-17.94 GB for photons and charged particles, respectively. CONCLUSION:The dosimetric accuracy, computational performance and open-source character of matRad encourages a future application of matRad for both educational and research purposes.

journal_name

Med Phys

journal_title

Medical physics

authors

Wieser HP,Cisternas E,Wahl N,Ulrich S,Stadler A,Mescher H,Müller LR,Klinge T,Gabrys H,Burigo L,Mairani A,Ecker S,Ackermann B,Ellerbrock M,Parodi K,Jäkel O,Bangert M

doi

10.1002/mp.12251

subject

Has Abstract

pub_date

2017-06-01 00:00:00

pages

2556-2568

issue

6

eissn

0094-2405

issn

2473-4209

journal_volume

44

pub_type

杂志文章
  • Active energy selective image detector for dual-energy computed radiography.

    abstract::A new energy selective detector for dual energy computed radiography has been developed that combines many of the advantages of x-ray tube voltage switching and single exposure double screen detectors. The new active detector utilizes electro-optical modulation of the response of the storage phosphor screens to allow ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.597831

    authors: Alvarez RE

    更新日期:1996-10-01 00:00:00

  • Response of an implantable MOSFET dosimeter to 192Ir HDR radiation.

    abstract::New in vivo dosimetry methods would be useful for clinical HDR brachytherapy. An implantable MOSFET Dose Verification System designed by Sicel Technologies, Inc. was examined for use with 192Ir HDR applications. This investigation demonstrated that varying the dose rate from 22 to 84 cGy/min did not change detector re...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3013574

    authors: Fagerstrom JM,Micka JA,DeWerd LA

    更新日期:2008-12-01 00:00:00

  • Validation of GPU-accelerated superposition-convolution dose computations for the Small Animal Radiation Research Platform.

    abstract:PURPOSE:The Small Animal Radiation Research Platform (SARRP) has been developed for conformal microirradiation with on-board cone beam CT (CBCT) guidance. The graphics processing unit (GPU)-accelerated Superposition-Convolution (SC) method for dose computation has been integrated into the treatment planning system (TPS...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12862

    authors: Cho N,Tsiamas P,Velarde E,Tryggestad E,Jacques R,Berbeco R,McNutt T,Kazanzides P,Wong J

    更新日期:2018-05-01 00:00:00

  • Experimental assessment of microwave ablation computational modeling with MR thermometry.

    abstract:PURPOSE:Computational models are widely used during the design and characterization of microwave ablation (MWA) devices, and have been proposed for pretreatment planning. Our objective was to assess three-dimensional (3D) transient temperature and ablation profiles predicted by MWA computational models with temperature...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.14318

    authors: Faridi P,Keselman P,Fallahi H,Prakash P

    更新日期:2020-09-01 00:00:00

  • Radiologic validation of a fast neutron multileaf collimator.

    abstract::Teletherapy with high linear energy transfer radiations (LET), perhaps more than with low LET types, requires careful beam collimation to limit effects to normal structures. Intensity modulated techniques may also hold promise in this regard. Accordingly, a remote computer-controlled, high-resolution multileaf collima...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2760026

    authors: Farr JB,Maughan RL,Yudelev M,Blosser E,Brandon J,Horste T,Forman JD

    更新日期:2007-09-01 00:00:00

  • The effect of angular spread on the intensity distribution of arbitrarily shaped electron beams.

    abstract::Knowledge of the relative intensity distribution at the patient's surface is essential for pencil beam calculations of three-dimensional dose distributions for arbitrarily shaped electron beams. To calculate the relative intensity distribution, the spatial spread resulting from angular spread is convolved with a two-d...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.596252

    authors: Mohan R,Chui CS,Fontenla D,Han K,Ballon D

    更新日期:1988-03-01 00:00:00

  • An accelerated threshold-based back-projection algorithm for compton camera image reconstruction.

    abstract:PURPOSE:Compton camera imaging (CCI) systems are currently under investigation for radiotherapy dose reconstruction and verification. The ability of such a system to provide real-time images during dose delivery will be limited by the computational speed of the image reconstruction algorithm. In this work, the authors ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3519873

    authors: Mundy DW,Herman MG

    更新日期:2011-01-01 00:00:00

  • Analysis and correction of imperfections in the image intensifier-TV-digitizer imaging chain.

    abstract::Image intensifier-television-video digitizer (IITVD) systems are commonly used for digital planar image acquisition in radiology. However, the well-known distortions inherent in these systems limit their utility in research and in some clinical applications where quantitatively correct images are required. Software co...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.596718

    authors: Boone JM,Seibert JA,Barrett WA,Blood EA

    更新日期:1991-03-01 00:00:00

  • Characterization of statistical prior image constrained compressed sensing (PICCS): II. Application to dose reduction.

    abstract:PURPOSE:The ionizing radiation imparted to patients during computed tomography exams is raising concerns. This paper studies the performance of a scheme called dose reduction using prior image constrained compressed sensing (DR-PICCS). The purpose of this study is to characterize the effects of a statistical model of x...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4773866

    authors: Lauzier PT,Chen GH

    更新日期:2013-02-01 00:00:00

  • GLAaS: an absolute dose calibration algorithm for an amorphous silicon portal imager. Applications to IMRT verifications.

    abstract::A new calibration algorithm (GLAaS) to derive absolute dose maps from images acquired with the Varian PV-aS500 electronic portal imager (based on amorphous silicon detectors) has been developed incorporating the dependence of detector response on primary and transmitted radiation and on field size. Detector calibratio...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2218314

    authors: Nicolini G,Fogliata A,Vanetti E,Clivio A,Cozzi L

    更新日期:2006-08-01 00:00:00

  • The x-ray time of flight method for investigation of ghosting in amorphous selenium-based flat panel medical x-ray imagers.

    abstract::Amorphous selenium (a-Se) based real-time flat-panel imagers (FPIs) are finding their way into the digital radiology department because they offer the practical advantages of digital x-ray imaging combined with an image quality that equals or outperforms that of conventional systems. The temporal imaging characteristi...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2042248

    authors: Rau AW,Bakueva L,Rowlands JA

    更新日期:2005-10-01 00:00:00

  • Noise, sampling, and the number of projections in cone-beam CT with a flat-panel detector.

    abstract:PURPOSE:To investigate the effect of the number of projection views on image noise in cone-beam CT (CBCT) with a flat-panel detector. METHODS:This fairly fundamental consideration in CBCT system design and operation was addressed experimentally (using a phantom presenting a uniform medium as well as statistically moti...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4875688

    authors: Zhao Z,Gang GJ,Siewerdsen JH

    更新日期:2014-06-01 00:00:00

  • Intensity modulated arc deliveries approximated by a large number of fixed gantry position sliding window dynamic multileaf collimator fields.

    abstract::The intensity modulated arc has been proposed as an alternative to tomotherapy. Treatment planing systems more typically model the conventional step and shoot or sliding window dynamic multileaf collimator (DMLC) deliveries, and may not support intensity modulated arc therapy (IMAT). As well, another potential drawbac...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.1508110

    authors: MacKenzie MA,Robinson DM

    更新日期:2002-10-01 00:00:00

  • SU-E-T-07: Edema Induced Changes in Tumor Cell Survival Fraction and Tumor Control Probability in 131Cs Permanent Prostate Implant Patients.

    abstract:PURPOSE:The purpose of this study is to estimate the effect of edema, developed during implant procedure, on tumor cell surviving fraction(SF) and tumor control probability(TCP) in the patients of prostate cancer who underwent 131 Cs permanent seed implants. METHODS:The impact of edema on SF and TCP, was calculated us...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4735061

    authors: Kehwar TS,Huq MS,Smith RP

    更新日期:2012-06-01 00:00:00

  • Contrast-detail analysis of image degradation due to lossy compression.

    abstract::A contrast-detail (CD) experiment was performed to study the effect of lossy compression on computed radiographic (CR) images. Digital CR images of a phantom were compressed by quantizing the full-frame discrete cosine transform and Huffman encoding the result. Since low-contrast detectability is directly linked to an...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.597489

    authors: Cook LT,Insana MF,McFadden MA,Hall TJ,Cox GG

    更新日期:1995-06-01 00:00:00

  • Quantitative study of focused ultrasound enhanced doxorubicin delivery to prostate tumor in vivo with MRI guidance.

    abstract:PURPOSE:The purpose of this study was to investigate the potential of MR-guided pulsed focused ultrasound (pFUS) for the enhancement of drug uptake in prostate tumors in vivo using doxorubicin (Dox). METHODS:An antitumor drug Dox, an orthotopic animal prostate tumor model using human prostate cancer, LNCaP cell line, ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4705346

    authors: Chen X,Cvetkovic D,Ma CM,Chen L

    更新日期:2012-05-01 00:00:00

  • Technical note: the calibration of 90Y-labeled SIR-Spheres using a nondestructive spectroscopic assay.

    abstract::90Y-labeled SIR-Spheres are currently used to treat patients with hepatic metastases secondary to colorectal adenocarcinoma. In general, the prescribed activity is based on empirical data collected during clinical trials. The activity of the source vial is labeled by the manufacturer as 3.0 GBq +/- 10% and is not inde...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2889621

    authors: Selwyn R,Micka J,DeWerd L,Nickles R,Thomadsen B

    更新日期:2008-04-01 00:00:00

  • Unexpectedly large dose rate dependent output from a linear accelerator.

    abstract::During our routine calibration of a Varian Clinac-20 linear accelerator, the absorbed dose for a fixed monitor unit (mu) was found to decrease with increasing dose rate. Between dose rates of 100 and 500 mu/min, there was up to 20% difference in absorbed dose for a 20-MeV electron beam. The cause of this problem was a...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.596192

    authors: Cheng PC,Kubo H

    更新日期:1988-09-01 00:00:00

  • On effective dose for radiotherapy based on doses to nontarget organs and tissues.

    abstract:PURPOSE:The National Council for Radiation Protection and Measurement (NCRP) published estimates for the collective population dose and the mean effective dose to the population of the United States from medical imaging procedures for 1980/1982 and for 2006. The earlier report ignored the effective dose from radiothera...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4906190

    authors: Uselmann AJ,Thomadsen BR

    更新日期:2015-02-01 00:00:00

  • Pretreatment patient-specific IMRT quality assurance: a correlation study between gamma index and patient clinical dose volume histogram.

    abstract:PURPOSE:The aim of this work is to investigate the predictive power of a common conventional intensity modulated radiation therapy (IMRT) quality assurance (QA) performance metric, the gamma passing rate (%GP), through the analysis of the sensitivity and of the correlation between %GP and different dose discrepancies b...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4767763

    authors: Stasi M,Bresciani S,Miranti A,Maggio A,Sapino V,Gabriele P

    更新日期:2012-12-01 00:00:00

  • A deep feature fusion methodology for breast cancer diagnosis demonstrated on three imaging modality datasets.

    abstract:BACKGROUND:Deep learning methods for radiomics/computer-aided diagnosis (CADx) are often prohibited by small datasets, long computation time, and the need for extensive image preprocessing. AIMS:We aim to develop a breast CADx methodology that addresses the aforementioned issues by exploiting the efficiency of pre-tra...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12453

    authors: Antropova N,Huynh BQ,Giger ML

    更新日期:2017-10-01 00:00:00

  • Selecting registration schemes in case of interstitial lung disease follow-up in CT.

    abstract:PURPOSE:Primary goal of this study is to select optimal registration schemes in the framework of interstitial lung disease (ILD) follow-up analysis in CT. METHODS:A set of 128 multiresolution schemes composed of multiresolution nonrigid and combinations of rigid and nonrigid registration schemes are evaluated, utilizi...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4923170

    authors: Vlachopoulos G,Korfiatis P,Skiadopoulos S,Kazantzi A,Kalogeropoulou C,Pratikakis I,Costaridou L

    更新日期:2015-08-01 00:00:00

  • Technical Note: Performance comparison of ultra-high-resolution scan modes of two clinical computed tomography systems.

    abstract:PURPOSE:Two ultra-high-resolution (UHR) computed tomography (CT) scanners are clinically available. One is achieved by a CT system with a 0.25 mm × 160 row detector (Detector-UHR), whereas the other is with a 0.6 mm × 32 row detector with in-plane comb filtering in a dual source CT (Comb-UHR). We compared radiation dos...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.13949

    authors: Kawashima H,Ichikawa K,Takata T,Nagata H,Hoshika M,Akagi N

    更新日期:2020-02-01 00:00:00

  • Determination of the initial beam parameters in Monte Carlo linac simulation.

    abstract::For Monte Carlo linac simulations and patient dose calculations, it is important to accurately determine the phase space parameters of the initial electron beam incident on the target. These parameters, such as mean energy and radial intensity distribution, have traditionally been determined by matching the calculated...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2168433

    authors: Aljarrah K,Sharp GC,Neicu T,Jiang SB

    更新日期:2006-04-01 00:00:00

  • Thermodynamic considerations of acrylic cement implant at the site of giant cell tumors of the bone.

    abstract::A discussion of the thermodynamic aspects of a relatively new treatment method for giant cell tumors of the bone is presented in this paper. The advantages of implanting methylmethacrylate acrylic bone cement into a curetted tumor site are briefly discussed and placed in perspective relative to more prevalent surgical...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.595902

    authors: Krishnan EC,Nelson C,Neff JR

    更新日期:1986-03-01 00:00:00

  • Dose and dose averaged LET comparison of ¹H, ⁴He, ⁶Li, ⁸Be, ¹⁰B, ¹²C, ¹⁴N, and ¹⁶O ion beams forming a spread-out Bragg peak.

    abstract:PURPOSE:Modern clinical accelerators are capable of producing ion beams from protons up to neon. This work compares the depth dose distribution and corresponding dose averaged linear energy transfer (LET) distribution, which is related to the biological effectiveness, for different ion beams (¹H, ⁴He, ⁶Li, ⁸Be, ¹⁰B, ¹²...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3662911

    authors: Kantemiris I,Karaiskos P,Papagiannis P,Angelopoulos A

    更新日期:2011-12-01 00:00:00

  • SU-E-T-555: Weighted Four-Dimensional IMRT Planning for Dynamic MLC Tracking Using a Practical and Simple Framework.

    abstract:PURPOSE:We present a simple, practical framework for truly 4D lung IMRT planning based on a weighted individual-phase optimization paradigm. This strategy is specifically developed for use in real-time tumor tracking delivery systems so as to utilize respiratory motion as an additional degree of freedom rather than a c...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4735644

    authors: Tachibana H,Cheung Y,Jain S,Sawant A

    更新日期:2012-06-01 00:00:00

  • Improved dynamic-programming-based algorithms for segmentation of masses in mammograms.

    abstract::In this paper, two new boundary tracing algorithms for segmentation of breast masses are presented. These new algorithms are based on the dynamic programming-based boundary tracing (DPBT) algorithm proposed in Timp and Karssemeijer, [S. Timp and N. Karssemeijer, Med. Phys. 31, 958-971 (2004)] The DPBT algorithm contai...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2791034

    authors: Rojas Domínguez A,Nandi AK

    更新日期:2007-11-01 00:00:00

  • Phantoms with 10BF3 detectors for boron neutron capture therapy applications.

    abstract::Two acrylic cube phantoms have been constructed for BNCT applications that allow the depth distribution of neutrons to be measured with miniature 10BF3 detectors in 0.5-cm steps beginning at 1-cm depth. Sizes and weights of the cubes are 14 cm, 3.230 kg, and 11 cm, 1.567 kg. Tests were made with the epithermal neutron...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.598354

    authors: Alburger DE,Raparia D,Zucker MS

    更新日期:1998-09-01 00:00:00

  • SU-E-I-108: Simple Approach to Determine Patient's Entrance Skin Exposure for Radiographic.

    abstract:PURPOSE:To evaluate a model which can predict patient's Entrance Skin Exposure (ESE) from height and weight and a single kVp mR/mAs measurement for radiographic. METHODS:Patient height, weight and post exposure mAs are obtained to assess ESE based on IAEA model. Patient's thickness is obtained by assuming a cylindrica...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4734825

    authors: Maricle S,Jeutang L,Prasad G

    更新日期:2012-06-01 00:00:00