 Research
 Open Access
 Published:
Characteristics of noncoplanar IMRT in the presence of targetembedded organs at risk
Radiation Oncology volume 10, Article number: 207 (2015)
Abstract
Background
The aim is to analyze characteristics and to study the potentials of noncoplanar intensity modulated radiation therapy (IMRT) techniques. The planning study applies to generalized organ at risk (OAR) – planning target volume (PTV) geometries.
Methods
The authors focus on OARs embedded in the PTV. The OAR shapes are spherically symmetric (A), cylindrical (B), and bended (C). Several IMRT techniques are used for the planning study: a) noncoplanar quasiisotropic; b) two sets of equidistant coplanar beams, half of beams incident in a plane perpendicular to the principal plane; c) coplanar equidistant (reference); d) coplanar plus one orthogonal beam. The number of beam directions varies from 9 to 16. The orientation of the beam sets is systematically changed; dose distributions resulting from optimal fluence are explored. A selection of plans is optimized with direct machine parameter optimization (DMPO) allowing 120 and 64 segments. The overall plan quality, PTV coverage, and OAR sparing are evaluated.
Results
For all fluence based techniques in cases A and C, plan quality increased considerably if more irradiation directions were used. For the cylindrically symmetric case B, however, only a weak beam number dependence was observed for the best beam set orientation, for which noncoplanar directions could be found where OAR and PTVprojections did not overlap. IMRT plans using quasiisotropical distributed noncoplanar beams showed stable results for all topologies A, B, C, as long as 16 beams were chosen; also the most unfavorable beam arrangement created results of similar quality as the optimally oriented coplanar configuration. For smaller number of beams or application in the trunk, a coplanar technique with additional orthogonal beam could be recommended. Techniques using 120 segments created by DMPO could qualitatively reproduce the fluence based results. However, for a reduced number of segments the beam number dependence declined or even reversed for the used planning system and the plan quality degraded substantially.
Conclusions
Topologies with targets encompassing sensitive OAR require sufficient number of beams of 15 or more. For the subgroup of topologies where beam incidences are possible which cover the whole PTV without direct OAR irradiation, the quality dependence on the number of beams is much less pronounced above 9 beams. However, these special noncoplanar beam directions have to be found. On the basis of this work the noncoplanar IMRT techniques can be chosen for further clinical planning studies.
Background
The aim of this work is to analyze characteristics of noncoplanar IMRT for planning target volumes (PTVs) encompassing an organ at risk (OAR).
Types of PTVOAR topologies
For intensity modulated radiotherapy, OARs surrounding a convex PTV lead to significantly different solutions than the reverse situation of a PTV encompassing an OAR: Only in the latter case can the OAR be completely spared, if the scatter is neglected (type I geometry, see Fig. 1). An ideal OAR sparing is achieved by the highly inhomogeneous fluence distribution peaking near the beam eye view projection of the OAR [1, 2]. The necessary fluence distribution can be numerically found for arbitrarily complicated PTVOAR geometries by the inverse planning procedures controlled by the set of dosevolume objectives [3]. For computed tomography (CT) image reconstruction, an incomplete scanning of the Radon space is known to produce aliasing artifacts in image reconstruction. To identify small structures, fine sampling and high resolution is necessary (see i.e. Kak p. 178 [4]). Similar mathematical formalism is also applicable to IMRTbased dose shaping. The generation of deliberately fine textured dose distributions will therefore follow the same rules: Similar to computed tomography, the minor distorted solutions need a sufficient number of projections (beam directions), best in an equidistant manner.
Note that one OAR or several OARs in a dense sequence distributed around a convex PTV cannot in principle be perfectly spared in the coplanar beam arrangement: there are always beam orientations critical for the dose saturation in the PTV, which irradiate through the OARs (type III geometry, Fig. 1). Exact solutions for the inverse planning problem do not exist for the positive definite solution space: negative fluences would be needed. The desired balance between the PTV coverage and the OAR sparing is reflected in the weights of dosevolume objectives used for the inverse planning and has to be controlled by the clinician. The optimal number of beams, their directions and fluence distributions depends on the properties of the OARs. A sparse distribution of OARs around PTV, which allows beam directions which perfectly spare the OAR (type II geometry, Fig. 1) leads to a problem of beam angle optimization (BAO). Often only few beams are needed for high quality plans, if appropriate directions are chosen. In real life, type I, II and III situations are mostly mixed. It should be borne in mind that if noncoplanar directions are allowed, a type I or type III problem can transform to type II. This work is mainly concentrated on type I and type II problems, a PTV encompassing an OAR.
Noncoplanar techniques vs. coplanar techniques
For type I topologies (PTV encompassing an OAR), it is not yet obvious whether the noncoplanar beam arrangement can improve OAR sparing, as the OAR can be perfectly spared also in the coplanar beam arrangements [1]. An increased sparing of healthy tissue is reported in the IMRT planning study for the simple spherical phantom and a noncoplanar beam arrangement [5]; an increase of the therapeutic width due to improved OAR sparing for the noncoplanar IMRT and a number of clinical cases is reported in [6]. Many compact (eyeballs, lenses, parotid glands, pituitary gland, etc.) or elongated (rectum, urethra, spinal cord, optic nerves) OARs exist that should benefit from noncoplanar IMRT. In the present work, the authors expand the scope of [5] and present the comparison of solutions for more complex phantom topologies and beam constellations. The purpose is to analyze the characteristics of noncoplanar beam incidence versus conventional coplanar IMRT for the type I geometry of PTV encompassing an OAR and to provide simple planning rules for IMRT/intensity modulated arc therapy (IMAT)/volumetric arc therapy (VMAT) planning studies with real patients that will follow. The clinical significance of the observed effects can be evaluated using appropriate biological models for various clinical situations and is beyond the scope of this work.
For type III topologies, the stereotactic treatments using noncoplanar beam arrangements are mostly used. Contrary to a coplanar beam arrangement, additional degrees of freedom of the noncoplanar beam arrangement would in most cases allow better OAR sparing without compromising PTV coverage. At the distances larger than a target diameter away from the target, the dose gradient approaches 1/r^{2} instead of 1/r behavior for noncoplanar vs. coplanar techniques, respectively. The price for the improved OAR sparing could be longer pathways in the body with additional dose load to the healthy tissue which potentially limits the use of noncoplanar techniques. This aspect of the problem cannot be discussed here, as the solution depends strongly on the actual biological and geometrical parameters.
Number of beams
The question of the number of beams has to be addressed early: what is the appropriate range of beam numbers to investigate? It is important to note that type I, type II, and type III geometries probably require different number of beams. E.g. for type III geometry, the dose load to surrounding organs can be spread over wider areas and “thinned out”. More available (noncoplanar) beam angles potentially allow better organ sparing, especially of serial organs. For type II geometry, an advantageous beam angle optimization (BAO) can find favorable beam arrangement: even five selected noncoplanar beams achieve the quality of nine equidistant coplanar beams; see e.g. Bangert et al. [7] who investigated type Ilike head and neck cases and observed the beam number dependencies that level off between 9 and 11 beams. A good example for a different behavior can be found in Rossi et al. [8]. The authors compare coplanar and noncoplanar techniques for stereotactic body radiotherapy (SBRT) prostate irradiation and observe clear beam number dependences above ten beams. This discrepancy can be understood as follows: Rossi’s prostate cases require clearly concave isodose lines towards rectum and embedded urethra, which we classify as type I geometry. Type I geometry is supposed to require more beams than type II or type III. The larger the OAR radius and smaller the PTVOAR distance (or the steeper the required dose gradient to the OAR), the more equidistant beams are required for successful dose painting, in analogy to inverse reconstruction methods in computer tomography. Equidistance was proposed by Schreibmann et al. as a class solution for the therapy of prostate carcinoma [9]. As noncoplanar beams are distributed over 2π solid angle (4π if beam entrance and exit are considered), the angular distance between the beams is inherently larger and the advantages of noncoplanar arrangements should become noticeable at least for the number of beams larger then 9–10. Therefore for this study we evaluated plans with 9–16 beams, 16 beams were chosen as an appropriate practical limit, see Results section.
Number of segments
The question of a sufficient number of segments is restricted to step and shoot IMRT. The number of necessary segments can be estimated based on the pioneering work on inverse planning by Brahme et al. [1]. They discovered the need of a substantial fluence increase near the BEV projection of the OAR to compensate the underdosage in the PTV without simultaneously to compromise the OAR sparing. For an important OAR surrounded by a PTV, the fluence profile featuring such “Brahmepeak” can be realized in practice by two leftsided, two rightsided, and a base segment (S0) [2]. The base segment was termed segment of the order 0, S0, segments which blend the OAR out were called segments of order 1, S1, narrow segments adjacent to the OAR were called segments of order 2, S2. According to this rough guess, up to five segments per beam are recommended.
Materials and methods
Materials
A virtual linac Elekta Synergy S with BeamModulator™ multileafcollimator (MLC; “S”), leaf width 4 mm, was commissioned in the therapy planning system (TPS) Pinnacle3™. The optimization engine was “Intensity Modulation™” module available in Pinnacle3™ with 80 steps for pure fluence distribution, and Step and Shoot Direct Machine Parameter Optimization (DMPO™) with two runs 40 steps each for DMPO, thereof ten steps for fluence optimization before sequencing. For DMPO, the total number of segments is given and it is up to a sequencer to distribute the segments among the beams. The resulting fluence distributions were investigated qualitatively.
The gradient based optimization algorithm is based on the composite objective value (COV) as objective function. It sums weighted relative quadratic deviations from the clinical dosevolume requirements. After the optimization procedure, the dose is normalized to the mean dose in the PTV (resulting in monitor unit (MU) corrections of less than 0.5 %).
Research topics
The purpose of this work was to investigate capabilities of IMRT with noncoplanar beams. Following parameters were varied for systematic studies:

OAR shape (Fig. 2):

A)
spherically symmetric (spherical),

B)
axially symmetric (cylindrical),

C)
mirror symmetric (bananashaped)

irradiation technique (Fig. 3):

a)
noncoplanar quasiisotropic,

b)
coplanar split in two orthogonal planes,

c)
coplanar equidistant,

d)
coplanar equidistant plus one additional orthogonal beam.
For fluence optimization, the relative orientation of the beams was maintained, and the orientation of the whole set of beams with respect to the PTVOAR arrangement was varied for the OAR shapes B and C to study the effect of the particular beam set orientation on the plan quality. Taking symmetries into account, up to 58 orientations were systematically investigated to allow an estimation of lowest and highest COV for each technique in any orientation. Dense beam axis positioning with small angles (3° or 5° for high and low number of beams, respectively) with respect to the OAR axis was especially examined to find the overall minimal COV.
OARs, PTVs, and irradiation techniques are described below in more detail. Plan quality is conjectured to depend strongly on:

number of beams, therefore for each plan 2–3 practical number of beams, between 9 and 16, was chosen: a: 10, 16; b, c: 9, 11, 15; d: 10, 12, 16; (always highly symmetric configurations)

number of segments, therefore for each beam arrangement optimal fluence distribution (no sequencing), “large” (120) and “small” (64) number of segments were considered.
Gantry and table angles for the 15 and 16 beam technique variants are presented in Table 1.
The following parameters were  not independently  varied for completeness (only for sequenced techniques):

