Search PubMed⌕ Search

Biomedical subjects

T R Mackie

Publications and source records attributed to T R Mackie.

At least 55 records · Page 3Linked to original sources
PubMed1998

Calculating dose distributions and wedge factors for photon treatment fields with dynamic wedges based on a convolution/superposition method.

A convolution/superposition based method was developed to calculate dose distributions and wedge factors in photon treatment fields generated by dynamic wedges. This algorithm used a dual source photon beam model that accounted for both primary photons from the target and secondary photons scattered from the machine head. The segmented treatment tables (STT) were used to calculate realistic photon fluence distributions in the wedged fields. The inclusion of the extra-focal photons resulted in more accurate dose calculation in high dose gradient regions, particularly in the beam penumbra. The wedge factors calculated using the convolution method were also compared to the measured data and showed good agreement within 0.5%. The wedge factor varied significantly with the field width along the moving jaw direction, but not along the static jaw or the depth direction. This variation was found to be determined by the ending position of the moving jaw, or the STT of the dynamic wedge. In conclusion, the convolution method proposed in this work can be used to accurately compute dose for a dynamic or an intensity modulated treatment based on the fluence modulation in the treatment field.

Humans↗
PubMed1998

Effects of changes in stopping-power ratios with field size on electron beam relative output factors.

Stopping-power ratios are a function of field size and vary with accelerators. To investigate how these variations affect relative output factor measurements made using ion chambers for electron beams, especially for small fields, (L/rho)air(water) is calculated using the Monte Carlo technique for different field sizes, beam energies, and accelerators and is compared to the data in TG-21 or TG-25, which are for mono-energetic broad beams. For very small field sizes defined by cutouts, if the change in (L/rho)air(water) with dmax is ignored (i.e., TG-25 is not carefully followed), there is an overestimate of relative output factors by up to 3%. Ignoring the field-size effect on stopping-power ratio adds an additional overestimate of up to one-half percent, and using mono-energetic stopping-power ratio data instead of realistic beam data gives another error, but in the opposite direction, of up to 0.7%. Due to the cancellation of these latter two errors, following TG-25 with (L/rho)air(water) data for broad mono-energetic beams will give the correct answer for the ROF measurement within 0.4% compared to using (L/rho)air(water) data for which the field-size effect is considered for realistic electron beams.

Biophysical Phenomena↗
PubMed1998

Synchroton-produced ultrasoft X rays: equivalent cell survival at the isoattenuating energies 273 eV and 860 eV.

In this paper we report on survival of Chinese hamster V79 and mouse C3H 10T1/2 cells after irradiation with synchrotron-produced 273 eV and 860 eV ultrasoft X rays. These two energies, which are available by multilayer monochromatization of the synchrotron output spectrum, exhibit equal attenuation within living cells. Such an isoattenuating energy pair allows the direct examination of how biological effectiveness varies with the energy of the ultrasoft X rays. In comparing survival results, we find similar biological effectiveness of these two energies for both the C3H 10T1/2 and the V79 cells. These results are not consistent with previous findings of increasing RBE with decreasing ultrasoft X-ray energies. In addition, after correcting for mean nuclear dose based on measurements of cell thickness obtained with confocal microscopy, we find no significant differences in survival between the two ultrasoft X-ray energies and 250 kVp X rays. These results suggest that RBE does not increase with decreasing energy of ultrasoft X rays between 860 eV and 273 eV. The possible impact of our results on past results for ultrasoft X rays is discussed.

Animals↗
PubMed1998

A comparison of cytotoxicity after whole- or partial-cell irradiation with synchrotron-produced ultrasoft X rays.

