Philip M. Evans

SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes

A software program, SpekCalc, is presented for the calculation of x-ray spectra from tungsten anode x-ray tubes. SpekCalc was designed primarily for use in a medical physics context, for both research and education purposes, but may also be of interest to those working with x-ray tubes in industry. Noteworthy is the particularly wide range of tube potentials (40-300 kVp) and anode angles (recommended: 6-30 degrees) that can be modelled: the program is therefore potentially of use to those working in superficial/orthovoltage radiotherapy, as well as diagnostic radiology. The utility is free to download and is based on a deterministic model of x-ray spectrum generation (Poludniowski 2007 Med. Phys. 34 2175). Filtration can be applied for seven materials (air, water, Be, Al, Cu, Sn and W). In this note SpekCalc is described and illustrative examples are shown. Predictions are compared to those of a state-of-the-art Monte Carlo code (BEAMnrc) and, where possible, to an alternative, widely-used, spectrum calculation program (IPEM78).

A cone-beam megavoltage CT scanner for treatment verification in conformal radiotherapy

PURPOSE A prototype scanner for large-volume megavoltage computed tomography (MVCT) in a clinical set-up is described. The ultimate aim is to improve treatment accuracy in conformal radiotherapy through patient set-up error reduction and transit dosimetry. MATERIALS AND METHODS The scanner consists of a custom-built 2D CsI(Tl) crystal array viewed by a lens and a CCD camera. Image acquisition is synchronized with radiation pulses. The 2D projections resulting from a single continuous 360 degrees gantry rotation are reconstructed using a cone-beam tomography algorithm. Prior to reconstruction, the raw projections are calibrated and corrected for centre of rotation movement and accelerator output fluctuation. The performance of the system has been evaluated by reconstructing projections of open fields, test objects and a humanoid phantom. RESULTS Hundreds of 2D projections can be acquired with a clinically-acceptable data collection time (about 2 min) and dose (approximately 40 cGy, with a possible four-fold reduction). A maximum density resolution of about 2% is achieved offering some soft tissue discrimination without using image enhancement tools. A spatial resolution of 2.5 mm is obtained. The reconstructed image intensity is linear with electron density over the range of interest. Coronal or sagittal slices through the 3D reconstruction of the humanoid phantom show a better delineation of structures than the corresponding portal images taken at the same orientation. CONCLUSIONS A similar image quality to our current single-slice MVCT scanner is achieved with the advantage of providing tens of tomographic slices for a single gantry rotation. This work demonstrates the feasibility of clinical cone-beam MVCT and indicates how this prototype can be improved.

The application of transit dosimetry to precision radiotherapy.

A method of using electronic portal imaging (EPI) for transit dosimetry is described. In this method, a portal image of the treatment field is first aligned with a digitally reconstructed radiograph (DRR) to geometrically relate the computed tomography (CT) scan, used to generate the DRR, with the EPI. Then the EPI is corrected for scatter within the patient to yield a map of primary fluence striking the detector. This is backprojected through the planning CT data set to yield a distribution of primary fluence within the patient. This distribution is then convolved with dose deposition kemels to yield a map of dose delivery within the patient. Such a distribution may be compared with the dose distribution resulting from the original treatment plan in order to evaluate the adequacy of the treatment. This method has been evaluated using a humanoid phantom. We find the transit dosimetry relative dose distribution when compared with film and thermoluminescent dosimeter (TLD) measurements and compared with our planning system to agree within 2% in the pelvic region of a humanoid phantom.

Calculation of x-ray spectra emerging from an x-ray tube. Part I. Electron penetration characteristics in x-ray targets

The penetration characteristics of electron beams into x-ray targets are investigated for incident electron kinetic energies in the range 50-150 keV. The frequency densities of electrons penetrating to a depth x in a target, with a fraction of initial kinetic energy, u, are calculated using Monte Carlo methods for beam energies of 50, 80, 100, 120 and 150 keV in a tungsten target. The frequency densities for 100 keV electrons in Al, Mo and Re targets are also calculated. A mixture of simple modeling with equations and interpolation from data is used to generalize the calculations in tungsten. Where possible, parameters derived from the Monte Carlo data are compared to experimental measurements. Previous electron transport approximations in the semiempirical models of other authors are discussed and related to this work. In particular, the crudity of the use of the Thomson-Whiddington law to describe electron penetration and energy loss is highlighted. The results presented here may be used towards calculating the target self-attenuation correction for bremsstrahlung photons emitted within a tungsten target.

