Antony Lomax

The calibration of CT Hounsfield units for radiotherapy treatment planning

Computer tomographic (CT) scans are used to correct for tissue inhomogeneities in radiotherapy treatment planning. In order to guarantee a precise treatment, it is important to obtain the relationship between CT Hounsfield units and electron densities (or proton stopping powers for proton radiotherapy), which is the basic input for radiotherapy planning systems which consider tissue heterogeneities. A method is described to determine improved CT calibrations for biological tissue (a stoichiometric calibration) based on measurements using tissue equivalent materials. The precision of this stoichiometric calibration and the more usual tissue substitute calibration is determined by a comparison of calculated proton radiographic images based on these calibrations and measured radiographs of a biological sample. It has been found that the stoichiometric calibration is more precise than the tissue substitute calibration.

The 200-MeV proton therapy project at the Paul Scherrer Institute: conceptual design and practical realization.

The new proton therapy facility is being assembled at the Paul Scherrer Institute (PSI). The beam delivered by the PSI sector cyclotron can be split and brought into a new hall where it is degraded from 590 MeV down to an energy in the range of 85-270 MeV. A new beam line following the degrader is used to clean the low-energetic beam in phase space and momentum band. The analyzed beam is then injected into a compact isocentric gantry, where it is applied to the patient using a new dynamic treatment modality, the so-called spot-scanning technique. This technique will permit full three-dimensional conformation of the dose to the target volume to be realized in a routine way without the need for individualized patient hardware like collimators and compensators. By combining the scanning of the focused pencil beam within the beam optics of the gantry and by mounting the patient table eccentrically on the gantry, the diameter of the rotating structure has been reduced to only 4 m. In the article the degrees of freedom available on the gantry to apply the beam to the patient (with two rotations for head treatments) are also discussed. The devices for the positioning of the patient on the gantry (x rays and proton radiography) and outside the treatment room (the patient transporter system and the modified mechanics of the computer tomograph unit) are briefly presented. The status of the facility and first experimental results are introduced for later reference.

Intensity modulation methods for proton radiotherapy

The characteristic Bragg peak of protons or heavy ions provides a good localization of dose in three dimensions. Through their ability to deliver laterally and distally shaped homogenous fields, protons have been shown to be a precise and practical method for delivering highly conformal radiotherapy. However, in an analogous manner to intensity modulation for photons, protons can be used to construct dose distributions through the application of many individually inhomogeneous fields, but with the localization of dose in the Bragg peak providing the possibility of modulating intensity within each field in two or three dimensions. We describe four different methods of intensity modulation for protons and describe how these have been implemented in an existing proton planning system. As a preliminary evaluation of the efficacy of these methods, each has been applied to an example case using a variety of field combinations. Dose-volume histogram analysis of the resulting dose distributions shows that when large numbers of fields are used, all techniques exhibit both good target homogeneity and sparing of neighbouring critical structures, with little difference between the four techniques being discerned. As the number of fields is decreased, however, only a full 3D modulation of individual Bragg peaks can preserve both target coverage and sparing of normal tissues. We conclude that the 3D method provides the greatest flexibility for constructing conformal doses in challenging situations, but that when large numbers of beam ports are available, little advantage may be gained from the additional modulation of intensity in depth.

Potential reduction of the incidence of radiation-induced second cancers by using proton beams in the treatment of pediatric tumors

PURPOSE To assess the potential influence of improved dose distribution with proton beams compared to conventional or intensity-modulated (IM) X-ray beams on the incidence of treatment-induced secondary cancers in pediatric oncology. METHODS AND MATERIALS Two children, one with a parameningeal rhabdomyosarcoma (RMS) and a second with a medulloblastoma, were used as models for the purpose of this study. After defining the target and critical structures, treatment plans were calculated and optimized, four for the RMS case (conventional X-ray, IM X-rays, protons, and IM protons) and three for the irradiation of the spinal axis in medulloblastoma (conventional X-ray, IM X-rays, protons). Secondary cancer incidence was estimated using a model based on Publication No. 60 of the International Commission on Radiologic Protection. This model allowed estimation of absolute risks of secondary cancer for each treatment plan based on dose-volume distributions for the nontarget organs. RESULTS Proton beams reduced the expected incidence of radiation-induced secondary cancers for the RMS patient by a factor of >or=2 and for the medulloblastoma case by a factor of 8 to 15 when compared with either IM or conventional X-ray plans. CONCLUSIONS The potential for a significant reduction in secondary cancers with pediatric cancers after using proton beams (forward planned or IM) in the treatment of RMS and MBD in children and adolescents represents an additional argument supporting the development of proton therapy for most radiotherapy indications in pediatric oncology.