We introduce a method which allows partial-volume irradiation of live cells using synchrotron-produced ultrasoft X rays and micro-fabricated irradiation masks. The masks were made by X-ray lithography at the University of Wisconsin Synchrotron Radiation Center, and they consist of 1.85-microm-wide stripes of gold 1.35 microm apart plated onto thin silicon nitride membranes. When placed adjacent to Mylar on which live cells are plated, these masks allow cells to be irradiated in a striped pattern with dimensions much smaller than the cell nuclei. Using 1340 eV synchrotron-produced X rays, we compare the survival of cells subjected to uniform irradiation and cells subjected to partial-volume irradiation. Our results show that, at equal mean dose to the nucleus (i.e. equal total energies deposited), survival is not statistically different for the two treatments over a wide range of doses. Thus imparting equal energies to smaller intranuclear volumes does not appear to enhance cell killing.

Animals↗
PubMed1997

Accurate characterization of Monte Carlo calculated electron beams for radiotherapy.

Monte Carlo studies of dose distributions in patients treated with radiotherapy electron beams would benefit from generalized models of clinical beams if such models introduce little error into the dose calculations. Methodology is presented for the design of beam models, including their evaluation in terms of how well they preserve the character of the clinical beam, and the effect of the beam models on the accuracy of dose distributions calculated with Monte Carlo. This methodology has been used to design beam models for electron beams from two linear accelerators, with either a scanned beam or a scattered beam. Monte Carlo simulations of the accelerator heads are done in which a record is kept of the particle phase-space, including the charge, energy, direction, and position of every particle that emerges from the treatment head, along with a tag regarding the details of the particle history. The character of the simulated beams are studied in detail and used to design various beam models from a simple point source to a sophisticated multiple-source model which treats particles from different parts of a linear accelerator as from different sub-sources. Dose distributions calculated using both the phase-space data and the multiple-source model agree within 2%, demonstrating that the model is adequate for the purpose of Monte Carlo treatment planning for the beams studied. Benefits of the beam models over phase-space data for dose calculation are shown to include shorter computation time in the treatment head simulation and a smaller disk space requirement, both of which impact on the clinical utility of Monte Carlo treatment planning.

Biophysical Phenomena↗
PubMed1997

An investigation of tomotherapy beam delivery.

Experimental simulations for tomotherapy beam delivery were performed using a computer-controlled phantom positioner, a cylindrical phantom, and a 6 MV x-ray slit beam. Both continuous helical beam and sequential segmented tomotherapy (SST) beam deliveries were evaluated. Beam junctioning problem due to couch indexing error or field width errors presented severe dose uniformity perturbations for SST, while the problem was minimized for helical beam delivery. Longitudinal breathing motions were experimentally simulated for helical and SST beam delivery. While motions reduced the dose uniformity perturbations for SST, small artifacts in dose uniformity can be introduced for helical beam delivery. With typical breath frequency and magnitude, for a slit beam of 2.0 cm width at 4 rpm, the dose uniformity perturbation was not significant. A running start/stop technique was implemented with helical beam delivery to sharpen the 20%-80% longitudinal dose fall-off from 1.5 to 0.5 cm. The latter was comparable to the corresponding dose penumbra of a conventional 6 MV 10 x 10 cm2 field. All together, helical beam delivery showed advantages over SST for tomotherapy beam delivery under similar delivery conditions.

Biophysical Phenomena↗
PubMed1997

Electron fluence correction factors for conversion of dose in plastic to dose in water.

