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

An uncertainty-incorporated method for fast beam angle selection in intensity-modulated proton therapy

Aim

We propose a novel metric called ψ - score to rank the Intensity Modulated Proton Therapy (IMPT) beams in the order of their optimality and robustness. The beams ranked based on this metric were accordingly chosen for IMPT optimization. The objective of this work is to study the effectiveness of the proposed method in various clinical cases.

Methods and Materials

We have used Pinnacle TPS (Philips Medical System V 16.2) for performing the optimization. To validate our approach, we have applied it in four clinical cases: Lung, Pancreas, Prostate+Node and Prostate. Basically, for all clinical cases, four set of plans were created using Multi field optimization (MFO) and Robust Optimization (RO) with same clinical objectives, namely (1) Conventional angle plan without Robust Optimization (CA Plan), (2) Suitable angle Plan without Robust Optimization (SA Plan), (3) Conventional angle plan with Robust Optimization (CA-RO Plan), (4) Suitable angle Plan with Robust Optimization (SA-RO Plan). Initial plan was generated with 20 equiangular beams starting from the gantry angle of 0°. In the corresponding SA Plan and SA-RO Plan, the beam angles were obtained using the guidance provided by ψ - score.

Results

All CA plans were compared against the SA plans in terms of Dose distribution, Dose volume histogram (DVH) and percentage of dose difference. The results obtained from the clinical cases indicate that the plan quality is considerably improved without significantly compromising the robustness when the beam angles are optimized using the proposed method. It takes approximately 10-15 min to find the suitable beam angles without Robust Optimization (RO), while it takes approximately 20-30 min to find the suitable beam angles with RO. However, the inclusion of RO in BAO did not result in a change in the final beam angles for anatomies other than lung.

Conclusion

The results obtained in different anatomic sites demonstrate the usefulness of our approach in improving the plan quality by determining optimal beam angles in IMPT.

A data-driven approach to optimal beam/arc angle selection for liver stereotactic body radiation therapy treatment planning

BACKGROUND

Stereotactic body radiation therapy (SBRT) for liver cancer has shown promising therapeutic effects. Effective treatment relies not only on the precise delivery provided by image-guided radiation therapy (IGRT) but also high dose gradient formed around the treatment volume to spare functional liver tissue, which is highly dependent on the beam/arc angle selection. In this study, we aim to develop a decision support model to learn human planner's beam navigation approach for beam angle/arc angle selection for liver SBRT.

METHODS

A total of 27 liver SBRT/HIGRT patients (10 IMRT, 17 VMAT/DCA) were included in this study. A dosimetric budget index was defined for each beam angle/control point considering dose penetration through the patient body and liver tissue. Optimal beam angle setting (beam angles for IMRT and start/terminal angle for VMAT/DCA) was determined by minimizing the loss function defined as the sum of total dosimetric budget index and beam span penalty function. Leave-one-out validation was exercised on all 27 cases while weighting coefficients in the loss function was tuned in nested cross validation. To compare the efficacy of the model, a model plan was generated using automatically generated beam setting while retaining the original optimization constraints in the clinical plan. Model plan was normalized to the same planning target volume (PTV) V100% as the clinical plans. Dosimetric endpoints including PTV D98%, D2%, liver V20Gy and total MU were compared between two plan groups. Wilcoxon Signed-Rank test was performed with the null hypothesis being that no difference exists between two plan groups.

RESULTS

Beam setting prediction was instantaneous. Mean PTV D98% was 91.3% and 91.3% (P=0.566), while mean PTV D2% was 107.9% and 108.1% (P=0.164) for clinical plan and model plan respectively. Liver V20Gy showed no significant difference (P=0.590) with 23.3% for clinical plan and 23.4% for the model plan. Total MU is comparable (P=0.256) between the clinical plan (avg. 2,389.6) and model plan (avg. 2,319.6).

CONCLUSIONS

