Email updates

Keep up to date with the latest news and content from Radiation Oncology and BioMed Central.

Open Access Highly Accessed Research

Scanned carbon beam irradiation of moving films: comparison of measured and calculated response

Christoph Bert1*, Daniel Richter12, Marco Durante12 and Eike Rietzel13

Author Affiliations

1 GSI Helmholtzzentrum für Schwerionenforschung, Planckstr. 1, 64291 Darmstadt, Germany

2 TU Darmstadt, Institut für Festkörperphysik, Hochschulstr. 6, 64289 Darmstadt, Germany

3 Siemens AG, Healthcare Sector, Imaging & Therapy, Particle Therapy, Hofmannstr. 26, 91052 Erlangen, Germany

For all author emails, please log on.

Radiation Oncology 2012, 7:55  doi:10.1186/1748-717X-7-55

The electronic version of this article is the complete one and can be found online at: http://www.ro-journal.com/content/7/1/55


Received:8 November 2011
Accepted:2 April 2012
Published:2 April 2012

© 2012 Bert et al; licensee 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.

Abstract

Background

Treatment of moving target volumes with scanned particle beams benefits from treatment planning that includes the time domain (4D). Part of 4D treatment planning is calculation of the expected result. These calculation codes should be verified against suitable measurements. We performed simulations and measurements to validate calculation of the film response in the presence of target motion.

Methods

All calculations were performed with GSI's treatment planning system TRiP. Interplay patterns between scanned particle beams and moving film detectors are very sensitive to slight deviations of the assumed motion parameters and therefore ideally suited to validate 4D calculations. In total, 14 film motion parameter combinations with lateral motion amplitudes of 8, 15, and 20 mm and 4 combinations for lateral motion including range changes were used. Experimental and calculated film responses were compared by relative difference, mean deviation in two regions-of-interest, as well as line profiles.

Results

Irradiations of stationary films resulted in a mean relative difference of -1.52% ± 2.06% of measured and calculated responses. In comparison to this reference result, measurements with translational film motion resulted in a mean difference of -0.92% ± 1.30%. In case of irradiations incorporating range changes with a stack of 5 films as detector the deviations increased to -6.4 ± 2.6% (-10.3 ± 9.0% if film in distal fall-off is included) in comparison to -3.6% ± 2.5% (-13.5% ± 19.9% including the distal film) for the stationary irradiation. Furthermore, the comparison of line profiles of 4D calculations and experimental data showed only slight deviations at the borders of the irradiated area. The comparisons of pure lateral motion were used to determine the number of motion states that are required for 4D calculations depending on the motion amplitude. 6 motion states per 10 mm motion amplitude are sufficient to calculate the film response in the presence of motion.

Conclusions

By comparison to experimental data, the 4D extension of GSI's treatment planning system TRiP has been successfully validated for film response calculations in the presence of target motion within the accuracy limitation given by film-based dosimetry.

Background

Radiotherapy of tumors influenced by organ motion benefits from dedicated treatment planning. Our treatment planning system for scanned particle beams (Treatment planning for Particles, TRiP) [1] was extended to include 4D functionality, i.e. time resolved dose calculation and treatment plan optimization [2]. The purpose of this contribution is to report measurements and corresponding film response calculations which were performed to confirm 4D calculations with TRiP. Such 4D calculations can, e.g., be used in treatment planning for patients with a moving tumor to check the robustness of the optimized 4D treatment plan by multiple calculations with simulated, time resolved, data of tumor motion and beam delivery.

Experiments were performed with a raster-scanned carbon beam [3] at GSI Helmholtz Centre for Heavy Ion Research in Darmstadt, Germany. Dose distributions delivered by scanned pencil beams are conformal to the target volume and allow for homogeneous coverage. But dose distributions are also sensitive to intrafractional organ motion because the combination of beam scanning and target motion represents a double dynamic system and typically results in interplay. Interplay causes deterioration of the dose distribution leading to hot and cold spots [4-6]. With techniques like gating [7], rescanning [4], or beam tracking [8,9], treatment of tumors that are subject to intrafractional organ motion like respiration with scanned particle beams can potentially be improved.