In radiation dosimetry protocols, plastic is allowed as a phantom material for the determination of absorbed dose to water in electron beams. The electron fluence correction factor is needed in conversion of dose measured in plastic to dose in water. There are large discrepancies among recommended values as well as measured values of electron fluence correction factors when polystyrene is used as a phantom material. Using the Monte Carlo technique, we have calculated electron fluence correction factors for incident clinical beam energies between 5 and 50 MeV as a function of depth for clear polystyrene, white polystyrene and PMMA phantom materials and compared the results with those recommended in protocols as well as experimental values from published data. In the Monte Carlo calculations, clinical beams are simulated using the EGS4 user-code BEAM for a variety of medical accelerators. The study shows that our calculated fluence correction factor, phi pw, is a function of depth and incident beam energy Eo with little dependence on other aspects of beam quality. However the phi pw values at dmax are indirectly influenced by the beam quality since they vary with depth and dmax also varies with the beam quality. Calculated phi pw values at dmax are in a range of 1.005-1.045 for a clear polystyrene phantom, 1.005-1.038 for a white polystyrene phantom and 0.996-1.016 for a PMMA phantom. Our values of phi pw are about 1-2% higher than those determined according to the AAPM TG-25 protocol at dmax for clear or white polystyrene. Our calculated values of phi pw also explain some of the variations of measured data because of its depth dependence. A simple formula is derived which gives the electron fluence correction factor phi pw as a function of R50 at dmax or at the depth of 0.6R50-0.1 for any clinical electron beam with energy between 5 and 25 MeV for three plastics: clear polystyrene, white polystyrene and PMMA. The study also makes a careful distinction between phi pw and the corresponding IAEA Code of Practice quantity, hm.

Biophysical Phenomena↗
PubMed1997

Calculating dose and output factors for wedged photon radiotherapy fields using a convolution/superposition method.

We have developed a convolution/superposition method to calculate dose distributions in photon treatment fields with beam modifiers such as physical wedges. The dose component due to wedge generated radiation was accounted for by using an extended phantom model, which integrated a wedge, an air gap, and a patient phantom as the calculation phantom. The inhomogeneities in the extended phantom and the effect of beam hardening by the wedge were both corrected for in the convolution dose calculation. The calculated dose was verified by Monte Carlo simulation of the same extended phantom. A new dual photon source model was also used in the convolution method to account for both primary photons from the target and extra-focal photons from the primary collimator and flattening filter. Thus, realistic photon energy fluence distributions in the extended phantom were used for the dose calculation. The calculated dose distributions and the wedge factors agreed with the measured data within 2% for a variety of treatment fields including asymmetric fields. Our results showed that the wedge-generated radiation could contribute a significant fraction of the total dose in patients. This dose component depends on a specific field configuration, thus wedge factor changes with photon energy, wedge angle, field size, depth, and patient phantom SSD. The variation of the wedge factor can be predicted accurately by our convolution approach with the extended phantom model, which allows for more accurate dose or monitor unit computation for photon fields with beam modifiers.

Algorithms↗
PubMed1997

Correcting kernel tilting and hardening in convolution/superposition dose calculations for clinical divergent and polychromatic photon beams.

To account for clinical divergent and polychromatic photon beams, we have developed kernel tilting and kernel hardening correction methods for convolution dose calculation algorithms. The new correction methods were validated by Monte Carlo simulation. The accuracy and computation time of the our kernel tilting and kernel hardening correction methods were also compared to the existing approaches including terma divergence correction, dose divergence correction methods, and the effective mean kernel method with no kernel hardening correction. Treatment fields of 10 x 10-40 x 40 cm2 (field size at source to axis distance (SAD)) with source to source distances (SSDs) of 60, 80, and 100 cm, and photon energies of 6, 10, and 18 MV have been studied. Our results showed that based on the relative dose errors at a depth of 15 cm along the central axis, the terma divergence correction may be used for fields smaller than 10 x 10 cm2 with a SSD larger than 80 cm; the dose divergence correction with an additional kernel hardening correction can reduce dose error and may be more applicable than the terma divergence correction. For both these methods, the dose error increased linearly with the depth in the phantom; the 90% isodose lines at the depth of 15 cm were shifted by about 2%-5% of the field width due to significant underestimation of the penumbra dose. The kernel hardening effect was less prominent than the kernel tilting effect for clinical photon beams. The dose error by using nonhardening corrected kernel is less than 2.0% at a depth of 15 cm along the central axis, yet it increased with a smaller field size and lower photon energy. The kernel hardening correction could be more important to compute dose in the fields with beam modifiers such as wedges when beam hardening is more significant. The kernel tilting correction and kernel hardening correction increased computation time by about 3 times, and 0.5-1 times, respectively. This can be justified by more accurate dose calculations for the majority of clinical treatments.