The delivery of intensity modulated radiotherapy to the breast using multiple static fields

BACKGROUND AND PURPOSE To develop a method of using a multileaf collimator (MLC) to deliver intensity modulated radiotherapy (IMRT) for tangential breast fields, using an MLC to deliver a set of multiple static fields (MSFs). MATERIALS AND METHODS An electronic portal imaging device (EPID) is used to obtain thickness maps of medial and lateral tangential breast fields. From these IMRT deliveries are designed to minimize the volume of breast above 105% of prescribed dose. The deliveries are universally-wedged beams augmented with a set of low dose shaped irradiations. Dosimetric and planning QA of this method has been compared with the standard, wedged treatment and the corresponding treatment using physical compensators. Several options for delivering the MSF treatment are presented. RESULTS The MSF technique was found to be superior to the standard technique (P value=0.002) and comparable with the compensated technique. Both IMRT methods reduced the volume of breast above 105% dose from a mean value of 12.0% of the total breast volume to approximately 2.8% of the total breast volume. CONCLUSIONS This MSF method may be used to reduce the high dose volume in tangential breast irradiation significantly. This may have consequences for long-term side effects, particularly cosmesis.

Prone versus supine positioning for whole and partial-breast radiotherapy: a comparison of non-target tissue dosimetry.

PURPOSE To compare non-target tissue (including left-anterior-descending coronary-artery (LAD)) dosimetry of prone versus supine whole (WBI) and partial-breast irradiation (PBI). METHODS AND MATERIALS Sixty-five post-lumpectomy breast cancer patients underwent CT-imaging supine and prone. On each dataset, the whole-breast clinical-target-volume (WB-CTV), partial-breast CTV (tumour-bed + 15 mm), ipsilateral-lung and chest-wall were outlined. Heart and LAD were outlined in left-sided cases (n=30). Tangential-field WBI and PBI plans were generated for each position. Mean LAD, heart, and ipsilateral-lung doses (x(mean)), maximum LAD (LAD(max)) doses, and the volume of chest-wall receiving 50 Gy (V(50Gy)) were compared. RESULTS Two-hundred and sixty plans were generated. Prone positioning reduced heart and LAD doses in 19/30 WBI cases (median reduction in LAD(mean)=6.2 Gy) and 7/30 PBI cases (median reduction in LAD(max)=29.3 Gy) (no difference in 4/30 cases). However, prone positioning increased cardiac doses in 8/30 WBI (median increase in LAD(mean)=9.5 Gy) and 19/30 PBI cases (median increase in LAD(max)=22.9 Gy) (no difference in 3/30 cases). WB-CTV>1000cm(3) was associated with improved cardiac dosimetry in the prone position for WBI (p=0.04) and PBI (p=0.04). Prone positioning reduced ipsilateral-lung(mean) in 65/65 WBI and 61/65 PBI cases, and chest-wall V(50Gy) in all WBI cases. PBI reduced normal-tissue doses compared to WBI in all cases, regardless of the treatment position. CONCLUSIONS In the context of tangential-field WBI and PBI, prone positioning is likely to benefit left-breast-affected women of larger breast volume, but to be detrimental in left-breast-affected women of smaller breast volume. Right-breast-affected women are likely to benefit from prone positioning regardless of breast volume.

Independent verification using portal imaging of intensity‐modulated beam delivery by the dynamic MLC technique

The use of intensity-modulated radiation fields in radiotherapy treatment has been shown to have the potential to deliver highly conformal dose distributions. One technique for delivering these intensity-modulated beams is a computerized dynamic multileaf collimator (MLC). A major current impediment to the development of dynamic MLC therapy is verification of these highly complex treatments. Electronic portal imaging is shown here to be a solution to this verification problem. Experimental results are presented showing that leaf penumbra measured with a portal imager can be used to accurately define the positions of moving leaves. The random error in these leaf positions is compared with mean leaf positions along each leaf bank and specified leaf positions at prescription control points to check mechanical performance. Individual leaves are also checked for systematic motion errors. All leaf positions are found to be well within the manufacturers specifications at all times. Finally, integral intensity images are presented that can be related to the dose distribution delivered. Portal imaging is shown to have the potential to become a valuable tool for the verification of dynamic MLC irradiation.

A linear array, scintillation crystal-photodiode detector for megavoltage imaging

