
INVITED PAPER 



Year : 2009  Volume
: 34
 Issue : 3  Page : 129132 

Optimization of beam angles for intensity modulated radiation therapy treatment planning using genetic algorithm on a distributed computing platform
Daryl P Nazareth, Stephen Brunner, Matthew D Jones, Harish K Malhotra, Mohammad Bakhtiari
Department of Radiation Medicine, Roswell Park Cancer Institute, Elm & Carlton Sts, Buffalo NY 14263, USA
Date of Submission  30Dec2008 
Date of Decision  13Mar2009 
Date of Acceptance  21Apr2009 
Date of Web Publication  3Aug2009 
Correspondence Address: Daryl P Nazareth Elm & Carlton Sts. Buffalo NY 14263 USA
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.54845
Abstract   
Planning intensity modulated radiation therapy (IMRT) treatment involves selection of several angle parameters as well as specification of structures and constraints employed in the optimization process. Including these parameters in the combinatorial search space vastly increases the computational burden, and therefore the parameter selection is normally performed manually by a clinician, based on clinical experience. We have investigated the use of a genetic algorithm (GA) and distributedcomputing platform to optimize the gantry angle parameters and provide insight into additional structures, which may be necessary, in the dose optimization process to produce optimal IMRT treatment plans. For an IMRT prostate patient, we produced the first generation of 40 samples, each of five gantry angles, by selecting from a uniform random distribution, subject to certain adjacency and opposition constraints. Dose optimization was performed by distributing the 40plan workload over several machines running a commercial treatment planning system. A score was assigned to each resulting plan, based on how well it satisfied clinicallyrelevant constraints. The second generation of 40 samples was produced by combining the highestscoring samples using techniques of crossover and mutation. The process was repeated until the sixth generation, and the results compared with a clinical (equallyspaced) gantry angle configuration. In the sixth generation, 34 of the 40 GA samples achieved better scores than the clinical plan, with the best plan showing an improvement of 84%. Moreover, the resulting configuration of beam angles tended to cluster toward the patient's sides, indicating where the inclusion of additional structures in the dose optimization process may avoid dose hot spots. Additional parameter selection in IMRT leads to a largescale computational problem. We have demonstrated that the GA combined with a distributedcomputing platform can be applied to optimize gantry angle selection within a reasonable amount of time.
Keywords: Intensity modulated radiation therapy, distributed computing, genetic algorithm, optimization
How to cite this article: Nazareth DP, Brunner S, Jones MD, Malhotra HK, Bakhtiari M. Optimization of beam angles for intensity modulated radiation therapy treatment planning using genetic algorithm on a distributed computing platform. J Med Phys 2009;34:12932 
How to cite this URL: Nazareth DP, Brunner S, Jones MD, Malhotra HK, Bakhtiari M. Optimization of beam angles for intensity modulated radiation therapy treatment planning using genetic algorithm on a distributed computing platform. J Med Phys [serial online] 2009 [cited 2019 Mar 26];34:12932. Available from: http://www.jmp.org.in/text.asp?2009/34/3/129/54845 
Introduction   
Intensity Modulated Radiation Therapy (IMRT) is a complex technology designed to deliver preciselymodulated and conformal radiation dose to a target. IMRT is employed clinically for several treatment sites, including the prostate, the head and neck region, and the brain.
IMRT planning requires the solution to the beamangle selection (BAS) problem, ^{ [1],[2],[3],[4],[5],[6],[7],[8]} which involves the selection of 510 angles from 360 possible gantry angles, subject to certain spacing and opposition constraints. This type of optimization problem is termed an integer programming (IP) problem, and involves a large number of binary variables. For example, in selecting 510 angles out of a set of 360, there are approximately 9 x 10^{ 19} possible candidates. In addition, in many clinics, the rotation angles of the treatment couch are also considered, further increasing the number of variables and possible solutions. Since no commercial software package exists which solves the BAS problem, clinicians generally select the gantry and couch angles manually based on clinical experience.
A second IMRT subproblem is that of dose optimization (DO). ^{ [9],[10],[11],[12],[13],[14],[15],[16],[17],[18],[19]} This problem involves optimizing beamlet weights in order to define an intensity map for each radiation field. The constraints involved are usually expressed as dosevolume histogram (DVH) constraints on targets and organs at risk. As noted above, no existing commercial technique attempts to solve the BAS and DO problems simultaneously. Other investigators have employed genetic algorithm approaches ^{ [20],[21]} to the BAS problem, but have not coupled it with a standard commercial dose optimization routine. As described below, we use the GA as an "outer loop" of our complete optimization technique.
The major impediment to solving the BAS, simultaneously with standard IMRT optimization, is the prohibitive size of the solution search space. However, the advent of highperformance computing and computer clusters has now made such attempts feasible. Our work involves the development of distributedcomputing tools to solve simultaneously the BAS and DO problems. These tools would represent significant improvement in current clinical IMRT treatmentplanning implementation.
Methods   
B1. Genetic Algorithm
Genetic algorithm^{ [22],[23],[24],[25],[26],[27] } (GA) is an adaptive heuristic search technique based on the biological principles of evolution and natural selection. It employs an initial random population which effectively samples the search space globally, along with a propagation technique to concentrate the search effort in promising regions. In this work, we will employ the GA to guide the search for optimal beam angles while simultaneously generating IMRT treatment plans.
The GA approach includes four steps: (1) producing the first generation of beamangle sample sets, (2) determining the "most fit" samples by means of a scoring function (3) combining the most promising samples using a crossover technique (4) applying a mutation scheme to a selection of the new generation of samples.
The search space consists of all candidate sets of gantry beam angles, subject to the following constraints: angles which differ by less than 30° are not permitted, and angles within 30° of opposition are not permitted. These are typical constraints employed in clinical IMRT treatment planning. The number of beams will be selected by the user, based on clinical considerations, although further development in this work would permit this parameter to be included in the search. As an example, we consider the treatment plan for a prostate patient in which the number of beams, N is equal to five. Producing the initial population of k samples is accomplished by random number generation.
For each of the k samples, an IMRT treatment plan is generated using the Eclipse software (Varian Medical Systems, Palo Alto, CA), with clinicallyrelevant DVH constraints. Currently, each plan is generated manually, as in standard clinical treatment planning; an automated method using a distributedcomputing approach is described in the ''Discussion and Conclusions'' section. Each plan then receives a score according to how well it satisfies the constraints; this enables the samples to be ranked. The best sample is cloned or copied, unchanged, to the next generation. The remaining k1 positions of the next generation are assigned to offspring of the best samples using the following method. For each pair of samples in the current generation, the average score is calculated. The best k1 pairs are placed into the mating pool . For each pair in this pool, a random binary crossover mask is generated and used to determine the beams angles of that pair's offspring. Each offspring then has a probability p of undergoing mutation , or alteration of one randomlyselected beam to a new random value.
The process is repeated with the new generation of k samples. After a number of generations, s , have been produced, the process is halted, and the sample with the best score is considered the solution to the BAS problem. In practice, the values of k , p , and s are selected based on the time and computational resources available, as discussed below. The GA process is illustrated in [Figure 1].
Evaluation
We analyzed an IMRT prostate case recently treated in our facility. The beam angles employed in the actual treatment were 45°, 105°, 180°, 255°, 315°. The DVH constraints for this case were those typically employed in our clinic. We used the GA algorithm described in Section A to produce each generation of beamangle samples. The plans were generated with the Eclipse software running on four separate machines.
After each IMRT plan was generated, it was normalized so that 95% of the PTV received the prescription dose of 81Gy, and was then evaluated using the following scoring function:
Where A_{ j} was the actual percentage volume receiving the constraint dose, C_{ j} was the percentage volume specified to receive the constraint dose, and w_{ i} was the relative weight assigned to the organ at risk (OAR). The sum over i included the bladder and rectum, and the sum over j included the constraints for each OAR. For this preliminary study, we used weights for the bladder and rectum of 1 and 2, respectively. Clearly, a lower score indicates a better treatment plan. The score for the clinical plan was 121.
We produced six generations of the GA, and in each case recorded the number of plans better than the clinical plan. We produced DVH figures to compare the best plan with the clinical plan.
Results   
In [Figure 2], we present the DVH's for the prostate, bladder, and rectum for the clinical and best GA plans. Clearly, the GA plan provides better sparing of the critical structures. [Figure 3] presents, for each generation of samples, the average and best plan scores.
Discussion and Conclusions   
The GA is an effective method of searching the space of gantry angles. We have shown that it can produce IMRT prostate plans with significantly better fitness scores than typical clinical plans. The GA implementation for the BASDO problem gains power when IMRT plans corresponding to the beamangle samples, for each iteration (see B1), are generated simultaneously. This requires distributing the samples over multiple IMRT workstations, or computing nodes , which is the subject of ongoing investigations. We have formed a collaboration with faculty members at the Center for Computational Research (CCR).^{ [28],[29]} Access to the vast computational resources of the CCR will permit us to implement the GA algorithm to run simultaneously on k machines, with k in the range of 25100 or more. The result is that the completion of each GA generation will require only approximately 30 minutes. This would enable 510 generations ( s =510) in the time currently required to produce one complete IMRT treatment plan.
References   
1.  Bortfeld T, Schlegel W. Optimization of beam orientations in radiation therapy: Some theoretical considerations. Phys Med Biol 1993;38:291304. 
2.  Das SK, Marks LB. Selection of coplanar or noncoplanar beams using threedimensional optimization based on maximum beam separation and minimized nontarget irradiation. Int J Radiat Oncol Biol Phys 1997;38:64355. 
3.  Djajaputra D, Wu Q, Wu Y, Mohan R. Algorithm and performance of a clinical IMRT beamangle optimization system. Phys Med Biol 2003;48:3191212. 
4.  Pugachev A, Li JG, Boyer AL, Hancock SL, Le QT, Donaldson SS, et al. Role of beam orientation optimization in intensitymodulated radiation therapy. Int J Radiat Oncol Biol Phys 2001;50:55160. 
5.  Rowbottom CG, Nutting CM, Webb S. Beamorientation optimization of intensitymodulated radiotherapy: clinical application to parotid gland tumours. Radiother Oncol 2001;59:16977. 
6.  Sailer SL, Rosenman JG, Symon JR, Cullip TJ, Chaney EL. The tetrad and hexad: Maximum beam separation as a starting point for noncoplanar 3D treatment planning: prostate cancer as a test case. Int J Radiat Oncol Biol Phys 1994;30:43946. 
7.  Sφderstrφm S, Brahme A. Which is the most suitable number of photon beam portals in coplanar radiation therapy? Int J Radiat Oncol Biol Phys 1995;33:1519. 
8.  D'Souza WD, Zhang HH, Nazareth DP, Shi L, Meyer RR. A nested partitions framework for beam angle optimization in intensitymodulated radiation therapy. Phys Med Biol 2008;53:3293307. 
9.  Chang SX, Cullip TJ, Rosenman JG, Halvorsen PH, Tepper JE. Dose optimization via indexdose gradient minimization. Med Phys 2002;29:113046. 
10.  Ezzell GA. Genetic and geometric optimization of threedimensional radiation therapy treatment planning. Med Phys 1996;23:293305. 
11.  Langer M, Morrill S, Brown R, Lee O, Lane R. A comparison of mixed integer programming and fast simulated annealing for optimizing beam weights in radiation therapy. Med Phys 1996;23:95764. 
12.  Lee EK, Fox T, Crocker I. Optimization of radiosurgery treatment planning via mixed integer programming. Med Phys 2000;27:9951004. 
13.  Llacer J. Inverse radiation treatment planning using the Dynamically Penalized Likelihood method. Med Phys 1997;24:175164. 
14.  Mageras GS, Mohan R. Application of fast simulated annealing to optimization of conformal radiation treatments. Med Phys 1993;20:63947. 
15.  Morrill SM, Lam KS, Lane RG, Langer M, Rosen II. Very fast simulated reannealing in radiation therapy treatment plan optimization. Int J Radiat Oncol Biol Phys 1995;31:17988. 
16.  Rosen II, Lane RG, Morrill SM, Belli JA. Treatment plan optimization using linear programming. Med Phys 1991;18:14152. 
17.  Spirou SV, Chui CS. A gradient inverse planning algorithm with dosevolume constraints. Med Phys 1998;25:32133. 
18.  Wu X, Zhu Y. A mixedencoding genetic algorithm with beam constraint for conformal radiotherapy treatment planning. Med Phys 2000;27:250816. 
19.  Langer M, Thai V, Papiez L. Improved leaf sequencing reduces segments or monitor units needed to deliver IMRT using multileaf collimators. Med Phys 2001;28:24508. 
20.  Hou Q, Wang J, Chen Y, Galvin JM. "Beam Orientation Optimization for IMRT by a Hybrid Method of the Genetic Algorithm and the Simulated Dynamics." Med Phys 2003;30:2360. 
21.  Lei J, Li Y. "An Approaching Genetic Algorithm for Automatic Beam Angle Selection in IMRT Planning." Comput Methods Programs Biomed 2009;93:25765. 
22.  Herrera F., Lozano, M.. Gradual distributed realcoded genetic algorithms. IEEE Transactions on Evolutionary Computation 2000;4:4363. 
23.  Voigt HM, Anheyer T. Modal mutations in evolutionary algorithms . Proceedings of the First IEEE Conference on Evolutionary Computation 1994;1:8892. 
24.  Chalermdamrichai V. A hybrid computerintelligent, userinteractive, process plan optimizer for fouraxis CNC turning centers. Industrial Engineering. Vol PhD. Madison: University of WisconsinMadison; 1999. 
25.  L. Shi SO, and N. Sun. Parallel algorithm for the traveling salesman problem. Computer and Operations Research 1999;26:371394. 
26.  L. Shi SOaQC. A new optimization framework for product design. Management Science 2001;47:168192. 
27.  Men LSaS. Optimal buffer allocation in production lines. IIE Transactions 2003;35:110. 
28.  M.D. Jones RY, C.P. Bhole. Hybrid MPIOpenMP programming for parallel OSEM PET reconstruction. IEEE Transactions on nuclear science 2006;53:27528. 
29.  Yao MDJaR. Parallel programming for OSEM reconstruction with MPT, OpenMP, and hybrid MPIOpenMP. IEEE Nuclear Science Syposium. Vol 5; 2004. p. 303642. 
[Figure 1], [Figure 2], [Figure 3]
This article has been cited by  1 
Integral dose investigation of noncoplanar treatment beam geometries in radiotherapy 

 Dan Nguyen,Peng Dong,Troy Long,Dan Ruan,Daniel A. Low,Edwin Romeijn,Ke Sheng   Medical Physics. 2014; 41(1): 011905   [Pubmed]  [DOI]   2 
