Clinical knowledge-based inverse treatment planning

Application and progress of artificial intelligence in radiation therapy dose prediction

Radiation therapy (RT) nowadays is a main treatment modality of cancer. To ensure the therapeutic efficacy of patients, accurate dose distribution is often required, which is a time-consuming and labor-intensive process. In addition, due to the differences in knowledge and experience among participants and diverse institutions, the predicted dose are often inconsistent. In last several decades, artificial intelligence (AI) has been applied in various aspects of RT, several products have been implemented in clinical practice and confirmed superiority. In this paper, we will review the research of AI in dose prediction, focusing on the progress in deep learning (DL).

Utilizing pre-determined beam orientation information in dose prediction by 3D fully-connected network for intensity modulated radiotherapy

Background

Although the effect of pre-determined beam orientation on dose distribution of intensity modulated radiotherapy (IMRT) has been well-documented, its impacts on dose prediction are less investigated. In this study, the direction map of beam orientation was incorporated into our proposed deep-learning network and utilized in dose prediction of IMRT plans consisting of multiple static fields.

Methods

The direction map was used to characterize the radiation path through region of interest along a beam orientation. Besides, the distance map was used to characterize the spatial distribution between organs at risk (OARs) and planning target volume (PTV). The input of prediction model consisted of CT image, mask image (for PTV and OARs), distance map, and direction map. The output of prediction model was the estimated dose distribution in three dimensions. A 3D fully-connected network composed of a down-sampling encoder and an up-sampling pyramid decoder was trained based on the calculated 3D dose distributions obtained from a treatment planning system. The voxel-level mean absolute error (MAE), dosimetric metrics, and dose-volume histogram were employed to assess the quality of the estimated dose distribution. Performance of the prediction model was evaluated in two aspects. First, the effectiveness of the new features, direction map, distance maps, and pyramid decoder on prediction accuracy of model were assessed. Second, the proposed model was compared with the other three published prediction models, 3D UNet, ResNet-anti-ResNet, U-ResNet-D for inter-model evaluation.

Results

The improvement of prediction accuracy was 0.38 with the input of direction map and 0.43 with the input of distance map. Our proposed model achieved the least MAE (3.97±1.42) compared with the other three models: (5.37±1.51) for ResNet-anti-ResNet, (4.45±1.52) for U-ResNet-D, and (4.53±1.72) for Unet-3D.

Conclusions

The preliminary result demonstrated that the prediction accuracy of the proposed model was higher than those of the other three state-of-the-art prediction models. The introduction of direction maps, distance map, and pyramid decoder can effectively improve the performance of the current deep-learning network-based prediction models.

Personalized automation of treatment planning in head-neck cancer: A step forward for quality in radiation therapy?

PURPOSE

To perform a comprehensive dosimetric and clinical evaluation of the new Pinnacle Personalized automated planning system for complex head-and-neck treatments.

METHODS

Fifteen consecutive head-neck patients were enrolled. Radiotherapy was prescribed using VMAT with simultaneous integrated boost strategy. Personalized planning integrates the Feasibility engine able to supply an "a priori" DVH prediction of the achievability of planning goals. Comparison between clinically accepted manually-generated (MP) and automated (AP) plans was performed using dose-volume histograms and a blinded clinical evaluation by two radiation oncologists. Planning time between MP and AP was compared. Dose accuracy was validated using the PTW Octavius-4D phantom together with the 1500 2D-array.

RESULTS

For similar targets coverage, AP plans reported less irradiation of healthy tissue, with significant dose reduction for spinal cord, brainstem and parotids. On average, the mean dose to parotids and maximal doses to spinal cord and brainstem were reduced by 13-15% (p < 0.001), 9% (p < 0.001) and 16% (p < 0.001), respectively. The integral dose was reduced by 16% (p < 0.001). The dose conformity for the three PTVs was significantly higher with AP plans (p < 0.001). The two oncologists chose AP plans in more than 80% of cases. Overall planning times were reduced to <30 min for automated optimization. All AP plans passed the 3%/2 mm γ-analysis by more than 95%.