The evidence driven beam setting model yielded similar plan quality as hand-crafted clinical plan. It is capable of capturing human's knowledge in beam selection decision making. This model could facilitate decision making for beam angle selection while eliminating lengthy trial-and-error process of adjusting beam setting during liver SBRT treatment planning.

A reinforcement learning application of a guided Monte Carlo Tree Search algorithm for beam orientation selection in radiation therapy

Current beam orientation optimization algorithms for radiotherapy, such as column generation (CG), are typically heuristic or greedy in nature because of the size of the combinatorial problem, which leads to suboptimal solutions. We propose a reinforcement learning strategy using a Monte Carlo Tree Search (MCTS) that can find a better beam orientation set in less time than CG. We utilize a reinforcement learning structure involving a supervised learning network to guide the MCTS and to explore the decision space of beam orientation selection problems. We previously trained a deep neural network (DNN) that takes in the patient anatomy, organ weights, and current beams, then approximates beam fitness values to indicate the next best beam to add. Here, we use this DNN to probabilistically guide the traversal of the branches of the Monte Carlo decision tree to add a new beam to the plan. To assess the feasibility of the algorithm, we used a test set of 13 prostate cancer patients, distinct from the 57 patients originally used to train and validate the DNN, to solve five-beam plans. To show the strength of the guided MCTS (GTS) compared to other search methods, we also provided the performances of Guided Search, Uniform Tree Search and Random Search algorithms. On average, GTS outperformed all the other methods. It found a better solution than CG in 237 s on average, compared to 360 s for CG, and outperformed all other methods in finding a solution with a lower objective function value in less than 1000 s. Using our GTS method, we could maintain planning target volume (PTV) coverage within 1% error similar to CG, while reducing the organ-at-risk mean dose for body, rectum, left and right femoral heads; the mean dose to bladder was 1% higher with GTS than with CG.

A fast deep learning approach for beam orientation optimization for prostate cancer treated with intensity-modulated radiation therapy

PURPOSE

Beam orientation selection, whether manual or protocol-based, is the current clinical standard in radiation therapy treatment planning, but it is tedious and can yield suboptimal results. Many algorithms have been designed to optimize beam orientation selection because of its impact on treatment plan quality, but these algorithms suffer from slow calculation of the dose influence matrices of all candidate beams. We propose a fast beam orientation selection method, based on deep learning neural networks (DNN), capable of developing a plan comparable to those developed by the state-of-the-art column generation (CG) method. Our model's novelty lies in its supervised learning structure (using CG to teach the network), DNN architecture, and ability to learn from anatomical features to predict dosimetrically suitable beam orientations without using dosimetric information from the candidate beams. This may save hours of computation.

METHODS

A supervised DNN is trained to mimic the CG algorithm, which iteratively chooses beam orientations one-by-one by calculating beam fitness values based on Karush-Kush-Tucker optimality conditions at each iteration. The DNN learns to predict these values. The dataset contains 70 prostate cancer patients - 50 training, 7 validation, and 13 test patients - to develop and test the model. Each patient's data contains 6 contours: PTV, body, bladder, rectum, and left and right femoral heads. Column generation was implemented with a GPU-based Chambolle-Pock algorithm, a first-order primal-dual proximal-class algorithm, to create 6270 plans. The DNN trained over 400 epochs, each with 2500 steps and a batch size of 1, using the Adam optimizer at a learning rate of 1 × 10-5 and a sixfold cross-validation technique.

RESULTS

