Pre-segmented 2-Step IMRT with subsequent direct machine parameter optimisation – a planning study
© Bratengeier et al; licensee BioMed Central Ltd. 2008
Received: 04 June 2008
Accepted: 06 November 2008
Published: 06 November 2008
Modern intensity modulated radiotherapy (IMRT) mostly uses iterative optimisation methods. The integration of machine parameters into the optimisation process of step and shoot leaf positions has been shown to be successful. For IMRT segmentation algorithms based on the analysis of the geometrical structure of the planning target volumes (PTV) and the organs at risk (OAR), the potential of such procedures has not yet been fully explored. In this work, 2-Step IMRT was combined with subsequent direct machine parameter optimisation (DMPO-Raysearch Laboratories, Sweden) to investigate this potential.
In a planning study DMPO on a commercial planning system was compared with manual primary 2-Step IMRT segment generation followed by DMPO optimisation. 15 clinical cases and the ESTRO Quasimodo phantom were employed. Both the same number of optimisation steps and the same set of objective values were used. The plans were compared with a clinical DMPO reference plan and a traditional IMRT plan based on fluence optimisation and consequent segmentation. The composite objective value (the weighted sum of quadratic deviations of the objective values and the related points in the dose volume histogram) was used as a measure for the plan quality. Additionally, a more extended set of parameters was used for the breast cases to compare the plans.
The plans with segments pre-defined with 2-Step IMRT were slightly superior to DMPO alone in the majority of cases. The composite objective value tended to be even lower for a smaller number of segments. The total number of monitor units was slightly higher than for the DMPO-plans. Traditional IMRT fluence optimisation with subsequent segmentation could not compete.
2-Step IMRT segmentation is suitable as starting point for further DMPO optimisation and, in general, results in less complex plans which are equal or superior to plans generated by DMPO alone.
Historically, inverse planning algorithms , derived from tomographic calculations, have played a more important role than "forward planning" techniques [2–4] for IMRT. In contrast to traditional "trial and error" methods for three-dimensional conformal radiotherapy (3DCRT), progressive "forward planning" techniques for IMRT analyse the geometrical constellation to create beam segments. The weights of these segments can then be optimised in the same fashion as in conventional inverse IMRT planning. A few hospitals have specialised in such "forward planning" and can compete with "inverse planning" techniques [5–8]. Their approach is presumed not to be as conformal and flexible as an inverse planning based one; often "forward planning" is used in terms of class solutions. However, algorithms such as 2-Step IMRT [9, 10] are utilized for a variety of tumour localizations [6, 11]. 2-Step IMRT is a segmentation technique, which creates segments, reflecting the shape of the tumour and OARs but also highly weighted narrow segments close to critical structures to obtain steep dose gradients . On the other hand, contemporary "inverse planning" uses iterative optimisation elements. The latest generation of planning systems integrates machine parameters in the iterative optimisation process, such as HYPERION  and direct aperture optimisation (DAO) [14–16]. The Pinnacle3® planning system refers to it as "direct machine parameter optimization" ("DMPO", Raysearch™ laboratories). Such algorithms are clearly superior to those with sequencing after a complete fluence optimization process . Fine-tuning of the segment parameters, such as with DMPO, could possibly also enhance "forward planning" techniques. The aim of this work is to explore whether a technique like 2-Step IMRT could not only compete with a former generation planning systems (fluence optimisation followed by subsequent segmentation) as shown before , but also compete with a planning system of the latest generation where segmentation is integrated in the optimisation procedure such as with DMPO. A further motivation for this planning study is to investigate the possibility of daily ad-hoc adaptation of IMRT-plans[10, 18] based on 2-Step IMRT. This could only be useful if the primary plan can concur with the results of up to date IMRT planning systems. Without an initial IMRT plan of sufficient quality, all adaptive efforts would be non-optimal. With this in mind, it should be noted that the adaptation process itself is not subject of this paper.
Planning procedure DMPO-25, DMPO-50
The IMRT optimization in Pinnacle with DMPO starts with a conformal beam of uniform intensity followed by four steps of fluence optimization. This is followed by a step that includes machine parameters: Leaf-positions and segment weights are varied within the limits of the linear accelerator. Two such DMPO plans were generated for each case used in this paper. The first with only 25 optimisation steps (DMPO-25), the second with additional 25 steps (DMPO-50), to be able to track the speed of convergence. A further increase of the number of steps does not necessarily lead to better results[19, 20].
For both DMPO-plans slightly higher segment numbers were allowed to give them a chance to generate better plans than in combination with 2-Step IMRT (see next section).
Planning procedure 2S-DMPO-25, 2S-DMPO-50
For 2-Step IMRT, the segments were manually shaped according to the 2-Step IMRT-algorithm, which is described in detail in the preceding publications[6, 9–11, 19]. Several classes (orders) of beam segments were defined. All are generated by use of the beams-eye-view (BEV)-projections of PTV and OARs
0th order "conformal" segments include the target volume neglecting any OAR
1st order "OAR-sparing" segments include the PTV excluding the OAR
2nd order "narrow" segments were defined adjacent to the OAR. Their definition is more complex and needs 3D information 
Planning procedure "IMRT"
For comparison, the "IMRT" algorithm of the planning system – detailed fluence optimisation followed by the segmentation – was used as benchmark.
This planning study utilized Pinnacle3®(Philips/ADAC laboratories) for all planning steps. 15 clinical cases were chosen: A: 4 prostates with simultaneous integrated boost (SIB), B: 4 head and neck (3 cases with SIB), C: 4 breast with parasternal lymph nodes or with protruding thoracic wall, D: 1 angiosarcoma of the scalp, E: 1 bone-metastasis partially surrounding the spine, F: 1 bone metastasis with SIB. Furthermore the ESTRO Quasimodo-Phantom (Q) was used in two ways: with 9 and 15 equidistant irradiation angles. This oval phantom (20 cm × 37 cm, length 16 cm) contains a horse-shoe-shaped PTV (max. diameter 16 cm) surrounding a circular OAR (diameter 5 cm) at a distance of 0.5 cm. The body outside the PTV (Body\PTV) is considered as organ at risk.
Planning parameters and optimisation goals
For all clinical cases except the breast tumours 9 equidistant fields were employed starting from 180° (prostate cases) or 0° (other cases). For the left breast beams at angles 240°, 250°, 265°, 285°, 310°, 345°, 10°, 30°, 45°, 55° were applied and for the right breast 305°, 315°, 330°, 350°, 15°, 50°, 75°, 95°, 110° 120°. One of the breast cases (C2) was supplied by the group of Fogliata and Cozzi (case #3) to be able to compare with international standards. All clinical cases were primarily planned and optimised by DMPO. The set of objective values was adapted according to the clinical requirements. The head and neck cases were based on RTOG 0522 but with more demanding specifications with respect to homogeneity in the planning target volumes. Prostate patient plans were planned as described in Guckenberger et al. All plans required a homogeneity expressed by a standard deviation of less than 3.5% in the central part of the PTV (which was defined as the PTV reduced by 5 mm) or the CTV. The breast case from Fogliata and Cozzi  was planned to fulfil the constraints cited herein as closely as possible (details see in the results section). Finally, a total of 11 up to 31 objectives for 6 – 14 considered volumes were created (see Additional file 1). Help volumes were typically defined as several concentric envelopes surrounding the PTVs, or additional structures to accentuate concave parts of the PTV, if no particular OAR was present. For all plans of a given case the identical set of objectives was used. Only objectives and no constraints were utilised for the planning procedure. This was to keep the optimisation process flexible. Mainly 4 points (2 MinDVH, 2 MaxDVH) of a DVH-curve were used to shape a PTV-DVH, 2 points to describe the shape of an OAR-curve (more, if the OAR was near to the PTV or overlapping, less, if the OAR was distant from it). The weights for the PTV were chosen between 10 and 100 to force the system to follow the objectives very tightly. The penalty of "uniform dose" as used in 10 of the cases for the (central) PTV aims at both directions, lower and higher doses. Here a weight of 1 was chosen. The objective weights of the organs at risk were mostly between 0.1 and 3, to pull the DVH-curve loosely to lower values. In summary, a dynamic range of 1000 was used for the objective weights.
If a clinical plan existed, it was used as the "reference plan". In all other cases the DMPO-25-plan served as reference (see Additional file 2: A1, C1b, C2, C3b, Q9, Q15)
Each of the optimisation runs took typically half an hour. The manual segment shape definition took between one and to two hours, depending on the complexity of the case (future computerisation should reduce this step to a couple of minutes or less).
All entities (A-F, Q)
Regarding the individual cases, in 16 of 22 situations (50 steps) and 12 of 22 (25 steps) the 2S-DMPO-plan was superior to the DMPO, taken the COV as a measure of quality as to be discussed later. It can be stated that 2S-DMPO plans were equivalent with DMPO (if not better than DMPO, for 50 optimisation steps). "IMRT", the classical fluence optimisation followed by the sequencing process (data not shown in Additional file 2, but summarized here), clearly cannot compete, neither with respect to the quality (mean relative COV: 1.92), nor regarding monitor units (relative mean value: 1.33) or segment numbers (mean number of segments:143-always more segments than with any other technique used in this work). The "IMRT" results should be seen as a reference for the comparison of 2S-DMPO and DMPO: Their differences are rather small in relation to "IMRT".
A closer look at the results shows differences for the groups of cases. In all groups except B and C the 2S-DMPO-50 COV was lower than that for DMPO-50. For the breast cases (C) three 2S-DMPO-plan were not advantageous. Therefore, the breast cases were scrutinized in more detail.
Breast case (C)
Adopting the same measures to cases C1, C3 and C4, 2S-DMPO-50 and DMPO-50 were better in 59 and 40 instances, respectively (110 parameters in total, not shown here in detail). However, also this result can also not identify an optimal method.
Quasimodo case (Q)
The Quasimodo  phantom constitutes the clearest geometrical constellation, which was challenging because the PTV partially surrounded the OAR with small distance between them. According Bohsung et al.  the dose is normalised to the mean dose in the PTV. The 95%-isodose should include 99% of the PTV-volume, whereas the 105%-isodose should encompass not more than 5% of the PTV-volume. The constraints concerning OAR and Body\PTV are given in Tab 4 (see Additional file 4). The quality index SD is defined simply as the sum of the differences of the constraints and the actual values if they violate the constraints. Ideally the SD is near to zero. The evaluation of the QUASIMODO-data in Tab. 4 shows notable advantage of the 2S-DMPO compared to DMPO. For 9 irradiation directions only much more segment numbers (75) of a DMPO plan lead to equivalent results as a 40 segment 2S-DMPO technique (SD = 6.1 vs. 6.8). For 15 gantry angles the desired quality index "0" could be achieved only by 2S-DMPO-50, indicating that the optimal conditions were almost met by 2-Step IMRT. Even 2-Step IMRT as segmentation alone with subsequent weight optimisation, without subsequent DMPO ("2S"), generated better results (SD = 4.6) than DMPO-25 with 25 iteration steps (SD = 4.9). As a secondary result, it should be noted that an increase of gantry angles – even with no increase in the number of segments (75 in both cases) – conceptually can enhance the quality of a plan.
Quality of 2S-DMPO
It should be noted that on average 20 "objectives" (see Additional file 1) were used as described in the materials section to shape the DVH of each plan; all of them were weighted quadratic deviations from a desired point of a DVH curve. Based on this, the CSV is well suited to assess the quality of an approach relative to a desired DVH, encouraging the authors to utilize it to analyze DMPO and 2S-DMPO. The parameter sets were primarily developed for pure DMPO-plans; so it is debatable whether they could be advantageous for DMPO. However, besides the breast cases that needed further consideration, most other case-types benefited from a 2-Step IMRT-pre-definition of the segments.
Furthermore the comparison with the published data of other groups in the Quasimodo-project  confirms the high quality of the 2S-DMPO-plans – without exception all 2S-DMPO plans would be placed in the leading group of all participants, the 15-field 2S-DMPO-50 result of this study was only achieved by an intensity modulated arc technique in the Quasimodo-group. Not even a prolongation of the optimisation process to 100 steps enabled the pure DMPO plan to match 2S-DMPO-50 (result not presented here). For the 9-field plan only a drastic increase of the segment number (from 45 to 75) led to a plan of similar quality as the 9-field 2S-DMPO. This suggests that the 2-Step IMRT pre-segmentation is proximate to the theoretical optimum which could not always be met by the pure DMPO.
For the breast data 2-Step IMRT was discussed critically in former work [9, 10]. Theoretical considerations  conjectured that for breast cases a third segment step could be favourable (that means i. e. splitting of the of 2nd order segments in a broader and a narrower part). The breast results of pure 2-Step IMRT  were the "worst" compared to all other cases. Nevertheless 2S-DMPO was comparable to pure DMPO in this study. A detailed view on the results in Tab. 3 (see Additional file 3) shows that both methods lead to similar results; especially when considering the wide range of results delivered by different working groups with different planning systems. With respect to all values, the results from this work are qualitatively better than most from the Fogliata study; in general they are in the leading group. The mean values of this study were better for 21 out of 22 parameters than the mean values over all participant systems. The results of this study closely resembled the results of the two variants of Pinnacle-applications "Pinn1" and "Pinn2" from the Fogliata-study . The "b" variants 2S-DMPO-50 DMPO-50 achieved better results in almost all categories besides the dose to the ipsilateral lung; however, the mean dose to both lungs was lower. The "c"-variants also tended to be better with respect to the lung doses and the minimum dose in the PTV with minor disadvantages with respect to hot spots in the PTV. The differences between DMPO and 2S-DMPO were much smaller within one set of optimisation parameters in contrast to the differences regarding other parameter sets or the Pinnacle3® plans presented by Fogliata et al.  in their study. The differences within all Pinnacle3® plans and the plans based on other planning systems from their study were even greater.
Plan quality and number of gantry angles
Do more gantry angles entail more segments? It can be shown, that increasing the number of gantry angles does not increase or even reduce the necessary number of segments: In a thought experiment an IMRT plan is optimized for a given configuration of gantry angles and number segments. The optimization is performed by minimizing an objective function (OF) of arbitrary type. Now, one or more additional gantry angles should be allowed, however, with the same number of segments in total. The system gains more degrees of freedom. Nevertheless, the primary plan is one possible solution, because all the new gantry angles can be rejected without increasing the value of the OF. Perhaps another configuration leads to a lower value for the OF and therefore to a better plan. In a next step it could be tried to reduce (!) the number of segments. That is just the situation as for the Quasimodo-case: For DMPO-50 and 9 gantry angles the 75 segment plan the quality score SD = 6.1 is much higher than for 15 gantry angles and 70 segments (SD = 2.7). For this reason few gantry angles are not automatically favourable.
Also, increasing the number of beams with constant number of segments need not be more time consuming. In all cases with equidistant gantry angles, a total gantry angle of 360° (n-1)/n must be covered (n, n': number of fields). Furthermore the time loss Δtg per starting or stopping-process of the gantry has to be considered. Comparing two techniques changing from n to n' fields, a time loss of t = 360° (1/n' - 1/n)/v + Δtg (n - n') results for a gantry velocity v. On the other hand, a time gain is to be taken into account, because during the gantry motion a segment can be formed, gaining ΔtS per such segment and overall gaining ΔtS (n' - n) comparing the two techniques. Typical values for two types of the accelerators in our department (values of the second one in brackets) are v = 360°/78 s (360°/135 s), Δtg = 4 s (2 s), Δtg = 8.5 s (4 s), resulting in a time gain of 11 s (the other linac type: 4 s time loss). for the 9-field technique in relation to a 5-field technique with the same number of segments, monitor units and does rate. Obviously the number of gantry angles must be discussed separately. Techniques using more gantry angles are not automatically slower and could be even advantageous with respect to time efficiency.
In summary, segment shape optimization with DMPO was nearly equivalent for both methods, the segment pre-definition by 2-Step IMRT and the pre-definition by inverse planning (with some advantages for 2S-DMPO). The 2-Step IMRT algorithm has been shown to be a well-suited method to pre-define segments, even for breast cases that were thought to be inappropriate in earlier studies for 2-Step IMRT alone. In all 22 cases and variants, 2-Step IMRT combined with DMPO was able to produce results which could compete with IMRT algorithms of the newest generation, i.e. DMPO alone. The quality of the plans were superior, when compared with segmentation after completion of the fluence optimisation.
One advantage of 2-Step IMRT is that it generates "intuitive" IMRT segments, reflecting the anatomy of the patient. In earlier work [9, 10] the author proposed 2-Step IMRT as a segmentation method with the potential for fast day-to-dayIMRT adaptation. This feature can now be investigated in further work.
The authors would like to thank L. Cozzi for providing the patient model C2 from  and the detailed planning results for this special case.
Also the authors thank J. Bohsung for providing the Quasimodo phantom .
The authors thank both for making the patient models available to the public:
- Bortfeld T, Burkelbach J, Boesecke R, Schlegel W: Methods of image reconstruction from projections applied to conformation radiotherapy. Phys Med Biol 1990,35(10):1423-1434. 10.1088/0031-9155/35/10/007View ArticlePubMedGoogle Scholar
- Bednarz G, Michalski D, Houser C, Huq MS, Xiao Y, Anne PR, Galvin JM: The use of mixed-integer programming for inverse treatment planning with pre-defined field segments. Phys Med Biol 2002, 47: 2235-2245. 10.1088/0031-9155/47/13/304View ArticlePubMedGoogle Scholar
- Claus F, De Gersem W, Vanhoutte I, Duthoy W, Remouchamps V, De Wagter C, De Neve W: Evaluation of a leaf position optimization tool for intensity modulated radiation therapy of head and neck cancer. Radiother Oncol 2001, 61: 281-286. 10.1016/S0167-8140(01)00441-8View ArticlePubMedGoogle Scholar
- Xiao Y, Galvin J, Hossain M, Valicenti R: An optimized forward-planning technique for intensity modulated radiation therapy. Med Phys 2000, 27: 2093-2099. 10.1118/1.1289255View ArticlePubMedGoogle Scholar
- Bos LJ, Schwarz M, Bar W, Alber M, Mijnheer BJ, Lebesque JV, Damen EM: Comparison between manual and automatic segment generation in step-and-shoot IMRT of prostate cancer. Med Phys 2004, 31: 122-130. 10.1118/1.1634481View ArticlePubMedGoogle Scholar
- Bratengeier K, Guckenberger M, Meyer J, Mueller G, Pfreundner L, Schwab F, Flentje M: A comparison between 2-Step IMRT and conventional IMRT planning. RadiotherOncol 2007, 84: 298-306. 10.1016/j.radonc.2007.06.018Google Scholar
- De Gersem W, Claus F, De Wagter C, De Neve W: An anatomy-based beam segmentation tool for intensity-modulated radiation therapy and its application to head-and-neck cancer. Int J Radiat Oncol Biol Phys 2001, 51: 849-859.View ArticlePubMedGoogle Scholar
- De Neve W, De Wagter C, De Jaeger K, Thienpont M, Colle C, Derycke S, Schelfhout J: Planning and delivering high doses to targets surrounding the spinal cord at the lower neck and upper mediastinal levels: static beam-segmentation technique executed with a multileaf collimator. RadiotherOncol 1996, 40: 271-279. 10.1016/0167-8140(96)01784-7Google Scholar
- Bratengeier K: 2-Step IMAT and 2-Step IMRT: A geometrical approach. Med Phys 2005, 32: 777-785. 10.1118/1.1861161View ArticlePubMedGoogle Scholar
- Bratengeier K: 2-Step IMAT and 2-Step IMRT in three dimensions. Med Phys 2005, 32: 3849-3861. 10.1118/1.2134928View ArticlePubMedGoogle Scholar
- Bratengeier K: IMRT Segment Predefinition using 2-Step IMRT. In Radiother Oncol. Volume 84. 8th Biennial ESTRO meeting on Physics and Radiation Technology for Clinical Radiotherapy; Lisboa; 2005:S190.
- Brahme A, Roos JE, Lax I: Solution of an integral equation encountered in rotation therapy. Phys Med Biol 1982,27(10):1221-1229. 10.1088/0031-9155/27/10/002View ArticlePubMedGoogle Scholar
- Alber M, Nusslin F: Optimization of intensity modulated radiotherapy under constraints for static and dynamic MLC delivery. Phys Med Biol 2001, 46: 3229-3239. 10.1088/0031-9155/46/12/311View ArticlePubMedGoogle Scholar
- van Asselen B, Schwarz M, van Vliet-Vroegindeweij C, Lebesque JV, Mijnheer BJ, Damen EM: Intensity-modulated radiotherapy of breast cancer using direct aperture optimization. Radiother Oncol 2006, 79: 162-169. 10.1016/j.radonc.2006.04.010View ArticlePubMedGoogle Scholar
- Cho PS, Marks RJ 2nd: Hardware-sensitive optimization for intensity modulated radiotherapy. Phys Med Biol 2000, 45: 429-440. 10.1088/0031-9155/45/2/312View ArticlePubMedGoogle Scholar
- Siebers JV, Lauterbach M, Keall PJ, Mohan R: Incorporating multi-leaf collimator leaf sequencing into iterative IMRT optimization. Med Phys 2002, 29: 952-959. 10.1118/1.1477230View ArticlePubMedGoogle Scholar
- Dobler B, Pohl F, Bogner L, Koelbl O: Comparison of direct machine parameter optimization versus fluence optimization with sequential sequencing in IMRT of hypopharyngeal carcinoma. Radiat Oncol 2007, 2: 33. 10.1186/1748-717X-2-33PubMed CentralView ArticlePubMedGoogle Scholar
- Bratengeier K: An IMRT adaptation concept for interfractional organ motion. In Radiother Oncol. Volume 76. 8th Biennial ESTRO meeting on Physics and Radiation Technology for Clinical Radiotherapy; Lisboa; 2005:S199.
- Bratengeier K, Meyer J, Schwab F, Vordermark D, Flentje M: Steep dose gradients for simultaneous integrated boost IMRT. Zeitschrift für Medizinische Physik, in press. accepted 23.7.2008
- Carlsson F, Forsgren A: Iterative regularization in intensity-modulated radiation therapy optimization. Med Phys 2006, 33: 225-234. 10.1118/1.2148918View ArticlePubMedGoogle Scholar
- Bohsung J, Gillis S, Arrans R, Bakai A, De Wagter C, Knoos T, Mijnheer BJ, Paiusco M, Perrin BA, Welleweerd H, Williams P: IMRT treatment planning: – a comparative inter-system and inter-centre planning exercise of the ESTRO QUASIMODO group. Radiother Oncol 2005, 76: 354-361. 10.1016/j.radonc.2005.08.003View ArticlePubMedGoogle Scholar
- Fogliata A, Nicolini G, Alber M, Asell M, Dobler B, El-Haddad M, Hardemark B, Jelen U, Kania A, Larsson M, Lohr F, Munger T, Negri E, Rodrigues C, Cozzi L: IMRT for breast. a planning study. Radiother Oncol 2005, 76: 300-310. 10.1016/j.radonc.2005.08.004View ArticlePubMedGoogle Scholar
- Guckenberger M, Flentje M: Intensity-Modulated Radiotherapy (IMRT) of Localized Prostate Cancer: A Review and Future Perspectives. Strahlenther Onkol 2007, 183: 57-62. 10.1007/s00066-007-1631-3View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.