4D MR Imaging of Respiratory Organ Motion and its Variability

This paper describes a method for 4D imaging, which is used to study respiratory organ motion, a key problem in various treatments. Whilst the commonly used imaging methods rely on simplified breathing patterns to acquire one breathing cycle, the proposed method was developed to study irregularities in organ motion during free breathing over tens of minutes. The method does not assume a constant breathing depth or even strict periodicity and does not depend on an external respiratory signal. Time-resolved 3D image sequences were reconstructed by retrospective stacking of dynamic 2D images using internal image-based sorting. The generic method is demonstrated for the liver and for the lung. Quantitative evaluations of the volume consistency show the advantages over one-dimensional measurements for image sorting. Dense deformation fields describing the respiratory motion were estimated from the reconstructed volumes using non-rigid 3D registration. All obtained motion fields showed variations in the range of minutes such as drifts and deformations, which changed both the exhalation position of the liver and the breathing pattern. The obtained motion data are used in proton therapy planning to evaluate dose delivery methodologies with respect to their motion sensitivity. Besides this application, the new possibilities of studying respiratory motion are valuable for other applications such as the evaluation of gating techniques with respect to residual motion.

Experimental characterization and physical modelling of the dose distribution of scanned proton pencil beams

In this paper we present the pencil beam dose model used for treatment planning at the PSI proton gantry, the only system presently applying proton therapy with a beam scanning technique. The scope of the paper is to give a general overview on the various components of the dose model, on the related measurements and on the practical parametrization of the results. The physical model estimates from first physical principles absolute dose normalized to the number of incident protons. The proton beam flux is measured in practice by plane-parallel ionization chambers (ICs) normalized to protons via Faraday-cup measurements. It is therefore possible to predict and deliver absolute dose directly from this model without other means. The dose predicted in this way agrees very well with the results obtained with ICs calibrated in a cobalt beam. Emphasis is given in this paper to the characterization of nuclear interaction effects, which play a significant role in the model and are the major source of uncertainty in the direct estimation of the absolute dose. Nuclear interactions attenuate the primary proton flux, they modify the shape of the depth-dose curve and produce a faint beam halo of secondary dose around the primary proton pencil beam in water. A very simple beam halo model has been developed and used at PSI to eliminate the systematic dependences of the dose observed as a function of the size of the target volume. We show typical results for the relative (using a CCD system) and absolute (using calibrated ICs) dosimetry, routinely applied for the verification of patient plans. With the dose model including the nuclear beam halo we can predict quite precisely the dose directly from treatment planning without renormalization measurements, independently of the dose, shape and size of the dose fields. This applies also to the complex non-homogeneous dose distributions required for the delivery of range-intensity-modulated proton therapy, a novel therapy technique developed at PSI.

Intensity modulated proton therapy and its sensitivity to treatment uncertainties 2: the potential effects of inter-fraction and inter-field motions

Simple tools for studying the effects of inter-fraction and inter-field motions on intensity modulated proton therapy (IMPT) plans have been developed, and have been applied to both 3D and distal edge tracking (DET) IMPT plans. For the inter-fraction motion, we have investigated the effects of misaligned density heterogeneities, whereas for the inter-field motion analysis, the effects of field misalignment on the plans have been assessed. Inter-fraction motion problems have been analysed using density differentiated error (DDE) distributions, which specifically show the additional problems resulting from misaligned density heterogeneities for proton plans. Likewise, for inter-field motion, we present methods for calculating motion differentiated error (MDE) distributions. DDE and MDE analysis of all plans demonstrate that the 3D approach is generally more robust to both inter-fraction and inter-field motions than the DET approach, but that strong in-field dose gradients can also adversely affect a plans robustness. An important additional conclusion is that, for certain IMPT plans, even inter-fraction errors cannot necessarily be compensated for by the use of a simple PTV margins, implying that more sophisticated tools need to be developed for uncertainty management and assessment for IMPT treatments at the treatment planning level.

Treatment planning and verification of proton therapy using spot scanning: Initial experiences

