A discrete differential evolution with local search particle swarm optimization to direct angle and aperture optimization in IMRT treatment planning problem

Intensity Modulated Radiation Therapy (IMRT) is a means of delivering radiation therapy where the intensity of the beam is varied within the treatment field. One of the main problems of IMRT, which becomes even more apparent as the complexity of the IMRT plan increases, is the amazing increase in the number of Monitor Units (MU) required to deliver a fractionated treatment. The difficulty with this increase in MU is its association with increased treatment times and a greater leakage of radiation from the MLCs increasing the total body dose and the risk of secondary cancers in patients. Many approaches have been taken to reduce the complexity. In this study, we used a direct aperture optimization (DAO) IMRT plan Direct aperture optimization IMRT is an IMRT method in which the aperture shapes and aperture weights are optimized simultaneously, and the MLC constraints and the number of segments are directly included in the optimization process. In DAO-IMRT planning, the planner specifies the planning objectives on the basis of the dose volume criteria for the target and critical structures as well as the number of beam segments to be delivered. The optimization only considers aperture shapes that satisfy the conditions set by the MLC. As a result, high-quality DAO-IMRT treatment plans can be generated using fewer segments (apertures) per beam. We investigated the dosimetric and technical differences in using DAO-IMRT. We also compared DAO-IMRT plans with standard beamlet-IMRT plans on phantom in terms of dosimetry, number of segments, and monitor units. In this study, we found that DAO-IMRT can achieve equal or better plans than standard beamlet-IMRT in dose optimization, secondary scattering dose and treatment time. The DAO-IMRT planning is proficient, and the plans can be delivered efficiently because small numbers of segments and monitor units are used.KeywordsIMRTnumber of segmentstreatment timemonitor unitssecondary scattering dose

Image Guided Radiation Therapy for Prostate IMRT: Rectum Volume Changes and Dosimetric Considerations

Currently, intensity-modulated radiation therapy (IMRT) is commonly used in radiotherapy clinics. However, designing a treatment plan with multiple beam angles depends on the experience of human planners, and is mostly achieved using a trial-and-error approach. It is preferrable but challenging to solve this issue automatically and mathematically using an optimization approach. The goal of this study is to develop a mixed-integer linear programming (MILP) approach for the beam angle optimization (BAO) of non-coplanar IMRT for liver cancer. MILP models for the BAO of both coplanar and non-coplanar IMRT treatment plans were developed. The beam angles of the IMRT plans were first selected by the MILP model built using mathematical optimization software. Next, the IMRT plans with the selected beam angles was created in a commercial treatment planning system. Finally, the fluence map and dose distribution of the IMRT plans were generated under pre-defined dose-volume constraints. The IMRT plans of 10 liver cancer patients previously treated at our institute were used to assessed the proposed MILP models. For each patient, both coplanar and non-coplanar IMRT plans with beam angles optimized by the MILP models were compared with the IMRT plan clinically approved by physicians. The MILP model-guided IMRT plans showed reduced doses for most of the organs at risk (OARs). Compared with the IMRT plans clinically approved by physicians, the doses for the spinal cord (28.5 vs. 36.1, P=0.001<0.05) and liver (27.6 vs. 29.1, P=0.005<0.05) decreased significantly in the IMRT plans with non-coplanar beams selected by the MILP models. The MILP model is an effective tool for the BAO in coplanar and non-coplanar IMRT treatment planning. It facilitates the automation of IMRT treatment planning for current high-precision radiotherapy.

Purpose: To evaluate the dosimetric characteristics of intensity-modulated radiotherapy (IMRT) treatment plan with beam angle optimization. Methods: Ten post-operation patients with cervical cancer were included in this analysis. Two IMRT plans using seven beams were designed in each patient. A standard coplanar equi-space beam angles were used in the first plan (plan 1), whereas the selection of beam angle was optimized by beam angle optimization algorithm in Varian Eclipse treatment planning system for the same number of beams in the second plan (plan 2). Two plans were designed for each patient with the same dose-volume constraints and prescription dose. All plans were normalized to the mean dose to PTV. The dose distribution in the target, the dose to the organs at risk and total MU were compared. Results: For conformity and homogeneity in PTV, no statistically differences were observed in the two plans. For the mean dose in bladder, plan 2 were significantly lower than plan 1(p<0.05). No statistically significant differences were observed between two plans for the mean doses in rectum, left and right femur heads. Compared with plan1, the average monitor units reduced 16% in plan 2. Conclusion: The IMRT plan based on beam angle optimization for cervical cancer could reduce the dose delivered to bladder and also reduce MU. Therefore there were some dosimetric advantages in the IMRT plan with beam angle optimization for cervical cancer.