The average and standard deviation of training, validation, and testing loss functions among the six folds were 0.62 ± 0.09%, 1.04 ± 0.06%, and 1.44 ± 0.11%, respectively. Using CG and supervised DNN, we generated two sets of plans for each scenario in the test set. The proposed method took at most 1.5 s to select a set of five beam orientations and 300 s to calculate the dose influence matrices for 5 beams and finally 20 s to solve the fluence map optimization (FMO). However, CG needed around 15 h to calculate the dose influence matrices of all beams and at least 400 s to solve both the beam orientation selection and FMO problems. The differences in the dose coverage of PTV between plans generated by CG and by DNN were 0.2%. The average dose differences received by organs at risk were between 1 and 6 percent: Bladder had the smallest average difference in dose received (0.956 ± 1.184%), then Rectum (2.44 ± 2.11%), Left Femoral Head (6.03 ± 5.86%), and Right Femoral Head (5.885 ± 5.515%). The dose received by Body had an average difference of 0.10 ± 0.1% between the generated treatment plans.

CONCLUSIONS

We developed a fast beam orientation selection method based on a DNN that selects beam orientations in seconds and is therefore suitable for clinical routines. In the training phase of the proposed method, the model learns the suitable beam orientations based on patients' anatomical features and omits time intensive calculations of dose influence matrices for all possible candidate beams. Solving the FMO to get the final treatment plan requires calculating dose influence matrices only for the selected beams.

Lung IMRT planning with automatic determination of beam angle configurations

Beam angle configuration is a major planning decision in intensity modulated radiation treatment (IMRT) that has a significant impact on dose distributions and thus quality of treatment, especially in complex planning cases such as those for lung cancer treatment. We propose a novel method to automatically determine beam configurations that incorporates noncoplanar beams. We then present a completely automated IMRT planning algorithm that combines the proposed method with a previously reported OAR DVH prediction model. Finally, we validate this completely automatic planning algorithm using a set of challenging lung IMRT cases. A beam efficiency index map is constructed to guide the selection of beam angles. This index takes into account both the dose contributions from individual beams and the combined effect of multiple beams by introducing a beam-spread term. The effect of the beam-spread term on plan quality was studied systematically and the weight of the term to balance PTV dose conformity against OAR avoidance was determined. For validation, complex lung cases with clinical IMRT plans that required the use of one or more noncoplanar beams were re-planned with the proposed automatic planning algorithm. Important dose metrics for the PTV and OARs in the automatic plans were compared with those of the clinical plans. The results are very encouraging. The PTV dose conformity and homogeneity in the automatic plans improved significantly. And all the dose metrics of the automatic plans, except the lung V_{5 Gy}, were statistically better than or comparable with those of the clinical plans. In conclusion, the automatic planning algorithm can incorporate non-coplanar beam configurations in challenging lung cases and can generate plans efficiently with quality closely approximating that of clinical plans.

Automatic learning-based beam angle selection for thoracic IMRT

PURPOSE

The treatment of thoracic cancer using external beam radiation requires an optimal selection of the radiation beam directions to ensure effective coverage of the target volume and to avoid unnecessary treatment of normal healthy tissues. Intensity modulated radiation therapy (IMRT) planning is a lengthy process, which requires the planner to iterate between choosing beam angles, specifying dose-volume objectives and executing IMRT optimization. In thorax treatment planning, where there are no class solutions for beam placement, beam angle selection is performed manually, based on the planner's clinical experience. The purpose of this work is to propose and study a computationally efficient framework that utilizes machine learning to automatically select treatment beam angles. Such a framework may be helpful for reducing the overall planning workload.

METHODS

The authors introduce an automated beam selection method, based on learning the relationships between beam angles and anatomical features. Using a large set of clinically approved IMRT plans, a random forest regression algorithm is trained to map a multitude of anatomical features into an individual beam score. An optimization scheme is then built to select and adjust the beam angles, considering the learned interbeam dependencies. The validity and quality of the automatically selected beams evaluated using the manually selected beams from the corresponding clinical plans as the ground truth.

RESULTS

