Improved Approximation Algorithms for Segment Minimization in Intensity Modulated Radiation Therapy

he segment minimization problem consists of finding the smallest set of integer matrices that sum to a given intensity matrix, such that each summand has only one non-zero value, and the non-zeroes in each row are consecutive. This has direct applications in intensity-modulated radiation therapy, an effective form of cancer treatment. We develop three approximation algorithms for matrices with arbitrarily many rows. Our first two algorithms improve the approximation factor from the previous best of $1+\log_2 h$ to (roughly) $3/2 \cdot (1+\log_3 h)$ and $11/6\cdot(1+\log_4{h})$, respectively, where $h$ is the largest entry in the intensity matrix. We illustrate the limitations of the specific approach used to obtain these two algorithms by proving a lower bound of $\frac{(2b-2)}{b}\cdot\log_b{h} + \frac{1}{b}$ on the approximation guarantee. Our third algorithm improves the approximation factor from $2 \cdot (\log D+1)$ to $24/13 \cdot (\log D+1)$, where $D$ is (roughly) the largest difference between consecutive elements of a row of the intensity matrix. Finally, experimentation with these algorithms shows that they perform well with respect to the optimum and outperform other approximation algorithms on 77% of the 122 test cases we consider, which include both real world and synthetic data.Comment: 18 page

A Hierachical Evolutionary Algorithm for Multiobjective Optimization in IMRT

Purpose: Current inverse planning methods for IMRT are limited because they are not designed to explore the trade-offs between the competing objectives between the tumor and normal tissues. Our goal was to develop an efficient multiobjective optimization algorithm that was flexible enough to handle any form of objective function and that resulted in a set of Pareto optimal plans. Methods: We developed a hierarchical evolutionary multiobjective algorithm designed to quickly generate a diverse Pareto optimal set of IMRT plans that meet all clinical constraints and reflect the trade-offs in the plans. The top level of the hierarchical algorithm is a multiobjective evolutionary algorithm (MOEA). The genes of the individuals generated in the MOEA are the parameters that define the penalty function minimized during an accelerated deterministic IMRT optimization that represents the bottom level of the hierarchy. The MOEA incorporates clinical criteria to restrict the search space through protocol objectives and then uses Pareto optimality among the fitness objectives to select individuals. Results: Acceleration techniques implemented on both levels of the hierarchical algorithm resulted in short, practical runtimes for optimizations. The MOEA improvements were evaluated for example prostate cases with one target and two OARs. The modified MOEA dominated 11.3% of plans using a standard genetic algorithm package. By implementing domination advantage and protocol objectives, small diverse populations of clinically acceptable plans that were only dominated 0.2% by the Pareto front could be generated in a fraction of an hour. Conclusions: Our MOEA produces a diverse Pareto optimal set of plans that meet all dosimetric protocol criteria in a feasible amount of time. It optimizes not only beamlet intensities but also objective function parameters on a patient-specific basis

Controlling Beam Complexity in Intensity Modulated Radiation Therapy.

External beam intensity modulated radiation therapy (IMRT) is a technique in which the spatial intensity of radiation from each beam direction can be modulated to provide superior conformality of dose to a tumor volume while sparing important normal tissues. A fundamental and potentially limiting feature of IMRT is the highly complex fields that can be created through inverse plan optimization. Highly modulated treatments are a large departure from conventional radiotherapy methods, are difficult to deliver accurately and efficiently, and can result in an undesirable increase in leakage dose being delivered to the patient. Longer deliveries may also increase the chance for patient motion during treatment and could potentially reduce the probability of controlling some tumors. The large intensity fluctuations observed in IMRT beams are often a result of the degeneracy of the optimization problem, and the types of optimization method and cost function used. This work demonstrates that beam complexity is a result of these two issues, and is dependent on the placement of dose evaluation points in the target and normal tissues. This research shows that (i) optimizing surfaces instead of discrete beamlet intensities to represent the beam can reduce the degrees of freedom in IMRT and results in much smoother beams at the expense of a slight increase in normal tissues, (ii) maximum beamlet intensity restrictions are useful for improved delivery efficiency, but may restrict the optimizer at low limits, and (iii) modulation penalties can be incorporated into the cost function to promote plan smoothness without sacrificing plan quality. Penalizing the overall plan modulation is an effective way to reduce modulation, but it falsely penalizes the desirable beam modulation as well as the undesirable modulation. To address this problem, diffusion principles are used to develop a spatially adaptive smoothing method that only penalizes the unnecessary beam modulation and can be used without degrading plan quality. This method is customizable to a variety of treatment scenarios. The clinical impact of reducing beam complexity is significant, as it can result in an improvement in delivery accuracy and efficiency, quicker optimization times, and increased robustness to point sampling and geometric uncertainty.Ph.D.Nuclear Engineering & Radiological SciencesUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/57648/2/mcoselmo_1.pd

The Possibilities and Dosimetric Limitations of MLC-Based Intensity-Modulated Radiotherapy Delivery and Optimization Techniques

