B Poppe

Two‐dimensional ionization chamber arrays for IMRT plan verification

In this paper we describe a concept for dosimetric treatment plan verification using two-dimensional ionization chamber arrays. Two different versions of the 2D-ARRAY (PTW-Freiburg, Germany) will be presented, a matrix of 16 x 16 chambers (chamber cross section 8 mm x 8 mm; the distance between chamber centers, 16 mm) and a matrix of 27 x 27 chambers (chamber cross section 5 mm x 5 mm; the distance between chamber centers is 10 mm). The two-dimensional response function of a single chamber is experimentally determined by scanning it with a slit beam. For dosimetric plan verification, the expected two-dimensional distribution of the array signals is calculated via convolution of the planned dose distribution, obtained from the treatment planning system, with the two-dimensional response function of a single chamber. By comparing the measured two-dimensional distribution of the array signals with the expected one, a distribution of deviations is obtained that can be subjected to verification criteria, such as the gamma index criterion. As an example, this verification method is discussed for one sequence of an IMRT plan. The error detection capability is demonstrated in a case study. Both versions of two-dimensional ionization chamber arrays, together with the developed treatment plan verification strategy, have been found to provide a suitable and easy-to-handle quality assurance instrument for IMRT.

DAVID - a translucent multi-wire transmission ionization chamber for in-vivo verification of IMRT and conformal irradiation techniques

Permanent in vivo verification of IMRT photon beam profiles by a radiation detector with spatial resolution, positioned on the radiation entrance side of the patient, has not been clinically available so far. In this work we present the DAVID system, which is able to perform this quality assurance measurement while the patient is treated. The DAVID system is a flat, multi-wire transmission-type ionization chamber, placed in the accessory holder of the linear accelerator and constructed from translucent materials in order not to interfere with the light field. Each detection wire of the chamber is positioned exactly in the projection line of a MLC leaf pair, and the signal of each wire is proportional to the line integral of the ionization density along this wire. Thereby, each measurement channel essentially presents the line integral of the ionization density over the opening width of the associated leaf pair. The sum of all wire signals is a measure of the dose-area product of the transmitted photon beam and of the total radiant energy administered to the patient. After the dosimetric verification of an IMRT plan, the values measured by the DAVID system are stored as reference values. During daily treatment the signals are re-measured and compared to the reference values. A warning is output if there is a deviation beyond a threshold. The error detection capability is a leaf position error of less than 1 mm for an isocentric 1 cm x 1 cm field, and of 1 mm for an isocentric 20 cm x 20 cm field.

The artefacts of radiochromic film dosimetry with flatbed scanners and their causation by light scattering from radiation-induced polymers.

Optical experiments and theoretical considerations have been undertaken in order to understand the causes of the orientation effect and the parabola effect, the artefacts impairing the desired light absorption measurement on radiochromic EBT3 films with flatbed scanners. EBT3 films exposed to doses up to 20.9 Gy were scanned with an Epson Expression 10000XL flatbed scanner in landscape and portrait orientation. The horizontally and vertically polarized light components of the scanner were determined, and another Epson Expression 10000XL flatbed scanner was disassembled to examine its optical components. The optical properties of exposed and unexposed EBT3 films were studied with incident polarized and unpolarized white light, and the transmitted red light was investigated for its polarization and scattering properties including the distribution of the scattering angles. Neutral density filters were studied for comparison. Guidance was sought from the theory of light scattering from rod-like macromolecular structures. The drastic dose-dependent variation of the transmitted total light current as function of the orientation of front and rear polarizers, interpreted by light scattering theory, shows that the radiation-induced polymerization of the monomers of EBT3 films produces light scattering oscillators preferably polarized at right angles with the coating direction of the film. The directional distribution of the scattered light is partly anisotropic, with a preferred scattering plane at right angles with the coating direction, indicating light scattering from stacks of coherently vibrating oscillators piled up along the monomer crystals. The polyester carrier film also participates in these effects. The orientation and parabola artefacts due to flatbed scanning of radiochromic films can be explained by the interaction of the polarization-dependent and anisotropic light scattering from exposed and unexposed EBT3 films with the quantitative difference between the scanners horizontally and vertically polarized light supply and with the limited directional acceptance of the scanners light recording system.

A new correction method serving to eliminate the parabola effect of flatbed scanners used in radiochromic film dosimetry