set of objectives which tend to achieve more PTV coverage or more OAR sparing: PTVconformal (15 objectives) or OARsparing (22 objectives), respectively,

dose grid 2 mm / 4 mm,

collimator angle 0°/20° (K00, K20, planned with Pinnacle3™ V 9.0 and 9.2 respectively).
The following combinations were chosen: K20_2mm_PTVconformal, K00_4mm_PTVconformal, K00_4mm_OARsparing.
For more extensive exploration of different orientations of the beam sets, only K00_4mm_PTVconformal was used.
PTVOAR topologies
The OAR shapes are motivated by the clinical cases mentioned in the introduction. Three types of OAR shapes shortly described above were embedded in a spherical PTV, 10 cm in diameter, centered in a ball shaped phantom, 20 cm in diameter. The outer shape is chosen to be as regular as possible, not to mask the effect of embedded organs. The OARs were inflated by 5 mm in all directions to define the cutout of the PTV. The following OAR shapes were considered (Fig. 2):

A)
ball of diameter 3 cm;

B)
cylinder of diameter 3 cm oriented parallel to the table axis with rounded ends, penetrating the PTV;

C)
bananashaped OAR of diameter 2 cm and bending radius 5 cm, penetrating the PTV.
In contrast to the cylindrical geometry of case B, for the cases A and C there is no gantry angle that would allow complete irradiation of the PTV combined with complete blocking of the OAR. It should be mentioned that A and C always lead to type I situations, also when switching from coplanar to noncoplanar techniques. Case B, however, becomes type IIlike for a noncoplanar technique: it is obvious to direct one of the beams along the cylinder axis for complete OAR blending. All elongated OARs were parallel to the table axis, similar to typical clinical objects (rectum, spinal cord).
For each PTVOAR topology from Fig. 2 several plans with 2–3 practical numbers of beams were created in each of the following categories (see Fig. 3 and Table 1):