The analysis included 149 clinically approved thoracic IMRT plans. For a randomly selected test subset of 27 plans, IMRT plans were generated using automatically selected beams and compared to the clinical plans. The comparison of the predicted and the clinical beam angles demonstrated a good average correspondence between the two (angular distance 16.8° ± 10°, correlation 0.75 ± 0.2). The dose distributions of the semiautomatic and clinical plans were equivalent in terms of primary target volume coverage and organ at risk sparing and were superior over plans produced with fixed sets of common beam angles. The great majority of the automatic plans (93%) were approved as clinically acceptable by three radiation therapy specialists.

CONCLUSIONS

The results demonstrated the feasibility of utilizing a learning-based approach for automatic selection of beam angles in thoracic IMRT planning. The proposed method may assist in reducing the manual planning workload, while sustaining plan quality.

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

Beam geometry selection using sequential beam addition

PURPOSE

The selection of optimal beam geometry has been of interest since the inception of conformal radiotherapy. The authors report on sequential beam addition, a simple beam geometry selection method, for intensity modulated radiation therapy.

METHODS

The sequential beam addition algorithm (SBA) requires definition of an objective function (score) and a set of candidate beam geometries (pool). In the first iteration, the optimal score is determined for each beam in the pool and the beam with the best score selected. In the next iteration, the optimal score is calculated for each beam remaining in the pool combined with the beam selected in the first iteration, and the best scoring beam is selected. The process is repeated until the desired number of beams is reached. The authors selected three treatment sites, breast, lung, and brain, and determined beam arrangements for up to 11 beams from a pool comprised of 25 equiangular transverse beams. For the brain, arrangements were additionally selected from a pool of 22 noncoplanar beams. Scores were determined for geometries comprised equiangular transverse beams (EQA), as well as two tangential beams for the breast case.

RESULTS

In all cases, SBA resulted in scores superior to EQA. The breast case had the strongest dependence on beam geometry, for which only the 7-beam EQA geometry had a score better than the two tangential beams, whereas all SBA geometries with more than two beams were superior. In the lung case, EQA and SBA scores monotonically improved with increasing number of beams; however, SBA required fewer beams to achieve scores equivalent to EQA. For the brain case, SBA with a coplanar pool was equivalent to EQA, while the noncoplanar pool resulted in slightly better scores; however, the dose-volume histograms demonstrated that the differences were not clinically significant.

CONCLUSIONS

For situations in which beam geometry has a significant effect on the objective function, SBA can identify arrangements equivalent to equiangular geometries but using fewer beams. Furthermore, SBA provides the value of the objective function as the number of beams is increased, allowing the planner to select the minimal beam number that achieves the clinical goals. The method is simple to implement and could readily be incorporated into an existing optimization system.

Nonisocentric treatment strategy for breast radiation therapy: a proof of concept study

PURPOSE

To propose a nonisocentric treatment strategy as a special form of station parameter optimized radiation therapy, to improve sparing of critical structures while preserving target coverage in breast radiation therapy.

METHODS AND MATERIALS

To minimize the volume of exposed lung and heart in breast irradiation, we propose a novel nonisocentric treatment scheme by strategically placing nonconverging beams with multiple isocenters. As its name suggests, the central axes of these beams do not intersect at a single isocenter as in conventional breast treatment planning. Rather, the isocenter locations and beam directions are carefully selected, in that each beam is only responsible for a certain subvolume of the target, so as to minimize the volume of irradiated normal tissue. When put together, the beams will provide an adequate coverage of the target and expose only a minimal amount of normal tissue to radiation. We apply the nonisocentric planning technique to 2 previously treated clinical cases (breast and chest wall).

RESULTS

The proposed nonisocentric technique substantially improved sparing of the ipsilateral lung. Compared with conventional isocentric plans using 2 tangential beams, the mean lung dose was reduced by 38% and 50% using the proposed technique, and the volume of the ipsilateral lung receiving ≥ 20 Gy was reduced by a factor of approximately 2 and 3 for the breast and chest wall cases, respectively. The improvement in lung sparing is even greater compared with volumetric modulated arc therapy.