A Treatment Planning System for Intensity-Modulated Small Animal Radiotherapy Realized Using an Existing Movable Rectangular Collimator

Purpose:To evaluate the dosimetric characteristics of intensity‐modulated radiotherapy (IMRT) treatment plan with beam angle optimization.Methods:Ten post‐operation patients with cervical cancer were included in this analysis. Two IMRT plans using seven beams were designed in each patient. A standard coplanar equi‐space beam angles were used in the first plan (plan 1), whereas the selection of beam angle was optimized by beam angle optimization algorithm in Varian Eclipse treatment planning system for the same number of beams in the second plan (plan 2). Two plans were designed for each patient with the same dose‐volume constraints and prescription dose. All plans were normalized to the mean dose to PTV. The dose distribution in the target, the dose to the organs at risk and total MU were compared.Results:For conformity and homogeneity in PTV, no statistically differences were observed in the two plans. For the mean dose in bladder, plan 2 were significantly lower than plan 1(p&lt;0.05). No statistically significant differences were observed between two plans for the mean doses in rectum, left and right femur heads. Compared with plan1, the average monitor units reduced 16% in plan 2.Conclusion:The IMRT plan based on beam angle optimization for cervical cancer could reduce the dose delivered to bladder and also reduce MU. Therefore there were some dosimetric advantages in the IMRT plan with beam angle optimization for cervical cancer.

We consider the problem of intensity-modulated radiation therapy (IMRT) treatment planning using direct aperture optimization. While this problem has been relatively well studied in recent years, most approaches employ a heuristic approach to the generation of apertures. In contrast, we use an exact approach that explicitly formulates the fluence map optimization (FMO) problem as a convex optimization problem in terms of all multileaf collimator (MLC) deliverable apertures and their associated intensities. However, the number of deliverable apertures, and therefore the number of decision variables and constraints in the new problem formulation, is typically enormous. To overcome this, we use an iterative approach that employs a subproblem whose optimal solution either provides a suitable aperture to add to a given pool of allowable apertures or concludes that the current solution is optimal. We are able to handle standard consecutiveness, interdigitation and connectedness constraints that may be imposed by the particular MLC system used, as well as jaws-only delivery. Our approach has the additional advantage that it can explicitly account for transmission of dose through the part of an aperture that is blocked by the MLC system, yielding a more precise assessment of the treatment plan than what is possible using a traditional beamlet-based FMO problem. Finally, we develop and test two stopping rules that can be used to identify treatment plans of high clinical quality that are deliverable very efficiently. Tests on clinical head-and-neck cancer cases showed the efficacy of our approach, yielding treatment plans comparable in quality to plans obtained by the traditional method with a reduction of more than 75% in the number of apertures and a reduction of more than 50% in beam-on time, with only a modest increase in computational effort. The results also show that delivery efficiency is very insensitive to the addition of traditional MLC constraints; however, jaws-only treatment requires about a doubling in beam-on time and number of apertures used. Finally, we showed the importance of accounting for transmission effects when assessing or, preferably, optimizing treatment plan quality.

Inverse planned intensity modulated radiation therapy (IMRT) has helped many centers implement highly conformal treatment planning with beamlet-based techniques. The many comparisons between IMRT and 3D conformal (3DCRT) plans, however, have been limited because most 3DCRT plans are forward-planned while IMRT plans utilize inverse planning, meaning both optimization and delivery techniques are different. This work avoids that problem by comparing 3D plans generated with a unique inverse planning method for 3DCRT called inverse-optimized 3D (IO-3D) conformal planning. Since IO-3D and the beamlet IMRT to which it is compared use the same optimization techniques, cost functions, and plan evaluation tools, direct comparisons between IMRT and simple, optimized IO-3D plans are possible. Though IO-3D has some similarity to direct aperture optimization (DAO), since it directly optimizes the apertures used, IO-3D is specifically designed for 3DCRT fields (i.e., 1-2 apertures per beam) rather than starting with IMRT-like modulation and then optimizing aperture shapes. The two algorithms are very different in design, implementation, and use. The goals of this work include using IO-3D to evaluate how close simple but optimized IO-3D plans come to nonconstrained beamlet IMRT, showing that optimization, rather than modulation, may be the most important aspect of IMRT (for some sites). The IO-3D dose calculation and optimization functionality is integrated in the in-house 3D planning/optimization system. New features include random point dose calculation distributions, costlet and cost function capabilities, fast dose volume histogram (DVH) and plan evaluation tools, optimization search strategies designed for IO-3D, and an improved, reimplemented edge/octree calculation algorithm. The IO-3D optimization, in distinction to DAO, is designed to optimize 3D conformal plans (one to two segments per beam) and optimizes MLC segment shapes and weights with various user-controllable search strategies which optimize plans without beamlet or pencil beam approximations. IO-3D allows comparisons of beamlet, multisegment, and conformal plans optimized using the same cost functions, dose points, and plan evaluation metrics, so quantitative comparisons are straightforward. Here, comparisons of IO-3D and beamlet IMRT techniques are presented for breast, brain, liver, and lung plans. IO-3D achieves high quality results comparable to beamlet IMRT, for many situations. Though the IO-3D plans have many fewer degrees of freedom for the optimization, this work finds that IO-3D plans with only one to two segments per beam are dosimetrically equivalent (or nearly so) to the beamlet IMRT plans, for several sites. IO-3D also reduces plan complexity significantly. Here, monitor units per fraction (MU/Fx) for IO-3D plans were 22%-68% less than that for the 1 cm × 1 cm beamlet IMRT plans and 72%-84% than the 0.5 cm × 0.5 cm beamlet IMRT plans. The unique IO-3D algorithm illustrates that inverse planning can achieve high quality 3D conformal plans equivalent (or nearly so) to unconstrained beamlet IMRT plans, for many sites. IO-3D thus provides the potential to optimize flat or few-segment 3DCRT plans, creating less complex optimized plans which are efficient and simple to deliver. The less complex IO-3D plans have operational advantages for scenarios including adaptive replanning, cases with interfraction and intrafraction motion, and pediatric patients.