CONCLUSION

Complex head-neck plans created using Personalized automated engine provided an overall increase of plan quality, in terms of dose conformity and sparing of normal tissues. The Feasibility module allowed OARs dose sparing well beyond the clinical objectives.

Automatic Planning for Nasopharyngeal Carcinoma Based on Progressive Optimization in RayStation Treatment Planning System

OBJECTIVE

To evaluate and quantify the planning performance of automatic planning (AP) with manual planning (MP) for nasopharyngeal carcinoma in the RayStation treatment planning system (TPS).

METHODS

A progressive and effective design method for AP of nasopharyngeal carcinoma was realized through automated scripts in this study. A total of 30 patients with nasopharyngeal carcinoma with initial treatment was enrolled. The target coverage, conformity index (CI), homogeneity index (HI), organs at risk sparing, and the efficiency of design and execution were compared between automatic and manual volumetric modulated arc therapy (VMAT) plans.

RESULTS

The results of the 2 design methods met the clinical dose requirement. The differences in D95 between the 2 groups in PTV1 and PTV2 showed statistical significance, and the MPs are higher than APs, but the difference in absolute dose was only 0.21% and 0.16%. The results showed that the conformity index of planning target volumes (PTV1, PTV2, PTVnd and PGTVnx+rpn [PGTVnx and PGTVrpn]), homogeneity index of PGTVnx+rpn, and HI of PTVnd in APs are better than that in MPs. For organs at risk, the APs are lower than the MPs, and the difference was statistically significant (P < .05). The manual operation time in APs was 83.21% less than that in MPs, and the computer processing time was 34.22% more.

CONCLUSION

IronPython language designed by RayStation TPS has clinical application value in the design of automatic radiotherapy plan for nasopharyngeal carcinoma. The dose distribution of tumor target and organs at risk in the APs was similar or better than those in the MPs. The time of manual operation in the plan design showed a sharp reduction, thus significantly improving the work efficiency in clinical application.

Impact of beam configuration on VMAT plan quality for Pinnacle3Auto-Planning for head and neck cases

Background

The purpose of this study was to compare automatically generated VMAT plans to find the superior beam configurations for Pinnacle³ Auto-Planning and share “best practices”.

Methods

VMAT plans for 20 patients with head and neck cancer were generated using Pinnacle³ Auto-Planning Module (Pinnacle³ Version 9.10) with different beam setup parameters. VMAT plans for single (V1) or double arc (V2) and partial or full gantry rotation were optimized. Beam configurations with different collimator positions were defined. Target coverage and sparing of organs at risk were evaluated based on scoring of an evaluation parameter set. Furthermore, dosimetric evaluation was performed based on the composite objective value (COV) and a new cross comparison method was applied using the COVs.

Results

The evaluation showed a superior plan quality for double arcs compared to one single arc or two single arcs for all cases. Plan quality was superior if a full gantry rotation was allowed during optimization for unilateral target volumes. A double arc technique with collimator setting of 15° was superior to a double arc with collimator 60° and a two single arcs with collimator setting of 15° and 345°.

Conclusion

The evaluation showed that double and full arcs are superior to single and partial arcs in terms of organs at risk sparing even for unilateral target volumes. The collimator position was found as an additional setup parameter, which can further improve the target coverage and sparing of organs at risk.

Automation in intensity modulated radiotherapy treatment planning-a review of recent innovations

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.

Isodose feature-preserving voxelization (IFPV) for radiation therapy treatment planning

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.

Machine Learning in Radiation Oncology: Opportunities, Requirements, and Needs

Machine learning (ML) has the potential to revolutionize the field of radiation oncology, but there is much work to be done. In this article, we approach the radiotherapy process from a workflow perspective, identifying specific areas where a data-centric approach using ML could improve the quality and efficiency of patient care. We highlight areas where ML has already been used, and identify areas where we should invest additional resources. We believe that this article can serve as a guide for both clinicians and researchers to start discussing issues that must be addressed in a timely manner.