Algorithms↗
PubMed1997

Analysis and convergence of the iterative convolution/superposition dose reconstruction technique for multiple treatment beams and tomotherapy.

An iterative convolution/superposition (C/S) algorithm has been created to reconstruct dose distributions in patients from exit dose measurements during a radiotherapy treatment. The method is based on an extended phantom which includes the patient CT representation and an electronic portal imaging device (EPID). The patient CT is assumed to be a true and rigid representation of the patient at the time of treatment. The C/S method computes the dose throughout the extended phantom which allows the exit dose to be predicted in the EPID. The process is then reversed to take the exit dose measurement and infer what the dose distribution must have been to produce the measured exit dose. The dose distribution is modeled without knowledge of the incident intensity distribution, and includes the effects of scatter in the computation. The iterative method begins by assuming that the exit primary energy fluence (PEF) is equal to the exit dose, the PEF is then backprojected through the extended phantom and superposed with the dose deposition kernel to determine a new prediction of the exit dose. The ratio of the computed PEF to exit dose is then multiplied by the measured exit dose image to produce a better representation of the exit PEF. Successive iterations then converge to the exit PEF image that would produce the measured exit dose image. Once convergence is established, the dose distribution is determined by backprojecting the exit PEF followed by superposition with the dose deposition kernel. The method is used to reconstruct the dose from a stimulated dynamic wedge and verified with film. Convergence and termination of the algorithm is then investigated with no noise and in the presence of noise. The method is then expanded to handle multiple treatment beams by separating the representation of the EPID from the patient or phantom representation in the computation process. Investigation of the effects of noise during the process of iterative dose reconstruction is necessary to understand the capabilities of the algorithm using exit dose images that may contain significant amounts of noise. The capability of the algorithm is evaluated for multiple field treatments to a cube phantom and a prostate patient CT representation in the presence of noise. The method is then used to simulate the dose reconstruction process for tomotherapy using 72 intensity-modulated fan beams. Dose reconstruction is shown to be capable of verifying the dose distributions in patients including multiple beams and dynamic collimation, provided the patient CT is known at the time of treatment.

Algorithms↗
PubMed1997

A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams.

A realistic model of photon beams generated by clinical linear accelerators has been incorporated in a convolution/superposition method to compute dose distributions in photon treatment fields. In this beam model, a primary photon source represents photons directly from the target, and an extra-focal photon source represents scattered photons from the primary collimator and the flattening filter. Monte Carlo simulation was used to study clinical linear accelerators producing photon beams. From the output of the Monte Carlo simulation, the fluence and spectral distributions of each photon component, as well as the geometrical characteristics of each photon source with respect to its distance to the isocenter and its source distribution, were analyzed. These quantities were used to reproduce realistic photon distributions in treatment fields, and thus to compute dose distributions using the convolution method. Our results showed that compared to the primary photon fluence, the extra-focal photon fluence from the primary collimator and the flattening filter was 11%-16% at the isocenter, among which 70% was contributed by the flattening filter. The variation of extra-focal photons in different treatment fields was predicted accurately by accounting for the finite size of the extra-focal source. Compared to measurements, dose distributions in photon treatment fields, including those of asymmetric jaw settings and at different SSDs were calculated accurately, particularly in the penumbral region, by using the convolution method with the new dual source photon beam model.

Computer Simulation↗
PubMed1997

Calculating output factors for photon beam radiotherapy using a convolution/superposition method based on a dual source photon beam model.