In a treatment planning system for intensity modulated radiation therapy (IMRT), the time sequence of multileaf collimator (MLC) settings are derived from an optimal fluence map as a postoptimization process using a software module called a "leaf sequencer." The dosimetric accuracy of the dynamic delivery depends on the functionality of the module and it is important to verify independently the correctness of the leaf sequences for each field of a patient treatment. This verification is unique to the IMRT treatment and has been done using radiographic film, electronic portal imaging device (EPID) or electronic imaging system (BIS). The measurement tests both the leaf sequencer and the dynamic multileaf collimator (MLC) delivery system, providing a reliable assurance of clinical IMRT treatment. However, this process is labor intensive and time consuming. In this paper, we propose to separate quality assurance (QA) of the leaf sequencer from the dynamic MLC delivery system. We describe a simple computer algorithm for the verification of the leaf sequences. The software reads in the leaf sequences and simulates the motion of the MLC leaves. The generated fluence map is then compared quantitatively with the reference map from the treatment planning system. A set of pre-defined QA indices is introduced to measure the "closeness" between the computed and the reference maps. The approach has been used to validate the CORVUS (NOMOS Co., Sewickley, PA) treatment plans. The results indicate that the proposed approach is robust and suitable to support the complex IMRT QA process.

MERT vs tangential photon beams vs IMRT for breast cancer: a comparative dosimetry study

Purpose: The main purpose of this study is to perform a treatment planning study on a synchronous bilateral non-small cell lung cancer case using three treatment modalities: uniform scanning proton therapy, RapidArc, and intensity modulated radiation therapy (IMRT). Methods : The maximum intensity projection (MIP) images obtained from the 4 dimensional-computed tomography (4DCT) scans were used for delineation of tumor volumes in the left and right lungs. The average 4D-CT was used for the treatment planning among all three modalities with identical patient contouring and treatment planning goal. A proton therapy plan was generated in XiO treatment planning system (TPS) using 2 fields for each target. For a comparative purpose, IMRT and RapidArc plans were generated in Eclipse TPS. Treatment plans were generated for a total dose of 74 CGE or Gy prescribed to each planning target volume (PTV) (left and right) with 2 CGE or Gy per fraction. In IMRT and RapidArc plans, normalization was done based on PTV coverage values in proton plans. Results: The mean PTV dose deviation from the prescription dose was lower in proton plan (within 3.4%), but higher in IMRT (6.5% to 11.3%) and RapidArc (3.8% to 11.5%) plans. Proton therapy produced lower mean dose to the total lung, heart, and esophagus when compared to IMRT and RapidArc. The relative volume of the total lung receiving 20, 10, and 5 CGE or Gy (V20, V10, and V5, respectively) were lower using proton therapy than using IMRT, with absolute differences of 9.71%, 22.88%, and 39.04%, respectively. The absolute differences in the V20, V10, and V5 between proton and RapidArc plans were 4.84%, 19.16%, and 36.8%, respectively, with proton therapy producing lower dosimetric values. Conclusion : Based on the results presented in this case study, uniform scanning proton therapy has a dosimetric advantage over both IMRT and RapidArc for a synchronous bi-lateral NSCLC, especially for the normal lung tissue, heart, and esophagus sparing. Further studies on a large group of patients with bi-lateral lung cancer are required to validate the dosimetric superiority of proton therapy over the IMRT and RapidArc. ------------------------------- Cite this article as: Rana S, Pokharel S, Zheng Y, Zhao L, Risalvato D, Vargas C, Cersonsky N. Treatment planning study comparing proton therapy, RapidArc and intensity modulated radiation therapy for a synchronous bilateral lung cancer case. Int J Cancer Ther Oncol 2014; 2 (2):020216. DOI: 10.14319/ijcto.0202.16

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