Automation of radiation treatment planning : Evaluation of head and neck cancer patient plans created by the Pinnacle3 scripting and Auto-Planning functions

BACKGROUND

Intensity-modulated radiotherapy (IMRT) techniques are now standard practice. IMRT or volumetric-modulated arc therapy (VMAT) allow treatment of the tumor while simultaneously sparing organs at risk. Nevertheless, treatment plan quality still depends on the physicist's individual skills, experiences, and personal preferences. It would therefore be advantageous to automate the planning process. This possibility is offered by the Pinnacle³ treatment planning system (Philips Healthcare, Hamburg, Germany) via its scripting language or Auto-Planning (AP) module.

MATERIALS AND METHODS

AP module results were compared to in-house scripts and manually optimized treatment plans for standard head and neck cancer plans. Multiple treatment parameters were scored to judge plan quality (100 points = optimum plan). Patients were initially planned manually by different physicists and re-planned using scripts or AP.

RESULTS AND DISCUSSION

Script-based head and neck plans achieved a mean of 67.0 points and were, on average, superior to manually created (59.1 points) and AP plans (62.3 points). Moreover, they are characterized by reproducibility and lower standard deviation of treatment parameters. Even less experienced staff are able to create at least a good starting point for further optimization in a short time. However, for particular plans, experienced planners perform even better than scripts or AP. Experienced-user input is needed when setting up scripts or AP templates for the first time. Moreover, some minor drawbacks exist, such as the increase of monitor units (+35.5% for scripted plans).

CONCLUSION

On average, automatically created plans are superior to manually created treatment plans. For particular plans, experienced physicists were able to perform better than scripts or AP; thus, the benefit is greatest when time is short or staff inexperienced.

Personalized treatment planning with a model of radiation therapy outcomes for use in multiobjective optimization of IMRT plans for prostate cancer

PURPOSE

To build a new treatment planning approach that extends beyond radiation transport and IMRT optimization by modeling the radiation therapy process and prognostic indicators for more outcome-focused decision making.

METHODS

An in-house treatment planning system was modified to include multiobjective inverse planning, a probabilistic outcome model, and a multi-attribute decision aid. A genetic algorithm generated a set of plans embodying trade-offs between the separate objectives. An influence diagram network modeled the radiation therapy process of prostate cancer using expert opinion, results of clinical trials, and published research. A Markov model calculated a quality adjusted life expectancy (QALE), which was the endpoint for ranking plans.

RESULTS

The Multiobjective Evolutionary Algorithm (MOEA) was designed to produce an approximation of the Pareto Front representing optimal tradeoffs for IMRT plans. Prognostic information from the dosimetrics of the plans, and from patient-specific clinical variables were combined by the influence diagram. QALEs were calculated for each plan for each set of patient characteristics. Sensitivity analyses were conducted to explore changes in outcomes for variations in patient characteristics and dosimetric variables. The model calculated life expectancies that were in agreement with an independent clinical study.

CONCLUSIONS

The radiation therapy model proposed has integrated a number of different physical, biological and clinical models into a more comprehensive model. It illustrates a number of the critical aspects of treatment planning that can be improved and represents a more detailed description of the therapy process. A Markov model was implemented to provide a stronger connection between dosimetric variables and clinical outcomes and could provide a practical, quantitative method for making difficult clinical decisions.

Simultaneous beam sampling and aperture shape optimization for SPORT

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 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.

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.

Automated population-based planning for whole brain radiation therapy