A realistic photon beam model based on Monte Carlo simulation of clinical linear accelerators was implemented in a convolution/superposition dose calculation algorithm. A primary and an extra-focal sources were used in this beam model to represent the direct photons from the target and the scattered photons from other head structures, respectively. The effect of the finite size of the extra-focal source was modeled by a convolution of the source fluence distribution with the collimator aperture function. Relative photon output in air (Sc) and in phantom (Scp) were computed using the convolution method with this new photon beam model. Our results showed that in a 10 MV photon beam, the Sc, Sp (phantom scatter factor), and Scp factors increased by 11%, 10%, and 22%, respectively, as the field size changed from 3 x 3 cm2 to 40 x 40 cm2. The variation of the Sc factor was contributed mostly by an increase of the extra-focal radiation with field size. The radiation backscattered into the monitor chamber inside the accelerator head affected the Sc by about 2% in the same field range. The output factors in elongated fields, asymmetric fields, and blocked fields were also investigated in this study. Our results showed that if the effect of the backscattered radiation was taken into account, output factors in these treatment fields can be predicted accurately by our convolution algorithm using the dual source photon beam model.

Computer Simulation↗
PubMed1996

Characteristics of secondary electrons produced by 6, 10 and 24 MV x-ray beams.

Megavoltage x-ray beams generated by linear accelerators (linacs) deliver their maximum dose a few centimetres below the treatment or phantom surface. This skin-sparing effect is degraded by the generation of secondary electrons as the x-ray beam passes to the patient or phantom. This work measures the characteristics of these electrons. A light-weight electromagnet was constructed that could be mounted in the block-tray position, 58 cm from the x-ray source of a Varian Clinac 2100C or 2500 linac. A field strength as high as 0.1500 T was generated, which was strong enough to sweep secondary electrons out of a 10 cm x 10 cm field. For 6, 10 and 24 MV x-ray beams, secondary contamination electrons produced 18, 38 and 65% of the surface dose, corresponding to 3, 5 and 12% of the maximum dose, respectively. A parameterized depth-dose curve for the contamination electrons was produced and was valid for all the x-ray energies studied.

Biophysical Phenomena↗
PubMed1996

Calculation of portal dose using the convolution/superposition method.

The convolution/superposition method was used to predict the dose throughout an extended volume, which includes a phantom and a portal imaging device. From the calculated dose volume, the dose delivered in the portal image plane was extracted and compared to a portal dose image. This comparison aids in verifying the beam configuration or patient setup after delivery of the radiation. The phantoms used to test the accuracy of this method include a solid water cube, a Nuclear Associates CT phantom, and an Alderson Rando thorax phantom. The dose distribution in the image plane was measured with film and an electronic portal imaging device in each case. The calculated portal dose images were within 4% of the measured images for most voxels in the central portion of the field for all of the extended volumes. The convolution/superposition method also enables the determination of the scatter and primary dose contributions using the particular dose deposition kernels for each contribution. The ratio of primary dose to total dose was used to extract the primary dose from the detected portal image, which enhances the megavoltage portal images by removing scatter blurring. By also predicting the primary energy fluence, we can find the ratio of computed primary energy fluence to total dose. Multiplying this ratio by the measured dose image estimates the relative primary energy fluence at the portal imager. The image of primary energy fluence possesses higher contrast and may be used for further quantitative image processing and dose modeling.

Biophysical Phenomena↗
PubMed1996

Modeling dose distributions from portal dose images using the convolution/superposition method.

