The perturbation effect in parallel-plate ionization chambers used for buildup measurements has been investigated. The fluence perturbation due to electrons emitted through the side walls are thoroughly investigated by measurements using film and extrapolation chambers and by calculations. The electron fluence varies both with side wall material and chamber geometry. In order to obtain a small perturbation effect, the chamber should have a large guard width compared to the electrode separation and the side walls should have as large an angle as possible with the central axis. The side wall should be of the same material as the rest of the chamber. The perturbation effect is also dependent on the electron contamination of the beam and angular distribution of the electron fluence. It is thus not possible to correct the perturbation effect in one parallel-plate chamber with fixed plate separation with correction factors obtained with extrapolation chambers of other dimensions. In order to make accurate surface dose measurements extrapolation chambers are therefore strongly recommended in favor of fixed parallel-plate chambers.

The output factor for a therapy photon beam is defined as the dose per monitor unit relative to the dose per monitor unit in a reference field. Convolution models for photon dose calculations yield the dose in units normalized to the incident energy fluence with phantom scatter intrinsically modeled. Output factors calculated with the convolution method as the dose per unit energy fluence relative to the calculated dose per unit energy fluence in a reference field could deviate as much as 5% if corrections are not made for perturbations due to treatment head scatter. Significant perturbations are particles backscattered from the collimators to the monitor and photons forward scattered from the filter and collimators in the treatment head. The forward scatter adds an "unmonitored" contribution to the total energy fluence of the beam. A model is developed that describes the field size dependence of these perturbations for conversion of output factors, calculated with the convolution method, to machine output factors as an integrated part in treatment planning. The necessary machine characteristics are derived from measurements of the output in air for a limited set of field sizes. The method has been tested using five different multileaf collimated irregular fields at 6 MV and for a large set of rectangular fields at 5, 6, and 18 MV and found to predict output factors with an accuracy better than 1%.

Three-dimensional dose planning systems employing accurate proton transport algorithms are essential for calculating absorbed dose distributions in proton therapy. In this paper, a pencil beam algorithm for the transport of protons in materials of interest for radiation therapy is developed. The Fermi-Eyges multiple-scattering theory is used to derive transport equations for calculating proton fluence and absorbed dose distributions. The multiple-scattering theory of Molière is used to predict mean square scattering angles and to develop an expression for calculating the root mean square (RMS) radial spread of a proton pencil beam, as a function of depth, in an arbitrary scattering material. A correction factor is suggested to account for the decrease in the radial spread at the end of the range due to range straggling. The effects of neglecting large-angle scattering events and the possibility of incorporating such events into the pencil beam algorithm are discussed. An energy scaling technique for determining the water-equivalent surface energy at a given depth in a heterogeneous scattering medium is developed. The water-equivalent energy, giving the same Molière scattering parameter B in water, is determined and the 1/e angle in water is scaled to the appropriate width in the scattering material. By using stored analytically or Monte Carlo calculated pencil beam distributions in water, the large-angle single-scattering events may be incorporated by approximating the scattering in an arbitrary material by the scattering in water for protons of the appropriate water-equivalent surface energy.

Determinations of the absorbed dose in a 170 MeV proton beam have been performed using seven ionization chambers of different types: five cylindrical (two FWT IC-18 and three NE-2571, of which one was modified to have the central electrode made of graphite) and two plane parallel (NACP-02 and Roos FK-6). The ionization was converted into absorbed dose in the proton beam according to the generalization of the formalism provided by the IAEA Code of Practice (TRS 277), which enables the use of the same equations for all kinds of beam used in radiotherapy. The absorbed dose obtained with the two IC-18 chambers, a chamber type commonly used as a reference in proton beams, was up to 1.5% lower than that obtained with the Farmer NE-2571 chamber, which was used as the reference in this work when calibration factors in terms of NK were used. To investigate this difference, experimental ND factors for six chambers (the two IC-18 chambers, the NACP-02, the FK-6 and two of the NE-2571 chambers) were determined in a high-energy electron beam. The procedure commonly recommended for plane parallel ion chambers was used for all the chambers, using the same reference chamber, a Farmer NE-2571. In the 170 MeV proton beam all the ND factors yielded consistent absorbed dose determinations within the estimated experimental uncertainties. This finding calls into question the value of the product kattkm for the IC-18 chamber given by the IAEA Code of Practice used in this comparison, and points at possible chamber to chamber variations that theoretical kattkm factors cannot predict. The investigations enabled the determination of the Pwall(60Co) factor of the Roos FK-6 plane parallel chamber, yielding 1.003 +/- 0.5%, and a correction for the effect of the aluminium central electrode of NE-2571 chambers in proton beams, equal to 1.003 +/- 0.4%. Two of the chambers (the plane parallel FK-6 and the modified cylindrical NE-2571) were provided with calibration factors in terms of absorbed dose to water, Nw, at the quality of 60Co by the Primary Standard Dosimetry Laboratory in Germany (PTB). Using the Nw formalism excellent agreement was found with the determination based on the experimental ND, giving support to the implementation of the NW procedure in therapeutic proton beams.