An imaging device has been developed to acquire images during external photon-beam radiotherapy treatments. It consists of a linear array of 128 zinc tungstate (ZnWO4) scintillation crystals each of which is individually optically coupled to a photodiode and associated electronics. The image is formed by scanning the linear array across the radiation field using a stepping motor under the control of a microcomputer. Image archive, display, and analysis are performed using a microVAX II computer. Results from a general theoretical analysis are presented before a detailed description of the particular detector construction. The mechanical design of the detector is such that the detector is automatically positioned to within a millimeter relative to the treatment source. This simplifies procedures for analyzing setup variations when comparing a treatment image to any other treatment, or planning, images. Image acquisition takes under 4 s with a contrast resolution of better than 1% at a spatial resolution of 2.5 mm in the object plane. The primary dose used to form these images is 0.55 cGy although the dose received by the patient will be closer to 25 cGy due to the linear scanning geometry and 3.8-s scan time that is used.

DOSE-VOLUME CONSTRAINTS TO REDUCE RECTAL SIDE EFFECTS FROM PROSTATE RADIOTHERAPY: EVIDENCE FROM MRC RT01 TRIAL ISRCTN 47772397

PURPOSE Radical radiotherapy for prostate cancer is effective but dose limited because of the proximity of normal tissues. Comprehensive dose-volume analysis of the incidence of clinically relevant late rectal toxicities could indicate how the dose to the rectum should be constrained. Previous emphasis has been on constraining the mid-to-high dose range (>/=50 Gy). Evidence is emerging that lower doses could also be important. METHODS AND MATERIALS Data from a large multicenter randomized trial were used to investigate the correlation between seven clinically relevant rectal toxicity endpoints (including patient- and clinician-reported outcomes) and an absolute 5% increase in the volume of rectum receiving the specified doses. The results were quantified using odds ratios. Rectal dose-volume constraints were applied retrospectively to investigate the association of constraints with the incidence of late rectal toxicity. RESULTS A statistically significant dose-volume response was observed for six of the seven endpoints for at least one of the dose levels tested in the range of 30-70 Gy. Statistically significant reductions in the incidence of these late rectal toxicities were observed for the group of patients whose treatment plans met specific proposed dose-volume constraints. The incidence of moderate/severe toxicity (any endpoint) decreased incrementally for patients whose treatment plans met increasing numbers of dose-volume constraints from the set of V30<or=80%, V40<or=65%, V50<or=55%, V60<or=40%, V65<or=30%, V70<or=15%, and V75<or=3%. CONCLUSION Considering the entire dose distribution to the rectum by applying dose-volume constraints such as those tested here in the present will reduce the incidence of late rectal toxicity.

Scattered radiation in portal images: A Monte Carlo simulation and a simple physical model

The scattered radiation in 6 MV radiotherapy portal images is analyzed. First, a quantity SPR* is studied, by means of Monte Carlo (MC) modeling. SPR* is defined as the ratio, on the central axis, of the signal due to scattered radiation to that due to the primary radiation. The detector model mimics a high-energy photon detector in the context of transit dosimetry. Second, a physical model of SPR* has been developed from first principles. For a cylindrical phantom, placed symmetrically about the isocenter, it predicts that SPR* depends on the area A at the isocenter of the circular field and the phantom thickness d as follows. SPR* = k0Ad(1 + k1d)(1 + k2A), where k0 = 0.0266(L1 + L2)2/(L1L2)2, k2 = - [L1(-2) + L2(-2) + (L1(-1) + L2(-1))2((2/3) + (3 kappa/2))]/2pi, L1 is the source-to-isocenter distance, L2 is the isocenter-to-detector distance, and kappa is the mean energy of the radiation beam (MeV/0.511). Constant k1, for which there is no simple expression, depends on L2. Comparison to the MC data shows that for 60 <or= L2 <or= 100 cm the dependence is weak and k1 approximately equal to 2 x 10(-3) cm-1. The root mean square (rms) agreement between the MC-derived values of SPR* and the physical model is better than 0.001 over a wide range of A and d values likely to be encountered in clinical practice for L2 >or= 50 cm. Third, experimental measurements of the scatter-to-primary ratio were obtained using our custom built imaging system mounted on a Philips SL25 linear accelerator. In the first experiment, A was varied from 40 to 400 cm2 with L1 = L2 = 100 cm with d = 20 cm. In the second experiment water depth d was varied from 0 to 28 cm with L1 = L2 = 100 cm and A = 200 cm2. The rms agreements between the MC data and the experiments were 0.0015 and 0.0045, respectively.