Treatment planning for whole‐brain radiation treatment is technically a simple process, but in practice it takes valuable clinical time of repetitive and tedious tasks. This report presents a method that automatically segments the relevant target and normal tissues, and creates a treatment plan in only a few minutes after patient simulation. Segmentation of target and critical structures is performed automatically through morphological operations on the soft tissue and was validated by comparing with manual clinical segmentation using the Dice coefficient and Hausdorff distance. The treatment plan is generated by searching a database of previous cases for patients with similar anatomy. In this search, each database case is ranked in terms of similarity using a customized metric designed for sensitivity by including only geometrical changes that affect the dose distribution. The database case with the best match is automatically modified to replace relevant patient info and isocenter position while maintaining original beam and MLC settings. Fifteen patients with marginally acceptable treatment plans were used to validate the method. In each of these cases the anatomy was accurately segmented, but the beams and MLC settings led to a suboptimal treatment plan by either underdosing the brain or excessively irradiating critical normal tissues. For each case, the anatomy was automatically segmented with the proposed method, and the automated and manual segmentations were then compared. The mean Dice coefficient was 0.97, with a standard deviation of 0.008 for the brain, 0.85±0.009 for the eyes, and 0.67±0.11 for the lens. The mean Euclidian distance was 0.13±0.13 mm for the brain, 0.27±0.31 for the eye, and 2.34±7.23 for the lens. Each case was then subsequently matched against a database of 70 validated treatment plans and the best matching plan (termed autoplanned), was compared retrospectively with the clinical plans in terms of brain coverage and maximum doses to critical structures. Maximum doses were reduced by a maximum of 8.37 Gy for the left eye (mean 2.08), 11.67 for the right eye (1.90) and, respectively, 25.44 (5.59) for the left lens and 24.40 (4.85) for the right lens. Time to generate the autoplan, including the segmentation, was 3−4 min. Automated database‐ based matching is an alternative to classical treatment planning that improves quality while providing a cost‐effective solution to planning through modifying previous validated plans to match a current patient's anatomy.

PACS number: 87.55.D, 87.55.tg, 87.57.nm

Dependence of achievable plan quality on treatment technique and planning goal refinement: a head-and-neck intensity modulated radiation therapy application

PURPOSE

To develop a practical workflow for retrospectively analyzing target and normal tissue dose-volume endpoints for various intensity modulated radiation therapy (IMRT) delivery techniques; to develop technique-specific planning goals to improve plan consistency and quality when feasible.

METHODS AND MATERIALS

A total of 165 consecutive head-and-neck patients from our patient registry were selected and retrospectively analyzed. All IMRT plans were generated using the same dose-volume guidelines for TomoTherapy (Tomo, Accuray), TrueBeam (TB, Varian) using fixed-field IMRT (TB_IMRT) or RAPIDARC (TB_RAPIDARC), or Siemens Oncor (Siemens_IMRT, Siemens). A MATLAB-based dose-volume extraction and analysis tool was developed to export dosimetric endpoints for each patient. With a fair stratification of patient cohort, the variation of achieved dosimetric endpoints was analyzed among different treatment techniques. Upon identification of statistically significant variations, technique-specific planning goals were derived from dynamically accumulated institutional data.

RESULTS

Retrospective analysis showed that although all techniques yielded comparable target coverage, the doses to the critical structures differed. The maximum cord doses were 34.1 ± 2.6, 42.7 ± 2.1, 43.3 ± 2.0, and 45.1 ± 1.6 Gy for Tomo, TB_IMRT, TB_RAPIDARC, and Siemens_IMRT plans, respectively. Analyses of variance showed significant differences for the maximum cord doses but no significant differences for other selected structures among the investigated IMRT delivery techniques. Subsequently, a refined technique-specific dose-volume guideline for maximum cord dose was derived at a confidence level of 95%. The dosimetric plans that failed the refined technique-specific planning goals were reoptimized according to the refined constraints. We observed better cord sparing with minimal variations for the target coverage and other organ at risk sparing for the Tomo cases, and higher parotid doses for C-arm linear accelerator-based IMRT and RAPIDARC plans.

CONCLUSION

Patient registry-based processes allowed easy and systematic dosimetric assessment of treatment plan quality and consistency. Our analysis revealed the dependence of certain dosimetric endpoints on the treatment techniques. Technique-specific refinement of planning goals may lead to improvement in plan consistency and plan quality.