a)
Noncoplanar quasiisotropic (NonCo): A quasiisotropic arrangement of 16 beams of triakis icosahedral symmetry (a combination of icosahedral and dodecahedral solids) [5] resulting in angular distances between the beams of 63°, 42° and 37°. It was compared to a dodecahedral combination of 10 fields alone. Both highly symmetric beam arrangements can be applied using only five different table angles (which is advantageous, since table rotations are timeconsuming).

b)
Field axes in two orthogonal planes each (2P): 9, 11, and 15 beams arranged in two equidistant coplanar sets, with half of beams applied perpendicular to the principal plane via table rotation to 90°.

c)
Coplanar equidistant (Co): classical coplanar technique with equidistant 9, 11 and 15 beams (reference).

d)
Coplanar equidistant supplemented by a single orthogonal beam along the table axis (Co + 1).
The number of beams was chosen to allow a high degree of symmetry. An entrance port of a beam never coincided with an exit port of another beam. All techniques cover maximally one hemisphere; the substitution of a beam by an opposed beam in all examined samples did not essentially influence the results. Except technique a, all others are suitable for irradiation of the trunk. For each combination of PTVOAR topology and beam arrangement, the dependence of the plan quality indices on the number of beams was investigated for pure fluence distribution, “large”, and “small” number of segments. Technique a (NonCo) was motivated by the question, whether many noncoplanar beam directions routinely used for standard stereotactic radiotherapy of type II and type III geometries are also useful in the reverse situation of the OAR surrounded by the target. The motivation behind the orthogonal beam sets of the technique b (2P) was to avoid the shadow zone of the cone beam geometries in the Radon space and to satisfy the TuySmith condition (see e.g. Buzug [10] pp. 368–378 or [11]). The technique c (Co) was chosen as a reference. Technique d (Co + 1) can be regarded as an intermediate technique between b and c.
Quality indicators
We evaluated plan quality of optimized plans in different ways. First, overall plan quality can be characterized by the composite objective value (COV) at the end of optimization. COV is the sum of all weighted objective values as described above. Objective values are defined as volumenormalized quadratic penalties determined for points in the dose volume histograms (DVHs). The lower the COV for a certain set of objectives, the better the beam/segment arrangement fulfills the dose requirements expressed through this set of objectives. Due to the high number of objectives (15 and 22, respectively), the deviation from the desired DVH shape is adequately represented by the residual COV at the end of optimization. COV is thus a plan quality indicator as long as all relevant DVH parameters are reflected in the objectives [12]. For the purpose of this study we used COV normalized to the best (lowest) COV among all techniques for a certain topology, nCOV.
Besides nCOV two indices were used to measure the overall plan quality in PTV and OAR:

quality index S _{ D }, sum of violations of dose requirements for PTV and OAR [13]. For the purpose of this study S _{ D } index considered violations of the following requirements: PTV V_{95%} ≥ 99 %; PTV V_{105%} ≤ 5 %; OAR V_{70%} ≤ 1 %; healthy tissue (whole phantom, PTV excluded) HT V_{80%} ≤ 15 %; HT V_{100%} ≤ 2 %; HT D_{max} ≤ 105 % . S _{ D } = 0 if all above conditions are fulfilled.