CONCLUSIONS

A nonisocentric implementation of station parameter optimized radiation therapy has been proposed for breast radiation therapy. The new treatment scheme overcomes the limitations of existing approaches and affords a useful tool for conformal breast radiation therapy, especially in cases with extreme chest wall curvature.

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

Improving IMRT delivery efficiency with reweighted L1-minimization for inverse planning

PURPOSE

This study presents an improved technique to further simplify the fluence-map in intensity modulated radiation therapy (IMRT) inverse planning, thereby reducing plan complexity and improving delivery efficiency, while maintaining the plan quality.

METHODS

First-order total-variation (TV) minimization (min.) based on L1-norm has been proposed to reduce the complexity of fluence-map in IMRT by generating sparse fluence-map variations. However, with stronger dose sparing to the critical structures, the inevitable increase in the fluence-map complexity can lead to inefficient dose delivery. Theoretically, L0-min. is the ideal solution for the sparse signal recovery problem, yet practically intractable due to its nonconvexity of the objective function. As an alternative, the authors use the iteratively reweighted L1-min. technique to incorporate the benefits of the L0-norm into the tractability of L1-min. The weight multiplied to each element is inversely related to the magnitude of the corresponding element, which is iteratively updated by the reweighting process. The proposed penalizing process combined with TV min. further improves sparsity in the fluence-map variations, hence ultimately enhancing the delivery efficiency. To validate the proposed method, this work compares three treatment plans obtained from quadratic min. (generally used in clinic IMRT), conventional TV min., and our proposed reweighted TV min. techniques, implemented by a large-scale L1-solver (template for first-order conic solver), for five patient clinical data. Criteria such as conformation number (CN), modulation index (MI), and estimated treatment time are employed to assess the relationship between the plan quality and delivery efficiency.

RESULTS

The proposed method yields simpler fluence-maps than the quadratic and conventional TV based techniques. To attain a given CN and dose sparing to the critical organs for 5 clinical cases, the proposed method reduces the number of segments by 10-15 and 30-35, relative to TV min. and quadratic min. based plans, while MIs decreases by about 20%-30% and 40%-60% over the plans by two existing techniques, respectively. With such conditions, the total treatment time of the plans obtained from our proposed method can be reduced by 12-30 s and 30-80 s mainly due to greatly shorter multileaf collimator (MLC) traveling time in IMRT step-and-shoot delivery.

CONCLUSIONS

The reweighted L1-minimization technique provides a promising solution to simplify the fluence-map variations in IMRT inverse planning. It improves the delivery efficiency by reducing the entire segments and treatment time, while maintaining the plan quality in terms of target conformity and critical structure sparing.

An adaptive planning strategy for station parameter optimized radiation therapy (SPORT): Segmentally boosted VMAT

PURPOSE

Conventional volumetric modulated arc therapy (VMAT) discretizes the angular space into equally spaced control points during planning and then optimizes the apertures and weights of the control points. The aperture at an angle in between two control points is obtained through interpolation. This approach tacitly ignores the differential need for intensity modulation of different angles. As such, multiple arcs are often required, which may oversample some angle(s) and undersample others. The purpose of this work is to develop a segmentally boosted VMAT scheme to eliminate the need for multiple arcs in VMAT treatment with improved dose distribution and∕or delivery efficiency.

METHODS

The essence of the new treatment scheme is how to identify the need of individual angles for intensity modulation and to provide the necessary beam intensity modulation for those beam angles that need it. We introduce a "demand metric" at each control point to decide which station or control points need intensity modulation. To boost the modulation at selected stations, additional segments are added in the vicinity of the selected stations. The added segments are then optimized together with the original set of station or control points as a whole. The authors apply the segmentally boosted planning technique to four previously treated clinical cases: two head and neck (HN) cases, one prostate case, and one liver case. The proposed planning technique is compared with conventional one-arc and two-arc VMAT.

