Abstract. An algorithm is described, based on ray-tracing and the beam's-eye-view, that exhaustively searches all permitted beam directions. The evaluation of ...
The British Journal of Radiology, 78 (2005), 316–327 DOI: 10.1259/bjr/13015047
E
2005 The British Institute of Radiology
Automatic selection of non-coplanar beam directions for three-dimensional conformal radiotherapy 1,2
J MEYER, PhD, 1S M HUMMEL, RT(T), CMD, 1P S CHO, PhD, 1M M AUSTIN-SEYMOUR, MD and M H PHILLIPS, PhD
1 1
University of Washington Medical Center, Department of Radiation Oncology, Box 356043, Seattle, WA 98195, USA and 2Universita¨t Wu¨rzburg, Klinik fu¨r Strahlentherapie, Josef-Schneider-Str. 11, D-97080 Wu¨rzburg, Germany
Abstract. An algorithm is described, based on ray-tracing and the beam’s-eye-view, that exhaustively searches all permitted beam directions. The evaluation of the search is based on a general cost function that can be adapted to the clinical objectives by means of parameters and weighting factors. The approach takes into account the constraints of the linear accelerator by discarding beam directions that are not permitted. A sensitivity analysis was carried out to determine appropriate parameters for different sized organs, and a prostate case was used to benchmark the approach. The algorithm was also applied to two clinical cases (brain and sinus) to test the benefits of the approach compared with manual angle selection. The time to perform a beam direction search was approximately 2 min for the coplanar and 12 min for the non-coplanar beam space. The angles obtained for the prostate case compared well with reports in the literature. For the brain case, the mean dose to the right and left optic nerves was reduced by 12% and 50%, respectively, whilst the target dose uniformity was improved. For the sinus case, the mean doses to the right and left parotid glands were reduced by 54% and 46%, respectively, to the right and left optic nerves by 37% and 62%, respectively, and to the optic chiasm by 39%, whilst the target dose uniformity was also improved. For the clinical cases the plans based on optimized beam directions were simpler and resulted in better sparing of critical structures compared with plans based on manual angle selection. The approach provides a practical alternative to elaborate and time consuming beam angle optimization schemes and is suitable for routine clinical usage.
The maximum radiation dose that can be deposited in a tumour is limited by the tolerance dose of the surrounding vital organs and critical structures [1]. Inherent objectives for a three-dimensional conformal radiation therapy (3DCRT) treatment plan are to use beam directions such that the irradiation of the organs at risk (OAR) is minimized so as to reduce normal tissue complications, and to ensure an adequate dose coverage of the planning target volume (PTV). This work concentrates on conventional 3D-CRT since it still constitutes most clinical treatments despite the increasingly widespread availability of intensity-modulated radiation therapy (IMRT). It was suggested by Webb [2] that only about 30% of clinical problems require the use of IMRT. The impetus of this work is to bridge the gap between 2D treatment planning and computer optimization of IMRT, since the efficacious use of 3D-CRT relies strongly on the judicious selection of non-coplanar beam angles. Rosen estimated in an industrial report in 1999 that only about one-third of facilities performed 3D-CRT regularly [3], but this number has probably increased significantly since then. The use of non-coplanar 3D-CRT beams, as opposed to coplanar arrangements, has the potential to produce superior plans due to its increased ability to avoid critical structures [4–6]. The beam angle orientation problem has been the subject of a number of studies. The different Received 16 April 2004 and in revised form 9 September 2004, accepted 29 November 2004. This work was partially supported by a grant from Elekta Oncology Systems.
316
approaches are based on geometric/volumetric methods [7–13], dose-based objectives [13–16], radiobiological objectives [5, 13] or IMRT specific objectives [17–19]. The geometric and volumetric methods for 3D-CRT are based on the assumption that the avoidance of the OAR is fundamental to the planning process. Further steps, such as optimization of beam weights [20, 21], beam portals [22, 23] or beam intensity modulation [2, 24], can be carried out independently during the next planning stage to achieve target dose homogeneity and to satisfy the dose constraints on the OAR. The work presented in this paper is also based on this principle. The approaches based on dose and radiobiological objectives as well as IMRT-specific algorithms, are more complex and some of them integrate the beam angle search into the actual plan optimization process [18, 19]. While these approaches have the most potential to find the optimal plan, they are computationally expensive since it is necessary to perform dose calculations during the optimization process. Furthermore, the standard deterministic optimization algorithms are not sufficient to prevent the algorithm becoming trapped in local minima [14, 25– 27]. The computational effort required to obtain optimal solutions, e.g. ranging from hours [18] up to days [6], makes some of these approaches prohibitive for routine clinical use, even when sophisticated methods are applied that reduce the computation time by an order of magnitude [28]. For coplanar plans, and only one or two OARs, it is fairly straightforward for an experienced dosimetrist to judge which beam directions are preferable. However, it is not feasible for a human planner to exploit the entire beam The British Journal of Radiology, April 2005
Beam direction search tool
space when multiple OARs of different size and radiosensitivity are present. Another difficulty stems from the inability to deliver isocentric beams from any direction owing to limitations relating to the dimensions of the treatment couch, linac design, patient dimensions and the location and size of the tumour. To represent accurately the permitted directions for each individual treatment requires precise modelling of the linac, the couch and the patient. This task is complex since the permitted beam space is a function of the relative location of the tumour [29–31]. Several attempts to automate the process of selecting beam directions have been described in the literature. Chen et al [8] and Myrianthopoulos et al [9] used the beam’seye-view (BEV) to obtain 2D plots of the percentage volume of OAR irradiated versus the gantry angle as an aid to select beam directions. Rowbottom et al [13] went a step further by combining the BEV volumetric approach with a simple dose model and biological considerations to obtain a single objective score for each beam direction. Pugachev and Xing [17] formulated an IMRT specific score function that was based on the maximum dose deliverable through a portal by means of an intensity modulated beam. The current work is similar to the work of the above authors. The BEV was used as a means to search for beam directions [8, 9], the cost function was introduced to obtain a score for each beam direction [13] and an exhaustive search was conducted to explore the non-coplanar beam space [17]. An exhaustive search circumvents the local minima problem and ensures that the results are reproducible independently of the starting conditions. In addition, the current work takes into account the couch-gantry-patient dimensions, evaluates the optimization by means of logical rules and uses a unique cost function for which optimal parameters were derived.
Methods and materials A beam direction search tool (BDST) was developed and implemented on the University of Washington treatment planning system (TPS) PRISM [32]. This tool evaluates a cost function, Ci for each beam direction, Wi(j,y), in the domain i51...t, where t is the total number of beam directions considered, and where each beam direction is defined by the gantry angle, j, and the couch angle, y. The BDST starts with both the gantry angle and couch angle set to zero, and then systematically works through all other permitted directions by increasing first the gantry angle and then the couch angle by multiples of a chosen angular increment, v. Beam directions that would cause the gantry to collide, either with the couch or with the patient, are excluded from the search. The beam space for the Elekta SL20 series linac (Elekta, Crawley, UK) was modelled for different treatment sites in combination with the Elekta PreciseTM couch and the female RandoTM phantom (The Phantom Laboratory, Salem, NY), by assuming a constant isocentre distance of 100 cm. An example of the beam space for tumours in the head is shown in Figure 1, where WP(j,y) denotes the permitted directions. Note that PRISM uses different angle conventions to the Elekta accelerator and therefore two axes are plotted in Figure 1. Each beam portal is shaped to match the projection of the PTV into the BEV and so only the rays represented by the solid circles in Figure 2a are considered. For each The British Journal of Radiology, April 2005
Figure 1. The area that is unshaded denotes the permitted beam space WP(j,y). The shaded region represents angles that result in a gantry collision and are, therefore, excluded from the beam search.
beam direction the BDST first computes distances along ray paths as illustrated in Figure 2b. For ray, j, these distances are the depth from the patient entry surface to the PTV, dj(PTV), the depth from the patient entry surface to each OAR in the beam, dj(OARk), where k is an index for the OAR in question, and the path length through the OAR, lj(OARk). By definition the central ray is the one that passes through the isocentre. The rays are projected onto a Cartesian grid located at a distance of 120 cm from the source (Figure 2) such that they are separated by equal increments, s, at the isocentre plane. The BDST uses the dj(PTV) values to determine the mean depth to the PTV entry surface for the ith beam direction, di(PTV), and similarly calculates a mean depth to each OARk entry surface, di(OARk). It also estimates the fractional volume of each OARk irradiated, vi(OARk), using the ray path lengths, ray separations s and a precalculation of the total OARk volume. The cost function was designed such that the lower its value the more preferable is the beam direction and it has two components that are added together. The first component, Ci(PTV), is based on the mean depth of the PTV, and the second, Ci(OAR), is a penalty term, which depends on the amount of OAR irradiation, thus Ci ~Ci (PTV)zCi (OAR)
ð1Þ
The first component term is computed using the equation Ci (PTV)~1 exp(kdi (PTV))
ð2Þ
where m is the linear attenuation coefficient in water for the proposed beam energy. The exponential term, with m set to a value of 0.0495 cm21 throughout this paper, corresponds to the attenuation of a 6 MV photon beam in water. Thus, the value of Ci(PTV) approaches zero when the PTV is at shallow depths and approaches unity for very large depths, so favouring those beam directions where the PTV is close to the beam entry surface. As several different OARs may be included in any particular beam direction, the second component is a weighted average of individual costs computed for each OAR in turn, where different weighting factors, wk, may 317
J Meyer, S M Hummel, P S Cho et al
Figure 2. Ray-tracing geometry: (a) beam’s-eye-view (BEV) illustrating the grid with uniform spacing s and (b) axial slice with planning target volume (PTV) and two organs at risk (OAR).
be assigned to give each different OAR a different relative importance. Each OAR has both a ‘‘depth cost’’ and a ‘‘volume cost’’. The depth cost for the kth OAR is calculated using the equation: Ci (di (OARk )~exp(kdi (OARk ))
ð3Þ
The same value of m is used as in Equation 1, but, unlike Equation 1, this gives higher costs (approaching unity) for the shallower depths where the relative dose to the OAR will be greater. If the OAR is completely outside the beam edge in the BEV the depth cost is set to zero. The volume cost for each OAR is set equal to the fractional volume irradiated, vi(OARk), which is already conveniently in the range zero to one. In combining the volume and depth costs a factor l (where 0,l,1) was introduced so as to control the relative importance of depth and volume, and to ensure that the total cost remains within the range 0 to 1; thus the cost associated with irradiating the kth OAR in the ith beam direction is quantified by: Ci (OARk )~jvi (OARk )z(1 j)Ci (di (OARk ))
ð4Þ
The angle optimization tool was designed so that the user specified the number of beam directions, nb, required for a plan and the minimum angle separation, Q, between selected angles. To further ensure beam diversity at least one of all nb beams had to be separated by an angle d, with 90 ˚¡d¡120 ˚ from the other beams that were selected. An empirical value of d590 ˚ was assumed throughout the work presented in this paper to ensure some dose uniformity in the PTV. 318
The optimization started by discarding all the beam directions that were not permitted within the specified gantry and couch ranges. The algorithm then iterated through the remaining angle combinations of j and y and created a rank list with the beam directions sorted according to the lowest cost. The direction that topped the list was selected as the first beam. The second direction was selected subject to the criterion that the absolute angle between the first and the second direction was greater or equal to Q. The value of Q depends on the total number of beams and the location of the tumour. However, experience has shown that a minimum separation of 40– 45 ˚ usually results in satisfactory plans. If the angle was not large enough, the algorithm proceeded to the next position in the rank list until it found an appropriate direction. After the second direction was selected, the algorithm searched for the third direction in the rank list that was separated by a minimum angle Q from any of the previously selected beam directions. The same procedure was carried out for all nb beam directions. For the selection of the last beam direction, the algorithm also examined whether any of the previously selected beam directions and the current candidate direction were separated by at least d, otherwise the angle was discarded and the algorithm proceeded to the next best direction until one meeting the criterion was found.
Parameter value selection Ray-tracing resolution The trade-off between computation time and a fine spatial resolution was explored by varying the value of The British Journal of Radiology, April 2005
Beam direction search tool
Figure 3. Test plan. (a) Axial slice and (b) lateral view (from left side).
s. Since the selection of s depends on the size of the organ, a simple test plan, shown in Figure 3, was produced with OAR of different realistic sizes, i.e. Ø(OAR1)510 mm, Ø(OAR2)520 mm, Ø(OAR3)535 mm and Ø(OAR4)550 mm. The axial length of each OAR was equal to its diameter, as illustrated in Figure 3b, and the diameter of the PTV was 25 mm. The diameter of the external contour was 120 mm. For the purpose of this sensitivity analysis, y was set to zero (reducing the search space to coplanar beam arrangements) the gantry range was set to js{0 ˚2360 ˚}, the angle resolution v510 ˚, l50.2, nb55, Q545 ˚ and the weighting factors were set to unity. The values of s investigated were 1 mm, 2 mm, 3.3 mm, 5 mm and 10 mm. The variation of the total cost (Equation 1) as a function of the ray-trace resolution was analysed. All the computation for this work was carried out on a PC running Linux with a CPU speed of 2.4 GHz.
Angle resolution Beam directions were sampled at uniform angle intervals v, as illustrated in Figure 3. For the search to be effective it is essential to keep v as small as possible. If the resolution is too coarse, potentially good beam directions could be missed. To determine an appropriate value for v, the beam space for the test plan in Figure 3 was searched for different values of v, namely 2 ˚, 5 ˚, 10 ˚, 15 ˚ and 30 ˚. The grid increment, s, was set to 3.3 mm and the other parameters remained as described in the previous paragraph.
Clinical examples The BDST was first applied to a prostate case and the resulting beam arrangement compared with clinically The British Journal of Radiology, April 2005
approved beam directions from the literature. The tool was then tested on two clinical cases with tumours in the head.
Prostate comparison The prostate is well suited for an initial comparison since the anatomy is comparable from patient to patient and coplanar beam arrangements are standard. In our institution, the gantry angles for standard 3D conformal prostate treatments are based on the recommendations by Pickett et al [33]. It is a coplanar 6-field conformal technique with opposing lateral beams and opposed oblique beams angled between 30 ˚ and 45 ˚ off the lateral beam position. The lateral beams penetrate through the femoral heads and the lateral oblique beams are placed such that they minimize the irradiation of the bladder, the rectum and the femoral heads. The dose to the femoral heads can be controlled by the relative beam weights of the opposed lateral beams. For the first run, the BDST was applied to an arbitrarily selected prostate case with the parameters for the search set to js{0 ˚2360 ˚}, ys{0 ˚}, s53.3 mm, v51 ˚, Q540 ˚, l50.2 and nb56. Only the rectum and the bladder were considered OAR and the corresponding weighting factors, wk, were set to 10 and 5, respectively. For the second run nb55 and v was increased to 5 ˚. In addition, the femoral heads were considered critical organs. The weighting factors applied were 10 for the rectum, 1 for the bladder and 2 for both femoral heads. These parameters were chosen following the findings by Rowbottom et al [13] who utilized a cost function with volumetric and biological elements to optimize the beam directions for prostate treatment with similar weighting factors. 319
J Meyer, S M Hummel, P S Cho et al
Clinical case studies Tumours in the head are especially suited for noncoplanar beam arrangements since numerous critical structures of different importance are to be spared and the permitted beam space is much larger since beams directed through the top of the head are an option. Two cases were studied to test the algorithm. Case I was a patient with a brain tumour and Case II was a patient with a sinus tumour. For each, an experienced senior dosimetrist had produced a 3D-CRT plan, which was used to treat the patient. These will be referred to as the conventional plans. The dosimetrist was confined to non-IMRT beams but could use any means available for 3D-CRT clinical routine planning. This included the use of non-coplanar beams, different energies, optimization of the MLC aperture, customized blocks, electron beams and wedges. These plans were compared, retrospectively, with plans based on beam directions obtained from the beam direction search, and these will be referred to as the optimized plans. The same dosimetrist produced the optimized plans using the same treatment planning tools as previously. Uniform margins of 5 mm for the brain case and 7 mm for the sinus case were applied, to the clinical target volume (CTV) in order to obtain the PTVs and a voxel size of 3.3 mm was used to calculate the final dose distribution. An additional margin of 5 mm was applied to the PTV for the purpose of the beam direction search. This was done to account for the fact that lateral electron transport takes place when the beams are delivered, and without this additional margin, beam directions for which the projection of an OAR is very close to the PTV, but does not overlap with it, might appear to be good. However, when the radiation is delivered the electron transport would result in a high dose to the OAR. The optimized plans were required to use an equal number of beams as the conventional plans. The beam weights and wedge angles were optimized manually for the first case study (brain) and re-optimized with simulated annealing for the second case study (sinus). The plans were compared by means of dose–volume histograms (DVH) and statistical evaluation based on the recommendations for reporting target volume and dose specifications for IMRT plans [34], which is based on ICRU50 [1]. Further evaluation makes use of the Webb and Nahum model [35] for tumour control probability (TCP), the LymanKutcher-Burman model [36, 37] for normal tissue complications probability (NTCP) and the radiation conformity index [38] to the 95% isodose surface (RCI95). The parameters for the NTCP model were based on data published by Emami et al [39].
Case study I: brain A 12-year-old paediatric patient with a cavernous meningioma was selected. The PTV was located inbetween the temporal lobes and overlapped approximately 25% of the right optic nerve (Figure 4). The patient had an acute onset of blindness in his right eye before radiotherapy. The left optic nerve was in close proximity to the PTV. The total dose prescribed to the PTV was 50.4 Gy. Since the tumour was impinging on the optic nerves it was important to reduce the dose to these structures, despite a tolerance dose of approximately 50 Gy, so as to further 320
Figure 4. Brain case: axial CT slice through a paediatric meningioma.
minimize the impact on the patient’s eyesight. This was considered important given the young age of the patient. The conventional plan consisted of three wedged beams. Two of the beams were opposed lateral coplanar 6 MV beams (W1(90 ˚,0 ˚) and W2(270 ˚,0 ˚)) and the third was a left superior oblique 18 MV beam (W3(100 ˚,300 ˚)). For the optimized plan the following parameters were applied for the beam direction search: js{0 ˚2360 ˚}, ys{0 ˚– 100 ˚,260 ˚–360 ˚}, s50.33, v510 ˚, Q545 ˚, l50.2, nb53 and the weighting factors wk for the OAR were set to 7, 5, 10 and 8, corresponding to the right and left eye, and the right and left optic nerve, respectively. A higher weighting was given to the OAR on the right lateral side compared with the corresponding counterparts on the left because of the overall location of the tumour and the proximity to the right optic nerve. The optic nerves were considered the most important structures, followed by the eyeballs. To allow a fair comparison, both plans were scaled to D95 so that the 95% isodose surface covered the entire PTV.
Case study II: sinus The second clinical case investigated was a large T4N0M0 sinus adenocarcinoma. A coronal slice through the centre of the tumour is shown in Figure 5. The volume of the PTV was approximately 128 cm3 and was close to the brain stem and the parotid glands. Superiorly, the PTV protruded between the eyeballs, which made it difficult to spare the optic nerves and the optic chiasm while at the same time delivering the prescribed dose of 70 Gy to the PTV. It was aimed for tolerance doses of 50 Gy for the optic nerves and the optic chiasm, 52 Gy for the brain stem and 10 Gy to the lens to minimize the risk of cataracts. A secondary priority was to prevent xerostomia by reducing the dose to the parotid glands. The British Journal of Radiology, April 2005
Beam direction search tool
The following parameters were used for the beam direction search for the optimized plan: js{0 ˚2360 ˚}, ys{0 ˚–100 ˚,260 ˚–360 ˚}, s53.3 mm, v510 ˚, Q545 ˚, nb56 and the weighting factors wk for the OAR were set to 7, 7, 5, 8, 9, 9, 10, 10, 10, 1 and 1 corresponding to the right and left eye, the spinal cord, the brain stem, the right and left optic nerve, the optic chiasm, the right and left lens and the right and left parotid gland, respectively.
Results Parameter value selection The results of the beam direction search for the test plan in Figure 3 are given as a function of the ray-trace resolution in Table 1 and as a function of the angular resolution in Table 2. With respect to both tables, the time to perform the search is given in column two and the top five angles in the rank list, and the selected angles are given in the third and fourth columns, respectively. Each individual cost term in Equation 1 was analysed for all of the test cases. Only the total cost is plotted as a function of the ray-trace resolution, s, in Figure 6a and as a function of the angle resolution, v, in Figure 6b.
Prostate comparison The time to perform the first run, for which the femoral heads were not considered OAR, was 112 s. The optimal gantry angles were found to be: 273 ˚, 96 ˚, 56 ˚, 313 ˚, 233 ˚ and 136 ˚ corresponding to the ranks #1, #3, #149, #158, #163 and #173. The first two angles were virtually opposed lateral fields rotated marginally, mainly to minimize the irradiation of the rectum, but also to reduce irradiation of the bladder. The rest of the beams were angled off the lateral beams by ¡40 ˚, which corresponds to the minimum angle separation Q. The second run lasted 24 s and included the femoral heads as critical structures. The following gantry angles were selected by the BDST: 300 ˚, 65 ˚, 115 ˚, 255 ˚ and 340 ˚ corresponding to ranks #1, #6, #7, #10 and #45, respectively. Figure 7 shows a comparison with the cost
Figure 5. Sinus case: coronal CT slice through tumour.
The conventional plan made use of six beams consisting of two anterior oblique 6 MV beams (W1(335 ˚,0 ˚) and W2(20˚,0 ˚)), a superior oblique 18 MV beam (W3(90˚,310 ˚)), a left lateral 6 MV beam (W4(90 ˚,0 ˚)), a right lateral 10 MV beam (W5(270 ˚,0 ˚)) and an anterior 6 MV beam (W6(0 ˚,0 ˚)). All but the anterior beam made use of a wedge and utilized customized beam blocks.
Table 1. The results of the beam direction search as a function of the ray-trace resolution s Ray-trace resolution, s (cm)
Computation time (s)
Top 5 ranked angles, j (degrees)
Selected angles, j (degrees)
0.1 0.25 0.33 0.5 1.0
67 11 7 3 1
330, 320, 320, 320, 320,
330, 320, 320, 320, 320,
240, 140, 260, 260, 260,
150, 330, 140, 140, 140,
250, 150, 330, 330, 330,
340 340 240 240 340
240, 140, 260, 260, 260,
150, 240, 140, 140, 140,
60, 180 170, 350 230, 350 230, 350 350, 230
Table 2. The results of the beam direction search as a function of the angle resolution v Angle resolution, v (degrees)
Computation time (s)
Top 5 ranked angles, j (degrees)
Selected angles, j (degrees)
2 5 10 15 30
34 14 7 5 2
320, 320, 320, 255, 330,
320, 320, 320, 255, 330,
The British Journal of Radiology, April 2005
318, 255, 260, 330, 240,
252, 260, 140, 240, 150,
254, 256 245, 140 330, 240 150, 315 60, 180
252, 255, 260, 330, 240,
140, 140, 140, 150, 150,
72, 350 350, 65 230, 350 225, 75 60, 180
321
J Meyer, S M Hummel, P S Cho et al
Figure 6. Cost as a function of (a) the ray-tracing resolution, s, and (b) as a function of the angle resolution, v, for the test case in Figure 3.
Figure 7. Cost function vs gantry-angle plot for prostate: (a) Graph reproduced from Rowbottom et al [13]. (b) Normalized cost corresponding to prostate case. The vertical lines indicate the selected gantry angles. Note that different cost functions were used and therefore the overall magnitude of the graphs is not comparable.
function vs gantry angle plot obtained by Rowbottom et al [13] (Figure 7a) and our results (Figure 7b).
Clinical case studies The times to perform the beam direction searches were 11 min for the brain case and 12 min for the sinus case. The normalized cost distributions over the permitted beam space are depicted in Figure 8a for the brain case and in Figure 8b for the sinus case. The beam space that was not permitted is depicted with a cost of unity (white). The numbers next to the diamonds denote the overall rankings. For the optimized treatment plan wedges were applied to all three beams for the brain case but no wedges were necessary for the sinus case. The beam energy for all optimized beams was 6 MV. DVHs for the brain case are shown in Figure 9, and the corresponding statistical evaluation is given in Table 3. 322
The graphs corresponding to the critical structures on the left side of the patient are depicted in black and those corresponding to the right side are in grey, and the conventional plan is plotted with dashed lines; whereas the optimized plan is plotted with solid lines. Both plans are within the dose specifications and do not violate the constraints. For the sinus case, DVHs are shown in Figure 10 for the PTV, the brain stem, the parotid glands, the optic nerves and the optic chiasm. The dose to the spinal cord was negligibly low for both the conventional and optimal plans. As the plans could not be normalized to the D95 of 70 Gy, due to the high dose to the optic chiasm for the conventional plan, they were scaled such that the maximum doses in the PTV were comparable. The statistical evaluation, including additional statistical information for the lenses, is given in Table 4. The British Journal of Radiology, April 2005
Beam direction search tool
Figure 8. Normalized cost distribution over the beam space WP(j,y). A circle corresponds to the conventional and a diamond to the optimized plan. The grey-scale bar represents the cost values. (a) Brain case and (b) sinus case.
Figure 9. Brain case: dose–volume histograms for the conventional plan (dashed) and the optimized plan (solid): (a) planning target volume (PTV) and (b) optic nerves. The vertical lines in (a) denote the ¡5% dose range.
Table 3. Case I. Brain: statistical evaluation Structure
Plan
Dose [Gy]
Discussion in [%]
Parameter value selection
Min. Mean Max. D95 V100 NTCP RCI PTV
Conv. Optim. Right optic Conv. nerve Optim. Left optic Conv. nerve Optim.
45.0 47.1 22.5 7.5 20.8 3.1
52.5 51.4 44.4 39.0 36.7 18.3
55.0 53.8 53.4 51.6 52.5 51.2
50.4 50.4 n.a. n.a. n.a. n.a.
95 95 n.a. n.a. n.a. n.a.
n.a. n.a. 0.26 0.07 0.01 ,0.01
0.27 0.23 n.a. n.a. n.a. n.a.
PTV, planning target volume; Conv., conventional; Optim., optimized; D95, dose that covers 95% of the PTV; V100, percentage of the PTV receiving the prescribed dose; NTCP, normal tissue complication probability; RCI, radiation conformity index; n.a., not applicable.
The British Journal of Radiology, April 2005
The weighting factors for this set-up were all set to unity so there was no preference for one OAR over another, and hence the cost for the OAR only reflects the relative penalty for its location and the partial volume irradiated. Owing to the normalization, the cost terms are also independent of the size of the OAR. Assuming that the finest ray-trace resolution and the smallest value of v result in the most accurate results, the anticipated tendency can be seen. The higher the value of s, the more deviation can be observed. This was most marked for the smallest OAR with a diameter of 10 mm. In Table 1 it can be seen that 323
J Meyer, S M Hummel, P S Cho et al
Figure 10. Sinus case: dose–volume histogram for the conventional plan (dashed) and the optimized plan (solid): (a) planning target volume (PTV) and brain stem, (b) parotid glands, (c) optic nerves and (d) optic chiasm. The vertical lines in (a) denote the ¡5% dose range. Table 4. Case II. Sinus: statistical evaluation Structure
Plan
Dose [Gy]
in [%]
Min. Mean Max. D95 V100 TCP NTCP PTV
Conv. Optim. Brain stem Conv. Optim. Right parotid Conv. Optim. Left parotid Conv. Optim. Right optic Conv. nerve Optim. Left optic Conv. nerve Optim. Optic chiasm Conv. Optim. Right lens Conv. Optim. Left lens Conv. Optim.
42.7 58.0 10.7 3.7 2.5 3.5 0.1 3.4 9.6 6.6 19.5 6.0 27.0 11.8 7.6 4.4 13.4 4.5
71.5 72.9 18.7 29.5 42.4 19.7 37.9 20.5 27.8 17.5 42.8 16.2 47.2 29.0 8.3 6.6 16.5 9.8
76.6 77.0 55.7 52.2 57.7 32.4 74.6 65.0 52.3 25.6 63.9 31.3 64.6 49.8 9.1 10.3 19.0 15.7
62.2 68.0 n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a.
83 89 n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a.
73 94 n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a.
n.a. n.a. 0.5 0.5 40 ,0.1 31 0.3 1.5 ,0.1 12 ,0.1 1.4 ,0.1 2 2 2 2
PTV, planning target volume; Conv., Conventional; Optim., optimized; D95, dose that covers 95% of the PTV; V100, percentage of the PTV receiving the prescribed dose; NTCP, normal tissue complication probability; RCI, radiation conformity index; n.a., not applicable.
324
the resolution 3.3 mm, 5 mm and 10 mm resulted in the same four beam directions to top the rank list with the first three of the selected beams being the same. However, reducing the resolution to 1 mm and 2.5 mm, the algorithm found beam directions with a lower cost. The selection of the value for the angle resolution v for the exhaustive search determines how efficiently the beam space is scanned. A value of 10 ˚ for v was acceptable; larger intervals could be used for larger OAR (.20 mm). Based on these results, values of s53.3 mm and v510 ˚ were selected for the clinical case studies due to the small size of the optic nerves. A value of Q545 ˚ was chosen for the minimum angle separation for the clinical cases with non-coplanar beam arrangements. A minimum angle separation is required to avoid overlapping beams within the body, and so prevents an unnecessarily high dose to the tissue located between the beam entrance surface and the tumour. In terms of dose homogeneity it is also desirable to have the beams maximally separated in space [11, 40]. The accumulation of all of the beams from one general direction makes it difficult to achieve target dose uniformity and difficult to conform the high dose region to the shape of the tumour. The British Journal of Radiology, April 2005
Beam direction search tool
Prostate comparison To validate the approach a comparison was drawn with our own prostate technique (first run) and with the extensive work that has been carried out on prostate planning at the Institute of Cancer Research and Royal Marsden Hospital [13] (second run). For the latter, Figure 7 illustrates that in both cases the gantry angles were selected in the valleys of the graph such that irradiation of the femoral heads was minimized. Based on clinical experience, e.g. symmetry of the beams leading to dose uniformity, Rowbottom et al [13] decided also to select an anterior beam (0 ˚) although the cost was higher than for other candidate gantry angles. The comparison of the cost vs gantry-angle plot revealed that the current methodology can reproduce both a standard 6-field beam arrangement as well as individualized beam angles that correspond well to the findings by other researchers.
Clinical case studies For each of the clinical test cases, the non-coplanar beam directions were checked by placing the treatment mask on the treatment couch, setting up the isocentre and positioning the couch and the gantry in the acquired positions. None of the beam directions violated any geometrical constraints imposed by the linac.
Case study I: brain In Figure 8a the locations of the angles that were selected manually for the conventional plan are superimposed (circles) on the cost distribution of the beam direction search. The existence of multiple local minima can be observed in the cost distribution, indicated by the low cost regions (dark). The non-coplanar beam with gantry angle 100 ˚ and couch angle 300 ˚ lies within the low cost area and shows that the dosimetrist was successful in approaching the best non-coplanar beam direction. The dosimetrist might have chosen a slightly different angle for reasons not directly considered in the beam direction search such as dose uniformity in the PTV. The two standard lateral coplanar beams with gantry angles 90 ˚ and 270 ˚ are not optimal with respect to the cost function used for the beam direction search. Comparing the DVH for the brain case in Figure 9 shows that the optimized plan is superior to the conventional plan for both the critical structures and the PTV. The PTV is more uniformly irradiated and the cold volume considerably reduced. This is also reflected in the statistical evaluation in Table 3. The NTCP values indicate that both plans do not overdose any critical structures and that the probability of complications is very low. However, any improvements are worthwhile, especially for a paediatric patient. In addition, the improved sparing of OAR without compromising the dose to the PTV is an important foundation for dose escalation, which can only be achieved safely if the OARs are not exposed to toxic dose levels. Case study II: sinus The cost distribution in Figure 8b shows that the BDST finds two non-coplanar beam directions (ranks #28, #77) of which #77 is located on the edge of the permitted beam space. It is a posterior beam that penetrates through the The British Journal of Radiology, April 2005
brain stem. In fact, it is this beam that enables better coverage of the PTV since it complements the other beams. It is also most likely that this beam direction would have not been explored during manual angle selection. The fact that the best beam with rank #1 is almost identical with one beam selected manually at gantry angle 335 ˚ highlights that the dosimetrist again made a good selection. One might wonder why there is no beam positioned in the low cost region on the bottom right and the top left of Figure 8b. After the first beam with rank #1 was selected all the beams in close proximity with a low rank were discarded if they were within a certain minimum angle Q. Furthermore, some of the beam angles are redundant. A simple example would be a vertex field that can be delivered using either Wi(90 ˚,270 ˚) or Wi(270 ˚,90 ˚). The DVH in Figure 10 shows that dose uniformity in the PTV was improved in the optimized plan. More importantly, the cold spot was greatly reduced which is reflected in the increase in the TCP value from 73% to 94%. The mean dose to the brainstem is higher for the optimized plan, which is due to the posterior beam. However, the maximum dose is lower, which is reflected in the similar NTCP values for both plans. All the other statistical values in Table 4 are clearly in favour of the optimized plan with greatly reduced mean and maximum doses for the parotid glands, the optic nerves and the optic chiasm. The most important difference between the two plans is the much improved sparing of both parotid glands. It has been reported that a mean parotid dose of more than 26 Gy slows salivary flow [41] and doses over 40 Gy may result in total dysfunction of salivary production [42]. The mean dose for the right parotid gland was reduced from 42 Gy to 20 Gy (improvement of 54%) for the optimized plan, which resulted in a reduction of the NTCP from 40% to ,0.1%. For the left parotid gland, the improvement in the NTCP was approximately 46%. This indicates that both parotid glands would be virtually complication free. The sixth beam selected (rank #77) went straight through the brain stem but avoided all the other critical structures. This selection was due to the design of the cost function and the choice of the weighting factors, which all had a similar magnitude. The weighting factors were based on common sense, but were chosen somewhat arbitrarily insofar as the magnitudes, as well as the relative difference could have been selected differently. Selecting weighting factors that are similar in magnitude can be considered as a soft constraint, and this is necessary if multiple objectives, in this case multiple OAR, are present. If one weighting factor is much higher than all the others, then this could result in a solution where only the one particular OAR would be spared. It should also be remembered that the individual cost terms are normalized to the volume of the respective organ, requiring the entire organ volume to be contoured. The beam weight for the beam that went straight through the brain stem was set relatively low by the beam weight optimization procedure. This is an interesting observation as it paves the way for the application of the BDST for IMRT planning [43]. If this particular beam were replaced by an IMRT beam it would have probably received a much higher overall weighting since IMRT has the ability to spare critical structures by modulating the intensity across the field. 325
J Meyer, S M Hummel, P S Cho et al
The main advantage of the BDST is its ability to search the permitted beam space within a time frame that is acceptable for clinical practice and to select beam directions in a rule-based manner. The approach is inferior to beam angle optimization techniques that integrate dose calculations, optimize wedge angles, etc. Nevertheless the results of this work have shown that the BDST is able to suggest beam directions that result in plans with greatly improved sparing of critical structures and is, therefore, a good compromise to more elaborate techniques. In this work, all the suggestions in terms of selected beam directions by the BDST were accepted. In clinical routine, however, the dosimetrist might only select some of the beam directions suggested by the BDST and combine them with standard beam directions or beam directions that have a high cost but when the portal is split into two subfields (e.g. one on the left and one on the right side of the spinal cord), or if the intensity is modulated, result in preferable beam directions. The BDST has shown that small changes in the beam directions (see clinical case I) can already result in considerable improvements in the reduction of dose to the normal tissues. It is noteworthy that none of the beams for the clinical case II made use of wedges or different energies, which highlights that the conventional plan was more complex and simpler plans might be possible. It would be of further interest to evaluate the trade-offs between the speed and the quality of a plan based on the current approach and more sophisticated techniques that include dose calculations during the optimization procedure. The cost maps obtained could be used to predetermine and separate the ‘‘good’’ beam directions from the ‘‘bad’’ ones. This could then be used as an input for more complex stochastic optimization procedures similar to the approach described by Pugachev et al [28] who used prior knowledge to reduce the optimization time.
Conclusion A general algorithm has been presented that explores the entire permitted beam space to determine beam directions that minimize the irradiation of healthy tissues and critical structures. The algorithm can be adjusted to different clinical goals and is a valuable tool for 3D-CRT treatment planning. The beam directions obtained for the prostate case compared well with standard beam directions and other work described in the literature. For the brain and the sinus cases the optimized plans were simpler and better than conventional plans. It is anticipated that the plans could be further improved if the intensity of the conformal beams could be modulated. Indeed, the fact that one beam penetrated through the brain stem in order to spare the other OAR makes it an interesting starting point for further research on the use of the BDST for IMRT, especially with regard to forward planned IMRT [44].
Acknowledgments The author would also like to thank Dr Homayon Parsai for his assistance with LaTeX during the preparation of this manuscript. 326
References 1. ICRU50. Prescribing, Recording and Reporting Photon Beam Therapy. Bethesda, MD: International Commission on Radiation Units and Measurements, 1993. 2. Webb S. Intensity-modulated radiation therapy. Bristol and Philadelphia, USA: Institute of Physics Publishing, 2001. 3. Rosen LR, quote in ‘Software and hardware advances bring greater accuracy, reduced toxicity to radiation therapy’, Technology Trends, May 1999. 4. Derycke S, Van Duyse B, De Gersem W, De Wagter C, De Neve W. Non-coplanar beam intensity modulation allows large dose escalation in stage III lung cancer. Radiother Oncol 1997;45:253–61. 5. So¨derstro¨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:151–9. 6. Pugachev A, Li JG, Boyer AL, Hancock SL, Le QT, Donaldson SS, et al. Role of beam orientation optimisation in intensity-modulated radiation therapy. Int J Radiat Oncol Biol Phys 2001;50:551–60. 7. Haas OC, Burnham KJ, Mills JA. Optimisation of beam orientation in radiotherapy using planar geometry. Phys Med Biol 1998;43:2179–93. 8. Chen GT, Spelbring DR, Pelizzari CA, Balter JM, Myrianthopoulos LC, Vijayakumar S, et al. The use of beam’s eye view volumetrics in the selection of non-coplanar radiation portals. Int J Radiat Oncol Biol Phys 1992;23:153–63. 9. Myrianthopoulos LC, Chen G, Vijayakumar S, Halpern H, Spelbring R, Pelizzari CA. Beam’s eye view volumetrics: an aid in rapid treatment plan development and evaluation. Int J Radiat Oncol Biol Phys 1992;23:367–75. 10. Ezzell GA. Genetic and geometric optimisation of threedimensional radiation therapy treatment planning. Med Phys 1996;23:293–305. 11. Das SK, Marks LB. Selection of coplanar or noncoplanar beams using three-dimensional optimisation based on maximum beam separation and minimized nontarget irradiation. Int J Radiat Oncol Biol Phys 1997;38:643–55. 12. Rowbottom CG, Webb S, Oldham M. Beam-orientation customisation using an artificial neural network. Phys Med Biol 1999;44:383–99. 13. Rowbottom CG, Webb S, Oldham M. Improvements in prostate radiotherapy from the customization of beam directions. Med Phys 1998;25:1171–9. 14. Bortfeld T, Schlegel W. Optimisation of beam orientations in radiation therapy: some theoretical considerations. Phys Med Biol 1993;38:291–304. 15. Rowbottom CG, Oldham M, Webb S. Constrained customization of non-coplanar beam orientations in radiotherapy of brain tumours. Phys Med Biol 1999;44:2251–62. 16. Pugachev AB, Boyer AL, Xing L. Beam orientation optimisation in intensity-modulated radiation treatment planning. Med Phys 2000;27:1238–45. 17. Pugachev AB, Xing L. Pseudo beam’s-eye-view as applied to beam orientation selection in intensity-modulated radiation therapy. Int J Radiat Oncol Biol Phys 2001;51:1361–70. 18. Djajaputra D, Wu Q, Wu Y, Mohan R. Algorithm and performance of a clinical IMRT beam-angle optimisation system. Phys Med Biol 2003;48:3191–212. 19. Li Y, Yao J, Yao D. Automatic beam angle selection in IMRT planning using genetic algorithm. Phys Med Biol 2004;49:1915–32. 20. Mohan R, Mageras GS, Baldwin B, Brewster LJ, Kutcher GJ, Leibel S, et al. Clinically relevant optimisation of 3-D conformal treatments. Med Phys 1992;19:933–44. 21. Wu X, Zhu Y. An optimisation method for importance factors and beam weights based on genetic algorithms for radiotherapy treatment planning. Phys Med Biol 2001;46: 1085–99.
The British Journal of Radiology, April 2005
Beam direction search tool 22. van Dalen S, Keijzer M, Huizenga H, Storchi PR. Optimisation of multileaf collimator settings for radiotherapy treatment planning. Phys Med Biol 2000;45:3615–25. 23. De Gersem W, Claus F, De Wagter C, Van Duyse B, De Neve W. Leaf position optimisation for step-and-shoot IMRT. Int J Radiat Oncol Biol Phys 2001;51:1371–88. 24. Bednarz G, Michalski D, Houser C, Huq MS, 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. 25. So¨derstro¨m S, Gustafsson A, Brahme A. The clinical value of different treatment objectives and degrees of freedom in radiation therapy optimisation. Radiother Oncol 1993;29: 148–63. 26. Stein J, Mohan R, Wang XH, Bortfeld T, Wu Q, Preiser K, et al. Number and orientations of beams in intensitymodulated radiation treatments. Med Phys 1997;24:149–60. 27. Alber M, Nuesslin F. An objective function for radiation treatment optimisation based on local biological measures. Phys Med Biol 1999;44:479–93. 28. Pugachev A, Xing L. Incorporating prior knowledge into beam orientation optimisation in IMRT. Int J Radiat Oncol Biol Phys 2002;54:1565–74. 29. Kessler ML, McShan DL, Fraass BA. A computer-controlled conformal radiotherapy system. III: graphical simulation and monitoring of treatment delivery. Int J Radiat Oncol Biol Phys 1995;33:1173–80. 30. Meyer J, Mills JA, Haas OC, Burnham KJ, Parvin EM. Accommodation of couch constraints for coplanar intensity modulated radiation therapy. Radiother Oncol 2001;61: 23–32. 31. Furhang EE, Hanley J, Chiu-Tsao ST, Toner S, Fan P, Gliedman P, et al. Clearance assurance for stereotactic radiosurgery and radiotherapy. Med Phys 2002;29:45–50. 32. Kalet IJ, Jacky JP, Austin-Seymour MM, Unger JM. Prism: a new approach to radiotherapy planning software. Int J Radiat Oncol Biol Phys 1996;36:451–61. 33. Pickett B, Roach 3rd M, Horine P, Verhey L, Phillips TL. Optimisation of the oblique angles in the treatment of prostate cancer during six-field conformal radiotherapy. Med Dosim 1994;19:237–54.
The British Journal of Radiology, April 2005
34. Intensity Modulated Radiation Therapy Collaborative Working Group. Intensity-modulated radiotherapy: current status and issues of interest. Int J Radiat Oncol Biol Phys 2001;51:880–914. 35. Webb S, Nahum AE. A model for calculating tumour control probability in radiotherapy including the effects of inhomogeneous distributions of dose and clonogenic cell density. Phys Med Biol 1993;38:653–66. 36. Lyman JT. Complication probability as assessed from dosevolume histograms. Radiat Res 1985;Suppl 8:S13–9. 37. Kutcher GJ, Burman C. Calculation of complication probability factors for non-uniform normal tissue irradiation: the effective volume method. Int J Radiat Oncol Biol Phys 1989;16:1623–30. 38. Kno¨o¨s T, Kristensen I, Nilsson P. Volumetric and dosimetric evaluation of radiation treatment plans: radiation conformity index. Int J Radiat Oncol Biol Phys 1998;42:1169–76. 39. Emami B, Brown JLA, Coia L, Goitein M, Munzenrider JE, Shank B, et al. Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys 1991;21:109–22. 40. 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:439–46. 41. Eisbruch A, Ten Haken RK, Kim HM, Marsh LH, Ship JA. Dose, volume, and function relationships in parotid salivary glands following conformal and intensity-modulated irradiation of head and neck cancer. Int J Radiat Oncol Biol Phys 1999;47:1458–60. 42. Makkonen TA, Nordman E. Estimation of long-term salivary gland damage induced by radiotherapy. Acta Oncol 1987;26:307–12. 43. Meyer J, Phillips MH. Beam direction optimisation for anatomy-based IMRT using beam’s-eye-view metric, Poster presentation at the UK Radiation Oncology Conference 2003, University of Bath, Bath, UK, 7–9 April 2003. 44. Xiao Y, Galvin L, Hossain M, Valicenti R. An optimised forward-planning technique for intensity modulated radiation therapy. Med Phys 2000;27:2093–9.
327