For a specific treatment plan, the interplay pattern of scanned beam and target motion is influenced by the target motion parameters (period, initial phase, amplitude) and the intensity controlled beam scanning process, especially scan speed which primarily depends on the particle extraction from the synchrotron [5]. 4D calculations incorporate these parameters to determine the interplay pattern or mitigated interplay pattern in case of beam tracking, gating, or rescanning. In a clinical setting, multiple calculations with varying motion and extraction parameters and/or mitigation techniques could be performed as part of treatment planning to estimate the range of possible outcomes. It is thus very important, that these calculations are valid. Their validity can be assessed by irradiating without motion mitigation, i.e. by producing interplay patterns. Validity is achieved if the calculated outcome is within the achievable experimental accuracy. We performed experiments with moving radiographic films as detector without using a motion mitigation technique. Radiographic films have a very high lateral spatial resolution and are thus ideally suited to detect interplay patterns. For measurements in a 3D irradiation field, multiple films can be stacked in depth alternated with spacers such as Lucite plates. In such 3D measurements the agreement of measured and calculated film blackening deviates between 10% and 30% due to the mixed radiation field of a carbon beam [10].

In order to allow a detailed study of the potential deviations between measurement and calculation we intentionally limited ourselves to measuring the film response and left-right target motion with depth-modulation by a stationary absorber positioned proximal of the target volume. Non-rigid motion that will be present in patients was not introduced on purpose. We designed the treatment plans in accordance with film (stacks) as target. These planar targets lead to irradiations with (multiple) pristine Bragg-peaks rather than treatment plans with a spread-out Bragg-peak (SOBP). In an SOBP irradiation, mitigation of interplay effects to some extent takes place intrinsically, e.g., due to rescanning effects in proximal target regions [11]. Thus, cause of deviations cannot be precisely studied. Since SOBPs are essential in clinical use, appropriate investigations have been performed by Gemmel et al. using Chinese hamster ovary cells to also incorporate the biological effect of carbon beams [12]. In addition, irradiations to a 3D geometry with an array of ionization chambers as detector are currently ongoing. Apart from comparing measured with calculated results that will validate the calculation we studied the impact of the number of motion states used for 4D calculation on the precision of the calculation in dependence on the motion amplitude.

Methods

Experimental methods

Setup

A schematic drawing of the setup is shown in Figure 1. The irradiations were performed in the patient treatment room at GSI. At isocenter, a periodically moving table was installed. The table moved radiographic films (Kodak X-Omat V, Kodak GmbH, Stuttgart) as detectors in left-right direction in beam's eye view. Motion amplitude, period, and initial motion states were controlled and recorded. In a first set of experiments, the effects of lateral target motion were investigated: Films were positioned free in air, no absorber material was used (Figure 1a). A second series of experiments was performed to assess range effects in the presence of target motion. An absorber with four different thicknesses was positioned stationary proximal to the motion table (Figure 1b). A stack of five films between absorber materials was used for these measurements (Figure 1b, c).

thumbnailFigure 1. Schematic drawing of the experimental setup. The particle beam was accelerated by a synchrotron and scanned over the moving radiographic film with the scanning magnets. a) For experiments with translational motion the films were positioned free in air moving left-right in beam's eye view. The motion was recorded in temporal correlation to the dose delivery by a calibrated CCD camera that monitored an infrared LED attached to the sliding table. b) In a second series of experiments an absorber (PMMA) with 4 different thicknesses that was positioned proximal to the moving table introduced range changes. As detector a stack of 5 films was used. Their positions relative to the Bragg peak is indicated in c).

To measure the motion of the moving table, a calibrated CCD camera monitored the movement of an IR-LED attached to the table. Data were acquired on a PC running LabVIEW (National Instruments Germany GmbH, LabVIEW, Version 7.1, Munich, Germany) with the Image Vision package (National Instruments Germany GmbH, Image Vision, Munich, Germany). The beam extraction pattern of the synchrotron was determined temporally correlated by measuring the irradiation time of each beam position during raster scanning.

Treatment plan and treatment delivery

Different treatment plans were used for the two experimental setups. For left-right motion without range modulation, a single iso-energy slice was used (carbon beam, 272 MeV/u). A quadratic area of 11 × 11 cm2 (50% iso-dose) was irradiated to ~0.9 Gy, raster points (7.7 mm FWHM pencil beams) were placed on a regular 2 mm grid (3013 beam positions), and scanned intensity-controlled [3] line-by-line with the faster scanning direction horizontally (x-axis, parallel to target motion). For the second setup including range changes, the treatment plan was based on a CT scan of the absorber and optimized such that the beam stops before the 5th film in the stack (see Figure 1c) for each absorber thickness. This setup results into four main beam-energies of 148, 155, 163, and 181 MeV/u for the four segments of different absorber thickness of 0, 5, 10, and 20 mm PMMA, respectively. The plan was optimized on a 2 mm grid (2179 beam positions), at 6.5 mm FWHM beam spot size to a lateral area of 72 × 72 mm2 (50% iso-dose).