The inverse planning of an intensity-modulated radiation therapy (IMRT) treatment requires decisions regarding the angles used for radiation incidence, even when arcs are used. The possibility of improving the quality of treatment plans by an optimized selection of the beam angle incidences-beam angle optimization (BAO)-is seldom done in clinical practice. The inclusion of noncoplanar beam incidences in an automated optimization routine is even more unusual. However, for some tumor sites, the advantage of considering noncoplanar beam incidences is well known. This paper presents the benefits of using a derivative-free multistart framework for the optimization of the noncoplanar BAO problem. Multistart methods combine a global strategy for sampling the search space with a local strategy for improving the sampled solutions. The proposed global strategy allows a thorough exploration of the continuous search space of the highly nonconvex BAO problem. To avoid local entrapment, a derivative-free method is used as local procedure. Additional advantages of the derivative-free method include the reduced number of function evaluations required to converge and the ability to use multithreaded computing. Twenty nasopharyngeal clinical cases were selected to test the proposed multistart framework. The planning target volumes included the primary tumor, the high and low risk lymph nodes. Organs-at-risk included the spinal cord, brainstem, optical nerves, chiasm, parotids, oral cavity, brain, thyroid, among others. For each case, a setup with seven equispaced beams was chosen and the resulting treatment plan, using a multicriteria optimization framework, was then compared against the coplanar and noncoplanar plans using the optimal beam setups obtained by the derivative-free multistart framework. The optimal noncoplanar beam setup obtained by the derivative-free multistart framework leads to high quality treatment plans with better target coverage and with improved organ sparing compared to treatment plans using equispaced or optimal coplanar beam angle setups. The noncoplanar treatment plans achieved, e.g., an average reduction in the mean dose of the oral cavity of 6.1 Gy and an average reduction in the maximum-dose of the brainstem of 7 Gy when compared to the equispaced treatment plans. The noncoplanar BAO problem is an extremely challenging multimodal optimization problem that can be successfully addressed through a thoughtful exploration of the continuous highly nonconvex BAO search space. The proposed framework is capable of calculating high quality treatment plans and thus can be an interesting alternative toward automated noncoplanar beam selection in IMRT treatment planning which is nowadays the natural trend in treatment planning.

Purpose:Hippocampus is one of the important organs which controls emotions, behaviors, movements the memorizing and learning ability. In the conventional head &amp; neck therapy position, it is difficult to perform the hippocampal‐sparing brain radiation therapy. The purpose of this study is to investigate optimal head angle which can save the hippocampal‐sparing and organ at risk (OAR) in conformal radiation therapy (CRT), Intensity modulation radiation therapy (IMRT) and helical tomotherapy (HT).Methods:Three types of radiation treatment plans, CRT, IMRT and Tomotherapy plans, were performed for 10 brain tumor patients. The image fusion between CT and MRI data were used in the contour due to the limited delineation of the target and OAR in the CT scan. The optimal condition plan was determined by comparing the dosimetric performance of the each plan with the use of various parameters which include three different techniques (CRT, IMRT, HT) and 4 angle (0, 15, 30, 40 degree). The each treatment plans of three different techniques were compared with the following parameters: conformity index (CI), homogeneity index (HI), target coverage, dose in the OARs, monitor units (MU), beam on time and the normal tissue complication probability (NTCP).Results:HI, CI and target coverage was most excellent in head angle 30 degree among all angle. When compared by modality, target coverage and CI showed good results in IMRT and TOMO than compared to the CRT. HI at the head angle 0 degrees is 1.137±0.17 (CRT), 1.085±0.09 (IMRT) and 1.077±0.06 (HT). HI at the head angle 30 degrees is 1.056±0.08 (CRT), 1.020±0.05 (IMRT) and 1.022±0.07 (HT).Conclusion:The results of our study show that when head angle tilted at 30 degree, target coverage, HI, CI were improved, and the dose delivered to OAR was reduced compared with conventional supine position in brain radiation therapy.This work was supported by the Radiation Technology R&amp;D program (No. 2013M2A2A7043498) and the Mid‐career Researcher Program (2014R1A2A1A10050270) through the National Research Foundation of Korea funded by the Ministry of Science, ICT&amp;Future Planning.

222. Optimization of pretreatment IMRT dose distribution verification: Bidimensional gaussian kernel implementation