Since the end of 1996, we have treated more than 160 patients at PSI using spot-scanned protons. The range of indications treated has been quite wide and includes, in the head region, base-of-skull sarcomas, low-grade gliomas, meningiomas, and para-nasal sinus tumors. In addition, we have treated bone sarcomas in the neck and trunk--mainly in the sacral area--as well as prostate cases and some soft tissue sarcomas. PTV volumes for our treated cases are in the range 20-4500 ml, indicating the flexibility of the spot scanning system for treating lesions of all types and sizes. The number of fields per applied plan ranges from between 1 and 4, with a mean of just under 3 beams per plan, and the number of fluence modulated Bragg peaks delivered per field has ranged from 200 to 45 000. With the current delivery rate of roughly 3000 Bragg peaks per minute, this translates into delivery times per field of between a few seconds to 20-25 min. Bragg peak weight analysis of these spots has shown that over all fields, only about 10% of delivered spots have a weight of more than 10% of the maximum in any given field, indicating that there is some scope for optimizing the number of spots delivered per field. Field specific dosimetry shows that these treatments can be delivered accurately and precisely to within +/-1 mm (1 SD) orthogonal to the field direction and to within 1.5 mm in range. With our current delivery system the mean widths of delivered pencil beams at the Bragg peak is about 8 mm (sigma) for all energies, indicating that this is an area where some improvements can be made. In addition, an analysis of the spot weights and energies of individual Bragg peaks shows a relatively broad spread of low and high weighted Bragg peaks over all energy steps, indicating that there is at best only a limited relationship between pencil beam weighting and depth of penetration. This latter observation may have some consequences when considering strategies for fast re-scanning on second generation scanning gantries.

Intensity modulated proton therapy: A clinical example

In this paper, we report on the clinical application of fully automated three-dimensional intensity modulated proton therapy, as applied to a 34-year-old patient presenting with a thoracic chordoma. Due to the anatomically challenging position of the lesion, a three-field technique was adopted in which fields incident through the lungs and heart, as well as beams directed directly at the spinal cord, could be avoided. A homogeneous target dose and sparing of the spinal cord was achieved through field patching and computer optimization of the 3D fluence of each field. Sensitivity of the resultant plan to delivery and calculational errors was determined through both the assessment of the potential effects of range and patient setup errors, and by the application of Monte Carlo dose calculation methods. Ionization chamber profile measurements and 2D dosimetry using a scintillator/CCD camera arrangement were performed to verify the calculated fields in water. Modeling of a 10% overshoot of proton range showed that the maximum dose to the spinal cord remained unchanged, but setup error analysis showed that dose homogeneity in the target volume could be sensitive to offsets in the AP direction. No significant difference between the MC and analytic dose calculations was found and the measured dosimetry for all fields was accurate to 3% for all measured points. Over the course of the treatment, a setup accuracy of +/-4 mm (2 s.d.) could be achieved, with a mean offset in the AP direction of 0.1 mm. Inhalation/exhalation CT scans indicated that organ motion in the region of the target volume was negligible. We conclude that 3D IMPT plans can be applied clinically and safely without modification to our existing delivery system. However, analysis of the calculated intensity matrices should be performed to assess the practicality, or otherwise, of the plan.

Intensity modulated proton therapy and its sensitivity to treatment uncertainties 1: the potential effects of calculational uncertainties

The effects of calculational uncertainties on 3D and distal edge tracking (DET) intensity modulated proton therapy (IMPT) treatment plans have been investigated. Dose calculation uncertainties have been assessed by comparing analytical and Monte Carlo dose calculations, and potential range uncertainties by recalculating plans with all CT values modified by +/-3%. Analysis of the volume of PTV agreeing to within +/-3% between the two calculations shows that the 3D approach provides significantly improved agreement (87.1 versus 80.3% of points for the 3D and DET approaches, respectively). For the DET approach, doses in the CTV have also been found to globally change by 5% as a result of 3% changes in CT value. When varying the intra-field gradients of the plans a similar trend is seen, but with the more complex plans also being found to be more sensitive to both uncertainties. In conclusion, the DET approach has been found to be relatively sensitive to the calculational errors investigated here. In contrast, the 3D approach appears to be quite robust, unless strong internal gradients are present. Nevertheless, the routine use of uncertainty analysis is advised when assessing all forms of IMPT plans.