Target motion was started in correlation with the irradiation. Before delivery of the dose pattern to a moving film, two cross-shaped or quadratic patterns were irradiated on the static film to define a precise 2D coordinate system. By means of the coordinate system the results of the experiments can be related to the results of the calculation.

Calculation of the expected film response

Calculation of the expected dose distribution in 4D, i.e. in the presence of organ motion was introduced shortly after 4DCT acquisitions became available [13,14]. The initial work focused on 3D conformal photon beam therapy and proton beam delivery with a passively shaped beam. For 4D calculation for scanned carbon ion beam therapy, TRiP was extended by Bert & Rietzel [2]. The basic principles of TRiP4D follow the ideas of Rietzel et al. and Keall et al. and are now also used to calculate the film response in the presence of motion. A brief summary of the underlying strategy follows.

Dose and film response calculations are based on:

i) the parameters of the treatment plan (beam position, beam energy, beam width, particles per beam position) and especially the scanning order of the beam positions,

ii) amplitude based motion states describing the time dependent target positions,

iii) translation vectors from all motion states to a reference motion state,

iv) the measured target motion data, and

v) scanning progress given by the measured irradiation time of each raster point.

The underlying principle is similar to 4DCT data acquisition: Based on (iv) target motion data and (v) scanning progress each rasterpoint of the delivered treatment plan (i) can be attributed to one of the motion states (ii). The motion amplitude detected with the CCD camera was used to determine the corresponding motion state. The peak-to-peak amplitude was divided into equally spaced motion states that correlate to specific positions of the target volume. Corresponding to each motion state, a sub-treatment plan containing the subset of raster points delivered during this motion state was constructed. Sub-treatment plans are used in combination with the corresponding motion state as quasi-static treatment plans and result in quasi-static sub-dose distributions (or sub-film responses). Before summation of all sub-dose distributions to derive the overall dose distribution, each of them is transformed to a reference motion state using the translation vectors (iii).

For calculation of film responses, the non-linear dose-response [15] of films has to be considered. The background-corrected relative blackening S/S0 of radiographic films can be described by

<a onClick="popup('http://www.ro-journal.com/content/7/1/55/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.ro-journal.com/content/7/1/55/mathml/M1">View MathML</a>

(1)

where S0 is the experimentally determined saturation level, m the slope of the blackening curve, and D the used dose-level. The summed relative blackening of the film Stot/S0 is not the sum of the transformed sub-film responses Si/S0 because the non-linear response Stot/S0 is determined as:

<a onClick="popup('http://www.ro-journal.com/content/7/1/55/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.ro-journal.com/content/7/1/55/mathml/M2">View MathML</a>

(2)

where NCT is the number of motion states and thus sub-film response distributions. The saturation level S0 is determined experimentally according to Spielberger et al. [10]. The slope m is particle and energy dependent and belongs to the base data which are calculated based on the slope of a 60Co blackening curve [16]. For these calculations, analysis of the 3D measurement data reported by Spielberger et al. results in an expected accuracy of 10-30% when comparing to measured data [10]. For the experiments with the film stack (Figure 1b), calculated data are plotted with an error bar of 10% to address the accuracy of the modeled film responses in a mixed radiation field.

For the experiments without absorbers (Figure 1a) fragmentation and energy degradation are negligible and m can be measured for the energy of the primary carbon beam. For this measurement, a single film is irradiated with several independent areas (e.g., multiple squares of 5 × 5 cm2) of different dose levels. For each area, S is determined with a densitometer and corrected by the background measured in unirradiated parts of the film. Using equation 1, m and S0 as well as their errors Δm and ΔS0 are then determined by a fitting algorithm. We determined the error of the calculated relative film response Stot/S0 (equation 2) by error propagation and show it as error bars on line profiles.