The use of intensity-modulated radiotherapy (IMRT) has increased extensively in the modern radiotherapy (RT) treatments over the past two decades. Radiation dose distributions can be delivered with higher conformality with IMRT when compared to the conventional 3D-conformal radiotherapy (3D-CRT). Higher conformality and target coverage increases the probability of tumour control and decreases the normal tissue complications. The primary goal of this work is to improve and evaluate the accuracy, efficiency and delivery techniques of RT treatments by using IMRT. This study evaluated the dosimetric limitations and possibilities of IMRT in small (treatments of head-and-neck, prostate and lung cancer) and large volumes (primitive neuroectodermal tumours). The dose coverage of target volumes and the sparing of critical organs were increased with IMRT when compared to 3D-CRT. The developed split field IMRT technique was found to be safe and accurate method in craniospinal irradiations. By using IMRT in simultaneous integrated boosting of biologically defined target volumes of localized prostate cancer high doses were achievable with only small increase in the treatment complexity. Biological plan optimization increased the probability of uncomplicated control on average by 28% when compared to standard IMRT delivery. Unfortunately IMRT carries also some drawbacks. In IMRT the beam modulation is realized by splitting a large radiation field to small apertures. The smaller the beam apertures are the larger the rebuild-up and rebuild-down effects are at the tissue interfaces. The limitations to use IMRT with small apertures in the treatments of small lung tumours were investigated with dosimetric film measurements. The results confirmed that the peripheral doses of the small lung tumours were decreased as the effective field size was decreased. The studied calculation algorithms were not able to model the dose deficiency of the tumours accurately. The use of small sliding window apertures of 2 mm and 4 mm decreased the tumour peripheral dose by 6% when compared to 3D-CRT treatment plan. A direct aperture based optimization (DABO) technique was examined as a solution to decrease the treatment complexity. The DABO IMRT technique was able to achieve treatment plans equivalent with the conventional IMRT fluence based optimization techniques in the concave head-and-neck target volumes. With DABO the effective field sizes were increased and the number of MUs was reduced with a factor of two. The optimality of a treatment plan and the therapeutic ratio can be further enhanced by using dose painting based on regional radiosensitivities imaged with functional imaging methods.Siirretty Doriast

Towards on-line plan adaptation of unified intensity-modulated arc therapy using a fast-direct aperture optimization algorithm

External beam radiotherapy (EBRT) plays a vital role in the treatment of cancer, with close to half of all cancer patients receiving EBRT at some point over their course of treatment. Although EBRT is a well-established form of treatment, there are a number of ways in which EBRT could still be improved in terms of quality and efficiency for treatment planning and radiation dose delivery. This thesis reports a series of improvements made to EBRT. First, we developed and evaluated a new treatment planning technique called unified intensity-modulated arc therapy (UIMAT) which combines the optimization and delivery of rotational volumetric modulated arc therapy (VMAT) and fixed-gantry intensity-modulated radiation therapy (IMRT). When retrospectively compared to clinical treatment plans using VMAT or IMRT alone, UIMAT plans reduced the dose to nearby critical structures by as much as 23% without compromising tumour volume coverage. The UIMAT plans were also more efficient to deliver. The reduction in normal tissue dose could help lower the probability of treatment-related toxicities, or alternatively could be used to improve tumour control probability, via dose escalation, while maintaining current dose limits for organs at risk. Second, we developed a new fast inverse direct aperture optimization (FIDAO) algorithm for IMRT, VMAT, and UIMAT treatment planning. FIDAO introduces modifications to the direct aperture optimization (DAO) process that help improve its computational efficiency. As demonstrated in several test cases, these modifications do not significantly impact the plan quality but reduced the DAO time by as much as 200-fold. If implemented with graphical processing units (GPUs), this project may allow for applications such as on-line treatment adaptation. Third, we investigated a method of acquiring tissue density information from cone-beam computed tomography (CBCT) datasets for on-line dose calculations, plan assessment, and potentially plan adaptation using FIDAO. This calibration technique accounts for patient-specific scattering conditions, demonstrated high dosimetric accuracy, and can be easily automated for on-line plan assessment. Collectively, these three projects will help reduce the normal tissue doses from EBRT, improve the planning and delivery efficiency, and pave the way for application like on-line plan assessment and adaptive radiotherapy in response to anatomical changes

Optimization Problems in Radiation Therapy Treatment Planning.

Radiation therapy is one of the most common methods used to treat many types of cancer. External beam radiation therapy and the models associated with developing a treatment plan for a patient are studied. External beams of radiation are used to deliver a highly complex so-called dose distribution to a patient that is designed to kill the cancer cells while sparing healthy organs and normal tissue. Treatment planning models and optimization are used to determine the delivery machine instructions necessary to produce a desirable dose distribution. These instructions make up a treatment plan. This thesis studies four problems in radiation therapy treatment plan optimization. First, treatment planners generate a plan with a number of competing treatment plan criteria. The relationship between criteria is not known a priori. A methodology is developed for physicians and treatment planners to efficiently navigate a clinically relevant region of the Pareto frontier generated by trading off these different criteria in an informed way. Second, the machine instructions for intensity modulated radiation therapy, a common treatment modality, consist of the locations of the external beams and the non-uniform intensity profiles delivered from each of these locations. These decisions are traditionally made with separate, sequential models. These decisions are integrated into a single model and propose a heuristic solution methodology. Third, volumetric modulated arc therapy (VMAT), a treatment modality where the beam travels in a coplanar arc around the patient while continuously delivering radiation, is a popular topic among optimizers studying treatment planning due to the difficult nature of the problem and the lack of a universally accepted treatment planning method. While current solution methodologies assume a predetermined coplanar path around the patient, that assumption is relaxed and the generation of a non-coplanar path is integrated into a VMAT planning algorithm. Fourth, not all patient information is available when developing a treatment plan pre-treatment. Some information, like a patient's sensitivity to radiation, can be realized during treatment through physiological tests. Methodologies of pre-treatment planning considering adaptation to new information are studied.PhDIndustrial and Operations EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/113366/1/troylong_1.pd

Approximated multileaf collimator field segmentation

In intensity-modulated radiation therapy the aim is to realize given intensity distributions as a superposition of differently shaped fields. Multileaf collimators are used for field shaping. This segmentation task leads to discrete optimization problems, that are considered in this dissertation. A variety of algorithms for exact and approximated segmentation, for different objective functions and various technical as well as dosimetric constraints are developed