Conformity number CN [14] for reference isodose D _{ RI } = 95 %; defined as CN := (TV _{ RI } ) ^{2} /(V _{ RI } TV)
for target volume TV, part of target volume surrounded by the reference isodoses RI, TV _{ RI }, and total volume V _{ RI } surrounded by the reference isodose.
Additionally, the necessary monitor units MU and the mean dose to the healthy tissue, HT D _{ mean }, were evaluated.
Results
Qualitative analysis of beam fluence
Figure 4 shows optimized fluence distributions for different combinations of PTVOAR topology and irradiation technique in standard orientation (see Fig. 3). Substructures of Fig. 4 can be viewed in more detail in the Additional file 1. The fluence Brahmepeak [1, 2, 5] adjacent to the OAR projection could clearly be seen as a dark shaded area parallel to the rotation axes of the gantry, see e.g. Fig. 4Ba. Technique b uses two orthogonal rotation axes, so for OAR B (see Fig. 4Bb), the fluence at table angle 90° included Brahme peaks also above and below the cylinder. For more isotropic irradiation (technique a) (Fig. 4Aa), Brahme peak surrounded the OAR [5]. For noncoplanar techniques a, b and d substructures of the fluence distribution reflected the symmetry of the beam arrangement. E.g. for Aa five table angles resulted in a C_{10}symmetrylike fluence structure. Co + 1 technique d (Fig. 4Ad) generated substructures in the coplanar fields: parallel stripes of reduced fluence tangential to the OAR along the direction of the noncoplanar beam. The substructures in the fluence of the noncoplanar field itself (Fig. 4Ad, lowest subfigure) reflected 15 equidistant coplanar fields.
We do not show fluences for different numbers of beams, in conclusion, more beams led to finer and less intense modulations; the maximal number of beams, 16, was chosen as a practical compromise. The noncoplanar technique a in the limiting case of infinitely many beams could be considered as a superposition of infinite number of coplanar irradiations around isotropical distributed rotation axes. For this case, all substructures except the Brahme peak wiped out. For a noncoplanar technique with limited number of beams substructures in the fluence distribution could not be avoided. Note, that Brahme peak required for noncoplanar IMRT is narrower than that for coplanar IMRT, as was observed in [5].
Plan quality for different beam sets and beam set orientations
For each irradiation technique a, b, c, d, the minimum (Min) and the maximum (Max) nCOV value were determined for all possible beam set orientations (Fig. 5).
Characteristics of the different topologies
For spherical symmetric geometry A, Min and Max were identical: there is no preferred beam set orientation for any irradiation technique. For cylindrical geometry B, Min and Max differed the most among all noncoplanar techniques. In addition, whereas for all geometries an obvious beam number dependence of Min was observed, Min for cylindrically symmetric geometry B depended only weakly on the number of beams for all noncoplanar techniques. The minimum of nCOV for geometry B always occurred when one of the beams was precisely incident along the cylinder axis. This is a pure type II situation: such a beam spares the OAR completely without compromising PTV coverage, the increased number of beams does not improve plan quality much. For the “bent” bananashaped geometry C the preferred beam arrangement corresponding to Min is less obvious, the number of beams plays more role: the dependence of Min on the number of beams is more pronounced.
Best choice of technique
For nontype II topologies the number of beams was more important than the choice of technique. Techniques a and d both provide good results. The quasiisotropic technique a was less sensitive with respect to the orientation of the beam set; however, this technique was very sensitive to changes of the beam number and demands 16 beams (or more). The Co + 1 technique d would be the most relevant technique, if the number of beams is restricted to about ten, see Figs. 5 and 6. Also, the 2Ptechnique b or Co + 1 d always increased the plan quality with respect to the coplanar situation at a comparable number of beams.
Plan quality for different numbers of segments
To obtain linac deliverable plans, the ideal fluences were sequenced into segments which were then run through DMPO algorithm. For each technique only one beam set orientation was chosen for sequencing.
As the quality of fluenceoptimized plans were largely independent from the combination of the collimator angle and the set of dose objectives (K20_PTVconformal, K00_PTVconformal, K00_OARsparing), only the mean values over three constellations are presented. The dependences of quality indices nCOV and S _{ D } on number of beams for chosen number of segments are shown in Fig. 6; they show the same trends. The CN index is not depicted, as its dependence was similar to that of nCOV and S _{ D }, but less pronounced. For all techniques with 120 segments the observable trends were qualitatively the same as for fluence based plans: both COV and S _{ D } indices decreased for increasing number of beams; the plan quality improved. The quality improvement was most pronounced for the quasiisotropic noncoplanar technique a. For the smaller numbers of beams and segments, the quality of the technique decreased. Note, that for 64 segments, a minimum S _{ D } occurred for the number of beams 11–12, i.e. at about five segments per beam as expected (see Background section). For the larger number of beams, i.e. smaller number of allowed segments per beam, the plan quality degraded. The authors suppose a problem with the DMPO sequencer in Pinnacle3™, which does not allow to distribute different segment orders to neighboring beams for decreased number of segments per beam as would be necessary ideally [15]. Such behavior is not expected for optimization algorithms like simulated annealing [16].
Monitor units and healthy tissue
We observed an increasing MU number for larger number of beams (Fig. 7b): always for 120 segments (typically around 3 % more MU for each additional beam), mostly for 64 segments. Plans using fewer segments needed less MU. However, the mean dose to the healthy tissue did not increase correspondingly (Fig. 7a). Smaller segment areas compensated the increased number of monitor units. HT D _{ mean } was almost independent from the beam arrangement for type I topologies A and C.
Discussion
The important observation is a clear dependence of the plan quality on the number of beams even above 9 to 11 beams, which is in agreement with results of other studies [8] for type I topologies like A and C. The effect of beam number is even more important than the choice of technique. Better dose painting results for the larger number of radiotherapy beams can be understood as an inverse analogy to the CT technique, where small angular steps between image acquisitions ensure small “reconstruction errors” (aliasing artifacts) [3, 17]. It remains to be investigated whether this effect is recognizable for the dense series of gantry angles of VMAT  provided the used optimization algorithm supports these potential abilities of VMAT. The fluenceoptimized twoplane technique b was always better than the beamnumberequivalent coplanar technique c. Whether this behavior can be traced back to a better coverage of the Radon space [10] could be tested using larger and smaller targets and OARs, leading to larger or smaller aperture angles of the beams and  corresponding  less or more coverage.
However, the situation in case B is quite different. Here the target and the OAR do not overlap in the eye view of the beam parallel to the common axis of OAR and PTV. So, PTV can be covered from this direction without OAR sparing to be compromised. This is a typical type II situation, where a good choice of a few beam directions can be more effective than an increased number of beams, and no additional beam can add serious benefits. Thus, in Fig. 6b, the quality indices nCOV and S _{ D } are nearly independent from the number of beams in the studied range for “well chosen” beam arrangement d. If such a “preferred” beam direction is not obvious, a beam angle optimization may be necessary.
The three types of PTVOAR topology of Fig. 1 help to systematize the choice of technique for daily work. This choice follows from the presented planning study and is in good agreement with previous studies.
In particular, Popple et al. [18] consider convex breast and mediastinal tumors which according to our classification fall into type II. The authors report that no more then 6–7 beams (11 beams for the most complicated head case) are needed for adequate PTV coverage/OAR sparing. Wang et al. [19] compare 5 beam BAO against 9 beam coplanar IMRT for two groups of paranasal sinus carcinoma cases, with one group characterized by the clear separation of PTV from OAR, typical for type II topology. This group profited most from 5 beam BAO, while the other group showed better results from conventional coplanar 9 beam IMRT treatment. Bangert et al. [7] report on a similar group of tumor localizations; a distinction like that of Wang is not made, however, the quality of their 11 coplanar beam plans is similar to that of optimized noncoplanar fivebeam technique. Note, that these results are in good agreement with our finding that a noncoplanar technique can be “caught up” by a conventional coplanar technique with large enough number of beams, if enough segments are allowed, see Fig. 6B. For the type III topology, Dong et al. [20] report optimal number of beams between 10 and 22 for an optimized quasiisotropic noncoplanar beam arrangement similar to our Q4π technique a. The high beam number probably reflects the need of steep dose gradients towards the surrounding tissue. However, type III was not the main focus of this work.
Perhaps due to principal optimization problems in BAO [21] and corresponding combinatorial problems of the concurring beam angle selection (BAS) method [7], only a limited number of sources is available. Mostly only low noncoplanar beam numbers were used (e.g. [22, 23]), below the range of our study.
For type I topology, a PTV partially enclosing an OAR, a sufficient number of beam directions (15 or more) is more important than the choice of technique. However, if such “good” OARsparing beam direction could be found, BAO becomes more important: for OAR structures that can be spared using selected beam directions (type II topology), a beam incidence along that direction is necessary; the total number of beams can then be reduced without considerable quality loss. If no “good” direction can be found and no BAO is available, for OARsensitive irradiations in the head the quasiisotropic technique with 16 beams in five table orientations is preferable; for Pinnacle3™ DMPO a sufficient number of segments has to be allowed for sequencing.
A technique with only one beam or half of beams applied perpendicular to the principal plane provides stable results of high quality for practical number of beams (9–16) and segments (under 120) potentially even in the case of trunk applications. The Co + 1 technique is recommended especially for Pinnacle DMPO if few segments are allowed. Note that only one additional table rotation is required for both techniques, and that the techniques can be partially applied as IMAT/VMAT. As more gantry angles improve target coverage and OAR sparing, a combination of noncoplanar technique with IMAT/VMAT seems promising.
Conclusions
The purpose of the present work is to study potentials of noncoplanar IMRT technique using schematic phantom cases which represent generalized OAR – PTV geometries. The work provides the basis for further planning studies with real clinical cases. The results, however, should not be applied to highly individual clinical cases without consideration.
In conclusion, if the extensive BAO is not available or meaningful, three types of PTVOAR topology introduced in the present study provide a good indicator for the choice between coplanar and noncoplanar beam arrangement and the choice of particular noncoplanar technique. The relative comparison of considered noncoplanar techniques to the reference coplanar technique is shown in Table 2. A good result (+ or ++) in a “Min” column means: If the best orientation of the beam set can be evaluated (using BAO), this technique is to be preferred with respect to the coplanar technique. A good result (+ or ++) in a “Max” column means: the most unfavourable result for the considered technique was better than the most unfavourable result for the coplanar reference technique; it could be preferable, if no BAO is available.
For pure type I problems (like A and C) in the head, ideally use a noncoplanar quasiisotropic technique a (Q4π) with 15 or more beams and many segments (100+) or sliding window technique. If less beams or less segments are allowed, it seems more practical to use d (Co + 1), the coplanar technique with additional orthogonal beam. For less segments also technique b (2P) with coplanar beams split regularly into two planes is an option. Techniques similar to b and d could also be applied in the trunk.
For type II topologies, also quasiisotropic noncoplanar technique a is recommended for the head, especially if no BAO can be performed or preferred beam directions are not obvious. For limited BAO computational time it could serve as starting constellation. A good choice is technique d, coplanar with one or several nearly orthogonal beams along potentially “good” axes which allow PTV coverage without OAR sparing to be compromised. The number of coplanar beams can then be reduced without appreciable quality loss. This technique can also be applied for trunk irradiations.
Type III topologies were not the topic of this work; however, for OARs encompassing a PTV a technique with minimum overlap between the beams should be preferred, to decrease maximum dose to the OAR, albeit at the expense of healthy tissue. This is again the case for the quasiisotropic noncoplanar technique a, which can also provide a good starting point for BAO (which is needed for local beam angle optimization [23]).
The future planning studies will concentrate on real clinical cases. The PTVOAR topologies will be grouped by their type (I, or type II, or III), the practical strategy for daily noncoplanar planning will be further developed, and the class solutions for standard clinical cases will be created. Planning studies are under way to test the presented recommendations.
Abbreviations
 2P:

Two orthogonal planes (planning technique)
 # beams:

Number of beams
 BAO:

Beam angle optimization
 BAS:

Beam angle selection
 CN:

Conformity number
 Co:

Coplanar equidistant (planning technique)
 Co + 1:

“Co” supplemented by a single orthogonal beam (planning technique)
 COV:

Composite objective value
 DMPO:

Direct machine parameter optimization
 DVH:

Dose volume histograms
 IMAT:

Intensity modulated arc therapy
 IMRT:

Intensity modulated radiation therapy
 MLC:

Multileafcollimator
 MU:

Monitor units
 nCOV:

Normalized composite objective value
 OAR:

Organ at risk
 PTV:

Planning target volume
 Q4π:

Noncoplanar quasiisotropic (planning technique)
 SBRT:

Stereotactic body radiotherapy
 S_{D} :

Quality index
 TPS:

Therapy planning system
 TV:

Target volume
 TV_{RI} :

Part of target volume surrounded by the reference isodose RI
 VMAT:

Volumetric arc therapy
 V_{RI} :

Total volume surrounded by the reference isodose RI
References
 1.
Brahme A, Roos JE, Lax I. Solution of an integral equation encountered in rotation therapy. Phys Med Biol. 1982;27:1221–9.
 2.
Bratengeier K. 2Step IMAT and 2Step IMRT: a geometrical approach. Med Phys. 2005;32:777–85.
 3.