PURPOSEnThe purpose of this study is the correction of the lateral scanner artifact, i.e., the effect that, on a large homogeneously exposed EBT3 film, a flatbed scanner measures different optical densities at different positions along the x axis, the axis parallel to the elongated light source. At constant dose, the measured optical density profiles along this axis have a parabolic shape with significant dose dependent curvature. Therefore, the effect is shortly called the parabola effect. The objective of the algorithm developed in this study is to correct for the parabola effect. Any optical density measured at given position x is transformed into the equivalent optical density c at the apex of the parabola and then converted into the corresponding dose via the calibration of c versus dose.nnnMETHODSnFor the present study EBT3 films and an Epson 10000XL scanner including transparency unit were used for the analysis of the parabola effect. The films were irradiated with 6 MV photons from an Elekta Synergy accelerator in a RW3 slab phantom. In order to quantify the effect, ten film pieces with doses graded from 0 to 20.9 Gy were sequentially scanned at eight positions along the x axis and at six positions along the z axis (the movement direction of the light source) both for the portrait and landscape film orientations. In order to test the effectiveness of the new correction algorithm, the dose profiles of an open square field and an IMRT plan were measured by EBT3 films and compared with ionization chamber and ionization chamber array measurement.nnnRESULTSnThe parabola effect has been numerically studied over the whole measuring field of the Epson 10000XL scanner for doses up to 20.9 Gy and for both film orientations. The presented algorithm transforms any optical density at position x into the equivalent optical density that would be measured at the same dose at the apex of the parabola. This correction method has been validated up to doses of 5.2 Gy all over the scanner bed with 2D dose distributions of an open square photon field and an IMRT distribution.nnnCONCLUSIONSnThe algorithm presented in this study quantifies and corrects the parabola effect of EBT3 films scanned in commonly used commercial flatbed scanners at doses up to 5.2 Gy. It is easy to implement, and no additional work steps are necessary in daily routine film dosimetry.

Experimental determination of the lateral dose response functions of detectors to be applied in the measurement of narrow photon-beam dose profiles.

This study aims at the experimental determination of the detector-specific 1D lateral dose response function K(x) and of its associated rotational symmetric counterpart K(r) for a set of high-resolution detectors presently used in narrow-beam photon dosimetry. A combination of slit-beam, radiochromic film, and deconvolution techniques served to accomplish this task for four detectors with diameters of their sensitive volumes ranging from 1 to 2.2 mm. The particular aim of the experiment was to examine the existence of significant negative portions of some of these response functions predicted by a recent Monte-Carlo-simulation (Looe et al 2015 Phys. Med. Biol. 60 6585-607). In a 6 MV photon slit beam formed by the Siemens Artiste collimation system and a 0.5 mm wide slit between 10 cm thick lead blocks serving as the tertiary collimator, the true cross-beam dose profile D(x) at 3 cm depth in a large water phantom was measured with radiochromic film EBT3, and the detector-affected cross-beam signal profiles M(x) were recorded with a silicon diode, a synthetic diamond detector, a miniaturized scintillation detector, and a small ionization chamber. For each detector, the deconvolution of the convolution integral M(x)u2009u2009=u2009u2009K(x)u2009u2009∗u2009u2009D(x) served to obtain its specific 1D lateral dose response function K(x), and K(r) was calculated from it. Fourier transformations and back transformations were performed using function approximations by weighted sums of Gaussian functions and their analytical transformation. The 1D lateral dose response functions K(x) of the four types of detectors and their associated rotational symmetric counterparts K(r) were obtained. Significant negative curve portions of K(x) and K(r) were observed in the case of the silicon diode and the diamond detector, confirming the Monte-Carlo-based prediction (Looe et al 2015 Phys. Med. Biol. 60 6585-607). They are typical for the perturbation of the secondary electron field by a detector with enhanced electron density compared with the surrounding water. In the cases of the scintillation detector and the small ionization chamber, the negative curve portions of K(x) practically vanish. It is planned to use the measured functions K(x) and K(r) to deconvolve clinical narrow-beam signal profiles and to correct the output factor values obtained with various high-resolution detectors.

Dosimetric characteristics of the novel 2D ionization chamber array OCTAVIUS Detector 1500.

