Source author record

David Craft

David Craft appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

Fast approximate delivery of fluence maps: the single map case

In this first paper of a two-paper series, we present a method for optimizing the dynamic delivery of fluence maps in radiation therapy. For a given fluence map and a given delivery time, the optimization of the leaf trajectories of a multi-leaf collimator to approximately form the given fluence map is a non-convex optimization problem. Its general solution has not been addressed in the literature, despite the fact that dynamic delivery of fluence maps has long been a common approach to intensity modulated radiation therapy. We model the leaf trajectory and dose rate optimization as a non-convex continuous optimization problem and solve it by an interior point method from randomly initialized feasible starting solutions. We demonstrate the method on a fluence map from a prostate case and a larger fluence map from a head-and-neck case. While useful for static beam IMRT delivery, our main motivation for this work is the extension to the case of sequential fluence map delivery, i.e. the case of VMAT, which is the topic of the second paper.

preprint2014arXiv

3D conformal planning using low segment multi-criteria IMRT optimization

Purpose: To evaluate automated multicriteria optimization (MCO)-- designed for intensity modulated radiation therapy (IMRT), but invoked with limited segmentation -- to efficiently produce high quality 3D conformal treatment (3D-CRT) plans. Methods: Ten patients previously planned with 3D-CRT were replanned with a low-segment inverse multicriteria optimized technique. The MCO-3D plans used the same number of beams, beam geometry and machine parameters of the corresponding 3D plans, but were limited to an energy of 6 MV. The MCO-3D plans were optimized using a fluence-based MCO IMRT algorithm and then, after MCO navigation, segmented with a low number of segments. The 3D and MCO-3D plans were compared by evaluating mean doses to individual organs at risk (OARs), mean doses to combined OARs, homogeneity indexes (HI), monitor units (MUs), physician preference, and qualitative assessments of planning time and plan customizability. Results: The MCO-3D plans significantly reduced the OAR mean doses and monitor units while maintaining good coverage and homogeneity of target volumes. MCO allows for more streamlined plan customization. All MCO-3D plans were preferred by physicians over their corresponding 3D plans. Conclusion: High quality 3D plans can be produced using IMRT optimization technology, resulting in automated field-in-field type plans with good monitor unit efficiency. Adopting this technology in a clinic could streamline treatment plan production.

preprint2014arXiv

Averaging VMAT treatment plans for multi-criteria navigation

The main approach to smooth Pareto surface navigation for radiation therapy multi-criteria treatment planning involves taking real-time averages of pre-computed treatment plans. In fluence-based treatment planning, fluence maps themselves can be averaged, which leads to the dose distributions being averaged due to the linear relationship between fluence and dose. This works for fluence-based photon plans and proton spot scanning plans. In this technical note, we show that two or more sliding window volumetric modulated arc therapy (VMAT) plans can be combined by averaging leaf positions in a certain way, and we demonstrate that the resulting dose distribution for the averaged plan is approximately the average of the dose distributions of the original plans. This leads to the ability to do Pareto surface navigation, i.e. interactive multi-criteria exploration of VMAT plan dosimetric tradeoffs.

preprint2014arXiv

Exploiting tumor shrinkage through temporal optimization of radiotherapy

In multi-stage radiotherapy, a patient is treated in several stages separated by weeks or months. This regimen has been motivated mostly by radiobiological considerations, but also provides an approach to reduce normal tissue dose by exploiting tumor shrinkage. The paper considers the optimal design of multi-stage treatments, motivated by the clinical management of large liver tumors for which normal liver dose constraints prohibit the administration of an ablative radiation dose in a single treatment. We introduce a dynamic tumor model that incorporates three factors: radiation induced cell kill, tumor shrinkage, and tumor cell repopulation. The design of multi-stage radiotherapy is formulated as a mathematical optimization problem in which the total dose to the liver is minimized, subject to delivering the prescribed dose to the tumor. Based on the model, we gain insight into the optimal administration of radiation over time, i.e. the optimal treatment gaps and dose levels. We analyze treatments consisting of two stages in detail. The analysis confirms the intuition that the second stage should be delivered just before the tumor size reaches a minimum and repopulation overcompensates shrinking. Furthermore, it was found that, for a large range of model parameters, approximately one third of the dose should be delivered in the first stage. The projected benefit of multi-stage treatments depends on model assumptions. However, the model predicts large liver dose reductions by more than a factor of two for plausible model parameters. The analysis of the tumor model suggests that substantial reduction in normal tissue dose can be achieved by exploiting tumor shrinkage via an optimal design of multi-stage treatments. This suggests taking a fresh look at multi-stage radiotherapy for selected disease sites where substantial tumor regression translates into reduced target volumes.

preprint2013arXiv

Guaranteed epsilon-optimal treatment plans with minimum number of beams for stereotactic body radiation therapy