This paper describes a dual-scattering-foil technique for flattening of radiotherapeutic charged particle beams. A theory for optimization of shapes and thicknesses of the scattering foils is presented. The result is a universal optimal secondary-scatterer profile, which can be adapted to any charged particle beam by a simple scaling procedure. The calculation of the mean square scattering angle of the beam after passing through the scattering foils is done using the generalized Fermi-Eyges model for charged particle transport. It is shown that the fluence profile in the plane of interest can be made flat to better than 1% inside a predefined beam radius provided the shaped secondary scatterer has the universal radial thickness profile. The thicknesses of the two foils are optimized to minimize the total energy loss. The theory has been tested experimentally in an 180 MeV clinical proton beam. The measured distributions agree well with the calculations.

The implementation of two algorithms for calculating dose distributions for radiation therapy treatment planning of intermediate energy proton beams is described. A pencil kernel algorithm and a depth penetration algorithm have been incorporated into a commercial three dimensional treatment planning system (Helax-TMS, Helax AB, Sweden) to allow conformal planning techniques using irregularly shaped fields, proton range modulation, range modification and dose calculation for non-coplanar beams. The pencil kernel algorithm is developed from the Fermi Eyges formalism and Molière multiple-scattering theory with range straggling corrections applied. The depth penetration algorithm is based on the energy loss in the continuous slowing down approximation with simple correction factors applied to the beam penumbra region and has been implemented for fast, interactive treatment planning. Modelling of the effects of air gaps and range modifying device thickness and position are implicit to both algorithms. Measured and calculated dose values are compared for a therapeutic proton beam in both homogeneous and heterogeneous phantoms of varying complexity. Both algorithms model the beam penumbra as a function of depth in a homogeneous phantom with acceptable accuracy. Results show that the pencil kernel algorithm is required for modelling the dose perturbation effects from scattering in heterogeneous media.

BackgroundSemi-automated segmentation using deformable registration of selected atlas cases consisting of expert segmented patient images has been proposed to facilitate the delineation of lymph node regions for three-dimensional conformal and intensity-modulated radiotherapy planning of head and neck and prostate tumours. Our aim is to investigate if fusion of multiple atlases will lead to clinical workload reductions and more accurate segmentation proposals compared to the use of a single atlas segmentation, due to a more complete representation of the anatomical variations.MethodsAtlases for lymph node regions were constructed using 11 head and neck patients and 15 prostate patients based on published recommendations for segmentations. A commercial registration software (Velocity AI) was used to create individual segmentations through deformable registration. Ten head and neck patients, and ten prostate patients, all different from the atlas patients, were randomly chosen for the study from retrospective data. Each patient was first delineated three times, (a) manually by a radiation oncologist, (b) automatically using a single atlas segmentation proposal from a chosen atlas and (c) automatically by fusing the atlas proposals from all cases in the database using the probabilistic weighting fusion algorithm. In a subsequent step a radiation oncologist corrected the segmentation proposals achieved from step (b) and (c) without using the result from method (a) as reference. The time spent for editing the segmentations was recorded separately for each method and for each individual structure. Finally, the Dice Similarity Coefficient and the volume of the structures were used to evaluate the similarity between the structures delineated with the different methods.ResultsFor the single atlas method, the time reduction compared to manual segmentation was 29% and 23% for head and neck and pelvis lymph nodes, respectively, while editing the fused atlas proposal resulted in time reductions of 49% and 34%. The average volume of the fused atlas proposals was only 74% of the manual segmentation for the head and neck cases and 82% for the prostate cases due to a blurring effect from the fusion process. After editing of the proposals the resulting volume differences were no longer statistically significant, although a slight influence by the proposals could be noticed since the average edited volume was still slightly smaller than the manual segmentation, 9% and 5%, respectively.ConclusionsSegmentation based on fusion of multiple atlases reduces the time needed for delineation of lymph node regions compared to the use of a single atlas segmentation. Even though the time saving is large, the quality of the segmentation is maintained compared to manual segmentation.

Dose painting prescriptions can be derived from retrospectively observed recurrence volumes spatial relation to pre-treatment FDG-PET image data. The ideal dose redistributions could significantly increase the TCP for patients with large tumor volumes and large spread in SUV from FDG-PET. The results yield a basis for prospective studies to determine the clinical value for dose painting of head and neck squamous cell carcinomas.