Bortfeld T, Burkelbach J, Boesecke R, Schlegel W. Methods of image reconstruction from projections applied to conformation radiotherapy. Phys Med Biol. 1990;35:1423–34.
 4.
Kak AC, M S: Principles of Computerized Tomographic Imaging. New York: IEEE Press; 1999/1988
 5.
Bratengeier K, Seubert B, Holubyev K, Schachner H. Considerations on IMRT for quasiisotropic noncoplanar irradiation. Phys Med Biol. 2012;57:7303–15.
 6.
Breedveld S, Storchi PR, Voet PW, Heijmen BJ. iCycle: Integrated, multicriterial beam angle, and profile optimization for generation of coplanar and noncoplanar IMRT plans. Med Phys. 2012;39:951–63.
 7.
Bangert M, Ziegenhein P, Oelfke U. Characterizing the combinatorial beam angle selection problem. Phys Med Biol. 2012;57:6707–23.
 8.
Rossi L, Breedveld S, Heijmen BJ, Voet PW, Lanconelli N, Aluwini S. On the beam direction search space in computerized noncoplanar beam angle optimization for IMRTprostate SBRT. Phys Med Biol. 2012;57:5441–58.
 9.
Schreibmann E, Xing L. Feasibility study of beam orientation classsolutions for prostate IMRT. Med Phys. 2004;31:2863–70.
 10.