Prior-knowledge treatment planning for volumetric arc therapy using feature-based database mining

Treatment planning for volumetric arc therapy (VMAT) is a lengthy process that requires many rounds of optimizations to obtain the best treatment settings and optimization constraints for a given patient's geometry. We propose a feature‐selection search engine that explores previously treated cases of similar anatomy, returning the optimal plan configurations and attainable DVH constraints. Using an institutional database of 83 previously treated cases of prostate carcinoma treated with volumetric‐modulated arc therapy, the search procedure first finds the optimal isocenter position with an optimization procedure, then ranks the anatomical similarity as the mean distance between targets. For the best matching plan, the planning information is reformatted to the DICOM format and imported into the treatment planning system to suggest isocenter, arc directions, MLC patterns, and optimization constraints that can be used as starting points in the optimization process. The approach was tested to create prospective treatment plans based on anatomical features that match previously treated cases from the institution database. By starting from a near‐optimal solution and using previous optimization constraints, the best matching test only required simple optimization steps to further decrease target inhomogeneity, ultimately reducing time spend by the therapist in planning arcs' directions and lengths.

PACS number: 87.55.D‐, 87.55.de

The use of a multiobjective evolutionary algorithm to increase flexibility in the search for better IMRT plans

PURPOSE

To evaluate how a more flexible and thorough multiobjective search of feasible IMRT plans affects performance in IMRT optimization.

METHODS

A multiobjective evolutionary algorithm (MOEA) was used as a tool to investigate how expanding the search space to include a wider range of penalty functions affects the quality of the set of IMRT plans produced. The MOEA uses a population of IMRT plans to generate new IMRT plans through deterministic minimization of recombined penalty functions that are weighted sums of multiple, tissue-specific objective functions. The quality of the generated plans are judged by an independent set of nonconvex, clinically relevant decision criteria, and all dominated plans are eliminated. As this process repeats itself, better plans are produced so that the population of IMRT plans will approach the Pareto front. Three different approaches were used to explore the effects of expanding the search space. First, the evolutionary algorithm used genetic optimization principles to search by simultaneously optimizing both the weights and tissue-specific dose parameters in penalty functions. Second, penalty function parameters were individually optimized for each voxel in all organs at risk (OARs) in the MOEA. Finally, a heuristic voxel-specific improvement (VSI) algorithm that can be used on any IMRT plan was developed that incrementally improves voxel-specific penalty function parameters for all structures (OARs and targets). Different approaches were compared using the concept of domination comparison applied to the sets of plans obtained by multiobjective optimization.

RESULTS

MOEA optimizations that simultaneously searched both importance weights and dose parameters generated sets of IMRT plans that were superior to sets of plans produced when either type of parameter was fixed for four example prostate plans. The amount of improvement increased with greater overlap between OARs and targets. Allowing the MOEA to search for voxel-specific penalty functions improved results for simple cases with three structures but did not improve results for a more complex case with seven structures. For this modification, the amount of improvement increased with less overlap between OARs and targets. The voxel-specific improvement algorithm improved results for all cases, and its clinical relevance was demonstrated in a complex prostate and a very complex head and neck case.

CONCLUSIONS

Using an evolutionary algorithm as a tool, it was found that allowing more flexibility in the search space enhanced performance. The two strategies of (a) varying the weights and reference doses in the objective function and (b) removing the constraint of equal penalties for all voxels in a structure both generated sets of plans that dominated sets of plans considered to be "Pareto optimal" within the conventional, more limited search space. When considering voxel-specific objectives, the very large search space can lead to convergence problems in the MOEA for complex cases, but this is not an issue for the VSI algorithm.

Toward truly optimal IMRT dose distribution: inverse planning with voxel-specific penalty

PURPOSE

To establish an inverse planning framework with adjustable voxel penalty for more conformal IMRT dose distribution as well as improved interactive controllability over the regional dose distribution of the resultant plan.

MATERIALS AND METHOD