A unified mixedinteger programming model for simultaneous fluence weight and aperture optimization in VMAT, Tomotherapy, and Cyberknife 

 Kerem Akartunali,Vicky MakHau,Thu Tran   Computers & Operations Research. 2014;   [Pubmed]  [DOI]   3 
Beam orientation in stereotactic radiosurgery using an artificial neural network 

 Agnieszka Skrobala,Julian Malicki   Radiotherapy and Oncology. 2014; 111(2): 296   [Pubmed]  [DOI]   4 
A surrogatebased metaheuristic global search method for beam angle selection in radiation treatment planning 

 H H Zhang,S Gao,W Chen,L Shi,W D D’Souza,R R Meyer   Physics in Medicine and Biology. 2013; 58(6): 1933   [Pubmed]  [DOI]   5 
A surrogatebased metaheuristic global search method for beam angle selection in radiation treatment planning 

 Zhang, H.H. and Gao, S. and Chen, W. and Shi, L. and DæSouza, W.D. and Meyer, R.R.   Physics in Medicine and Biology. 2013; 58(6): 19331946   [Pubmed]   6 
Comparison of beam angle selection strategies for intracranial IMRT 

 Bangert, M. and Ziegenhein, P. and Oelfke, U.   Medical Physics. 2013; 40(1)   [Pubmed]   7 
