Modelling of output factors for conformal megavoltage X-ray beams.

A method is described for calculating the output from conformally shaped megavoltage X-ray beams. The model has been developed for Varian accelerators but is shown to work for accelerators from another manufacturer. The use of dynamic wedging and both static and dynamic multileaf collimated beams are included in the model. For any linear accelerator, the data required are a set of measured output factors for square beams, an in-air profile and a limited number of readily available parameters defining the geometry of the head of the accelerator. The three components of the output, namely primary, head scatter and phantom scatter are modelled and calculated individually for any point in a beam. An optimization procedure is developed that automatically determines the eight parameters required to model an accelerator in order for these calculations to be performed. The performance of the method is demonstrated for shaped beams using asymmetric and multileaf collimation, both with and without wedging, and for a range of beam energies. The model has been incorporated into a computer program that is used clinically.

Planning of beam intensity modulation using an advanced 3D dose calculation algorithm and a simulated annealing method.

PURPOSE: The aim of this work was to develop a fast inverse planning algorithm that will calculate optimum beam intensity distributions and beam shapes, and to incorporate the algorithm into a three-dimensional CT planning system. METHOD: The algorithm is based on the technique of simulated annealing and produces beam intensity distributions that could in principle be implemented clinically, either by the use of compensators or dynamic multileaf collimation. Dose distributions are calculated using a voxel beam model based on a spherical co-ordinate system, and transformations are given allowing the dose to be determined at any point within the patient. The dose calculation algorithm calculates primary and scattered dose separately from a knowledge of tissue/air ratios and differential scatter/air ratios, and both are corrected for the presence of heterogeneities in three dimensions. Specific attention is given to the execution time of the algorithm, and the methods developed allow satisfactory results to be achieved in calculation times which are sufficiently fast to be used interactively in the planning system. Several objective functions have been developed and can be selected in a simple manner by the user. In general, these attempt to achieve a uniform dose within the target while limiting the dose to organs at risk, either by upper dose limits or by specifying constraints on their dose volume histograms. RESULTS: The beam intensity distributions produced from the optimization have been used automatically by the forward planning system to produce three-dimensional dose distributions, and the results obtained in a number of clinical situations are presented. CONCLUSIONS: The inverse planning algorithm developed has been successfully incorporated into a three-dimensional planning system and is capable of producing beam intensity modulated distributions for clinical implementation. The execution time of the algorithm is sufficiently fast to be used as an optimization tool in an interactive forward planning system.

The use of compensators to optimise the three dimensional dose distribution in radiotherapy of the intact breast.

BACKGROUND AND PURPOSE: Dose heterogeneity in tangential breast irradiation has been shown to be as high as 20% and may lead to problems in local control and cosmesis. In this study, dose heterogeneity in three dimensions (3D) in the breast irradiated with wedged tangential beams is assessed and the improvement which can be made by the use of individualised two dimensional (2D) compensators is established. The compensation required is calculated in two ways: (I) by an iterative technique giving a uniform dose on a plane through the isocentre normal to the central axis of each beam, and (II) by inverse planning using an optimisation technique based on simulated annealing. MATERIALS AND METHODS: A total of 17 patients with histologically proven T0-3, N0, N1, M0 breast cancer undergoing breast irradiation following wide local excision, were CT scanned using contiguous 1 cm slices from approximately 2 cm superior to 2 cm inferior of the irradiated volume. The dose distributions are determined using a 3D algorithm that calculates primary and scatter dose separately using a differential scatter air ratio method and corrects both for the presence of heterogeneities. The iterative technique achieves a dose variation of better than 0.5% on the plane through the isocentre with compensation on both beams. Compensation for the lateral beam only is calculated using the optimisation technique in order to minimise the scatter dose to the contralateral breast. The optimisation algorithm minimises the dose variance over the target and sets upper dose limits for the lung and the remainder of the irradiated volume. RESULTS: For the group of patients the average dose heterogeneity in 3D using wedges is 12% (range 8-17%), which reduces to 8% (5-16%) using compensation on a plane and to 5% (4-7%) using the optimisation technique. CONCLUSIONS: Inverse planning is normally used for complex radiotherapy techniques but when applied to tangential breast irradiation, can reduce the dose heterogeneity through the breast as a whole to as little as 4%, with potential benefits in local control and cosmesis.