RESULTS

The proposed segmentally boosted VMAT technique achieves better critical structure sparing than one-arc VMAT with similar or better target coverage in all four clinical cases. The segmentally boosted VMAT also outperforms two-arc VMAT for the two complicated HN cases, yet with ∼30% reduction in the machine monitor units (MUs) relative to two-arc VMAT, which leads to less leakage∕scatter dose to the patient and can potentially translate into faster dose delivery. For the less challenging prostate and liver cases, similar critical structure sparing as the two-arc VMAT plans was obtained using the segmentally boosted VMAT. The benefit for the two simpler cases is the reduction of MUs and improvement of treatment delivery efficiency.

CONCLUSIONS

Segmentally boosted VMAT achieves better dose conformality and∕or reduced MUs through effective consideration of the need of individual beam angles for intensity modulation. Elimination of the need for multiple arcs in rotational arc therapy while improving the dose distribution should lead to improved workflow and treatment efficacy, thus may have significant implication to radiation oncology practice.

Bridging the gap between IMRT and VMAT: dense angularly sampled and sparse intensity modulated radiation therapy

PURPOSE

To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency.

METHODS

In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency.

RESULTS

Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.

CONCLUSION

DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

An algorithm for fast beam angle selection in intensity modulated radiotherapy

PURPOSE

This article aims to introduce a novel algorithm for fast beam angle selection in intensity modulated radiotherapy (IMRT).

METHODS

The algorithm models the optimization problem as a beam angle ranking problem and chooses suitable beam angles according to their rank. A new parameter called "beam intensity profile perturbation score (BIPPS)" is used for ranking the beam angles. The BIPPS-based beam angle ranking implicitly accounts for the dose-volume effects of the involved structures. A simulated phantom case with obvious optimal beam angles is used to verify the validity of the presented technique. In addition, the efficiency of the algorithm was examined in three clinical cases (prostate, pancreas, and head and neck) in terms of DVH and dose distribution. In all cases, the judgment of the algorithm's efficiency was based on the comparison between plans with equidistant beams (equal-angle-plan) and plans with beams obtained using the algorithm (suitable-angle-plan).

RESULTS

It is observed from the study that the beam angle ranking function over BIPPS instantly picks up a suitable set of beam angles for a specific case. It takes only about 15 min for choosing the suitable beam angles even for the most complicated cases. The DVHs and dose distributions confirm that the proposed algorithm can efficiently reduce the mean or maximum dose to OARs, while guaranteeing the target coverage and dose uniformity. On the average, about 17% reduction in the mean dose to critical organs, such as rectum, bladder, kidneys and parotids, is observed. Also, about 12% (averaged) reduction in the maximum dose to critical organs (spinal cord) is observed in the clinical cases presented in this study.

CONCLUSIONS

This study demonstrates that the algorithm can be effectively applied to IMRT scenarios to get fast and case specific beam angle configurations.

A methodology for selecting the beam arrangement to reduce the intensity-modulated radiation therapy (IMRT) dose to the SPECT-defined functioning lung

Macroaggregated albumin single-photon emission computed tomography (MAA-SPECT) provides a map of the spatial distribution of lung perfusion. Our previous work developed a methodology to use SPECT guidance to reduce the dose to the functional lung in IMRT planning. This study aims to investigate the role of beam arrangement on both low and high doses in the functional lung. In our previous work, nine-beam IMRT plans were generated with and without SPECT guidance and compared for five patients. For the current study, the dose-function histogram (DFH) contribution for each of the nine beams for each patient was calculated. Four beams were chosen based on orientation and DFH contributions to create a SPECT-guided plan that spared the functional lung and maintained target coverage. Four-beam SPECT-guided IMRT plans reduced the F(20) and F(30) values by (16.5 +/- 6.8)% and (6.1 +/- 9.2)%, respectively, when compared to nine-beam conventional IMRT plans. Moreover, the SPECT-4F Plan reduces F(5) and F(13) for all patients by (11.0 +/- 8.2)% and (6.1 +/- 3.6)%, respectively, compared to the SPECT Plan. Using fewer beams in IMRT planning may reduce the amount of functional lung that receives 5 and 13 Gy, a factor that has recently been associated with radiation pneumonitis.