In the proposed coarse-to-fine planning scheme, a conventional inverse planning with organ specific parameters is first performed. The voxel penalty scheme is then "switched on" by allowing the prescription dose to change on an individual voxel scale according to the deviation of the actual voxel dose from the ideally desired dose. The rationale here is intuitive: when the dose at a voxel does not meet its ideal dose, it simply implies that this voxel is not competitive enough when compared with the ones that have met their planning goal. In this case, increasing the penalty of the voxel by varying the prescription can boost its competitiveness and thus improve its dose. After the prescription adjustment, the plan is re-optimized. The dose adjustment/re-optimization procedure is repeated until the resultant dose distribution cannot be improved anymore. The prescription adjustment on a finer scale can be accomplished either automatically or manually. In the latter case, the regions/voxels where a dose improvement is needed are selected visually, unlike in the automatic case where the selection is done purely based on the difference of the actual dose at a given voxel and its ideal prescription. The performance of the proposed method is evaluated using a head and neck and a prostate case.

RESULTS

An inverse planning framework with the voxel-specific penalty is established. By adjusting voxel prescriptions iteratively to boost the region where large mismatch between the actual calculated and desired doses occurs, substantial improvements can be achieved in the final dose distribution. The proposed method is applied to a head and neck case and a prostate case. For the former case, a significant reduction in the maximum dose to the brainstem is achieved while the PTV dose coverage is greatly improved. The doses to other organs at risk are also reduced, ranging from 10% to 30%. For the prostate case, the use of the voxel penalty scheme also results in vast improvements to the final dose distribution. The PTV experiences improved dose uniformity and the mean dose to the rectum and bladder is reduced by as much as 15%.

CONCLUSION

Introduction of the spatially non-uniform and adjustable prescription provides room for further improvements of currently achievable dose distributions and equips the planner with an effective tool to modify IMRT dose distributions interactively. The technique is easily implementable in any existing inverse planning platform, which should facilitate clinical IMRT planning process and, in future, off-line/on-line adaptive IMRT.

Search for IMRT inverse plans with piecewise constant fluence maps using compressed sensing techniques

An intensity-modulated radiation therapy (IMRT) field is composed of a series of segmented beams. It is practically important to reduce the number of segments while maintaining the conformality of the final dose distribution. In this article, the authors quantify the complexity of an IMRT fluence map by introducing the concept of sparsity of fluence maps and formulate the inverse planning problem into a framework of compressing sensing. In this approach, the treatment planning is modeled as a multiobjective optimization problem, with one objective on the dose performance and the other on the sparsity of the resultant fluence maps. A Pareto frontier is calculated, and the achieved dose distributions associated with the Pareto efficient points are evaluated using clinical acceptance criteria. The clinically acceptable dose distribution with the smallest number of segments is chosen as the final solution. The method is demonstrated in the application of fixed-gantry IMRT on a prostate patient. The result shows that the total number of segments is greatly reduced while a satisfactory dose distribution is still achieved. With the focus on the sparsity of the optimal solution, the proposed method is distinct from the existing beamlet- or segment-based optimization algorithms.

Comparison of intensity modulated x-ray therapy and intensity modulated proton therapy for selective subvolume boosting: a phantom study

Selective subvolume boosting can theoretically improve tumour control probability while maintaining normal tissue complication probabilities similar to those of uniform dose distributions. In this work the abilities of intensity-modulated x-ray therapy (IMXT) and intensity-modulated proton therapy (IMPT) to deliver boosts to multiple subvolumes of varying size and proximities are compared in a thorough phantom study. IMXT plans were created using the step-and-shoot (IMXT-SAS) and helical tomotherapy (IMXT-HT) methods. IMPT plans were created with the spot scanning (IMPT-SS) and distal gradient tracking (IMPT-DGT) methods. IMPT-DGT is a generalization of the distal edge tracking method designed to reduce the number of proton beam spots required to deliver non-uniform dose distributions relative to IMPT-SS. The IMPT methods were delivered over both 180 degrees and 360 degrees arcs. The IMXT-SAS and IMPT-SS methods optimally satisfied the non-uniform dose prescriptions the least and the most, respectively. The IMPT delivery methods reduced the normal tissue integral dose by a factor of about 2 relative to the IMXT delivery methods, regardless of the delivery arc. The IMPT-DGT method reduced the number of proton beam spots by a factor of about 3 relative to the IMPT-SS method.