To transform the sub-dose distributions of the different motion states to the reference motion state, a one-dimensional translation vector was sufficient. To determine the required number of motion states in dependence on the motion amplitude, film response distributions were calculated for 2, 3, ..., 30, 31 motion states for each experimental data set.

Data analysis

Film response

The irradiated films were developed (developing machine M35 with developer DX31 and fixer FX31, Kodak GmbH, Stuttgart) and digitized (FIPS Plus LS75, PTW, Freiburg, Germany) in a 1 mm grid as reported by Spielberger et al. [10]. The experimental error of the relative blackening S/S0 was determined in error propagation with ΔS = 0.02 as stated in the manual of the densitometer. ΔS0 was determined experimentally as described in 2.2. To allow comparison between measured and calculated film responses, each image pair was aligned according to its cross-shaped 2D coordinate system including translations and rotation.

Quantitative comparisons were performed as follows:

i) For each pixel the difference between experimental and calculated film response was calculated and plotted as [(S/S0)experiment-(S/S0)calculated]/(S/S0)experiment.

ii) Horizontal and vertical line profiles were extracted and plotted for experimental image, calculated image, and difference image including error bars as described above.

iii) The mean difference (± standard deviation, SD) was calculated based on a difference image. This analysis was performed for two regions-of-interest (ROI): The larger ROI, labeled ROIirradiated, was defined by all pixels with a calculated film response (S/S0)calculated > 0.1; the smaller ROI, labeled ROItarget, included the area of the stationary irradiation only (dashed line in Figure 2a for the pure translational motion).

thumbnailFigure 2. Results of all parameter combinations. Measured (A-P) and calculated (a-p) film response distributions for the different parameter combinations (amplitude, initial phase, period). In A) the dashed square indicates the small region-of-interest used in the analysis (ROItarget).

Determination of the number of motion states

For each experimental data set without absorbers (Figure 1a), film response distributions were calculated using 2-31 motion states, i.e. 3D "time-frames" in the 4DCT. The required number of motion states was analyzed qualitatively by visual inspection of the resulting saturation curve. To determine the required number of motion states for a given amplitude quantitatively, the data were fitted with a standard saturation function: d(p) = a · exp(-b · p) + c with d relative film response difference and p number of motion states. The fitting parameters a, b, and c represent the initial difference, the saturation rate, and the saturation level of the film response difference at a quasi-continuous number of motion states, respectively. The required number of motion states per amplitude pr was defined at 5% deviation from saturation level c. The required number of motion states pr was calculated for each measurement. Mean pr values per motion amplitude as well as individual pr values with error bars that result from error propagation using the fitting uncertainties of a, b, and c are reported.

Results

Film response calibration

For the experiments without absorbers, analysis of the film response in dependence on the deposited dose resulted in a slope m of the blackening curve of mC-12, E = 272 MeV/u = 0.33/Gy with a fitting error of Δm = 0.0125/Gy. The film response in saturation S0 was 4.0 with a fitting error of ΔS0 = 0.04. These values were used in all subsequent film response calculations.

Lateral target motion

Stationary reference irradiation

For assessment of the film response calculation accuracy for lateral target motion, the setup displayed in Figure 1a was used. Figure 3 shows the results for irradiating a stationary film as a reference for irradiations with target motion. The subfigures show the measured and the calculated film responses, the difference image between calculation and measurement, as well as a horizontal and a vertical line profile. The differences between calculated and measured film responses are -1.52% ± 2.06% inside the ROItarget and 2.3% ± 4.92% inside the ROIirradiated. The largest deviations are observed at the lateral borders in the dose gradients. The profiles demonstrate that these effects are not due to mispositioning.

thumbnailFigure 3. Film response for stationary irradiation. Results for the stationary film that were used as reference. The image shows the 2D experimental film response distribution, the calculated distribution, the relative difference between experimental and calculated film response, as well as profiles along the dashed lines in the difference distribution.

Overview of experiments with lateral motion

Fourteen motion parameter combinations were studied. Amplitudes ranged from 8-20 mm, periods from 3-7 s, and initial motion starting phases from 0-270°. Figure 2 displays experimental as well as calculated film responses (interplay patterns) for the studied motion parameter combinations. Visual inspection yields good agreement between calculation and experiment. Details of the statistical analysis within the two ROIs are provided in Table 1.

Table 1. Deviation analysis - single film, lateral motion

Required number of motion states