A multidisciplinary approach to improving the quality of tangential chest wall and breast irradiation for carcinoma of the breast.

Physicists, radiotherapists and radiographers have worked together to enhance the quality and accuracy of radiotherapy for tangential irradiation of the chest wall and breast. Each stage of the process has been reassessed and improved. A technique has been developed on the treatment simulator which determines the required beam directions and sizes in a straightforward manner. A computerised tomography facility has also been developed on the simulator and this provides one or more slices for planning through the treatment volume, thus allowing accurate determination of external contour and lung position with the patient in the treatment position. The beam edge entry points and the isocentre position can be seen from radiation opaque markers placed on the skin surface, allowing accurate reproduction of the treatment beam positions as set on the simulator. A photon beam algorithm that corrects for the changes in scatter dose in a 3-dimensional (3-D) inhomogeneous situation has been developed and applied to tangential chest wall irradiation. This has shown large differences (up to 10%) in dose compared to a conventional 2-dimensional algorithm. The changes in dose distribution due to the accurate determination of lung position have also been investigated. A method has been developed of measuring the volume of ipsilateral lung irradiated. Dose-area histograms are used to evaluate the fraction of the lung area irradiated in the central slice, and an estimate of the volume irradiated can be made using a beam's eye radiograph and the simulator CT image.

A technique for the fast calculation of three-dimensional photon dose distributions using the superposition model.

Techniques for reducing computation time in 3D photon dose calculations are addressed with specific emphasis given to the convolution/superposition approach. A single polyenergetic superposition model calculating absorbed dose per incident photon fluence (Gy cm2) was developed in terms of TERMA and a total energy deposition kernel (a total point spread function). A novel approach was devised for reducing calculation time. The method, named the CF method, was based on the use of a conventional, fast model (here a modified power-law method was used) for the generation of 3D dose distributions on a fine dose matrix. Superposition calculations were carried out on a coarse matrix and calculation speed was increased simply by reducing the number of calculations. A set of correction factors was derived on the coarse grid from the ratio of the dose values from superposition to those from the conventional algorithm. These were interpolated onto the fine matrix and used to modify the dose calculation from the conventional algorithm. The method was tested in a worst-case example where large dose gradients were present and in a clinically relevant irradiation geometry. It is shown that the time required for the generation of a 3D matrix with superposition can be reduced by at least a factor of 100 with no significant loss in accuracy.

A new technique for radiotherapy planning using quadratic programming.

Previous attempts at optimization of radiotherapy planning are described and criticized, and consideration of these attempts has resulted in the development of a new technique using quadratic programming. Uniformity of tumour dose is selected as the most important feature of any plan, and this is achieved by minimizing the variance of the dose to preselected points within the tumour. The dose to vulnerable regions can be constrained not to exceed a given percentage of the mean tumour dose. Optimization of field weight and field type is possible. The operation of the system is described and some typical results are given.

Rapid depth dose determination by a computer-controlled dosimetry system.

A comprehensive depth dose measuring system controlled by a small computer has been designed to enable full depth dose information for a linear accelerator to be obtained within one day. An ionization microchamber is moved rapidly under computer control in a horizontal plane to preselected points within a water phantom by two stepping motors to a spatial resolution of 0-025 mm. The ionization current pulses are converted to DC voltages by an interface providing complete charge integration and yet having a rapid response to drastic changes in dose level. These are fed to the computer for analysis and subsequent conversion to depth dose data. The chamber path follows a preselected fan-shaped grid after automatic determination of the X-ray beam central axis. The combination of rapid chamber motion with halts for charge integration and computer sampling has resulted in fast data acquisition independent of measurement response times. This information, which is stored on disc, can be directly read by a treatment planning program. The overall accuracy of measurement lies within 0-5% of the true depth dose at any point. Central axis and transverse profiles can be assessed at the time of measurement through visualization on the computer storage oscilloscope.

A beam model for three-dimensional radiotherapy planning.