Stereotactic body radiotherapy (SBRT) is characterized by delivering a high amount of dose in a short period of time. In SBRT the dose is delivered using open fields (e.g., beam's-eye-view) known as "apertures". Mathematical methods can be used for optimizing treatment planning for delivery of sufficient dose to the cancerous cells while keeping the dose to surrounding organs at risk (OARs) minimal. Two important elements of a treatment plan are quality and delivery time. Quality of a plan is measured based on the target coverage and dose to OARs. Delivery time heavily depends on the number of beams used in the plan since the setup times for different beam directions constitute a large portion of the delivery time. Therefore the ideal plan, in which all potential beams can be used simultaneously, will be associated with a long impractical delivery time. We use the dose to OARs in the ideal plan to find the plan with the minimum number of beams which is guaranteed to be epsilon-optimal (i.e., a predetermined maximum deviation from the ideal plan is guaranteed). Since the treatment plan optimization is inherently a multicriteria optimization problem, the planner can navigate the ideal dose distribution Pareto surface and select a plan of desired target coverage versus OARs sparing, and then use the proposed technique to reduce the number of beams while guaranteeing epsilon-optimality. We use mixed integer programming (MIP) for optimization. To reduce the computation time for the resultant MIP, we use two heuristics: a beam elimination scheme and a family of heuristic cuts, known as "neighbor cuts", based on the concept of "adjacent beams". We show the effectiveness of the proposed solution technique on two clinical cases, a liver and a lung case. Based on our technique we propose an algorithm for fast generation of epsilon-optimal plans.

preprint2013arXiv

Multi-criteria optimization methods in radiation therapy planning: a review of technologies and directions

We review the field of multi-criteria optimization for radiation therapy treatment planning. Special attention is given to the technique known as Pareto surface navigation, which allows physicians and treatment planners to interactively navigate through treatment planning options to get an understanding of the tradeoffs (dose to the target versus over-dosing of important nearby organs) involved in each patient's plan. We also describe goal programming and prioritized optimization, two other methods designed to handle multiple conflicting objectives. Issues related to nonconvexities, both in terms of dosimetric goals and the fact that the mapping from controllable hardware parameters to patient doses is usually nonconvex, are discussed at length since nonconvexities have a large impact on practical solution techniques for Pareto surface construction and navigation. A general planning strategy is recommended which handles the issue of nonconvexity by first finding an ideal Pareto surface with radiation delivered from many preset angles. This can be cast as a convex optimization problem. Once a high quality solution is selected from the Pareto surface, a sparse version (which can mean fewer beams, fewer segments, less leaf travel for arc therapy techniques, etc.) is obtained using an appropriate sparsification heuristic. We end by discussing issues of efficiency regarding the planning and the delivery of radiation therapy.

preprint2013arXiv

Two Effective Heuristics for Beam Angle Optimization in Radiation Therapy

In radiation therapy, mathematical methods have been used for optimizing treatment planning for delivery of sufficient dose to the cancerous cells while keeping the dose to critical surrounding structures minimal. This optimization problem can be modeled using mixed integer programming (MIP) whose solution gives the optimal beam orientation as well as optimal beam intensity. The challenge, however, is the computation time for this large scale MIP. We propose and investigate two novel heuristic approaches to reduce the computation time considerably while attaining high-quality solutions. We introduce a family of heuristic cuts based on the concept of 'adjacent beams' and a beam elimination scheme based on the contribution of each beam to deliver the dose to the tumor in the ideal plan in which all potential beams can be used simultaneously. We show the effectiveness of these heuristics for intensity modulated radiation therapy (IMRT) and stereotactic body radiation therapy (SBRT) on a clinical liver case.

preprint2012arXiv

Deliverable navigation for multicriteria step and shoot IMRT treatment planning

We consider Pareto surface based multi-criteria optimization for step and shoot IMRT planning. By analyzing two navigation algorithms, we show both theoretically and in practice that the number of plans needed to form convex combinations of plans during navigation can be kept small (much less than the theoretical maximum number needed in general, which is equal to the number of objectives for on-surface Pareto navigation). Therefore a workable approach for directly deliverable navigation in this setting is to segment the underlying Pareto surface plans and then enforce the mild restriction that only a small number of these plans are active at any time during plan navigation, thus limiting the total number of segments used in the final plan.

preprint2012arXiv

Optimal partial-arcs in VMAT treatment planning

Purpose: To improve the delivery efficiency of VMAT by extending the recently published VMAT treatment planning algorithm vmerge to automatically generate optimal partial-arc plans. Methods and materials: A high-quality initial plan is created by solving a convex multicriteria optimization problem using 180 equi-spaced beams. This initial plan is used to form a set of dose constraints, and a set of partial-arc plans is created by searching the space of all possible partial-arc plans that satisfy these constraints. For each partial-arc, an iterative fluence map merging and sequencing algorithm (vmerge) is used to improve the delivery efficiency. Merging continues as long as the dose quality is maintained above a user-defined threshold. The final plan is selected as the partial arc with the lowest treatment time. The complete algorithm is called pmerge. Results: Partial-arc plans are created using pmerge for a lung, liver and prostate case, with final treatment times of 127, 245 and 147 seconds. Treatment times using full arcs with vmerge are 211, 357 and 178 seconds. Dose quality is maintained across the initial, vmerge, and pmerge plans to within 5% of the mean doses to the critical organs-at-risk and with target coverage above 98%. Additionally, we find that the angular distribution of fluence in the initial plans is predictive of the start and end angles of the optimal partial-arc. Conclusions: The pmerge algorithm is an extension to vmerge that automatically finds the partial-arc plan that minimizes the treatment time. VMAT delivery efficiency can be improved by employing partial-arcs without compromising dose quality. Partial arcs are most applicable to cases with non-centralized targets, where the time savings is greatest.

preprint2011arXiv

A Dynamic Programming Approach to Adaptive Fractionation

We conduct a theoretical study of various solution methods for the adaptive fractionation problem. The two messages of this paper are: (i) dynamic programming (DP) is a useful framework for adaptive radiation therapy, particularly adaptive fractionation, because it allows us to assess how close to optimal different methods are, and (ii) heuristic methods proposed in this paper are near-optimal, and therefore, can be used to evaluate the best possible benefit of using an adaptive fraction size. The essence of adaptive fractionation is to increase the fraction size when the tumor and organ-at-risk (OAR) are far apart (a "favorable" anatomy) and to decrease the fraction size when they are close together. Given that a fixed prescribed dose must be delivered to the tumor over the course of the treatment, such an approach results in a lower cumulative dose to the OAR when compared to that resulting from standard fractionation. We first establish a benchmark by using the DP algorithm to solve the problem exactly. In this case, we characterize the structure of an optimal policy, which provides guidance for our choice of heuristics. We develop two intuitive, numerically near-optimal heuristic policies, which could be used for more complex, high-dimensional problems. Furthermore, one of the heuristics requires only a statistic of the motion probability distribution, making it a reasonable method for use in a realistic setting. Numerically, we find that the amount of decrease in dose to the OAR can vary significantly (5 - 85%) depending on the amount of motion in the anatomy, the number of fractions, and the range of fraction sizes allowed. In general, the decrease in dose to the OAR is more pronounced when: (i) we have a high probability of large tumor-OAR distances, (ii) we use many fractions (as in a hyper-fractionated setting), and (iii) we allow large daily fraction size deviations.

preprint2011arXiv

Including robustness in multi-criteria optimization for intensity-modulated proton therapy

We present a method to include robustness into a multi-criteria optimization (MCO) framework for intensity-modulated proton therapy (IMPT). The approach allows one to simultaneously explore the trade-off between different objectives as well as the trade-off between robustness and nominal plan quality. In MCO, a database of plans each emphasizing different treatment planning objectives, is pre-computed to approximate the Pareto surface. An IMPT treatment plan that strikes the best balance between the different objectives can be selected by navigating on the Pareto surface. In our approach, robustness is integrated into MCO by adding robustified objectives and constraints to the MCO problem. Uncertainties of the robust problem are modeled by pre-calculated dose-influence matrices for a nominal scenario and a number of pre-defined error scenarios. A robustified objective represents the worst objective function value that can be realized for any of the error scenarios. The optimization method is based on a linear projection solver and is capable of handling large problem sizes resulting from a fine dose grid resolution, many scenarios, and a large number of proton pencil beams. A base-of-skull case is used to demonstrate the robust optimization method. It is demonstrated that the robust optimizationmethod reduces the sensitivity of the treatment plan to setup and range errors to a degree that is not achieved by a safety margin approach. A chordoma case is analyzed in more detail to demonstrate the involved tradeoffs between target underdose and brainstem sparing as well as robustness and nominal plan quality. The latter illustrates the advantage of MCO in the context of robust planning. For all cases examined, the robust optimization for each Pareto optimal plan takes less than 5 min on a standard computer, making a computationally friendly interface possible to the planner.

preprint2011arXiv

Multicriteria VMAT optimization

Purpose: To make the planning of volumetric modulated arc therapy (VMAT) faster and to explore the tradeoffs between planning objectives and delivery efficiency. Methods: A convex multicriteria dose optimization problem is solved for an angular grid of 180 equi-spaced beams. This allows the planner to navigate the ideal dose distribution Pareto surface and select a plan of desired target coverage versus organ at risk sparing. The selected plan is then made VMAT deliverable by a fluence map merging and sequencing algorithm, which combines neighboring fluence maps based on a similarity score and then delivers the merged maps together, simplifying delivery. Successive merges are made as long as the dose distribution quality is maintained. The complete algorithm is called VMERGE. Results: VMERGE is applied to three cases: a prostate, a pancreas, and a brain. In each case, the selected Pareto-optimal plan is matched almost exactly with the VMAT merging routine, resulting in a high quality plan delivered with a single arc in less than five minutes on average. VMERGE offers significant improvements over existing VMAT algorithms. The first is the multicriteria planning aspect, which greatly speeds up planning time and allows the user to select the plan which represents the most desirable compromise between target coverage and organ at risk sparing. The second is the user-chosen epsilon-optimality guarantee of the final VMAT plan. Finally, the user can explore the tradeoff between delivery time and plan quality, which is a fundamental aspect of VMAT that cannot be easily investigated with current commercial planning systems.