Buzug TM. Computed tomography  from photon statistics to modern conebeam CT. Berlin, Heidelberg: Springer; 2008.
 11.
Lin WT. Computationally efficient cone beam CT reconstruction algorithm using circleandline orbit. In: Medical imaging physics of medical imaging. 1999. p. 933.
 12.
Bratengeier K, Meyer J, Flentje M. Presegmented 2Step IMRT with subsequent direct machine parameter optimisation  a planning study. Radiat Oncol. 2008;3:38.
 13.
Bohsung J, Gillis S, Arrans R, Bakai A, De Wagter C, Knoos T, et al. IMRT treatment planning: a comparative intersystem and intercentre planning exercise of the ESTRO QUASIMODO group. Radiother Oncol. 2005;76:354–61.
 14.
van’t Riet A, Mak AC, Moerland MA, Elders LH, van der Zee W. A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. Int J Radiat Oncol Biol Phys. 1997;37:731–6.
 15.
Bratengeier K, Gainey MB, Flentje M. Fast IMRT by increasing the beam number and reducing the number of segments. Radiat Oncol. 2011;6:170.
 16.
Webb S. Optimisation of conformal radiotherapy dose distributions by simulated annealing [published erratum appears in Phys Med Biol 1990 Feb;35(2):297]. Phys Med Biol. 1989;34:1349–70.
 17.
Holmes T, Mackie TR. A filtered backprojection dose calculation method for inverse treatment planning. Med Phys. 1994;21:303–13.
 18.
Popple RA, Brezovich IA, Fiveash JB. Beam geometry selection using sequential beam addition. Med Phys. 2014;41:051713.
 19.
Wang X, Zhang X, Dong L, Liu H, Gillin M, Ahamad A, et al. Effectiveness of noncoplanar IMRT planning using a parallelized multiresolution beam angle optimization method for paranasal sinus carcinoma. Int J Radiat Oncol Biol Phys. 2005;63:594–601.
 20.
Dong P, Lee P, Ruan D, Long T, Romeijn E, Yang Y, et al. 4pi noncoplanar liver SBRT: a novel delivery technique. Int J Radiat Oncol Biol Phys. 2013;85:1360–6.
 21.
Meedt G, Alber M, Nusslin F. Noncoplanar beam direction optimization for intensitymodulated radiotherapy. Phys Med Biol. 2003;48:2999–3019.
 22.
Aleman DM, Romeijn HE, Dempsey JF. A response surface approach to beam orientation optimization in intensitymodulated radiation therapy treatment planning. INFORMS J Comput. 2009;21:62–76.
 23.
Craft D. Local beam angle optimization with linear programming and gradient search. Phys Med Biol. 2007;52:N127–35.
Acknowledgements
The publication of this work was funded by the German Research Foundation (DFG) and the University of Wuerzburg in the funding programme “Open Access Publishing”.
The authors thank B. Seubert for assistance and the Philips Corp. for supplying a research version of the Pinnacle3™ planning system.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
KB was responsible for the primary concept and the design of the study; he performed the calculations and drafted the manuscript. KH critically reviewed the study and revised the manuscript. Both authors read and approved the final manuscript.
Additional file
Additional file 1:
Typical fluences for noncoplanar and coplanar techniques. A: spherical OAR; B: cylindrical OAR; C: bananashaped OAR; a: Quasiisotropic technique Q4π (16 beams), b: beams arranged in two orthogonal planes 2P (15 beams); c: coplanar technique Co (15 beams); d: coplanar + 1 noncoplanar beam Co + 1 (16 beams); o: table angle = 0°^ gantry angle = 0°; grey thick bordered box: table angle = 90°; black thin bordered box: table angle ≠ 90° ^ ≠ 0°. (PNG 10614 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Bratengeier, K., Holubyev, K. Characteristics of noncoplanar IMRT in the presence of targetembedded organs at risk. Radiat Oncol 10, 207 (2015). https://doi.org/10.1186/s1301401504945
Received:
Accepted:
Published:
Keywords
 Planning Target Volume
 Intensity Modulate Radiation Therapy
 Monitor Unit
 Plan Quality
 Planning Target Volume Coverage