|Year : 2011 | Volume
| Issue : 2 | Page : 85-94
A hybrid algorithm for instant optimization of beam weights in anatomy-based intensity modulated radiotherapy: A performance evaluation study
Ranganathan Vaitheeswaran1, VK Sathiya Narayanan2, Janhavi R Bhangle2, Amit Nirhali2, Namita Kumar2, Sumit Basu2, Vikram Maiya2
1 Siemens Ltd., HealthCare Sector, Pune, India
2 Department of Radiation Oncology, Ruby Hall Clinic, Pune, India
|Date of Submission||16-Jul-2010|
|Date of Decision||14-Aug-2010|
|Date of Acceptance||09-Sep-2010|
|Date of Web Publication||18-Apr-2011|
Healthcare Sector, Siemens Ltd., 403A, ICC Trade Tower, Tower B, Senapati Bapat Road, Pune - 411 016
Source of Support: None, Conflict of Interest: None
| Abstract|| |
The study aims to introduce a hybrid optimization algorithm for anatomy-based intensity modulated radiotherapy (AB-IMRT). Our proposal is that by integrating an exact optimization algorithm with a heuristic optimization algorithm, the advantages of both the algorithms can be combined, which will lead to an efficient global optimizer solving the problem at a very fast rate. Our hybrid approach combines Gaussian elimination algorithm (exact optimizer) with fast simulated annealing algorithm (a heuristic global optimizer) for the optimization of beam weights in AB-IMRT. The algorithm has been implemented using MATLAB software. The optimization efficiency of the hybrid algorithm is clarified by (i) analysis of the numerical characteristics of the algorithm and (ii) analysis of the clinical capabilities of the algorithm. The numerical and clinical characteristics of the hybrid algorithm are compared with Gaussian elimination method (GEM) and fast simulated annealing (FSA). The numerical characteristics include convergence, consistency, number of iterations and overall optimization speed, which were analyzed for the respective cases of 8 patients. The clinical capabilities of the hybrid algorithm are demonstrated in cases of (a) prostate and (b) brain. The analyses reveal that (i) the convergence speed of the hybrid algorithm is approximately three times higher than that of FSA algorithm; (ii) the convergence (percentage reduction in the cost function) in hybrid algorithm is about 20% improved as compared to that in GEM algorithm; (iii) the hybrid algorithm is capable of producing relatively better treatment plans in terms of Conformity Index (CI) [~ 2% - 5% improvement] and Homogeneity Index (HI) [~ 4% - 10% improvement] as compared to GEM and FSA algorithms; (iv) the sparing of organs at risk in hybrid algorithm-based plans is better than that in GEM-based plans and comparable to that in FSA-based plans; and (v) the beam weights resulting from the hybrid algorithm are about 20% smoother than those obtained in GEM and FSA algorithms. In summary, the study demonstrates that hybrid algorithms can be effectively used for fast optimization of beam weights in AB-IMRT.
Keywords: Anatomy-based IMRT, hybrid algorithm, intensity modulated radiotherapy, optimization, fast simulated annealing
|How to cite this article:|
Vaitheeswaran R, Sathiya Narayanan V K, Bhangle JR, Nirhali A, Kumar N, Basu S, Maiya V. A hybrid algorithm for instant optimization of beam weights in anatomy-based intensity modulated radiotherapy: A performance evaluation study. J Med Phys 2011;36:85-94
|How to cite this URL:|
Vaitheeswaran R, Sathiya Narayanan V K, Bhangle JR, Nirhali A, Kumar N, Basu S, Maiya V. A hybrid algorithm for instant optimization of beam weights in anatomy-based intensity modulated radiotherapy: A performance evaluation study. J Med Phys [serial online] 2011 [cited 2019 Jun 19];36:85-94. Available from: http://www.jmp.org.in/text.asp?2011/36/2/85/79693
| Introduction|| |
Recently, there has been a growing interest in aperture-based inverse planning (ABIP) for IMRT, as ABIP can significantly reduce the number of segments and monitor units. , This is accomplished without loss of dose coverage for the targets and with sparing of nearby critical structures. Also, IMRT plans with pre-defined anatomy-based MLC fields, known as anatomy-based IMRT (AB-IMRT), could be considered to reduce both the treatment complexity and verification burden. ,, The optimization of the beam weights in AB-IMRT was addressed by many investigators using different methods. ,,,, In general,the heuristic methods such as simulated annealing (SA) and genetic algorithms (GAs) are capable of escaping local optima and thus able to arrive at a global optimum. 
The simulated annealing method simulates the slow cooling of a sample to find low-energy states. This technique has been applied to problems in radiotherapy, especially in IMRT. ,,, Recently several enhancements of simulated annealing method have been developed, such as parallel tempering approach.  In general, the method of simulated annealing can provide well-acceptable results in IMRT optimization as compared to any other optimization algorithms, mainly due to its ability to escape from the local optima.  However, if time is a critical factor, simulated annealing method may deliver suboptimal solutions as it employs a random search technique. 
On the other hand, several very efficient exact optimization algorithms have been developed in recent years. ,, These algorithms can now be applied to some problems of IMRT as the system sizes which can be treated are now much larger than those being treated 10 years ago. The advantage of using such exact optimization algorithms is that they take very less time as compared to iterative and heuristic algorithms. However, applying of such non-iterative methods may produce suboptimal solutions in some situations, like in those where they can get trapped into the possible local minima.
In this work, a simple and efficient optimization algorithm for AB-IMRT, called "hybrid algorithm," is introduced in response to the drawbacks mentioned above. Our proposal is that by integrating an exact optimization algorithm with a heuristic optimization algorithm, the advantages of both the algorithms can be integrated into the created hybrid algorithm, which will lead to an efficient global optimizer solving the problem at a very fast rate. Our hybrid approach combines Gaussian elimination method (GEM) (exact optimizer) and fast simulated annealing (FSA) algorithm (a heuristic global optimizer) for the optimization of beam weights in AB-IMRT. The numerical and the clinical characteristics of the hybrid algorithm are compared with those of the GEM and FSA algorithms. In the numerical analysis, the numerical characteristics of the hybrid algorithm, such as convergence, convergence rate, consistency, number of trails and overall optimization speed, were analyzed for 8 patients in comparison with those of the GEM and FSA algorithms. The clinical capabilities of the hybrid algorithm were demonstrated in (i) prostate and (ii) brain cases.
| Material and Methods|| |
Anatomy-based intensity modulated radiotherapy
We used a simple anatomy-based segmentation method ,, for manually generated anatomy-based MLC fields in AB-IMRT. More details on how the anatomy-based fields are generated can be seen in our recent publication.  [Figure 1] shows an example of a set of anatomy-based MLC fields for a particular beam angle.
|Figure 1: a: An example of a set of anatomy-based MLC fields for a particular beam angle. In the give beam angle, blocking of spinal cord present within the BEV of the target volume results in the subsequent fields as shown in Figure 1a. Also few beams directly passing through spinal cord avoiding the rest of the target volume is used in AB-IMRT plans as shown in Figure 1b in order to produce the desired dose distribution|
Figure 1: b: An example of a set of anatomy-based MLC fields for a particular beam angle. In the give beam angle, blocking of spinal cord present within the BEV of the target volume results in the subsequent fields as shown in Figure 1a. Also few beams directly passing through spinal cord avoiding the rest of the target volume is used in AB-IMRT plans as shown in Figure 1b in order to produce the desired dose distribution
Click here to view
The proposed hybrid approach combines GEM algorithm (an exact approach) and FSA algorithm (a heuristic approach) for the optimization of beam weights in AB-IMRT using a quadratic dose-based cost function. In linear algebra, GEM is a powerful algorithm that can be used to determine the exact solutions of a system of linear equations. ,,, Recently, the use of Gaussian elimination algorithm for optimizing beam weights in AB-IMRT has been demonstrated.  In our sequential optimization approach, the initial approximate solutions (beam weights) are obtained using GEM. These initial solutions are, in turn, fed into FSA algorithm for further optimization. We used MATLAB software package, which incorporates GEM and SA algorithms in its optimization tool box. In the simulated annealing (SA) optimization module, we used a fast-cooling scheme to speed up the annealing process (FSA scheme). The processes involved in the hybrid algorithm are given in the flow chart [Figure 2].
|Figure 2: A flow chart of the proposed hybrid optimization algorithm for anatomy-based inverse planning|
Click here to view
Patient contours are first generated in the CMS-XiO® (4.3.1) treatment planning system, and the same system is used for dose calculation as well. In this planning system, the dose is calculated using a fast convolution superposition algorithm.  A dose grid size of 3 × 3 × 3 mm was used throughout the study. The slice thickness of the CT images used for planning purpose was 3 mm.
We have tested the hybrid algorithm in patient treatments. To date, many patients have been successfully delivered AB-IMRT plans optimized using the hybrid algorithm. The AB-IMRT treatment is delivered using Siemens ONCOR Impression Plus Linear Accelerator. This linac has a facility called automatic field sequencing (AFS) with which it is possible to deliver non-IMRT treatments involving several beams and beam segments very quickly without the operator interventions. This AFS facility makes the AB-IMRT treatment delivery very fast and efficient.
This analysis will demonstrate how the mathematical properties of the hybrid algorithm could be exploited in the optimization of single-criterion functions for AB-IMRT. To perform the analysis, different data sets (A, B, C, D, E, F, G and H) were generated that belonged to different patient cases (HandN 1, HandN 2, Brain 1, Brain 2, Abdomen 1, Abdomen 2, Pelvic 1 and Pelvic 2, respectively), representing a typical AB-IMRT planning situation. Each patient data set is represented using a quadratic dose-based cost function. The patient data sets used for the numerical analysis comprised the following:
The number of gantry angles, number of apertures per beam angle and their shapes were adapted to the anatomy of the given case. [Table 1] gives the summary of the AB-IMRT plans for these cases used in the numerical analysis. Then, an attempt was made to minimize the cost functions by optimizing beam weights for each data set. The optimizations for the above-mentioned patient data sets were performed using (a) hybrid algorithm (FSA+GEM), (b) GEM and (c) FSA in order to understand the numerical abilities of the hybrid approach.
- Beam parameters (gantry angles, number of apertures and their shapes)
- Optimization parameters (dose constraints and penalties)
- User-defined dose-control points (voxel samples)
|Table 1: Summary of AB-IMRT plan details for the cases used in the numerical analysis|
Click here to view
In this section, we have presented a detailed account of the clinical performance of the hybrid algorithm done for 2 patient cases (prostate and brain) in comparison with GEM and FSA algorithms. The patient cases are chosen on the basis that there is a considerable geometric and dosimetric complexity involved in the planning and three dimensional conformal radiotherapy (3D CRT) is apparently not capable of producing the desired dose distribution in the cases taken for the study. In order to simplify the process of optimization, we have systematically sampled a number of dose control points in the regions of planning target volume (PTV), normal tissues and in the surrounding regions (to control spillage) only for which the dose will be calculated during optimization. A sampling density (r) of 1 point/cm 3 (approximately) was used in both the patient cases.Also,a differential tissue penalty scheme was usedin both the cases in order to prioritize the goals. The PTV coverage was given the highest penalty, and the dose to the OARs was given relatively lower penalties. Moreover, the plan quality obtained using the three different algorithms were analyzed in terms of Conformity Index and Homogeneity Index. Here, the Conformity Index (CI)  is defined as
CI = 1 + V n /V t
where V n is the volume of normal tissue receiving the prescribed dose, and V t is the volume of the target receiving the prescribed dose.
The Homogeneity Index (HI)  is defined as
HI = [D max− D min ]/D prescribed
where D min (dose to 2% of the PTV), D max (dose to 98% of the PTV) and D prescribed are the minimum, maximum and prescribed doses, respectively.
| Results|| |
The numerical analysis shows [Figure 3], [Figure 4] and [Figure 5] that the convergence speed increases significantly for the hybrid algorithm as compared to the fast simulated algorithm due to the inclusion of Gaussian elimination method. Also, the consistency of the hybrid algorithm is not affected by the inclusion of Gaussian elimination method algorithm with the FSA algorithm. Most importantly, the number of iterations required to optimize the AB-IMRT problem is dramatically reduced for hybrid algorithm as compared to FSA algorithm.
|Figure 3: A plot of the final cost function for the patient data sets A, B, C, D, E, F, G and H in GEM, FSA and hybrid algorithms|
Click here to view
|Figure 4: A plot of the number of iterations taken for hybrid and FSA algorithms in patient data sets A, B, C, D, E, F, G and H|
Click here to view
|Figure 5: A plot of the time taken for optimization in hybrid and FSA algorithms for patient data sets A, B, C, D, E, F, G and H|
Click here to view
The convergence is about 20% improved in hybrid algorithm as compared to that in GEM algorithm and is comparable to that in FSA algorithm in most of the cases. In our study, the convergence was measured by the percentage reduction in the cost function for a particular case. In other words, the more the percentage reduction in the cost function, the better is the convergence of the optimization algorithm. But the percentage reduction in the cost function value does not necessarily mean an equal amount of improvement in the required outcomes. However, the dose distributions obtained in the 8 patient cases used for the numerical analysis indicate that the plan with lower final cost function value appears to have apparently better dose distribution as compared to the plan with higher final cost function value. This feature is common for all the 8 patients cases included in the numerical analysis and 2 patient cases included in the clinical analysis. This observation confirms that in a single-criterion optimization, the reduction in the cost function can be used as an approximate indication of the corresponding improvement in the dose distribution. The numerical analysis also points out that an algorithm having relatively good convergence characteristics will obviously lead to better dose distribution.
The hybrid algorithm was used to generate an AB-IMRT plan for a prostate case. In this case, we planned a dose of 57.6 Gy for the gross disease in a single phase comprising 37 fractions. The gross disease was delineated using CT images and was defined as clinical target volume (CTV). The planning target volume (PTV) was drawn with a 3-mm margin to the CTV. The volume of PTV was 160 cc. The geometry of the PTV was very complex as it was overlapping on the nearby rectum and bladder volumes. Moreover, we wanted to restrict the dose to bladder and rectum as much as possible. Because of the geometric and dosimetric complexities, we considered to execute AB-IMRT plan instead of 3DCRT plan for this case.
The summary of the treatment goals for this case is given in [Table 2]. Six 6-MV beams were used with 4 apertures per beam, resulting in a total of 32 beam segments. The OARs included in this case were the rectum, bladder and the two femoral heads. In order to compare the performance of the hybrid algorithm with that of GEM and FSA algorithms, respectively, we planned the same case using GEM and FSA algorithms. All plans were normalized based on the dose volume histograms (DVHs) such that 95% of the target volume was covered by the prescribed dose.
|Table 2: Summary of treatment goals for the AB-IMRT plans for prostate and brain cases|
Click here to view
[Table 3] gives the overall summary of the results obtained using hybrid, GEM and FSA algorithms.The axial dose distribution and DVHs obtained using hybrid, GEM and FSA algorithms are shown in [Figure 6]a and b, [Figure 7]a and b. The reduction in cost function in terms of number of iterations for hybrid and FSA algorithms is shown in [Figure 8].The Conformity Index (CI) for hybrid, GEM and FSA algorithms is 1.44, 1.49 and 1.47, respectively. The CI values are corresponding to the 95% dose coverage. The Homogeneity Index (HI) for hybrid, GEM and FSA algorithms is 0.25, 0.28 and 0.26, respectively. [Figure 9] shows the comparison of beam weights obtained using hybrid, GEM and FSA algorithms.
|Figure 6: A comparison of final dose distributions in an axial slice for the prostate case obtained in (a) hybrid algorithm, (b) GEM algorithm and (c) FSA algorithm. The thick red line shows the planning target volume (PTV)|
Click here to view
|Figure 7: A DVH comparison for the prostate case between (a) hybrid and GEM algorithms, the solid lines denoting GEM-based plan and the dotted lines denoting hybrid algorithm-based plan; and (b) hybrid and FSA algorithms, the solid lines denoting FSA algorithm-based plan and dotted lines denoting hybrid algorithm-based plan|
Click here to view
|Figure 8: A plot of the reduction in cost function with the number of iterations for the prostate case in hybrid and FSA algorithms|
Click here to view
|Figure 9: A comparison of beam weights obtained for the prostate case using hybrid, GEM and FSA algorithms|
Click here to view
|Table 3: Summary of dose-volume indices obtained for prostate and brain cases in hybrid, gaussian elimination method and fast simulated annealing algorithms|
Click here to view
The case is a typical brain lesion of volume 308 cc (PTV). A total dose of 50.4 Gy in 25 fractions was prescribed for the PTV in this case. The sensitive normal structures such as brainstem, optic nerves and eyes were close to the target volume. Especially, we wanted to restrict the dose to the portion of the brainstem volume which was not overlapping on the target volume (rest of the brainstem), while maintaining good dose coverage to the PTV. Hence it was decided to go for AB-IMRT plan instead of 3DCRT. The summary of the treatment goals for this case is given in [Table 2]. Seven 6-MV beams were used with 3 apertures per beam, resulting in a total of 29 beam segments. The direct exposure to right and left eyes was avoided in the plan. The OARs included were the brainstem, right and left eye, and optic nerves. All plans were normalized based on the DVHs such that 95% of the target volume was covered by the prescribed dose.
[Table 3] gives the overall summary of the results obtained using hybrid, GEM and FSA algorithms.The axial dose distribution and DVHs obtained using hybrid, GEM and FSA algorithms are shown in [Figure 10]a and b, [Figure 11]a and b. The reduction in cost function in terms of number of iterations for hybrid and FSA algorithms is shown in [Figure 12].The Conformity Index (CI) for hybrid, GEM and FSA algorithms is 1.37, 1.43 and 1.41, respectively. The Homogeneity Index (HI) for hybrid, GEM and FSA algorithms is 0.38, 0.41 and 0.41, respectively. [Figure 13] shows the comparison of beam weights obtained using hybrid, GEM and FSA algorithms.
|Figure 10: A comparison of final dose distributions in an axial slice for the brain case obtained in a) hybrid algorithm, b) GEM algorithm and c) FSA algorithm. The thick red line shows the planning target volume (PTV)|
Click here to view
|Figure 11: A DVH comparison for the brain case between (a) hybrid and GEM algorithms, the solid lines denoting GEM-based plan and the dotted lines denoting hybrid algorithm-based plan; and (b) hybrid and FSA algorithms, the solid lines denoting FSA-based plan and the dotted lines denoting hybrid algorithm-based plan|
Click here to view
|Figure 12: A plot of the reduction in cost function with the number of iterations for the brain case in hybrid and FSA algorithms|
Click here to view
|Figure 13: A comparison of beam weights obtained for the brain case using hybrid, GEM and FSA algorithms|
Click here to view
| Discussion|| |
First, the comparison of the numerical capabilities of the hybrid algorithm with those of the GEM and FSA algorithms as shown in [Figure 3], [Figure 4] and [Figure 5] gives a clear indication that the proposed strategy (hybrid algorithm) gives a better result in terms of convergence (as compared to GEM) and convergence rate (as compared to FSA). The optimization using the hybrid algorithm is almost three times faster than that obtained using the FSA algorithm.
The dose distribution and DVH comparisons demonstrate that the plan obtained using hybrid algorithm offers better PTV dose conformity and dose homogeneity as compared to the plans obtained using GEM and FSA algorithms. On an average, one can observe about 2% to 5% improvement in dose conformity and 4% to 10% improvement in dose homogeneity in plans optimized using hybrid algorithm as compared to GEM- and FSA-based plans. From [Figure 7] and [Figure 11] and [Table 3], one can observe that the OAR-sparing is improved in the hybrid algorithm-based plans as compared to GEM-based plans. However, the plans obtained using FSA algorithm offer better OAR-sparing as compared to both GEM-based and hybrid algorithm-based plans. An impressive advantage with hybrid algorithm is the huge reduction in the number of iterations as compared to FSA algorithm [Figure 14]a, as a consequence of which the optimization speed is considerably improved with hybrid algorithm [Figure 14]b in prostate and brain cases.
|Figure 14: A graphical representation of the comparison of (a) the number of iterations taken for hybrid and FSA algorithms for prostate and brain cases, (b) the time taken for hybrid and FSA algorithms for prostate and brain cases|
Click here to view
The plan quality obtained with hybrid algorithm is improved as compared to that obtained with GEM-based algorithm, which is due to two reasons: First, the GEM component in the hybrid algorithm gets the initial set of solutions, which when fed into the FSA component drives the heuristic process towards the goal in an effective way. Secondly, the simulated annealing component in the hybrid algorithm escapes the solution from the possible local minima during the optimization. The general observation is that the dose conformation to the tumor is very good in both the cases; and furthermore, the dose gradient in the region proximal to the OAR is steep, ensuring a good OAR protection.
Moreover,[Figure 9] and [Figure 13] indicate that the hybrid algorithm is able to generate relatively smooth beam weights for the anatomy segments as compared to GEM and FSA algorithms. The smoothness of a set of beam weights obtained using an algorithm was measured using the standard deviation (SD) value of the beam weights. The larger the SD value, the greater is the variation or fluctuation in the beam weights. The SD values for hybrid, GEM and FSA algorithms were 14.1, 21.6 and 16.5, respectively, for the prostate case. Similarly, for the brain case, the SD values were 9.7, 18.6 and 12, respectively, for hybrid, GEM and FSA algorithms. Therefore, on an average, the beam weights resulting from the hybrid algorithm are about 20% smoother than those obtained from GEM and FSA algorithms. It is well known that smooth beam weights can translate into an efficient treatment delivery. ,,, Henceit is always desirable to achieve smooth beam weights. However there is no considerable difference in the total monitor units (MUs) obtained in plans created using GEM, FSA and hybrid algorithms as shown in [Figure 15].
|Figure 15: A comparison of monitor units in prostate and brain cases with respect to the optimization algorithm|
Click here to view
There is an extensivedebate about whether the cost functions really possess local minima,and whether they are sufficiently close to the global minimum so that the finding of global minimum becomes insignificant.  However, with the available techniques to handle local minima, along withmodern fast computers, they can be invoked even in the absenceof conclusive proof of the existence of local minima.  The simulatedannealing technique generally guarantees global optimum in the final outcome, however with a compromise of optimization speed, which diminishes the real importance of such algorithms in a clinic. In this given situation, a methodology that helps speeding up the optimization process without degrading the quality of final solutions will be very useful in a clinic.
| Conclusion|| |
A hybrid optimization algorithm for anatomy-based IMRT (AB-IMRT) is introduced, which integrates an exact solver (GEM) with a global optimizer (FSA) in order to get better solutions at faster rate. For the cases presented, the implemented hybrid optimizer was able to produce treatment plans comparable to FSA-based plans in terms of target coverage and OAR-sparing with a remarkable improvement in the optimization speed (about three times faster than FSA algorithm). We believe that such an improvement in optimization speed will lead to highly efficient workflow in AB-IMRT optimization, which in turn can be helpful to produce better treatment plans.
| Acknowledgment|| |
The authors would like to thank Siemens Ltd., India, for providing access to MATLAB software, without which this research work would not have been possible.
| References|| |
|1.||De Gersem W, Claus F, De Wagter C, W De Neve. An anatomybased segmentation tool for intensity-modulated radiation therapy and its application to head-and-neck cancer. Int J Radiat Oncol Biol Phys 2001;51:849-59. |
|2.||Shepard D, Earl M, Li X, Naqvi S, Yu C. Direct aperture optimization: A turnkey solution for step-and-shoot IMRT. Med Phys 2002;29:1007-18. |
|3.||Mageras G, Mohan R. Application of fast simulated annealing to optimization of conformal radiation treatments. Med Phys 1993;20:639-47. |
|4.||Xiao Y, Galvin J, Hossain M, Valicenti R. An optimized forward planning technique for intensity modulated radiation therapy. Med Phys 2000;27:2093-9. |
|5.||Zhang G, Jiang Z, Shepard D, Zhang B, Yu C. Direct aperture optimization of breast IMRT and the dosimetric impact of respiration motion. Phys Med Biol 2006;51:357-69. |
|6.||Bednarz G , Michalski D, Houser C, Saiful Huq M, Xiao Y, Anne PR, et al. The use of mixed-integer programming for inverse treatment planning with pre-defined field segments. Phys Med Biol 2002;47:2235-45. |
|7.||Bortfeld T, Urkelbach J, Boesecke R, Schlegel W. Methods of image reconstruction from projections applied to conformation therapy. Phys Med Biol 1990;35:1423-34. |
|8.||Cotrutz C, Xing L. Segment-based dose optimization using a genetic algorithm. Phys Med Biol 2003;48:2987-98. |
|9.||Yaohang L, Protopopescub VA, Nikita A, Xinyu Z, Andrey G. Hybrid parallel tempering and simulated annealing method. Appl Math Comput 2009;212:216-28. |
|10.||Rieger H. Application of exact combinatorial optimization algorithms to the physics of disordered systems. Comput Phys Commun 2002;147:702-6. |
|11.||Dhaenens C, Lemesrea J, Talbi EG. K-PPM: A new exact method to solve multi-objective combinatorial optimization problems. Eur J Oper Res 2010;200:45-53. |
|12.||Dumitrescu I, Stützle T. Usage of exact algorithms to enhance stochastic local search algorithms. Ann Inf Syst 2009;10:103-34. |
|13.||Ranganathan V, Sathiya Narayanan VK, Bhangle JR, Gupta KK, Basu S, Maiya V, et al. Performance evaluation of an algorithm for fast optimization of beam weights in anatomy-based IMRT. J Med Phys 2010;35:104-12. |
|14.||Isaacson E, Keller HB. Analysis of Numerical Methods. New York: Willey Pub.; 1966. p. 220. |
|15.||Westlake JR.A Handbook of Numerical Matrix Inversion and Solution of Linear Equations. New York: Willey Pub.; 1968. p. 312. |
|16.||Chawla M. A parallel Gaussian elimination method for general linear systems. Int J Comput Math 1992;42:1029-265. |
|17.||Grimes RG. Solving Systems of Large Dense Linear Equations. J Supercomput 1988;1:291-9. |
|18.|| Bedford JL. Speed versus accuracy in a fast convolution photon dose calculation for conformal radiotherapy. Phys Med Biol 2002;47:3475-84. |
|19.||Wiggenraad J, Petoukhova L, Versluis L, Santvoort PC. Stereotactic Radiotherapy of intracranial tumors: A comparison of intensity-modulated radiotherapy and dynamic conformal arc.Int J Radiation Oncology Biol Phys 2009;74:1018-26. |
|20.||Wu Q, Mohan R, Morris M, Lauve A, Schmidt-Ullrich R. Simultaneous integrated boost intensity-modulated radiotherapy for locally advanced head-and-neck squamous cell carcinomas. I: Dosimetric results. Int J Radiat Oncol Biol Phys 2003;56:573-85. |
|21.||Mohan R, Arnfield M, Tong S, Wu Q, Siebers J. The impact of fluctuations in intensity patterns on the number of monitor units and the quality and accuracy of intensity modulated radiotherapy. Med Phys 2000;27:1226-37. |
|22.||Spirou SV, Fournier-Bidoz N, Yang J, Chui CS, Ling CC. Smoothing intensity modulated beam profiles to improve the efficiency of delivery. Med Phys 2001;28:2105-12. |
|23.||Ma L. Smoothing intensity modulated treatment delivery under hardware constraints. Med Phys 2002;29:2937-45. |
|24.||Webb S. A simple method to control aspects of fluence modulation in IMRT planning. Phys Med Biol 2001;46:187-95. |
|25.||Llacer J, Deasy JO, Bortfeld TR, Solberg TD, Promberger C. Absence of multiple local minima in intensity modulated optimization with dose-volume constraints. Phys Med Biol 2003;48:183-210. |
|26.||Webb S. The physical basis of IMRT and inverse planning.Br J Radiol 2003;76:678-89. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10], [Figure 11], [Figure 12], [Figure 13], [Figure 14], [Figure 15]
[Table 1], [Table 2], [Table 3]
|This article has been cited by|
||An experimental investigation on the effect of beam angle optimization on the reduction of beam numbers in IMRT of head and neck tumors
| || Sathiya Narayanan, V.K., Vaitheeswaran, R., Bhangle, J.R., Basu, S., Maiya, V., Zade, B. |
| ||Journal of Applied Clinical Medical Physics. 2012; 13(4): 36-43 |