The dependence of the mean relative difference between measurement and calculation inside the ROIirradiated on the number of motion states used for calculations is plotted in Figure 4a grouped for the amplitudes of 8, 15, and 20 mm. All curves plateau above film response calculations with 10-15 motion states. Figure 4b shows the result of the pr calculation, i.e. the required number of motion states to reach 95% of the plateau value. Despite the comparably large fluctuations at 20 mm amplitude, there is a linear relationship between amplitude a and the required number of motion states: pr = 0.58 · a - 0.05. For each 10 mm motion amplitude, ~6 motion states are thus sufficient for dose calculations. If the desired level of agreement is increased to 98% (99%) of the saturation value, the required number of motion states per 10 mm target motion increases to 7 (8).

thumbnailFigure 4. Influence of the number of motion states. a) Influence of the number of motion states on the relative difference between calculated and experimental film response distribution. Data are grouped according to the motion amplitude. The dashed lines are added to guide the eye. b) The required number of motion states pr was determined for each measured parameter combination at 8, 15, and 20 mm motion amplitude a. A linear fit based on the mean of each measured motion amplitude yielded pr = 0.58 · a - 0.05.

Validation of calculated film response

To exclude influences of the number of motion states pr on the presented results, calculations for the following analyses were performed with 20 motion states. The results of difference image analysis (mean relative difference) in the two ROIs are given in Table 1 for all 14 image pairs. Agreement between experimental and calculated film response is higher in the smaller ROItarget (-0.92% ± 1.30%) than in the larger ROIirradiated (-4.52% ± 2.39%). In general, agreement is lower for the four distributions at 20 mm motion amplitude.

Figures 5 and 6 show examples for a motion parameter combination with good agreement between measurement and calculation in the statistical analysis (Figure 5) and a parameter combination with lower agreement (Figure 6). Both images show the experimental and calculated film response distributions, the relative differences to calculations, as well as horizontal and vertical line profiles.

thumbnailFigure 5. Specific example: good agreement. Example of good agreement between experimental and calculated film response (Figure 2e). Apart from the distributions of the experimental response, the calculated response, and the relative difference, horizontal and vertical profiles are shown along the dashed lines indicated in the distribution of the relative difference. Details of the analysis are reported in Table 1. δ indicates the vertical shift of the line profiles in order to show all of them in the same figure.

thumbnailFigure 6. Specific example: less good agreement. Example of less good agreement between experimental and calculated film response (Figure 2o). Apart from the distributions of the experimental response, the calculated response, and the relative difference, horizontal and vertical profiles are shown along the dashed lines indicated in the distribution of the relative difference. Details of the analysis are reported in Table 1. δ indicates the vertical shift of the line profiles in order to show all of them in the same figure.

For both examples, the line profiles indicate that the calculation reproduced the measured interplay pattern well. While in Figure 5 the deviations are negligible, the example with lower agreement in Figure 6 reveals deviations mainly at the transitions from irradiated to unirradiated areas and in the fine-structure of the interplay pattern (e.g. the "cross" centered at x = 125 mm/y = 125 mm or the area 50 mm < y < 75 mm).

Lateral target motion including range changes

For assessment of combined lateral and longitudinal target motion, four different experiments each with a stack of 5 films were performed in comparison to a stationary experiment with the same setup. Motion parameters of the individual experiments are listed in Table 2.

Table 2. Motion parameters of film stack experiments

Figures 7 and 8 show the results of these measurements that were performed with the setup shown in Figure 1b. Figure 7 shows the results for film 3 of all measurements (stationary reference and four different motion parameter combinations). The direct comparison of calculated and experimental film response distributions shows good agreement. For the right-top part of the pattern (bolus: 20 mm, see Figure 1c), there are deviations of > 10% between experimental data and calculation even for irradiation of a stationary film (setup in Figure 1a). For experiments with target motion, the calculations reproduced the distinct interplay patterns. Comparison of the two horizontal profiles yields similar results but the profiles also reveal that positions of some interplay-induced peaks slightly deviate for some of the experiments (e.g. NO2 and NO5, y = 55 mm).

thumbnailFigure 7. Film stack measurements - all parameter combinations. Results of the film stack measurements that included ranges. Displayed are the data for film 3 for the stationary irradiation as well as the 4 irradiations that included lateral motion. The distributions show experimental response, calculated response, their relative difference, as well as profiles along the horizontal lines indicated in the experimental distributions. Details of the analysis are reported in Table 3.