Formulating adaptive radiation therapy (ART) treatment planning into a closed-loop control framework

While ART has been studied for years, the specific quantitative implementation details have not. In order for this new scheme of radiation therapy (RT) to reach its potential, an effective ART treatment planning strategy capable of taking into account the dose delivery history and the patient's on-treatment geometric model must be in place. This paper performs a theoretical study of dynamic closed-loop control algorithms for ART and compares their utility with data from phantom and clinical cases. We developed two classes of algorithms: those Adapting to Changing Geometry and those Adapting to Geometry and Delivered Dose. The former class takes into account organ deformations found just before treatment. The latter class optimizes the dose distribution accumulated over the entire course of treatment by adapting at each fraction, not only to the information just before treatment about organ deformations but also to the dose delivery history. We showcase two algorithms in the class of those Adapting to Geometry and Delivered Dose. A comparison of the approaches indicates that certain closed-loop ART algorithms may significantly improve the current practice. We anticipate that improvements in imaging, dose verification and reporting will further increase the importance of adaptive algorithms.

IMRT treatment planning based on prioritizing prescription goals

Determining the 'best' optimization parameters in IMRT planning is typically a time-consuming trial-and-error process with no unambiguous termination point. Recently we and others proposed a goal-programming approach which better captures the desired prioritization of dosimetric goals. Here, individual prescription goals are addressed stepwise in their order of priority. In the first step, only the highest order goals are considered (target coverage and dose-limiting normal structures). In subsequent steps, the achievements of the previous steps are turned into hard constraints and lower priority goals are optimized, in turn, subject to higher priority constraints. So-called 'slip' factors were introduced to allow for slight, clinically acceptable violations of the constraints. Focusing on head and neck cases, we present several examples for this planning technique. The main advantages of the new optimization method are (i) its ability to generate plans that meet the clinical goals, as well as possible, without tuning any weighting factors or dose-volume constraints, and (ii) the ability to conveniently include more terms such as fluence map smoothness. Lower level goals can be optimized to the achievable limit without compromising higher order goals. The prioritized prescription-goal planning method allows for a more intuitive and human-time-efficient way of dealing with conflicting goals compared to the conventional trial-and-error method of varying weighting factors and dose-volume constraints.

Overview of image-guided radiation therapy

Radiation therapy has gone through a series of revolutions in the last few decades and it is now possible to produce highly conformal radiation dose distribution by using techniques such as intensity-modulated radiation therapy (IMRT). The improved dose conformity and steep dose gradients have necessitated enhanced patient localization and beam targeting techniques for radiotherapy treatments. Components affecting the reproducibility of target position during and between subsequent fractions of radiation therapy include the displacement of internal organs between fractions and internal organ motion within a fraction. Image-guided radiation therapy (IGRT) uses advanced imaging technology to better define the tumor target and is the key to reducing and ultimately eliminating the uncertainties. The purpose of this article is to summarize recent advancements in IGRT and discussed various practical issues related to the implementation of the new imaging techniques available to radiation oncology community. We introduce various new IGRT concepts and approaches, and hope to provide the reader with a comprehensive understanding of the emerging clinical IGRT technologies. Some important research topics will also be addressed.

Dose–volume based ranking of incident beam direction and its utility in facilitating IMRT beam placement

PURPOSE

Beam orientation optimization in intensity-modulated radiation therapy (IMRT) is computationally intensive, and various single beam ranking techniques have been proposed to reduce the search space. Up to this point, none of the existing ranking techniques considers the clinically important dose-volume effects of the involved structures, which may lead to clinically irrelevant angular ranking. The purpose of this work is to develop a clinically sensible angular ranking model with incorporation of dose-volume effects and to show its utility for IMRT beam placement.