Future Issues in Highly Conformal Radiotherapy for Head and Neck Cancer

Improving the conformity of the radiation dose to targets in the head and neck promises reduced toxicity and, in some cases, potentially improved local-regional tumor control. Intensity-modulated radiotherapy (IMRT) is a method that allows highly conformal delivery of radiotherapy. In recent years, its use has spread rapidly in both academic and community radiation oncology facilities. The use of IMRT has raised multiple issues related to target definition, optimal treatment delivery methods, and the need to account for anatomic changes occurring during therapy. Some of these issues are reviewed in this article.

Direct-aperture optimization applied to selection of beam orientations in intensity-modulated radiation therapy

Direct-aperture optimization (DAO) was applied to iterative beam-orientation selection in intensity-modulated radiation therapy (IMRT), so as to ensure a realistic segmental treatment plan at each iteration. Nested optimization engines dealt separately with gantry angles, couch angles, collimator angles, segment shapes, segment weights and wedge angles. Each optimization engine performed a random search with successively narrowing step sizes. For optimization of segment shapes, the filtered backprojection (FBP) method was first used to determine desired fluence, the fluence map was segmented, and then constrained direct-aperture optimization was used thereafter. Segment shapes were fully optimized when a beam angle was perturbed, and minimally re-optimized otherwise. The algorithm was compared with a previously reported method using FBP alone at each orientation iteration. An example case consisting of a cylindrical phantom with a hemi-annular planning target volume (PTV) showed that for three-field plans, the method performed better than when using FBP alone, but for five or more fields, neither method provided much benefit over equally spaced beams. For a prostate case, improved bladder sparing was achieved through the use of the new algorithm. A plan for partial scalp treatment showed slightly improved PTV coverage and lower irradiated volume of brain with the new method compared to FBP alone. It is concluded that, although the method is computationally intensive and not suitable for searching large unconstrained regions of beam space, it can be used effectively in conjunction with prior class solutions to provide individually optimized IMRT treatment plans.

Novel lung IMRT planning algorithms with nonuniform dose delivery strategy to account for respiratory motion

To effectively deliver radiation dose to lung tumors, respiratory motion has to be considered in treatment planning. In this paper we first present a new lung IMRT planning algorithm, referred as the dose shaping (DS) method, that shapes the dose distribution according to the probability distribution of the tumor over the breathing cycle to account for respiratory motion. In IMRT planning a dose-based convolution method was generally adopted to compensate for random organ motion by performing 4-D dose calculations using a tumor motion probability density function. We modified the CON-DOSE method to a dose volume histogram based convolution method (CON-DVH) that allows nonuniform dose distribution to account for respiratory motion. We implemented the two new planning algorithms on an in-house IMRT planning system that uses the Eclipse (Varian, Palo Alto, CA) planning workstation as the dose calculation engine. The new algorithms were compared with (1) the conventional margin extension approach in which margin is generated based on the extreme positions of the tumor, (2) the dose-based convolution method, and (3) gating with 3 mm residual motion. Dose volume histogram, tumor control probability, normal tissue complication probability, and mean lung dose were calculated and used to evaluate the relative performance of these approaches at the end-exhale phase of the respiratory cycle. We recruited six patients in our treatment planning study. The study demonstrated that the two new methods could significantly reduce the ipsilateral normal lung dose and outperformed the margin extension method and the dose-based convolution method. Compared with the gated approach that has the best performance in the low dose region, the two methods we proposed have similar potential to escalate tumor dose, but could be more efficient because dose is delivered continuously.