thumbnailFigure 8. Film stack measurements - experiment NO2. Results of the film stack measurements that included ranges. Displayed are the data of all films for experiment NO2 (see Table 3 and Figure 7). The distributions show experimental response, calculated response, their relative difference, as well as profiles along the horizontal lines indicated in the experimental distributions. Details of the analysis are reported in Table 3. Note that the profiles for film 5 are scaled by a factor of 4.

In Figure 8, the responses of all 5 films in the stack for measurement NO2 of Figure 7 are shown. The calculations reproduce the experimental film responses well. Interference patterns as well as height of the peaks that is dominated by the depth of the films in the stack are modeled within the experimental errors with the exemption of some peak positions.

Analysis of the difference images shows that the deviations between experiments and calculations are within 10% for most of the irradiated area (green regions in Figure 8). Comparable to the stationary setup, the top-right part of the irradiations with target motion indicates slightly higher deviations. The differences in some peak positions results in blue or red bands in the difference images that are also prominent at some transition areas from irradiated to unirradiated regions as for the results for lateral motion only (section 3.2).

Table 3 shows the statistical data of the pixel-by-pixel comparison between measured and calculated data. The mean relative deviation in ROIirradiated is -11.1% (SD 11.7%) for ROItarget the deviations are -10.3% (SD 9%) in comparison to -13.1% (SD 20.9%) and -13.5% (SD 19.9%) for the stationary irradiation. If film #5 that is positioned in the distal fall-off region is omitted from the analysis the values change to -6.0% (SD 2.5%)/-6.4% (SD 2.6%) in comparison to -2.7% (SD 3.3%)/-3.6% (SD 2.5%) for the stationary irradiation for ROIirradiated/ROItarget.

Table 3. Deviation analysis - film stack

Discussion

Experiments with moving radiographic films were performed with the intention to validate the calculations of our 4D treatment planning system in the presence of translational target motion. For non-rigid motion patterns the calculation routines will be identical but non-rigid registration routines need to be incorporated and validated which was not part of this study. Calculation routines are in principle identical for rescanning [4], beam gating [7,17,18], and also beam tracking [8,9,19]. Since the calculation of film responses relates on the calculation of the delivered dose, the presented results therefore form the basis of 4D dose calculation precision for motion mitigated irradiation schemes. A full treatment plan validation needs to incorporate optimization parameters that depend on the motion mitigation technique and should be complemented with, e.g., ionization chamber data that have a higher accuracy than radiographic films at a lower spatial resolution. The validation of 4D optimization parameters such as compensation vectors for beam tracking, the size of the gating window in gated irradiations, or the number of required rescans as well as aspects of biological treatment planning have not been investigated.

Irradiations of moving detectors without motion mitigation were chosen for experimental validation because resulting dose distributions provide prominent interplay patterns and are very sensitive to slight parameter changes and thus also to slight mismatches in the calculation. To measure interplay patterns accurately, a detector with high spatial resolution is required. Radiographic films are ideally suited for this purpose despite their dependence on the irradiation field and their limited accuracy. For mixed irradiation fields that are, e.g., given in our experimental setup (Figure 1b) for combined lateral and longitudinal motion, analysis of the data reported by Spielberger et al. yields an agreement between measurement and calculation of ~10% in the target regions and within 30% at the border of the irradiation field for stationary irradiations [10]. To exclude the dependency on the composition of the irradiation field, we also performed experiments in which the films were irradiated with a single particle energy in the plateau region of the Bragg peak (setup shown in Figure 1a). This irradiation scheme allows a calibration of the detector as described in section 2.2.1., which reduced the uncertainty at least in the central parts of the irradiation field (stationary irradiation in Figure 3). For the borders, i.e. the transitions from irradiated to unirradiated areas, the precision of the detector is mainly determined by the beam shape. Imperfections such as deviations from the nominal Gaussian beam shape lead to deviations that were visible as the blue colored "margin" (rel. difference > 10%) in Figure 3, 5, and 6. In the central part of the target volume these deviations from the nominal beam profile lead to stripe patterns in the case of stationary beam delivery (Figure 3) with an acceptable dosimetric impact because neighboring beam spots have a sufficiently large overlap. In case of a moving detector, interplay results in a decreased overlap and the dose contribution of a single beam position has a higher impact. In the experiments with the mixed irradiation field, the necessity of a precise beam shape is in general less prominent since there is additional dose contribution from neighboring iso-energy slices. However, also Spielberger et al. reported increased deviations of up to 30% at the borders of the irradiation field [10]. Due to these slight beam shape distortions but also due to the uncertainties in the extraction rate and motion trajectory measurements the deviations in calculated film response distributions in the presence of motion are expected to be higher than for irradiations of stationary detectors.