METHODS AND MATERIALS

The general consideration in constructing this angular ranking function is that a beamlet/beam is preferable if it can deliver a higher dose to the target without exceeding the tolerance of the sensitive structures located on the path of the beamlet/beam. In the previously proposed dose-based approach, the beamlets are treated independently and, to compute the maximally deliverable dose to the target volume, the intensity of each beamlet is pushed to its maximum intensity without considering the values of other beamlets. When volumetric structures are involved, the complication arises from the fact that there are numerous dose distributions corresponding to the same dose-volume tolerance. In this situation, the beamlets are not independent and an optimization algorithm is required to find the intensity profile that delivers the maximum target dose while satisfying the volumetric constraints. In this study, the behavior of a volumetric organ was modeled by using the equivalent uniform dose (EUD). A constrained sequential quadratic programming algorithm (CFSQP) was used to find the beam profile that delivers the maximum dose to the target volume without violating the EUD constraint or constraints. To assess the utility of the proposed technique, we planned a head-and-neck and abdominal case with and without the guidance of the angular ranking information. The qualities of the two types of IMRT plans were compared quantitatively.

RESULTS

An effective angular ranking model with consideration of volumetric effect has been developed. It is shown that the previously reported dose-based angular ranking represents a special case of the general formalism proposed here. Application of the technique to a abdominal and a head-and-neck IMRT case indicated that the proposed technique is capable of producing clinically sensible angular ranking. In both cases, we found that the IMRT plans obtained under the guidance of EUD-based angular ranking were improved in comparison with that obtained using the conventional uniformly spaced beams.

CONCLUSIONS

The EUD-based function is a general approach for angular ranking and allows us to identify the potentially good and bad angles for clinically complicated cases. The ranking can be used either as a guidance to facilitate the manual beam placement or as prior information to speed up the computer search for the optimal beam configuration. Thus the proposed technique should have positive clinical impact in facilitating the IMRT planning process.

Towards biologically conformal radiation therapy (BCRT): Selective IMRT dose escalation under the guidance of spatial biology distribution

It is well known that the spatial biology distribution (e.g., clonogen density, radiosensitivity, tumor proliferation rate, functional importance) in most tumors and sensitive structures is heterogeneous. Recent progress in biological imaging is making the mapping of this distribution increasingly possible. The purpose of this work is to establish a theoretical framework to quantitatively incorporate the spatial biology data into intensity modulated radiation therapy (IMRT) inverse planning. In order to implement this, we first derive a general formula for determining the desired dose to each tumor voxel for a known biology distribution of the tumor based on a linear-quadratic model. The desired target dose distribution is then used as the prescription for inverse planning. An objective function with the voxel-dependent prescription is constructed with incorporation of the nonuniform dose prescription. The functional unit density distribution in a sensitive structure is also considered phenomenologically when constructing the objective function. Two cases with different hypothetical biology distributions are used to illustrate the new inverse planning formalism. For comparison, treatments with a few uniform dose prescriptions and a simultaneous integrated boost are also planned. The biological indices, tumor control probability (TCP) and normal tissue complication probability (NTCP), are calculated for both types of plans and the superiority of the proposed technique over the conventional dose escalation scheme is demonstrated. Our calculations revealed that it is technically feasible to produce deliberately nonuniform dose distributions with consideration of biological information. Compared with the conventional dose escalation schemes, the new technique is capable of generating biologically conformal IMRT plans that significantly improve the TCP while reducing or keeping the NTCPs at their current levels. Biologically conformal radiation therapy (BCRT) incorporates patient-specific biological information and provides an outstanding opportunity for us to truly individualize radiation treatment. The proposed formalism lays a technical foundation for BCRT and allows us to maximally exploit the technical capacity of IMRT to more intelligently escalate the radiation dose.