A comprehensive beam model for three-dimensional radiotherapy planning is described. Input of patient data is either from an arbitrary spaced set of computed tomography (CT) scans throughout the treatment volume or from outlines traced on a digitizer. Any gantry, head or table rotation is allowed as well as any beam portal shape including internal blocking. Account is taken of beam modification due to the presence of the beam flattening filter, wedges and compensators. Both portal image and beam's eye view displays are readily available. A dose calculation algorithm is presented which uses a differential scatter-air ratio method to calculate the scatter dose in the presence of inhomogeneities. The beam model and algorithm form the basis of an interactive three-dimensional planning system running on a UNIX workstation, and performance times for dose calculation are given.

A generic computer program for checking photon beam dose calculations.

A computer program for checking photon external beam dose calculations is described. It provides a check of the treatment planning calculation by the use of machine data that are independent of that used in the initial calculation. Data required to specify any radiation beam are easily and rapidly set up and the program works for all the major manufacturers' machines. The user interface uses an interactive screen for machine data input and also for dose checking, where dose calculation is performed in real time as data are entered. The dose resulting from the planned field parameters is calculated and compared with the prescribed dose with a warning provided if the agreement is outside a set range (+/-5%). Its purpose is to act as a final quality assurance check to ensure that no significant errors occur in the monitor unit calculation.

The use of an image processing system in radiotherapy simulation.

Image quality produced in radiotherapy simulation is often so poor in certain anatomical sites (e.g. the pelvis) that it can adversely affect clinical judgement with regard to treatment set-up. In treatment planning both the field arrangement and the availability of accurate anatomical information are important, and should ideally be provided at the time of simulation by methods based on the principle of computerised tomography. In this work an image processing system has been used as a link between the video output from the simulator and the treatment planning computer in order to implement a CT facility on the simulator, as well as allowing real-time image enhancement to be investigated as a method of improving picture quality. This approach also enables a new method of simulation to be considered which reduces patient dose and prolongs the useful working life of the X-ray tube.

The reconstruction and dose-rate calculation of radium needle implants.

A comprehensive computer program has been developed for carrying out the various tasks associated with radium needle implants. Data input for the program are obtained by digitizing an orthogonal pair of radiographs with an analogue tracing system. Corresponding needles in the two radiographs are automatically located, and the implant may be rotated to any desired position in three dimensional space, at the same time presenting a view of the rotation on a storage oscilloscope. Automatic determination of the mean plane of a planar implant is possible, and the dose-rate distribution may be calculated in any specified plane. A Paterson--Parker exposure rate calculation may be performed with a biological dose and treatment time correction incorporated. The program is written in Fortran and runs on a small interactive computer system.

A comprehensive radiotherapy planning system implemented in Fortran on a small interactive computer.

A suite of Fortran programs for carrying out the various dose computational aspects of radiotherapy has been developed on an enhanced RAD8 computing system. The hardward configuration of the computer is described and the main features of the programs in the suite are discussed. The include: (a) beam data input systems for use with linear accelerators and cobalt units; (b) static and rotational teletherapy planning, with computer optimization in the static planning; (c) irregular field calculations with isodose visualization; (d) interstitial calculations including routines which will reconstruct a radium needle implant in three dimensions in addition to presenting the isodose distribution in any desired plane. The problems of implementing the programs on another computer system are discussed.

Stage IB cervical carcinoma: a clinical audit.

All patients with FIGO Stage IB cervical cancer registered with the Department of Clinical Oncology at the Western General Hospital, Edinburgh, during the 6 years from 1979 to 1984 have been reviewed, as part of a continuing programme of clinical audit. Of the 140 patients with Stage IB disease, 68 (49%) were treated by primary surgery of whom 44 (31%) also received adjuvant radiotherapy. Radical radiotherapy was the definitive treatment for 69 patients (49%). Three patients (2%) were not treated with curative intent. The crude 5-year survival rate for all cases was 72% and the cause-specific 5-year survival rate was 78%. Local tumour control at 5 years was 72%. There was no significant difference in outcome between the surgically treated and irradiated groups of patients. Age, histology and nodal status did not influence outcome. Irradiated patients with bulky tumours fared significantly worse than the other patients who received radical radiotherapy. Multivariate analysis of all patients revealed no significant independent prognostic variables. Primary surgery appears to confer no benefit over radical radiotherapy in terms of either survival or local control. Treatment-related late bladder and bowel morbidity was, however, significantly worse in irradiated patients.