With increasing motion amplitude interplay effects increase and therefore areas with transitions from irradiated to unirradiated increase. This explains the larger deviations between measured and calculated film response distributions for 20 mm motion amplitude. The effect is also visible by comparing the results of the lateral target motion for the large ROI that includes the complete film response distribution and the small ROI which encompasses only the central region of the films. Agreement is better within the small ROI: maximum mean deviation -2.94% and a maximum SD of 17.02% in comparison to a maximum mean deviation of 9.84% (SD 68.98%) in the large ROI. For the lateral motion including range changes the deviations between the two ROIs is less pronounced.

A comparison between experiments in the plateau region of the Bragg peak with only lateral motion and the film stack experiments in the Bragg peak region that included range changes shows an increased deviation between calculation and experiment for the film stack measurements. This is most likely due to the additional sensitivity of the calculations to the exact description of the experimental setup. Already in the stationary experiment, the deviations are increased especially for film #5 that is positioned in the distal fall-off region of the Bragg peak which is most sensitive to position deviations due to the sharp gradient. Similar arguments are valid for films 1-4 since we did not irradiate a spread-out Bragg peak but a pristine peak. The most likely reasons for position deviations are uncertainties in the bolus thicknesses and the range precision in the treatment planning process.

The determined number of 6 required motion states per 10 mm motion amplitude is only valid for the distal slice of an extended target volume. In more proximal regions intrinsic averaging effects, namely the dose contributions from irradiation of more distal slices in proximal slices, lead to reduced sensitivity to target motion [18]. If 6 states per 10 mm motion amplitude are sufficient for 4D calculation as well as optimization, standard 4DCT protocols with ~10 reconstructed motion states will be sufficient for typical intra-fractional motion amplitudes of, e.g., lung tumors [20,21].

Conclusion

GSI's 4D treatment planning system TRiP [1,2] was validated by comparing calculated film responses to experimental data for translational motion geometries. Experimental and calculated film responses for lateral motion without range changes show good agreement of -0.92 ± 1.3% in an region of interest covering the target area in comparison to < 10% reported by [10] for stationary targets. If range changes are introduced, the agreement is still given with a deviation of -6.4 ± 2.6% (-10.3 ± 9.0% if film #5 in distal fall-off is included) in comparison to -3.6% (-13.5% including film #5) for the stationary irradiation. By performing calculations with varying number of motion states 6 motion states per 10 mm motion amplitude were determined to be sufficient.

Competing interests

GSI Biophysics currently receives funding from Siemens AG, Healthcare Sector, Imaging & Therapy, Particle Therapy. ER currently is an employee of Siemens AG, but most of the work in the presented manuscript was performed prior he entered Siemens AG.

Authors' contributions

CB, ER: design of study, CB data acquisition and majority of data analysis, all: data analysis and manuscript generation. All authors read and approved the final manuscript.

Acknowledgements

The authors thank Prof. Dr. Dr. h.c. Gerhard Kraft for fruitful discussions and advices regarding this project. We further acknowledge the help of Drs. Nikolaus Kurz and Dieter Schardt in setting up the data acquisition system.