Post-treatment dose verification refers to the process of reconstructing delivered dose distributions internal to a patient from information obtained during the treatment. The exit dose is commonly used to describe the dose beyond the exit surface of the patient from a megavoltage photon beam. Portal imaging provides a method of determining the dose in a plane distal to a patient from a megavoltage therapeutic beam. This exit dose enables reconstruction of the dose distribution from external beam radiation throughout the patient utilizing the convolution/superposition method and an extended phantom. An iterative convolution/superposition algorithm has been created to reconstruct dose distributions in patients from exit dose measurements during a radiotherapy treatment. The method is based on an extended phantom that includes the patient CT representation and an electronic portal imaging device (EPID). The convolution/superposition method computes the dose throughout the extended phantom, which allows the portal dose image to be predicted in the EPID. The process is then reversed to take the portal dose measurement and infer what the dose distribution must have been to produce the measured portal dose. The dose distribution is modeled without knowledge of the incident intensity distribution, and includes the effects of scatter in the computation. The iterative method begins by assuming that the primary energy fluence (PEF) at the portal image plane is equal to the portal dose image, the PEF is then back-projected through the extended phantom and convolved with the dose deposition kernel to determine a new prediction of the portal dose image. The image of the ratio of the computed PEF to the computed portal dose is then multiplied by the measured portal dose image to produce a better representation of the PEF. Successive iterations of this process then converge to the exiting PEF image that would produce the measured portal dose image. Once convergence is established, the dose distribution is determined by back-projecting the PEF and convolving with the dose deposition kernel. The method is accurate, provided the patient representation during treatment is known. The method was used on three phantoms with a photon energy of 6 MV to verify convergence and accuracy of the algorithm. The reconstructed dose volumes agree to within 3% of the forward computation dose volumes. Furthermore, this technique assumes no prior knowledge of the incident fluence and therefore may better represent the dose actually delivered.

Algorithms↗
PubMed1995

An iterative filtered backprojection inverse treatment planning algorithm for tomotherapy.

PURPOSE: An inverse treatment planning algorithm for tomotherapy is described. METHODS AND MATERIALS: The algorithm iteratively computes a set of nonnegative beam intensity profiles that minimizes the least-square residual dose defined in the target and selected normal tissue regions of interest. At each iteration the residual dose distribution is transformed into a set of residual beam profiles using an inversion method derived from filtered backprojection image reconstruction theory. These "residual" profiles are used to correct the current beam profile estimates resulting in new profile estimates. Adaptive filtering is incorporated into the inversion model so that the gross structure of the dose distribution is optimized during initial iterations of the algorithm, and the fine structure corresponding to edges is obtained at later iterations. A three dimensional, kernel based, convolutions/superposition dose model is used to compute dose during each iteration. RESULTS: Two clinically relevant treatment planning examples are presented illustrating the use of the algorithm for planning conformal radiotherapy of the breast and the prostate. Solutions are generally achieved in 10-20 iterations requiring about 20 h of CPU time using a midrange workstation. The majority of the calculation time is spent on the three-dimensional dose calculation. CONCLUSIONS: The inverse treatment planning algorithm is a useful research tool for exploring the potential of tomotherapy for conformal radiotherapy. Further work is needed to (a) achieve clinically acceptable computation times; (b) verify the algorithm using multileaf collimator technology; and (c) extend the method to biological objectives.

Algorithms↗
PubMed1995

MMC--a high-performance Monte Carlo code for electron beam treatment planning.

The macro Monte Carlo (MMC) method has been developed to improve the speed of traditional Monte Carlo (MC) high-energy electron transport calculations without loss in accuracy. The MMC algorithm uses results derived from conventional MC simulations of electron transport through macroscopic spheres of various radii and consisting of a variety of media. Based on these results, electrons are transported in macroscopic steps through the absorber. The absorber geometry is represented by a three-dimensional (3D) density matrix, typically derived from computer tomographic (CT) data. Energy lost by the electrons along their paths through the absorber is scored in a 3D dose matrix. Transport of secondary electrons and bremsstrahlung photons is taken into account. Major modifications of the original implementation of the MMC algorithm have resulted in an improved version of the code, resolving earlier problems with electron transport across interfaces of different materials, and running at a substantially higher speed. Furthermore, the code has been integrated into a clinical 3D treatment planning system. MMC results are in good agreement with results from conventional MC codes and are obtained with a speed gain of about one order of magnitude for clinically relevant irradiation situations. Calculation times to obtain a relative statistical accuracy of 2% per dose grid voxel for small electron field sizes are short enough to be routinely useful in radiotherapy clinics on present day affordable workstation computers. Considering speed, accuracy and memory requirements, MMC is a promising alternative to currently available electron dose planning algorithms.

Algorithms↗