1
|
Zarepisheh M, Hong L, Zhou Y, Huang Q, Yang J, Jhanwar G, Pham HD, Dursun P, Zhang P, Hunt MA, Mageras GS, Yang JT, Yamada Y, Deasy JO. Automated and Clinically Optimal Treatment Planning for Cancer Radiotherapy. INFORMS JOURNAL ON APPLIED ANALYTICS 2022; 52:69-89. [PMID: 35847768 PMCID: PMC9284667 DOI: 10.1287/inte.2021.1095] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Each year, approximately 18 million new cancer cases are diagnosed worldwide, and about half must be treated with radiotherapy. A successful treatment requires treatment planning with the customization of penetrating radiation beams to sterilize cancerous cells without harming nearby normal organs and tissues. This process currently involves extensive manual tuning of parameters by an expert planner, making it a time-consuming and labor-intensive process, with quality and immediacy of critical care dependent on the planner's expertise. To improve the speed, quality, and availability of this highly specialized care, Memorial Sloan Kettering Cancer Center developed and applied advanced optimization tools to this problem (e.g., using hierarchical constrained optimization, convex approximations, and Lagrangian methods). This resulted in both a greatly improved radiotherapy treatment planning process and the generation of reliable and consistent high-quality plans that reflect clinical priorities. These improved techniques have been the foundation of high-quality treatments and have positively impacted over 4,000 patients to date, including numerous patients in severe pain and in urgent need of treatment who might have otherwise required longer hospital stays or undergone unnecessary surgery to control the progression of their disease. We expect that the wide distribution of the system we developed will ultimately impact patient care more broadly, including in resource-constrained countries.
Collapse
Affiliation(s)
- Masoud Zarepisheh
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Linda Hong
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Ying Zhou
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Qijie Huang
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Jie Yang
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Gourav Jhanwar
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Hai D Pham
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Pinar Dursun
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Pengpeng Zhang
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Margie A Hunt
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Gig S Mageras
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| | - Jonathan T Yang
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York
| | - Yoshiya Yamada
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York
| | - Joseph O Deasy
- Departments of Medical Physics, Memorial Sloan Kettering Cancer Center, New York
| |
Collapse
|
2
|
Huang C, Yang Y, Panjwani N, Boyd S, Xing L. Pareto Optimal Projection Search (POPS): Automated Radiation Therapy Treatment Planning by Direct Search of the Pareto Surface. IEEE Trans Biomed Eng 2021; 68:2907-2917. [PMID: 33523802 PMCID: PMC8526351 DOI: 10.1109/tbme.2021.3055822] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
OBJECTIVE Radiation therapy treatment planning is a time-consuming, iterative process with potentially high inter-planner variability. Fully automated treatment planning processes could reduce a planner's active treatment planning time and remove inter-planner variability, with the potential to tremendously improve patient turnover and quality of care. In developing fully automated algorithms for treatment planning, we have two main objectives: to produce plans that are 1) Pareto optimal and 2) clinically acceptable. Here, we propose the Pareto optimal projection search (POPS) algorithm, which provides a general framework for directly searching the Pareto front. METHODS Our POPS algorithm is a novel automated planning method that combines two main search processes: 1) gradient-free search in the decision variable space and 2) projection of decision variables to the Pareto front using the bisection method. We demonstrate the performance of POPS by comparing with clinical treatment plans. As one possible quantitative measure of treatment plan quality, we construct a clinical acceptability scoring function (SF) modified from the previously developed general evaluation metric (GEM). RESULTS On a dataset of 21 prostate cases collected as part of clinical workflow, our proposed POPS algorithm produces Pareto optimal plans that are clinically acceptable in regards to dose conformity, dose homogeneity, and sparing of organs-at-risk. CONCLUSION Our proposed POPS algorithm provides a general framework for fully automated treatment planning that achieves clinically acceptable dosimetric quality without requiring active planning from human planners. SIGNIFICANCE Our fully automated POPS algorithm addresses many key limitations of other automated planning approaches, and we anticipate that it will substantially improve treatment planning workflow.
Collapse
|
3
|
Huang C, Yang Y, Xing L. Fully automated noncoplanar radiation therapy treatment planning. Med Phys 2021; 48:7439-7449. [PMID: 34519064 DOI: 10.1002/mp.15223] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Revised: 08/23/2021] [Accepted: 08/30/2021] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To perform fully automated noncoplanar (NC) treatment planning, we propose a method called NC-POPS to produce NC plans using the Pareto optimal projection search (POPS) algorithm. METHODS NC radiation therapy treatment planning has the potential to improve dosimetric quality as compared to traditional coplanar techniques. Likewise, automated treatment planning algorithms can reduce a planner's active treatment planning time and remove inter-planner variability. Our NC-POPS algorithm extends the original POPS algorithm to the NC setting with potential applications to both intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT). The proposed algorithm consists of two main parts: (1) NC beam angle optimization (BAO) and (2) fully automated inverse planning using the POPS algorithm. RESULTS We evaluate the performance of NC-POPS by comparing between various NC and coplanar configurations. To evaluate plan quality, we compute the homogeneity index (HI), conformity index (CI), and dose-volume histogram statistics for various organs-at-risk (OARs). As compared to the evaluated coplanar baseline methods, the proposed NC-POPS method achieves significantly better OAR sparing, comparable or better dose conformity, and similar dose homogeneity. CONCLUSIONS Our proposed NC-POPS algorithm provides a modular approach for fully automated treatment planning of NC IMRT cases with the potential to substantially improve treatment planning workflow and plan quality.
Collapse
Affiliation(s)
- Charles Huang
- Department of Bioengineering, Stanford University, Stanford, California, USA
| | - Yong Yang
- Department of Radiation Oncology, Stanford University, Stanford, California, USA
| | - Lei Xing
- Department of Radiation Oncology, Stanford University, Stanford, California, USA
| |
Collapse
|
4
|
Wu VW, Epelman MA, Pasupathy KS, Sir MY, Deufel CL. A new optimization algorithm for HDR brachytherapy that improves DVH-based planning: Truncated Conditional Value-at-Risk (TCVaR). Biomed Phys Eng Express 2020; 6. [PMID: 35102005 DOI: 10.1088/2057-1976/abb4bc] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Accepted: 09/02/2020] [Indexed: 11/12/2022]
Abstract
Purpose:To introduce a new optimization algorithm that improves DVH results and is designed for the type of heterogeneous dose distributions that occur in brachytherapy.Methods:The new optimization algorithm is based on a prior mathematical approach that uses mean doses of the DVH metric tails. The prior mean dose approach is referred to as conditional value-at-risk (CVaR), and unfortunately produces noticeably worse DVH metric results than gradient-based approaches. We have improved upon the CVaR approach, using the so-called Truncated CVaR (TCVaR), by excluding the hottest or coldest voxels in the structure from the calculations of the mean dose of the tail. Our approach applies an iterative sequence of convex approximations to improve the selection of the excluded voxels. Data Envelopment Analysis was used to quantify the sensitivity of TCVaR results to parameter choice and to compare the quality of a library of 256 TCVaR plans created for each of prostate, breast, and cervix treatment sites with commercially-generated plans.Results:In terms of traditional DVH metrics, TCVaR outperformed CVaR and the improvements increased monotonically as more iterations were used to identify and exclude the hottest/coldest voxels from the optimization problem. TCVaR also outperformed the Eclipse-Brachyvision TPS, with an improvement in PTVD95% (for equivalent organ-at-risk doses) of up to 5% (prostate), 3% (breast), and 1% (cervix).Conclusions:A novel optimization algorithm for HDR treatment planning produced plans with superior DVH metrics compared with a prior convex optimization algorithm as well as Eclipse-Brachyvision. The algorithm is computationally efficient and has potential applications as a primary optimization algorithm or quality assurance for existing optimization approaches.
Collapse
Affiliation(s)
- Victor W Wu
- Department of Industrial and Operations Engineering, University of Michigan, Ann Arbor, MI 48109, United States of America.,Department of Radiation Oncology, Mayo Clinic, Rochester, MN 55905, United States of America
| | - Marina A Epelman
- Department of Industrial and Operations Engineering, University of Michigan, Ann Arbor, MI 48109, United States of America
| | - Kalyan S Pasupathy
- Department of Health Sciences Research, Mayo Clinic, Rochester, MN 55905, United States of America.,Kern Center for the Science of Health Care Delivery, Mayo Clinic, Rochester, MN 55905, United States of America
| | - Mustafa Y Sir
- Department of Health Sciences Research, Mayo Clinic, Rochester, MN 55905, United States of America.,Kern Center for the Science of Health Care Delivery, Mayo Clinic, Rochester, MN 55905, United States of America
| | - Christopher L Deufel
- Department of Radiation Oncology, Mayo Clinic, Rochester, MN 55905, United States of America
| |
Collapse
|
5
|
Hong L, Zhou Y, Yang J, Mechalakos JG, Hunt MA, Mageras GS, Yang J, Yamada J, Deasy JO, Zarepisheh M. Clinical Experience of Automated SBRT Paraspinal and Other Metastatic Tumor Planning With Constrained Hierarchical Optimization. Adv Radiat Oncol 2019; 5:1042-1050. [PMID: 33083666 PMCID: PMC7557131 DOI: 10.1016/j.adro.2019.11.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Revised: 11/12/2019] [Accepted: 11/26/2019] [Indexed: 11/18/2022] Open
Abstract
Purpose We report on the clinical performance of a fully automated approach to treatment planning based on a Pareto optimal, constrained hierarchical optimization algorithm, named Expedited Constrained Hierarchical Optimization (ECHO). Methods and materials From April 2017 to October 2018, ECHO produced 640 treated plans for 523 patients who underwent stereotactic body radiation therapy (RT) for paraspinal and other metastatic tumors. A total of 182 plans were for 24 Gy in a single fraction, 387 plans were for 27 Gy in 3 fractions, and the remainder were for other prescriptions or fractionations. Of the plans, 84.5% were for paraspinal tumors, with 69, 302, and 170 in the cervical, thoracic, and lumbosacral spine, respectively. For each case, after contouring, a template plan using 9 intensity modulated RT fields based on disease site and tumor location was sent to ECHO through an application program interface plug-in from the treatment planning system. ECHO returned a plan that satisfied all critical structure hard constraints with optimal target volume coverage and the lowest achievable normal tissue doses. Upon ECHO completion, the planner received an e-mail indicating the plan was ready for review. The plan was accepted if all clinical criteria were met. Otherwise, a limited number of parameters could be adjusted for another ECHO run. Results The median planning target volume size was 84.3 cm3 (range, 6.9-633.2). The median time to produce 1 ECHO plan was 63.5 minutes (range, 11-340 minutes) and was largely dependent on the field sizes. Of the cases, 79.7% required 1 run to produce a clinically accepted plan, 13.3% required 1 additional run with minimal parameter adjustments, and 7.0% required ≥2 additional runs with significant parameter modifications. All plans met or bettered the institutional clinical criteria. Conclusions We successfully implemented automated stereotactic body RT paraspinal and other metastatic tumors planning. ECHO produced high-quality plans, improved planning efficiency and robustness, and enabled expedited treatment planning at our clinic.
Collapse
Affiliation(s)
- Linda Hong
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Ying Zhou
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Jie Yang
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - James G Mechalakos
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Margie A Hunt
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Gig S Mageras
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Jonathan Yang
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Josh Yamada
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Joseph O Deasy
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Masoud Zarepisheh
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York
| |
Collapse
|
6
|
Zarepisheh M, Hong L, Zhou Y, Oh JH, Mechalakos JG, Hunt MA, Mageras GS, Deasy JO. Automated intensity modulated treatment planning: The expedited constrained hierarchical optimization (ECHO) system. Med Phys 2019; 46:2944-2954. [PMID: 31055858 DOI: 10.1002/mp.13572] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 04/25/2019] [Accepted: 04/25/2019] [Indexed: 11/09/2022] Open
Abstract
PURPOSE To develop and implement a fully automated approach to intensity modulated radiation therapy (IMRT) treatment planning. METHOD The optimization algorithm is developed based on a hierarchical constrained optimization technique and is referred internally at our institution as expedited constrained hierarchical optimization (ECHO). Beamlet contributions to regions-of-interest are precomputed and captured in the influence matrix. Planning goals are of two classes: hard constraints that are strictly enforced from the first step (e.g., maximum dose to spinal cord), and desirable goals that are sequentially introduced in three constrained optimization problems (better planning target volume (PTV) coverage, lower organ at risk (OAR) doses, and smoother fluence map). After solving the optimization problems using external commercial optimization engines, the optimal fluence map is imported into an FDA-approved treatment planning system (TPS) for leaf sequencing and accurate full dose calculation. The dose-discrepancy between the optimization and TPS dose calculation is then calculated and incorporated into optimization by a novel dose correction loop technique using Lagrange multipliers. The correction loop incorporates the leaf sequencing and scattering effects into optimization to improve the plan quality and reduce the calculation time. The resultant optimal fluence map is again imported into TPS for leaf sequencing and final dose calculation for plan evaluation and delivery. The workflow is automated using application program interface (API) scripting, requiring user interaction solely to prepare the contours and beam arrangement prior to launching the ECHO plug-in from the TPS. For each site, parameters and objective functions are chosen to represent clinical priorities. The first site chosen for clinical implementation was metastatic paraspinal lesions treated with stereotactic body radiotherapy (SBRT). As a first step, 75 ECHO paraspinal plans were generated retrospectively and compared with clinically treated plans generated by planners using VMAT (volumetric modulated arc therapy) with 4 to 6 partial arcs. Subsequently, clinical deployment began in April, 2017. RESULTS In retrospective study, ECHO plans were found to be dosimetrically superior with respect to tumor coverage, plan conformity, and OAR sparing. For example, the average PTV D95%, cord and esophagus max doses, and Paddick Conformity Index were improved, respectively, by 1%, 6%, 14%, and 15%, at a negligible 3% cost of the average skin D10cc dose. CONCLUSION Hierarchical constrained optimization is a powerful and flexible tool for automated IMRT treatment planning. The dosimetric correction step accurately accounts for detailed dosimetric multileaf collimator and scattering effects. The system produces high-quality, Pareto optimal plans and avoids the time-consuming trial-and-error planning process.
Collapse
Affiliation(s)
- Masoud Zarepisheh
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Linda Hong
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Ying Zhou
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Jung Hun Oh
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - James G Mechalakos
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Margie A Hunt
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Gig S Mageras
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Joseph O Deasy
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| |
Collapse
|
7
|
Hussein M, Heijmen BJM, Verellen D, Nisbet A. Automation in intensity modulated radiotherapy treatment planning-a review of recent innovations. Br J Radiol 2018; 91:20180270. [PMID: 30074813 DOI: 10.1259/bjr.20180270] [Citation(s) in RCA: 142] [Impact Index Per Article: 23.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Radiotherapy treatment planning of complex radiotherapy techniques, such as intensity modulated radiotherapy and volumetric modulated arc therapy, is a resource-intensive process requiring a high level of treatment planner intervention to ensure high plan quality. This can lead to variability in the quality of treatment plans and the efficiency in which plans are produced, depending on the skills and experience of the operator and available planning time. Within the last few years, there has been significant progress in the research and development of intensity modulated radiotherapy treatment planning approaches with automation support, with most commercial manufacturers now offering some form of solution. There is a rapidly growing number of research articles published in the scientific literature on the topic. This paper critically reviews the body of publications up to April 2018. The review describes the different types of automation algorithms, including the advantages and current limitations. Also included is a discussion on the potential issues with routine clinical implementation of such software, and highlights areas for future research.
Collapse
Affiliation(s)
- Mohammad Hussein
- 1 Metrology for Medical Physics Centre, National Physical Laboratory , Teddington , UK
| | - Ben J M Heijmen
- 2 Division of Medical Physics, Erasmus MC Cancer Institute , Rotterdam , The Netherlands
| | - Dirk Verellen
- 3 Faculty of Medicine and Pharmacy, Vrije Universiteit Brussel (VUB) , Brussels , Belgium.,4 Radiotherapy Department, Iridium Kankernetwerk , Antwerp , Belgium
| | - Andrew Nisbet
- 5 Department of Medical Physics, Royal Surrey County Hospital NHS Foundation Trust , Guildford , UK.,6 Department of Physics, University of Surrey , Guildford , UK
| |
Collapse
|
8
|
Lin YH, Hong LX, Hunt MA, Berry SL. Use of a constrained hierarchical optimization dataset enhances knowledge-based planning as a quality assurance tool for prostate bed irradiation. Med Phys 2018; 45:4364-4369. [PMID: 30168160 DOI: 10.1002/mp.13163] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2017] [Revised: 05/15/2018] [Accepted: 08/21/2018] [Indexed: 11/10/2022] Open
Abstract
PURPOSE To investigate whether building a knowledge-based planning (KBP) model with prostate bed plans constructed from constrained hierarchical optimization (CHO) would result in more efficient model construction with more consistent output than a model built using plans from a traditional, trial-and-error-based optimization (TEO) technique. METHODS Three KBP models were constructed from plans from subsets of 58 post-prostatectomy patients treated with intensity-modulated radiation therapy. TEO54 was built from 54 TEO plans, selected to represent typical clinical variations in target and organ-at-risk sizes and shapes. CHO30 and TEO30 were built from the same 30 patients populated with CHO and TEO plans, respectively. The three models were each applied to a new set of 18 patient scans and dose-volume histogram estimates (DVHEs) were generated for rectal and bladder walls and compared for each patient. RESULTS CHO30 resulted in a significantly tighter range in DVHEs (P < 0.01) for both the rectal and bladder walls compared with either of the TEO models, indicating less uncertainty in the dose estimation. Plans resulting from KBP optimization using each model were very similar. CONCLUSION Populating a KBP model with CHO data resulted in a high quality model. Since CHO plans can be generated automatically offline in a process that necessitates little to no user interaction, a CHO-KBP model can quickly adapt to changes in plan evaluation criteria or planning techniques without the need to wait to accrue sufficient numbers of clinical TEO plans. This may facilitate the use of KBP approaches for initial or ongoing quality assurance procedures and plan quality audits.
Collapse
Affiliation(s)
- Yen Hwa Lin
- Graduate School of Biomedical and Health Sciences, Hiroshima University, 1-2-3, Kasumi, Minami-ku, Hiroshima, 734-8551, Japan
| | - Linda X Hong
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY, 10065, USA
| | - Margie A Hunt
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY, 10065, USA
| | - Sean L Berry
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY, 10065, USA
| |
Collapse
|
9
|
Liu H, Xing L. Isodose feature-preserving voxelization (IFPV) for radiation therapy treatment planning. Med Phys 2018; 45:3321-3329. [PMID: 29772065 DOI: 10.1002/mp.12977] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Revised: 04/27/2018] [Accepted: 05/07/2018] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Inverse planning involves iterative optimization of a large number of parameters and is known to be a labor-intensive procedure. To reduce the scale of computation and improve characterization of isodose plan, this paper presents an isodose feature-preserving voxelization (IFPV) framework for radiation therapy applications and demonstrates an implementation of inverse planning in the IFPV domain. METHODS A dose distribution in IFPV scheme is characterized by partitioning the voxels into subgroups according to their geometric and dosimetric values. Computationally, the isodose feature-preserving (IFP) clustering combines the conventional voxels that are spatially and dosimetrically close into physically meaningful clusters. A K-means algorithm and support vector machine (SVM) runs sequentially to group the voxels into IFP clusters. The former generates initial clusters according to the geometric and dosimetric information of the voxels and SVM is invoked to improve the connectivity of the IFP clusters. To illustrate the utility of the formalism, an inverse planning framework in the IFPV domain is implemented, and the resultant plans of three prostate IMRT and one head-and-neck cases are compared quantitatively with that obtained using conventional inverse planning technique. RESULTS The IFPV generates models with significant dimensionality reduction without compromising the spatial resolution seen in traditional downsampling schemes. The implementation of inverse planning in IFPV domain is demonstrated. In addition to the improved computational efficiency, it is found that, for the cases studied here, the IFPV-domain inverse planning yields better treatment plans than that of DVH-based planning, primarily because of more effective use of both geometric and dose information of the system during plan optimization. CONCLUSIONS The proposed IFPV provides a low parametric representation of isodose plan without compromising the essential characteristics of the plan, thus providing a practically valuable framework for various applications in radiation therapy.
Collapse
Affiliation(s)
- Hongcheng Liu
- Department of Industrial and Systems Engineering, University of Florida, Gainesville, FL, 32611-6595, USA
| | - Lei Xing
- Department of Radiation Oncology, Stanford University School of Medicine, Stanford, CA, 94305-5847, USA
| |
Collapse
|
10
|
Liu H, Dong P, Xing L. Using measurable dosimetric quantities to characterize the inter-structural tradeoff in inverse planning. Phys Med Biol 2017; 62:6804-6821. [PMID: 28447959 DOI: 10.1088/1361-6560/aa6fcb] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Traditional inverse planning relies on the use of weighting factors to balance the conflicting requirements of different structures. Manual trial-and-error determination of weighting factors has long been recognized as a time-consuming part of treatment planning. The purpose of this work is to develop an inverse planning framework that parameterizes the dosimetric tradeoff among the structures with physically meaningful quantities to simplify the search for clinically sensible plans. In this formalism, instead of using weighting factors, the permissible variation range of the prescription dose or dose volume histogram (DVH) of the involved structures are used to characterize the 'importance' of the structures. The inverse planning is then formulated into a convex feasibility problem, called the dosimetric variation-controlled model (DVCM), whose goal is to generate plans with dosimetric or DVH variations of the structures consistent with the pre-specified values. For simplicity, the dosimetric variation range for a structure is extracted from a library of previous cases which possess similar anatomy and prescription. A two-phase procedure (TPP) is designed to solve the model. The first phase identifies a physically feasible plan to satisfy the prescribed dosimetric variation, and the second phase automatically improves the plan in case there is room for further improvement. The proposed technique is applied to plan two prostate cases and two head-and-neck cases and the results are compared with those obtained using a conventional CVaR approach and with a moment-based optimization scheme. Our results show that the strategy is able to generate clinically sensible plans with little trial and error. In all cases, the TPP generates a very competitive plan as compared to those obtained using the alternative approaches. Particularly, in the planning of one of the head-and-neck cases, the TPP leads to a non-trivial improvement in the resultant dose distribution-the fractional volumes receiving a dose above 20 Gy for the spinal cord are reduced by more than 40% when compared to the alternative schemes, while maintaining the same PTV coverage. With physically more meaningful modeling of the inter-structural tradeoff, the reported technique enables us to substantially reduce the need for trial-and-error adjustment of the model parameters. The new formalism also opens new opportunities for incorporating prior knowledge to facilitate the treatment planning process.
Collapse
Affiliation(s)
- Hongcheng Liu
- Department of Radiation Oncology, Stanford University, 875 Blake Wilbur Drive, Stanford, CA 94305-5847, United States of America
| | | | | |
Collapse
|
11
|
Engberg L, Forsgren A, Eriksson K, Hårdemark B. Explicit optimization of plan quality measures in intensity-modulated radiation therapy treatment planning. Med Phys 2017; 44:2045-2053. [DOI: 10.1002/mp.12146] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2016] [Revised: 01/25/2017] [Accepted: 01/25/2017] [Indexed: 11/08/2022] Open
Affiliation(s)
- Lovisa Engberg
- Optimization and Systems Theory, Department of Mathematics; KTH Royal Institute of Technology; Stockholm SE-100 44 Sweden
| | - Anders Forsgren
- Optimization and Systems Theory, Department of Mathematics; KTH Royal Institute of Technology; Stockholm SE-100 44 Sweden
| | - Kjell Eriksson
- RaySearch Laboratories; Sveavägen 44 Stockholm SE-103 65 Sweden
| | - Björn Hårdemark
- RaySearch Laboratories; Sveavägen 44 Stockholm SE-103 65 Sweden
| |
Collapse
|
12
|
McGeachy P, Villarreal-Barajas JE, Zinchenko Y, Khan R. Modulated photon radiotherapy (XMRT): an algorithm for the simultaneous optimization of photon beamlet energy and intensity in external beam radiotherapy (EBRT) planning. Phys Med Biol 2016; 61:1476-98. [PMID: 26808280 DOI: 10.1088/0031-9155/61/4/1476] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
This is a proof of principle study on an algorithm for optimizing external beam radiotherapy in terms of both photon beamlet energy and fluence. This simultaneous beamlet energy and fluence optimization is denoted modulated photon radiotherapy (XMRT). XMRT is compared with single-energy intensity modulated radiotherapy (IMRT) for five clinically relevant test geometries to determine whether treating beamlet energy as a decision variable improves the dose distributions. All test geometries were modelled in a cylindrical water phantom. XMRT optimized the fluence for 6 and 18 MV beamlets while IMRT optimized with only 6 MV and only 18 MV. CERR (computational environment for radiotherapy research) was used to calculate the dose deposition matrices and the resulting dose for XMRT and IMRT solutions. Solutions were compared via their dose volume histograms and dose metrics, such as the mean, maximum, and minimum doses for each structure. The homogeneity index (HI) and conformity number (CN) were calculated to assess the quality of the target dose coverage. Complexity of the resulting fluence maps was minimized using the sum of positive gradients technique. The results showed XMRT's ability to improve healthy-organ dose reduction while yielding comparable coverage of the target relative to IMRT for all geometries. All three energy-optimization approaches yielded similar HI and CNs for all geometries, as well as a similar degree of fluence map complexity. The dose reduction provided by XMRT was demonstrated by the relative decrease in the dose metrics for the majority of the organs at risk (OARs) in all geometries. Largest reductions ranged between 5% to 10% in the mean dose to OARs for two of the geometries when compared with both single-energy IMRT schemes. XMRT has shown potential dosimetric benefits through improved OAR sparing by allowing beam energy to act as a degree of freedom in the EBRT optimization process.
Collapse
Affiliation(s)
- Philip McGeachy
- Department of Physics and Astronomy, University of Calgary, Calgary, AB T2N 1N4, Canada. Department of Medical Physics, Tom Baker Cancer Centre, Calgary, AB T2N 4N2, Canada
| | | | | | | |
Collapse
|
13
|
Zarepisheh M, Uribe-Sanchez AF, Li N, Jia X, Jiang SB. A multicriteria framework with voxel-dependent parameters for radiotherapy treatment plan optimization. Med Phys 2014; 41:041705. [PMID: 24694125 DOI: 10.1118/1.4866886] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To establish a new mathematical framework for radiotherapy treatment optimization with voxel-dependent optimization parameters. METHODS In the treatment plan optimization problem for radiotherapy, a clinically acceptable plan is usually generated by an optimization process with weighting factors or reference doses adjusted for a set of the objective functions associated to the organs. Recent discoveries indicate that adjusting parameters associated with each voxel may lead to better plan quality. However, it is still unclear regarding the mathematical reasons behind it. Furthermore, questions about the objective function selection and parameter adjustment to assure Pareto optimality as well as the relationship between the optimal solutions obtained from the organ-based and voxel-based models remain unanswered. To answer these questions, the authors establish in this work a new mathematical framework equipped with two theorems. RESULTS The new framework clarifies the different consequences of adjusting organ-dependent and voxel-dependent parameters for the treatment plan optimization of radiation therapy, as well as the impact of using different objective functions on plan qualities and Pareto surfaces. The main discoveries are threefold: (1) While in the organ-based model the selection of the objective function has an impact on the quality of the optimized plans, this is no longer an issue for the voxel-based model since the Pareto surface is independent of the objective function selection and the entire Pareto surface could be generated as long as the objective function satisfies certain mathematical conditions; (2) All Pareto solutions generated by the organ-based model with different objective functions are parts of a unique Pareto surface generated by the voxel-based model with any appropriate objective function; (3) A much larger Pareto surface is explored by adjusting voxel-dependent parameters than by adjusting organ-dependent parameters, possibly allowing for the generation of plans with better trade-offs among different clinical objectives. CONCLUSIONS The authors have developed a mathematical framework for radiotherapy treatment optimization using voxel-based parameters. The authors can improve the plan quality by adjusting voxel-based weighting factors and exploring the unique and large Pareto surface which include all the Pareto surfaces that can be generated by organ-based model using different objective functions.
Collapse
Affiliation(s)
- Masoud Zarepisheh
- Center for Advanced Radiotherapy Technologies and Department of Radiation Medicine and Applied Sciences, University of California San Diego, La Jolla, California 92037-0843
| | - Andres F Uribe-Sanchez
- Center for Advanced Radiotherapy Technologies and Department of Radiation Medicine and Applied Sciences, University of California San Diego, La Jolla, California 92037-0843
| | - Nan Li
- Center for Advanced Radiotherapy Technologies and Department of Radiation Medicine and Applied Sciences, University of California San Diego, La Jolla, California 92037-0843
| | - Xun Jia
- Center for Advanced Radiotherapy Technologies and Department of Radiation Medicine and Applied Sciences, University of California San Diego, La Jolla, California 92037-0843
| | - Steve B Jiang
- Center for Advanced Radiotherapy Technologies and Department of Radiation Medicine and Applied Sciences, University of California San Diego, La Jolla, California 92037-0843
| |
Collapse
|
14
|
Kalantzis G, Apte A. A novel reduced-order prioritized optimization method for radiation therapy treatment planning. IEEE Trans Biomed Eng 2014; 61:1062-70. [PMID: 24658231 DOI: 10.1109/tbme.2013.2293779] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
In this study, a novel reduced order prioritized algorithm is presented for optimization in radiation therapy treatment planning. The proposed method consists of three stages. In the first stage, the intensity space was sampled by solving a series of unconstrained optimization problems. The objective function of the first stage is expressed as a scalarized weighted sum of partial objectives for the target and organ at risk. Latin hypercube sampling was utilized to define the weights for each run of the unconstrained optimizations. In the second stage, principal component analysis is applied to the solutions determined in the first stage to identify the major eigen modes in the intensities space, significantly reducing the number of independent variables. In the third stage, treatment planning goals/objectives are prioritized, and the problem is solved in the reduced order space. After each objective is optimized, that objective function is converted into a constraint for the lower-priority objectives. In the current formulation, a slip factor is used to relax the hard constraints for planning target volume (PTV) coverage. The applicability of the proposed method is demonstrated for one prostate and one lung intensity-modulated radiation therapy treatment plan. Upon completion of the sequential prioritized optimization, the mean dose at the rectum and bladder was reduced by 21.3% and 22.4%, respectively. Additionally, we investigated the effect of the slip factor 's' on PTV coverage and we found minimal degradation of the tumor dose (∼4%). Finally, the speed up factors upon the dimensionality reduction were as high as 49.9 without compromising the quality of the results.
Collapse
|
15
|
Long T, Matuszak M, Feng M, Fraass BA, Ten Haken RK, Romeijn HE. Sensitivity analysis for lexicographic ordering in radiation therapy treatment planning. Med Phys 2012; 39:3445-55. [PMID: 22755724 DOI: 10.1118/1.4720218] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To introduce a method to efficiently identify and calculate meaningful tradeoffs between criteria in an interactive IMRT treatment planning procedure. The method provides a systematic approach to developing high-quality radiation therapy treatment plans. METHODS Treatment planners consider numerous dosimetric criteria of varying importance that, when optimized simultaneously through multicriteria optimization, yield a Pareto frontier which represents the set of Pareto-optimal treatment plans. However, generating and navigating this frontier is a time-consuming, nontrivial process. A lexicographic ordering (LO) approach to IMRT uses a physician's criteria preferences to partition the treatment planning decisions into a multistage treatment planning model. Because the relative importance of criteria optimized in the different stages may not necessarily constitute a strict prioritization, the authors introduce an interactive process, sensitivity analysis in lexicographic ordering (SALO), to allow the treatment planner control over the relative sequential-stage tradeoffs. By allowing this flexibility within a structured process, SALO implicitly restricts attention to and allows exploration of a subset of the Pareto efficient frontier that the physicians have deemed most important. RESULTS Improvements to treatment plans over a LO approach were found by implementing the SALO procedure on a brain case and a prostate case. In each stage, a physician assessed the tradeoff between previous stage and current stage criteria. The SALO method provided critical tradeoff information through curves approximating the relationship between criteria, which allowed the physician to determine the most desirable treatment plan. CONCLUSIONS The SALO procedure provides treatment planners with a directed, systematic process to treatment plan selection. By following a physician's prioritization, the treatment planner can avoid wasting effort considering clinically inferior treatment plans. The planner is guided by criteria importance, but given the information necessary to accurately adjust the relative importance at each stage. Through these attributes, the SALO procedure delivers an approach well balanced between efficiency and flexibility.
Collapse
Affiliation(s)
- T Long
- Department of Industrial and Operations Engineering, University of Michigan, 1205 Beal Avenue, Ann Arbor, MI 48109-2117, USA
| | | | | | | | | | | |
Collapse
|
16
|
Abstract
Despite many studies over the last 3 decades that have attempted to explicitly quantify the decision-making process for radiotherapy treatment plan evaluation, judgments of an individual plan's degree of quality are still largely subjective and can show inter- and intra-practitioner variability even if the clinical treatment goals are the same. Several factors conspire to confound the full quantification of treatment plan quality, including uncertainties in dose response of cancerous and normal tissue, the rapid pace of new technology adoption, and the human component of treatment planning. However, new developments in clinical informatics and automation are lowering the bar for developing and implementing quantitative metrics into the treatment planning process. This review discusses general strategies for using quantitative metrics in the treatment planning process and presents a case study in intensity-modulated radiation therapy planning whereby control was established on a variable system via such techniques.
Collapse
Affiliation(s)
- Kevin L Moore
- Department of Radiation Oncology, Washington University School of Medicine, St Louis, MO 63110, USA.
| | | | | | | |
Collapse
|
17
|
Prioritized optimization in intensity modulated proton therapy. Z Med Phys 2012; 22:21-8. [DOI: 10.1016/j.zemedi.2011.05.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2011] [Revised: 05/06/2011] [Accepted: 05/11/2011] [Indexed: 11/20/2022]
|
18
|
Craft DL, Hong TS, Shih HA, Bortfeld TR. Improved planning time and plan quality through multicriteria optimization for intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys 2011; 82:e83-90. [PMID: 21300448 DOI: 10.1016/j.ijrobp.2010.12.007] [Citation(s) in RCA: 132] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2010] [Revised: 11/30/2010] [Accepted: 12/07/2010] [Indexed: 10/18/2022]
Abstract
PURPOSE To test whether multicriteria optimization (MCO) can reduce treatment planning time and improve plan quality in intensity-modulated radiotherapy (IMRT). METHODS AND MATERIALS Ten IMRT patients (5 with glioblastoma and 5 with locally advanced pancreatic cancers) were logged during the standard treatment planning procedure currently in use at Massachusetts General Hospital (MGH). Planning durations and other relevant planning information were recorded. In parallel, the patients were planned using an MCO planning system, and similar planning time data were collected. The patients were treated with the standard plan, but each MCO plan was also approved by the physicians. Plans were then blindly reviewed 3 weeks after planning by the treating physician. RESULTS In all cases, the treatment planning time was vastly shorter for the MCO planning (average MCO treatment planning time was 12 min; average standard planning time was 135 min). The physician involvement time in the planning process increased from an average of 4.8 min for the standard process to 8.6 min for the MCO process. In all cases, the MCO plan was blindly identified as the superior plan. CONCLUSIONS This provides the first concrete evidence that MCO-based planning is superior in terms of both planning efficiency and dose distribution quality compared with the current trial and error-based IMRT planning approach.
Collapse
Affiliation(s)
- David L Craft
- Department of Radiation Oncology, Massachusetts General Hospital, Harvard Medical School, Boston, MA 02114, USA.
| | | | | | | |
Collapse
|
19
|
Experience-based quality control of clinical intensity-modulated radiotherapy planning. Int J Radiat Oncol Biol Phys 2011; 81:545-51. [PMID: 21277097 DOI: 10.1016/j.ijrobp.2010.11.030] [Citation(s) in RCA: 240] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2010] [Revised: 10/12/2010] [Accepted: 11/16/2010] [Indexed: 12/22/2022]
Abstract
PURPOSE To incorporate a quality control tool, according to previous planning experience and patient-specific anatomic information, into the intensity-modulated radiotherapy (IMRT) plan generation process and to determine whether the tool improved treatment plan quality. METHODS AND MATERIALS A retrospective study of 42 IMRT plans demonstrated a correlation between the fraction of organs at risk (OARs) overlapping the planning target volume and the mean dose. This yielded a model, predicted dose = prescription dose (0.2 + 0.8 [1 - exp(-3 overlapping planning target volume/volume of OAR)]), that predicted the achievable mean doses according to the planning target volume overlap/volume of OAR and the prescription dose. The model was incorporated into the planning process by way of a user-executable script that reported the predicted dose for any OAR. The script was introduced to clinicians engaged in IMRT planning and deployed thereafter. The script's effect was evaluated by tracking δ = (mean dose-predicted dose)/predicted dose, the fraction by which the mean dose exceeded the model. RESULTS All OARs under investigation (rectum and bladder in prostate cancer; parotid glands, esophagus, and larynx in head-and-neck cancer) exhibited both smaller δ and reduced variability after script implementation. These effects were substantial for the parotid glands, for which the previous δ = 0.28 ± 0.24 was reduced to δ = 0.13 ± 0.10. The clinical relevance was most evident in the subset of cases in which the parotid glands were potentially salvageable (predicted dose <30 Gy). Before script implementation, an average of 30.1 Gy was delivered to the salvageable cases, with an average predicted dose of 20.3 Gy. After implementation, an average of 18.7 Gy was delivered to salvageable cases, with an average predicted dose of 17.2 Gy. In the prostate cases, the rectum model excess was reduced from δ = 0.28 ± 0.20 to δ = 0.07 ± 0.15. On surveying dosimetrists at the end of the study, most reported that the script both improved their IMRT planning (8 of 10) and increased their efficiency (6 of 10). CONCLUSIONS This tool proved successful in increasing normal tissue sparing and reducing interclinician variability, providing effective quality control of the IMRT plan development process.
Collapse
|
20
|
Zinchenko Y, Craig T, Keller H, Terlaky T, Sharpe M. Controlling the dose distribution with gEUD-type constraints within the convex radiotherapy optimization framework. Phys Med Biol 2008; 53:3231-50. [DOI: 10.1088/0031-9155/53/12/011] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|