PURPOSEnThe dosimetric properties of the OCTAVIUS Detector 1500 (OD1500) ionization chamber array (PTW-Freiburg, Freiburg, Germany) have been investigated. A comparative study was carried out with the OCTAVIUS Detector 729 and OCTAVIUS Detector 1000 SRS arrays.nnnMETHODSnThe OD1500 array is an air vented ionization chamber array with 1405 detectors in a 27 × 27 cm(2) measurement area arranged in a checkerboard pattern with a chamber-to-chamber distance of 10 mm in each row. A sampling step width of 5 mm can be achieved by merging two measurements shifted by 5 mm, thus fulfilling the Nyquist theorem for intensity modulated dose distributions. The stability, linearity, and dose per pulse dependence were investigated using a Semiflex 31013 chamber (PTW-Freiburg, Freiburg, Germany) as a reference detector. The effective depth of measurement was determined by measuring TPR curves with the array and a Roos chamber type 31004 (PTW-Freiburg, Freiburg, Germany). Comparative output factor measurements were performed with the array, the Semiflex 31010 ionization chamber and the Diode 60012 (both PTW-Freiburg, Freiburg, Germany). The energy dependence of the OD1500 was measured by comparing the arrays readings to those of a Semiflex 31010 ionization chamber for varying mean photon energies at the depth of measurement, applying to the Semiflex chamber readings the correction factor kNR for nonreference conditions. The Gaussian lateral dose response function of a single array detector was determined by searching the convolution kernel suitable to convert the slit beam profiles measured with a Diode 60012 into those measured with the arrays central chamber. An intensity modulated dose distribution measured with the array was verified by comparing a OD1500 measurement to TPS calculations and film measurements.nnnRESULTSnThe stability and interchamber sensitivity variation of the OD1500 array were within ±0.2% and ±0.58%, respectively. Dose linearity was within 1% over the range from 5 to 1000 MU. The effective point of measurement of the OD1500 for dose measurements in RW3 phantoms was determined to be (8.7 ± 0.2) mm below its front surface. Output factors showed deviations below 1% for field sizes exceeding 4 × 4 cm(2). The dose per pulse dependence was smaller than 0.4% for doses per pulse from 0.2 to 1 mGy. The energy dependence of the array did not exceed ±0.9%. The parameter σ of the Gaussian lateral dose response function was determined as σ6MV = (2.07 ± 0.02) mm for 6 MV and σ15MV = (2.09 ± 0.02) mm for 15 MV. An IMRT verification showed passing rates well above 90% for a local 3 mm/3% criterion.nnnCONCLUSIONSnThe OD1500 arrays dosimetric properties showed the applicability of the array for clinical dosimetry with the possibility to increase the spatial sampling frequency and the coverage of a dose distribution with the sensitive areas of ionization chambers by merging two measurements.

Magnetic fields are causing small, but significant changes of the radiochromic EBT3 film response to 6 MV photons

The optical density (OD) of EBT3 radiochromic films (Ashland Specialty Ingredients, Bridgewater, NJ, USA) exposed to absorbed doses to water up to Du2009u2009=u2009u200920 Gy in magnetic fields of Bu2009u2009=u2009u20090.35 and 1.42 T was measured in the three colour channels of an Epson Expression 10000XL flatbed scanner. A 7u2009cm wide water phantom with fixed film holder was placed between the pole shoes of a constant-current electromagnet with variable field strength and was irradiated by a 6 MV photon beam whose axis was directed at right angles with the field lines. The doses at the film position at water depth 5u2009cm were measured with a calibrated ionization chamber when the magnet was switched off and were converted to the doses in presence of the magnetic field via the monitor units and by a Monte Carlo-calculated correction accounting for the slight change of the depth dose curves in magnetic fields. In the presence of the 0.35 and 1.42 T fields small negative changes of the OD values at given absorbed doses to water occurred and just significantly exceeded the uncertainty margin given by the stochastic and the uncorrected systematic deviations. This change can be described by au2009u2009+2.1% change of the dose values needed to produce a given optical density in the presence of a 1.42 T field. The thereby modified OD versus D function remained unchanged irrespective of whether the original short film side-the preference direction of the monomer crystals of the film-was directed parallel or orthogonal to the magnetic field. The orientation effect, the difference between the optical densities measured in the portrait or landscape film positions on the scanner bed caused by the reflection of polarised light in the scanners mirror system, remained unaltered after EBT3 film exposure in magnetic fields. An independent optical bench investigation of EBT3 films exposed to doses of 10 and 20 Gy at 0.35 and 1.42 T showed that the direction of the electric vector of polarised light experiencing the largest transmission through EBT3 films remained unaltered after film exposure in the magnetic fields. The observed small modification of the OD versus D curve of the radiochromic film EBT3 in the range up to 20 Gy and 1.42 T, hardly exceeding the experimental uncertainty margin, numerically confirms other recent studies on EBT3 film. A stronger magnetic field effect had been observed with the previous product EBT2 exposed to 60Co gamma radiation at 0.35 T.