Characterizing the combinatorial beam angle selection problem 

 Bangert, M. and Ziegenhein, P. and Oelfke, U.   Physics in Medicine and Biology. 2012; 57(20): 67076723   [Pubmed]   8 
Uncertainty incorporated beam angle optimization for IMPT treatment planning 

 Cao, W. and Lim, G.J. and Lee, A. and Li, Y. and Liu, W. and Ronald Zhu, X. and Zhang, X.   Medical Physics. 2012; 39(8): 52485256   [Pubmed]   9 
Uncertainty incorporated beam angle optimization for IMPT treatment planning 

 Wenhua Cao,Gino J. Lim,Andrew Lee,Yupeng Li,Wei Liu,X. Ronald Zhu,Xiaodong Zhang   Medical Physics. 2012; 39(8): 5248   [Pubmed]  [DOI]   10 
Characterizing the combinatorial beam angle selection problem 

 Mark Bangert,Peter Ziegenhein,Uwe Oelfke   Physics in Medicine and Biology. 2012; 57(20): 6707   [Pubmed]  [DOI]   11 
Introducing multiple treatment planbased comparison to investigate the performance of gantry angle optimisation (GAO) in IMRT for head and neck cancer 

 Maria Thor, Hunor Benedek, Tommy Knöös, Per Engström, Claus F. Behrens, Anna Karlsson Hauer, David Sjöström, Crister Ceberg   Acta Oncologica. 2012; : 1   [VIEW]  [DOI]   12 
A twophase method for selecting IMRT treatment beam angles: BranchandPrune and local neighborhood search 

 Gino J. Lim, Wenhua Cao   European Journal of Operational Research. 2011;   [VIEW]  [DOI]   13 
A selfadaptive casebased reasoning system for dose planning in prostate cancer radiotherapy 

 Nishikant Mishra, Sanja Petrovic, Santhanam Sundar   Medical Physics. 2011; 38(12): 6528   [VIEW]  [DOI]   14 
Computational enhancements to fluence map optimization for total marrow irradiation using IMRT 

 D.M. Aleman, V.V. Mišić, M.B. Sharpe   Computers & Operations Research. 2011;   [VIEW]  [DOI]   15 
An algorithm for fast beam angle selection in intensity modulated radiotherapy 

 R. Vaitheeswaran, V. K. Sathiya Narayanan, Janhavi R. Bhangle, Amit Nirhali, Namitha Kumar, Sumit Basu, Vikram Maiya   Medical Physics. 2010; 37(12): 6443   [VIEW]  [DOI]  