References

  1. Krämer M, Jäkel O, Haberer T, Kraft G, Schardt D, Weber U: Treatment planning for heavy-ion radiotherapy: physical beam model and dose optimization.

    Phys Med Biol 2000, 45:3299-3317. PubMed Abstract | Publisher Full Text OpenURL

  2. Bert C, Rietzel E: 4D treatment planning for scanned ion beams.

    Radiat Oncol 2007, 2:24. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  3. Haberer T, Becher W, Schardt D, Kraft G: Magnetic scanning system for heavy ion therapy.

    Nucl Instrum Meth A 1993, 330:296-305. Publisher Full Text OpenURL

  4. Phillips MH, Pedroni E, Blattmann H, Boehringer T, Coray A, Scheib S: Effects of respiratory motion on dose uniformity with a charged particle scanning method.

    Phys Med Biol 1992, 37:223-233. PubMed Abstract | Publisher Full Text OpenURL

  5. Bert C, Grözinger SO, Rietzel E: Quantification of interplay effects of scanned particle beams and moving targets.

    Phys Med Biol 2008, 53:2253-2265. PubMed Abstract | Publisher Full Text OpenURL

  6. Furukawa T, Inaniwa T, Sato S, Shirai T, Mori S, Takeshita E, et al.: Moving target irradiation with fast rescanning and gating in particle therapy.

    Med Phys 2010, 37:4874-4879. PubMed Abstract | Publisher Full Text OpenURL

  7. Minohara S, Kanai T, Endo M, Noda K, Kanazawa M: Respiratory gated irradiation system for heavy-ion radiotherapy.

    Int J Radiat Oncol 2000, 47:1097-1103. Publisher Full Text OpenURL

  8. Keall PJ, Kini VR, Vedam SS, Mohan R: Motion adaptive x-ray therapy: a feasibility study.

    Phys Med Biol 2001, 46:1-10. PubMed Abstract | Publisher Full Text OpenURL

  9. Bert C, Saito N, Schmidt A, Chaudhri N, Schardt D, Rietzel E: Target motion tracking with a scanned particle beam.

    Med Phys 2007, 34:4768-4771. PubMed Abstract | Publisher Full Text OpenURL

  10. Spielberger B, Krämer M, Scholz M, Kraft G: Three-dimensional dose verification in complex particle radiaion fields based on X-ray films.

    Nucl Instrum Meth B 2003, 209:277-282. OpenURL

  11. Rietzel E, Bert C: Respiratory motion management in particle therapy.

    Med Phys 2010, 37:449-460. PubMed Abstract | Publisher Full Text OpenURL

  12. Gemmel A, Rietzel E, Kraft G, Durante M, Bert C: Calculation and experimental verification of the RBE-weighted dose for scanned ion beams in the presence of target motion.

    Phys Med Biol 2011, 56:7337-7351. PubMed Abstract | Publisher Full Text OpenURL

  13. Rietzel E, Chen GTY, Choi NC, Willet CG: Four-dimensional image-based treatment planning: Target volume segmentation and dose calculation in the presence of respiratory motion.

    International Journal of Radiation Oncology*Biology*Physics 2005, 61:1535-1550. Publisher Full Text OpenURL

  14. Keall PJ: 4-dimensional computed tomography imaging and treatment planning.

    Semin Radiat Oncol 2004, 14:81-90. PubMed Abstract | Publisher Full Text OpenURL

  15. Spielberger B, Scholz M, Krämer M, Kraft G: Calculation of the x-ray film response to heavy charged particle irradiation.

    Phys Med Biol 2002, 47:4107-4120. PubMed Abstract | Publisher Full Text OpenURL

  16. Spielberger B, Scholz M, Krämer M, Kraft G: Experimental investigations of the response of films to heavy-ion irradiation.

    Phys Med Biol 2001, 46:2889-2897. PubMed Abstract | Publisher Full Text OpenURL

  17. Kubo HD, Hill BC: Respiration gated radiotherapy treatment: a technical study.

    Phys Med Biol 1996, 41:83-91. PubMed Abstract | Publisher Full Text OpenURL

  18. Bert C, Gemmel A, Saito N, Rietzel E: Gated irradiation with scanned particle beams.

    Int J Radiat Oncol 2009, 73:1270-1275. Publisher Full Text OpenURL

  19. Grözinger SO, Li Q, Rietzel E, Haberer T, Kraft G: 3D online compensation of target motion with scanned particle beam.

    Radiother Oncol 2004, 73(Suppl 2):S77-S79. PubMed Abstract OpenURL

  20. Erridge SC, Seppenwoolde Y, Muller SH, Jaeger KD, Belderbos JSA, Boersma LJ, et al.: Portal imaging to assess set-up errors, tumor motion and tumor shrinkage during conformal radiotherapy of non-small cell lung cancer.

    Radiother Oncol 2003, 66:75-85. PubMed Abstract | Publisher Full Text OpenURL

  21. Seppenwoolde Y, Shirato H, Kitamura K, Shimizu S, Lebesque JV, Miyasaka K: Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy.

    Int J Radiat Oncol 2002, 53:822-834. Publisher Full Text OpenURL