Changes of the optical characteristics of radiochromic films in the transition from EBT3 to EBT-XD films

A new type of radiochromic film, the EBT-XD film, has been introduced with the aim to reduce the orientation effect and the lateral response artifact occurring in the use of radiochromic films together with flatbed scanners. The task of the present study is to quantify the changes of optical characteristics involved with the transition from the well-known EBT3 films to the new EBT-XD films, using the optical bench arrangement already applied by Schoenfeld et al (2014 Phys. Med. Biol. 59 3575-97). Largely reduced polarization effects and the almost complete loss of the anisotropy of the scattered light produced in a radiation-exposed film have been observed. The Rayleigh-Debye-Gans theory is used to understand these optical changes as arising from the reduced length-to-width ratio of the LiPCDA polymer crystals in the active layer of the EBT-XD film. The effect of these changes on the flatbed scanning artifacts will be shortly addressed, but treated in more detail in a further paper.

WE-D-AUD B-02: On the Influences of the Detector Size and Sampling Frequency On IMRT Verifications with 2D Arrays

Purpose: The influences of single detector size and sampling frequency of detector arrays on IMRT verification are discussed. Method and Materials: The 2D‐ARRAY Type 10024 (PTW‐Freiburg, single chamber cross section 5×5 mm2, center‐to‐center distance between chambers 10 mm) is analyzed as an example. Due to scattering effects at the ridges between the chambers, the Full‐Width‐At‐Half‐Maximum of the lateral detector response function is approximately 7 mm. By shifting the array in 3 steps of 5 mm, re‐measuring the dose distribution and arranging the data in a matrix, the sampling rate of 0.1 mm−1 can be increased to 0.2 mm−1 (sampling distance of 5mm). Results: The measurement process with detector arrays can be described as a two step process: 1. Convolution of the dose distribution with the response function of a single chamber of the array. 2. Sampling of the convolved dose distribution with the chosen sampling rate. Step 1 results in a small deviation of measured and real doses in the region of steep dose gradients. A mathematical model is introduced to estimate the deviation by consideration of detector size and initial dose penumbra. For the chosen array, the deviation in the region of clinical relevant doses is shown to be approximately 1mm. In the Fourier space, step 2 leads to a periodic replication of the Fourier Transform of the convolved dose distribution in intervals of the sampling frequency 0.2 mm−1. For various IMRT distributions we can show that the maximum spatial frequency does not exceed 0.1mm−1. According means that the sampling frequency is sufficient. Conclusion: With a sampling distance of 5 mm, the measurement of typical IMRT dose distributions with the 2D‐ARRAY complies with the Nyquist criterion. The results can be generalized to other arrays to analyze the limits of applicability for IMRT verification measurements.

SU‐FF‐T‐21: A Method to Increase the Resolution of IMRT Plan Verification with a Two‐Dimensional Ionisation Chamber Array

Purpose: We have developed a method for high resolution dosimetric plan verification with two‐dimensional ionization chamber arrays. Methods and Materials: The used 2D‐Array (PTW Freiburg, Germany) contains a matrix of 27 × 27 ionization chambers, each with an entrance window of 5 mm × 5 mm and with 5 mm wide ridges between the chambers. For an IMRT plan verification, the calculated dose distribution of the patient is exported to a CT of a phantom containing the 2D‐Array. The finite size of the chambers is accounted for by convolving the calculated dose distribution with the lateral transfer function of a single chamber. Considering the Nyquist theorem the chamber‐to‐chamber distance of 10 mm permits to resolve spatial frequencies up to 0,5/cm. Results: The resolution of the system can be doubled to 1,0/cm by shifting the array 5 mm in both x and y‐direction and repeating the measurement. By this the chambers will be positioned where the ridges have been in the first measurement. In most IMRT sequences the sizes of field elements are usually not smaller than 1 cm × 1cm, therefore the verification of these techniques can be achieved with sufficient resolution. Since lateral dose gradients are much higher in segmental than in dynamic IMRT techniques, the sampling distance may be increased for the latter. Conclusion: As most planning systems offer the possibility of limiting the field size to areas larger then 1 cm × 1 cm, a single measurement will mostly be sufficient for the verification. Clinical examples show the wide and easy applicability of the described methods. Conflict of Interest: The method was developed in cooperation with PTW‐Freiburg, Germany.