Search

Objective. Reducing plan complexity in intensity modulated radiation therapy (IMRT) to ensure dosimetric accuracy and delivery efficiency of the radiation treatment plans. We propose a novel approach by representing the beamlet intensities using an incomplete wavelet basis that explicitly excludes fluctuating intensity maps from the decision space (explicit hard constraint). This technique provides a built-in wavelet-induced smoothness that improves both dosimetric plan quality and delivery efficiency. Approach. The beamlet intensity maps need to be especially smooth in the leaf travel direction (referred to as the X-direction). We treat the intensity map of each beam as a 2D image and represent it using the wavelets corresponding to low-frequency changes in the X-direction (i.e. approximation and horizontal). The absence of wavelets corresponding to high-frequency changes (i.e. vertical and diagonal) induces built-in smoothness. We still utilize a regularization term in the objective function to promote smoothness in the Y-direction (perpendicular to the X-direction) and further possible smoothness in the X-direction. This technique has been tested on three patient cases of different disease sites (paraspinal, lung, prostate) and all final evaluations and comparisons have been performed on an FDA-approved commercial treatment planning system (Varian EclipseTM). Main results. Wavelet-induced smoothness reduced monitor units by about 10%, 45%, and 14% for paraspinal, lung, and prostate cases, respectively. It also improved organ at risk sparing, especially on the complex paraspinal case where it resulted in about 7%, 13%, and 14% less mean dose to esophagus, lung, and cord, respectively. Moreover, built-in wavelet-induced smoothness desensitizes the results to changing the weight associated to the regularization term, and thereby mitigates the weight fine-tuning difficulty. Significance. Fluence smoothness is often achieved by smoothing the beamlet intensity maps using a proper regularization term in the objective function aiming at disincentivizing fluctuation in the beamlet intensities (implicit soft constraint). This work reports a novel application of wavelets in imposing an explicit smoothness hard constraint in the search space using an incomplete wavelet basis. This idea has been successfully applied to exclude complex and clinically irrelevant radiation plans from the search space. The code and pertained models along with a sample dataset are released on our LowDimRT GitHub (https://github.com/PortPy-Project/LowDimRT).

Purpose We present a method for fully automated generation of high quality robust proton treatment plans using hierarchical optimization. To fill the gap between the two common extreme robust optimization approaches, that is, stochastic and worst-case, a robust optimization approach based on the p-norm function is used whereby a single parameter, p , can be used to control the level of robustness in an intuitive way. Methods A fully automated approach to treatment planning using Expedited Constrained Hierarchical Optimization (ECHO) is implemented in our clinic for photon treatments. ECHO strictly enforces critical (inviolable) clinical criteria as hard constraints and improves the desirable clinical criteria sequentially, as much as is feasible. We extend our in-house developed ECHO codes for proton therapy and integrate it with a new approach for robust optimization. Multiple scenarios accounting for both setup and range uncertainties are included (13scenarios), and the maximum/mean/dose-volume constraints on organs-at-risk (OARs) and target are fulfilled in all scenarios. We combine the objective functions of the individual scenarios using the p-norm function. The p-norm with a parameter p = 1 or p = ∞ result in the stochastic or the worst-case approach, respectively; an intermediate robustness level is obtained by employing p -values in-between. While the worst-case approach only focuses on the worst-case scenario(s), the p-norm approach with a large p value ( p ≈ 20 ) resembles the worst-case approach without completely neglecting other scenarios. The proposed approach is evaluated on three head-and-neck (HN) patients and one water phantom with different parameters, p ∈ 1 , 2 , 5 , 10 , 20 . The results are compared against the stochastic approach (p-norm approach with p = 1 ) and the worst-case approach, as well as the nonrobust approach (optimized solely on the nominal scenario). Results The proposed algorithm successfully generates automated robust proton plans on all cases. As opposed to the nonrobust plans, the robust plans have narrower dose volume histogram (DVH) bands across all 13 scenarios, and meet all hard constraints (i.e., maximum/mean/dose-volume constraints) on OARs and the target for all scenarios. The spread in the objective function values is largest for the stochastic approach ( p = 1 ) and decreases with increasing p toward the worst-case approach. Compared to the worst-case approach, the p-norm approach results in DVH bands for clinical target volume (CTV) which are closer to the prescription dose at a negligible cost in the DVH for the worst scenario, thereby improving the overall plan quality. On average, going from the worst-case approach to the p-norm approach with p = 20 , the median objective function value across all the scenarios is improved by 15% while the objective function value for the worst scenario is only degraded by 3%. Conclusion An automated treatment planning approach for proton therapy is developed, including robustness, dose-volume constraints, and the ability to control the robustness level using the p-norm parameter p , to fit the priorities deemed most important.

PURPOSE To present a fully automated treatment planning process for proton therapy including beam angle selection using a novel Bayesian optimization approach and previously developed constrained hierarchical fluence optimization method. METHODS We adapted our in-house automated intensity modulated radiation therapy (IMRT) treatment planning system, which is based on constrained hierarchical optimization and referred to as ECHO (expedited constrained hierarchical optimization), for proton therapy. To couple this to beam angle selection, we propose using a novel Bayesian approach. By integrating ECHO with this Bayesian beam selection approach, we obtain a fully automated treatment planning framework including beam angle selection. Bayesian optimization is a global optimization technique which only needs to search a small fraction of the search space for slowly varying objective functions (i.e., smooth functions). Expedited constrained hierarchical optimization is run for some initial beam angle candidates and the resultant treatment plan for each beam configuration is rated using a clinically relevant treatment score function. Bayesian optimization iteratively predicts the treatment score for not-yet-evaluated candidates to find the best candidate to be optimized next with ECHO. We tested this technique on five head-and-neck (HN) patients with two coplanar beams. In addition, tests were performed with two noncoplanar and three coplanar beams for two patients. RESULTS For the two coplanar configurations, the Bayesian optimization found the optimal beam configuration after running ECHO for, at most, 4% of all potential configurations (23 iterations) for all patients (range: 2%-4%). Compared with the beam configurations chosen by the planner, the optimal configurations reduced the mandible maximum dose by 6.6 Gy and high dose to the unspecified normal tissues by 3.8 Gy, on average. For the two noncoplanar and three coplanar beam configurations, the algorithm converged after 45 iterations (examining

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.

Purpose Dose-volume constraints (DVCs) continue to be common features in intensity-modulated radiation therapy (IMRT) prescriptions, but they are non-convex and difficult to incorporate. We propose computationally efficient methods to incorporate dose-volume constraints (DVCs) into automated IMRT planning. Methods We propose a two-phase approach: in phase-1, we solve a convex approximation with DVCs. Although this convex approximation does not guarantee DVC satisfaction, it provides crucial initial information about voxels likely to receive doses below DVC thresholds. Subsequently, phase-2 solves an optimization problem with maximum dose constraints imposed on those subthreshold voxels. We further categorize DVCs into hard- and soft-DVCs, where hard-DVCs are strictly enforced by the optimization and soft-DVCs are encouraged in the objective function. We tested this approach in our automated treatment planning system which is based on hierarchical constrained optimization. Performance is demonstrated on a series of paraspinal, lung, oligometastasis, and prostate cases as well as a small paraspinal case for which we can computationally afford to obtain a ground-truth by solving a non-convex optimization problem. Results The proposed algorithm successfully meets all the hard-DVCs while increasing the overall computational time of the baseline planning process (without DVCs) by 20%, 10%, and 11% for paraspinal, oligometastasis, and prostate cases, respectively. For a soft-DVC applied to the lung case, the dose-volume histogram curve moves toward the desired direction and the computational time is increased by 11%. For a low-resolution paraspinal case, the ground-truth solution process using mixed-integer programming methods required 15 h while the proposed algorithm converges in only 2 min with a proximal solution. Conclusions A computationally tractable algorithm to handle hard- and soft-DVCs is developed which is capable of satisfying DVCs without any parameter tweaking. Although the algorithm is demonstrated in our in-house developed automated treatment planning system, it can potentially be used in any constrained optimization framework.

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.

Purpose: Station parameter optimized radiation therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital linear accelerators, in which the station parameters of a delivery system, such as aperture shape and weight, couch position/angle, gantry/collimator angle, can be optimized simultaneously. SPORT promises to deliver remarkable radiation dose distributions in an efficient manner, yet there exists no optimization algorithm for its implementation. The purpose of this work is to develop an algorithm to simultaneously optimize the beam sampling and aperture shapes. Methods: The authors build a mathematical model with the fundamental station point parameters as the decision variables. To solve the resulting large-scale optimization problem, the authors devise an effective algorithm by integrating three advanced optimization techniques: column generation, subgradient method, and pattern search. Column generation adds the most beneficial stations sequentially until the plan quality improvement saturates and provides a good starting point for the subsequent optimization. It also adds the new stations during the algorithm if beneficial. For each update resulted from column generation, the subgradient method improves the selected stations locally by reshaping the apertures and updating the beam angles toward a descent subgradient direction. The algorithm continues to improve the selected stations locallymore » and globally by a pattern search algorithm to explore the part of search space not reachable by the subgradient method. By combining these three techniques together, all plausible combinations of station parameters are searched efficiently to yield the optimal solution. Results: A SPORT optimization framework with seamlessly integration of three complementary algorithms, column generation, subgradient method, and pattern search, was established. The proposed technique was applied to two previously treated clinical cases: a head and neck and a prostate case. It significantly improved the target conformality and at the same time critical structure sparing compared with conventional intensity modulated radiation therapy (IMRT). In the head and neck case, for example, the average PTV coverage D99% for two PTVs, cord and brainstem max doses, and right parotid gland mean dose were improved, respectively, by about 7%, 37%, 12%, and 16%. Conclusions: The proposed method automatically determines the number of the stations required to generate a satisfactory plan and optimizes simultaneously the involved station parameters, leading to improved quality of the resultant treatment plans as compared with the conventional IMRT plans.« less

Purpose:Adaptive radiation therapy(ART) can reduce normal tissue toxicity and/or improve tumorcontrol through treatment adaptations based on the current patient anatomy. Developing an efficient re‐planning algorithm is an important step for ART. For the re‐planning process, manual trial‐and‐error approach to fine‐tune planning parameters is a time‐ consuming and resource‐requiring task. However, prior information in the original plan, such as dose distribution or dose‐volume histogram (DVH) can be employed to facilitate the re‐planning. The goal of this work is to develop a re‐planning algorithm that automates voxel weighting factor adjustments to generate a plan with close, or possibly better, DVH curves compared with original plan. Methods: Our algorithm iterates the following two loops. The inner loop is the traditional fluence map optimization, in which we optimize a quadratic objective function penalizing the deviation of the dose received by each voxel from its prescribed or threshold dose with a set of fixed weighting factors. In outer loop, the weighting factors in the objective function for each voxel are heuristically adjusted according to the deviation of the DVH curves in the calculated plan from those in the original plan. The process is repeated until the result converges, or the maximum iteration step is reached. The whole algorithm is implemented on GPU for high efficiency. Results: We have tested our algorithm on three 8‐field head‐and‐neck cases. Compared with the DVH curves in original plan, the DVH curves in the resulting plan using our algorithm with 30 iterations are almost better for every structure. We can finish the re‐optimization process around 30 seconds. Conclusions: Adjusting voxel‐based weighting factors automatically, by comparing the DVH curves, seems to be a promising approach to avoid the tedious trial‐and‐error scheme for ART re‐planning. This work is supported by Varian Medical Systems through a Master Research Agreement.

This article deals with the scalarization techniques in multi-objective optimization problems (MOPs). In this article, the nonlinear version of the weighted sum method is considered as the scalarization problem, and the relation of the set of properly efficient points and the set of scalarization solutions (optimal solutions of the scalarized problem) is studied. In particular, conditions are determined under which the optimal solutions of the scalarized problem are properly efficient points of the MOP, and the conditions under which the properly efficient points could be characterized in terms of the scalarization solutions are presented.

Intensity-modulated radiation therapy (IMRT) currently plays an important role in radiotherapy, but its treatment plan quality can vary significantly among institutions and planners. Treatment plan quality control (QC) is a necessary component for individual clinics to ensure that patients receive treatments with high therapeutic gain ratios. The voxel-weighting factor-based plan re-optimization mechanism has been proved able to explore a larger Pareto surface (solution domain) and therefore increase the possibility of finding an optimal treatment plan. In this study, we incorporated additional modules into an in-house developed voxel weighting factor-based re-optimization algorithm, which was enhanced as a highly automated and accurate IMRT plan QC tool (TPS-QC tool). After importing an under-assessment plan, the TPS-QC tool was able to generate a QC report within 2 minutes. This QC report contains the plan quality determination as well as information supporting the determination. Finally, the IMRT plan quality can be controlled by approving quality-passed plans and replacing quality-failed plans using the TPS-QC tool. The feasibility and accuracy of the proposed TPS-QC tool were evaluated using 25 clinically approved cervical cancer patient IMRT plans and 5 manually created poor-quality IMRT plans. The results showed high consistency between the QC report quality determinations and the actual plan quality. In the 25 clinically approved cases that the TPS-QC tool identified as passed, a greater difference could be observed for dosimetric endpoints for organs at risk (OAR) than for planning target volume (PTV), implying that better dose sparing could be achieved in OAR than in PTV. In addition, the dose-volume histogram (DVH) curves of the TPS-QC tool re-optimized plans satisfied the dosimetric criteria more frequently than did the under-assessment plans. In addition, the criteria for unsatisfied dosimetric endpoints in the 5 poor-quality plans could typically be satisfied when the TPS-QC tool generated re-optimized plans without sacrificing other dosimetric endpoints. In addition to its feasibility and accuracy, the proposed TPS-QC tool is also user-friendly and easy to operate, both of which are necessary characteristics for clinical use.

Station parameter optimized radiation therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital linear accelerators, in which the station parameters of a delivery system, such as aperture shape and weight, couch position/angle, gantry/collimator angle, can be optimized simultaneously. SPORT promises to deliver remarkable radiation dose distributions in an efficient manner, yet there exists no optimization algorithm for its implementation. The purpose of this work is to develop an algorithm to simultaneously optimize the beam sampling and aperture shapes.The authors build a mathematical model with the fundamental station point parameters as the decision variables. To solve the resulting large-scale optimization problem, the authors devise an effective algorithm by integrating three advanced optimization techniques: column generation, subgradient method, and pattern search. Column generation adds the most beneficial stations sequentially until the plan quality improvement saturates and provides a good starting point for the subsequent optimization. It also adds the new stations during the algorithm if beneficial. For each update resulted from column generation, the subgradient method improves the selected stations locally by reshaping the apertures and updating the beam angles toward a descent subgradient direction. The algorithm continues to improve the selected stations locally and globally by a pattern search algorithm to explore the part of search space not reachable by the subgradient method. By combining these three techniques together, all plausible combinations of station parameters are searched efficiently to yield the optimal solution.A SPORT optimization framework with seamlessly integration of three complementary algorithms, column generation, subgradient method, and pattern search, was established. The proposed technique was applied to two previously treated clinical cases: a head and neck and a prostate case. It significantly improved the target conformality and at the same time critical structure sparing compared with conventional intensity modulated radiation therapy (IMRT). In the head and neck case, for example, the average PTV coverage D99% for two PTVs, cord and brainstem max doses, and right parotid gland mean dose were improved, respectively, by about 7%, 37%, 12%, and 16%.The proposed method automatically determines the number of the stations required to generate a satisfactory plan and optimizes simultaneously the involved station parameters, leading to improved quality of the resultant treatment plans as compared with the conventional IMRT plans.

Stereotactic body radiotherapy (SBRT) shows promise in unresectable pancreatic cancer, though this treatment modality has high rates of normal tissue toxicity. This study explores the dosimetric utility of daily adaptive re-planning with pancreas SBRT. We used a previously developed supercomputing online re-planning environment (SCORE) to re-plan 10 patients with pancreas SBRT. Tumor and normal tissue contours were deformed from treatment planning computed tomographies (CTs) and transferred to daily cone-beam CT (CBCT) scans before re-optimizing each daily treatment plan. We compared the intended radiation dose, the actual radiation dose, and the optimized radiation dose for the pancreas tumor planning target volume (PTV) and the duodenum. Treatment re-optimization improved coverage of the PTV and reduced dose to the duodenum. Within the PTV, the actual hot spot (volume receiving 110% of the prescription dose) decreased from 4.5% to 0.5% after daily adaptive re-planning. Within the duodenum, the volume receiving the prescription dose decreased from 0.9% to 0.3% after re-planning. It is noteworthy that variation in the amount of air within a patient׳s stomach substantially changed dose to the PTV. Adaptive re-planning with pancreas SBRT has the ability to improve dose to the tumor and decrease dose to the nearby duodenum, thereby reducing the risk of toxicity.

To develop a novel algorithm that incorporates prior treatment knowledge into intensity modulated radiation therapy optimization to facilitate automatic treatment planning and adaptive radiotherapy (ART) replanning.The algorithm automatically creates a treatment plan guided by the DVH curves of a reference plan that contains information on the clinician-approved dose-volume trade-offs among different targets/organs and among different portions of a DVH curve for an organ. In ART, the reference plan is the initial plan for the same patient, while for automatic treatment planning the reference plan is selected from a library of clinically approved and delivered plans of previously treated patients with similar medical conditions and geometry. The proposed algorithm employs a voxel-based optimization model and navigates the large voxel-based Pareto surface. The voxel weights are iteratively adjusted to approach a plan that is similar to the reference plan in terms of the DVHs. If the reference plan is feasible but not Pareto optimal, the algorithm generates a Pareto optimal plan with the DVHs better than the reference ones. If the reference plan is too restricting for the new geometry, the algorithm generates a Pareto plan with DVHs close to the reference ones. In both cases, the new plans have similar DVH trade-offs as the reference plans.The algorithm was tested using three patient cases and found to be able to automatically adjust the voxel-weighting factors in order to generate a Pareto plan with similar DVH trade-offs as the reference plan. The algorithm has also been implemented on a GPU for high efficiency.A novel prior-knowledge-based optimization algorithm has been developed that automatically adjust the voxel weights and generate a clinical optimal plan at high efficiency. It is found that the new algorithm can significantly improve the plan quality and planning efficiency in ART replanning and automatic treatment planning.

Adaptive radiation therapy (ART) can reduce normal tissue toxicity and/or improve tumor control through treatment adaptations based on the current patient anatomy. Developing an efficient and effective re-planning algorithm is an important step toward the clinical realization of ART. For the re-planning process, manual trial-and-error approach to fine-tune planning parameters is time-consuming and is usually considered unpractical, especially for online ART. It is desirable to automate this step to yield a plan of acceptable quality with minimal interventions. In ART, prior information in the original plan is available, such as dose?volume histogram (DVH), which can be employed to facilitate the automatic re-planning process. The goal of this work is to develop an automatic re-planning algorithm to generate a plan with similar, or possibly better, DVH curves compared with the clinically delivered original plan. Specifically, our algorithm iterates the following two loops. An inner loop is the traditional fluence map optimization, in which we optimize a quadratic objective function penalizing the deviation of the dose received by each voxel from its prescribed or threshold dose with a set of fixed voxel weighting factors. In outer loop, the voxel weighting factors in the objective function are adjusted according to the deviation of the current DVH curves from those in the original plan. The process is repeated until the DVH curves are acceptable or maximum iteration step is reached. The whole algorithm is implemented on GPU for high efficiency. The feasibility of our algorithm has been demonstrated with three head-and-neck cancer IMRT cases, each having an initial planning CT scan and another treatment CT scan acquired in the middle of treatment course. Compared with the DVH curves in the original plan, the DVH curves in the resulting plan using our algorithm with 30 iterations are better for almost all structures. The re-optimization process takes about 30?s using our in-house optimization engine.

Purpose: To propose a new parallel-friendly optimization algorithm, based on alternating direction method of multipliers (ADMM), for solving optimization problems emerging in treatment planning and imaging Methods: ADMM, emerging as a powerful tool for distributed optimization, is employed here for IMRT treatment planning. We modify the existing ADMM by integrating that with two optimization techniques known as Barzilai-Borwein gradient method and line search. We apply original and integrated ADMM, with various penalty parameter values, on three IMRT treatment planning cases and we compare their performance. Each algorithm terminates once it found a solution within the relative error of 1E-4 from the optimal solution. We also compare the performance of original and integrated ADMM against the commercial optimization software CPLEX. Results: Our experiments on three different cancer sites (prostate, head and neck, and GYN) indicate that integrated ADMM is much faster than original ADMM (5.8 times faster on average). Moreover, while the performance of original ADMM heavily depends on the selected penalty parameter, integrated ADMM performs very robust against the penalty parameter. Compared to the commercial software CPLEX, integrated ADMM turns out to be around 4–7 times faster. Conclusion: Integrated ADMM is a very efficient optimization algorithm to cope with the large-scale optimization problems. It is a parallel-friendly optimization algorithm and provides an idea platform for cloud-based treatment planning and imaging.

Purpose: Stereotactic body radiotherapy (SBRT) shows promise in unresectable pancreatic cancer with high rates of local control compared to conventionally fractionated radiotherapy. One critical limitation with pancreas SBRT arises from toxicity in the nearby duodenum, which is a dynamic and distensible organ that changes shape daily. The purpose of this pilot study was to evaluate the potential gains of daily adaptive re‐planning with SBRT in pancreatic cancer. Methods: Pancreas SBRT patients were re‐planned using our previously developed GPU‐based supercomputing online re‐planning environment (SCORE), which automatically and quickly generates plans from daily cone‐beam CT (CBCT) data with the original plan and DVH as a reference. For each patient, the tumor volume and normal structure contours were deformed from the original planning CT to each daily CBCT. The deformed contours were reviewed and edited by a physician, and then the plan was re‐optimized to account for daily changes in anatomy. The daily adaptive re‐planned dose was compared to the actual dose received for the target volume and nearby duodenum. Results: Fourteen daily CBCT scans were used to generate daily adaptive radiation plans. For all patients, the average target coverage for the PTV improved, and the cold spot decreased with daily adaptive re‐planning. The cold spot (minimum dose to 1% of PTV) decreased on average from 25.1% to 4.6% (p=0.002, paired t‐test). With daily adaptive re‐planning the hot spot in the duodenum (volume that received 100% of dose) decreased from 2.8% to 1.7% (p=0.09). Doses to other organs were clinically insignificant, and did not significantly change with daily adaptive radiotherapy. Conclusion: This pilot study shows that daily adaptive re‐planning with pancreas SBRT has the potential to improve target coverage, and reduce dose to normal organs. Further study into the technical implementation of daily adaptive re‐planning is ongoing. Varian Medical Systems through a Master Research Agreement

Purpose: To discover the geometric features that best identify an anatomically‐similar reference patient from a library of previous treatments to predict achievable dose distributions for automatic IMRT treatment planning of a new patient. Methods: Using a library of previously‐treated patients, two similarity functions are defined in the domain of the dose distributions and geometrical features that match a patient with their most similar member from the cohort. The goal is to attain the same patient matches from these similarity functions, so that when a new patient without dosimetric information is input into the geometric similarity function a similar patient from the library can be identified with good dosimetric predictive power for the new patient. The dosimetric similarity function is defined by a weighted sum of the difference area between two DVHs for all organs and PTV, with the weights tuned to reflect clinical priorities. The generalized geometric similarity function was defined as a nonlinear combination of organ/PTV volume, organ/PTV overlapping volume, and Dice index between a pair of two organs/PTVs/intersections. A genetic‐algorithm‐based learning technique was then employed to find the optimal combination of geometrical features to predict patient similarity. Results: For a database of 10 prostate cancer treated patients, a combination of bladder volume and (PTV, rectum) overlapping Dice index achieved 80% accuracy compared to the DVH similarity function. Validation on a library of 15 patients reduced the accuracy to 60%, implying that other geometric features could be needed. Conclusion: We have developed a method using machine learning to identify the best match of a new patient to the prior patients to guide automatic IMRT planning. Initial results from this feasibility study have modest predictive power but future work will focus in the incorporation of other geometric features, e.g. mutual information and overlap‐volume‐histograms, to increase selection accuracy.

Purpose: In treatment plan optimization for Intensity Modulated Radiation Therapy (IMRT), after a plan is initially developed by a dosimetrist, the attending physician evaluates its quality and often would like to improve it. As opposed to having the dosimetrist implement the improvements, it is desirable to have the physician directly and efficiently modify the plan for a more streamlined and effective workflow. In this project, we developed an interactive optimization system for physicians to conveniently and efficiently fine-tune iso-dose curves. Methods: An interactive interface is developed under C++/Qt. The physician first examines iso-dose lines. S/he then picks an iso-dose curve to be improved and drags it to a more desired configuration using a computer mouse or touchpad. Once the mouse is released, a voxel-based optimization engine is launched. The weighting factors corresponding to voxels between the iso-dose lines before and after the dragging are modified. The underlying algorithm then takes these factors as input to re-optimize the plan in near real-time on a GPU platform, yielding a new plan best matching the physician's desire. The re-optimized DVHs and iso-dose curves are then updated for the next iteration of modifications. This process is repeated until a physician satisfactory plan is achieved.more » Results: We have tested this system for a series of IMRT plans. Results indicate that our system provides the physicians an intuitive and efficient tool to edit the iso-dose curves according to their preference. The input information is used to guide plan re-optimization, which is achieved in near real-time using our GPU-based optimization engine. Typically, a satisfactory plan can be developed by a physician in a few minutes using this tool. Conclusion: With our system, physicians are able to manipulate iso-dose curves according to their preferences. Preliminary results demonstrate the feasibility and effectiveness of this tool.« less