Search PubMed⌕ Search

Biomedical subjects

T R Mackie

Publications and source records attributed to T R Mackie.

At least 19 recordsLinked to original sources
PubMed2000

Iterative approaches to dose optimization in tomotherapy.

This paper will present the results of an investigation into three iterative approaches to inverse treatment planning. These techniques have been examined in the hope of developing an optimization algorithm suitable for the large-scale problems that are encountered in tomotherapy. The three iterative techniques are referred to as the ratio method, iterative least-squares minimization and the maximum-likelihood estimator. Our results indicate that each of these techniques can serve as a useful tool in tomotherapy optimization. As compared with other mathematical programming techniques, the iterative approaches can reduce both memory demands and time requirements. In this paper, the results from small- and large-scale optimizations will be analysed. It will also be demonstrated that the flexibility of the iterative techniques can be greatly enhanced through the use of dose-volume histogram based penalty functions and/or through the use of weighting factors assigned to each region of the patient. Finally, results will be presented from an investigation into the stability of the iterative techniques.

Algorithms↗
PubMed2000

Megavoltage CT image reconstruction during tomotherapy treatments.

An integrated tomotherapy system allows for improved radiotherapy verification by enabling the collection of megavoltage computed tomography (MVCT) images before or after treatment delivery. In this investigation, the possibility of collecting MV tomographic data and reconstructing images during a tomotherapy treatment is examined. By overcoming difficulties with the normalization of modulated treatment data and with the incompleteness of treatment data, it is possible to use data collected during tomotherapeutic treatments for MVCT reconstruction. The benefits of these techniques include potential increases in patient throughput, reductions in imaging dose, visualization of the patient in the treatment position and improvements in image contrast.

Animals↗
PubMed2000

Calibration of a tomotherapeutic MVCT system.

Megavoltage CT provides the ability to image the patient before, during or after a radiotherapy treatment. This allows one to verify not only the placement of a patient's external boundary, but also the locations of internal anatomy. In addition, the reconstructed MVCT values are potentially useful for treatment planning inhomogeneity corrections and dose reconstruction. To this end, dosimetric calibration of the University of Wisconsin Tomotherapy Benchtop MVCT system was investigated. It was found that MVCT values correlate extremely well with electron density and that unlike kilovoltage CT, this correlation is well maintained for higher atomic number materials. Improvements of the order of 1% in the dosimetric calculations of high atomic number materials should be possible by deriving input images from MVCT as opposed to kVCT, and calibrating in terms of electron density, as opposed to physical density.

Calibration↗
PubMed2000

Megavoltage CT imaging as a by-product of multileaf collimator leakage.

In addition to their potential for the delivery of highly conformal radiation therapy treatments, tomotherapeutic treatments also feature increased potential for verification. For example, megavoltage CT allows one to use the megavoltage linac to generate tomographic images of the patient in the treatment position. This is typically done before or after radiation therapy treatments. However, it is also possible to collect MVCT images entirely during the treatment itself. This process utilizes the leakage radiation through the closed leaves of the Nomos MIMiC MLC, along with slight inefficiencies in treatment delivery, to generate MVCT images during treatment that require neither additional time nor dose. The image quality is limited, yet sufficient to see a patient's external boundary, density differences over 8% for 25.0 mm objects and resolutions of 3.0 mm for high-contrast objects. Such images can potentially be viewed during treatment, used to flag additional CT immediately after the treatment and provide a representation of the patient's exact position during treatment for use with dose reconstruction.

Humans↗
PubMed2000

Modeling photon output caused by backscattered radiation into the monitor chamber from collimator jaws using a Monte Carlo technique.

Dose per monitor unit in photon fields generated by clinical linear accelerators can be affected by the backscattered radiation into the monitor chamber from collimator jaws. Thus, it is necessary to account for the backscattered radiation in computing monitor unit setting for a treatment field. In this work, we investigated effects of the backscatter from collimator jaws based on Monte Carlo simulations of a clinical linear accelerator. The backscattered radiation scored within the monitor chamber was identified as originating either from the upper jaws (Y jaws), or from the lower jaws (X jaws). From the results of Monte Carlo simulations, ratios of the monitor-chamber-scored dose caused by the backscatter to the dose caused by the forward radiation, R(x,y), were modeled as functions of the individual X and Y jaw positions. The amount of the backscattered radiation for any field setting was then computed as a compound contribution from both the X and Y jaws. The dose ratios of R(x,y) were then used to calculate the change in photon output caused by the backscatter, Scb(x,y). Results of these calculations were compared with available measured data based on counting the electron pulses or charge from the electron target of an accelerator. Data from this study showed that the backscattered radiation contributes approximately 3% to the monitor-chamber-scored dose. A majority of the backscattered radiation comes from the upper jaws, which are located closer to the monitor chamber. The amount of the backscatter decreases approximately in a linear fashion with the jaw opening. This results in about a 2% increase of photon output from a 10 cm x 10 cm field to a 40 cm x 40 cm field. The off-axis location of the jaw opening does not have a significant effect on the magnitude of the backscatter. The backscatter effect is significant for monitor chambers using kapton windows, particularly for treatment fields using moving jaws. Applying the backscatter correction improves the accuracy of monitor-unit calculation using a model-based dose calculation algorithm such as the convolution method.

Algorithms↗
PubMed1999

Tomotherapy.

Tomotherapy is delivery of intensity-modulated, rotational radiation therapy using a fan-beam delivery. The NOMOS (Sewickley, PA) Peacock system is an example of sequential (or serial) tomotherapy that uses a fast-moving, actuator-driven multileaf collimator attached to a conventional C-arm gantry to modulate the beam intensity. In helical tomotherapy, the patient is continuously translated through a ring gantry as the fan beam rotates. The beam delivery geometry is similar to that of helical computed tomography (CT) and requires the use of slip rings to transmit power and data. A ring gantry provides a stable and accurate platform to perform tomographic verification using an unmodulated megavoltage beam. Moreover, megavoltage tomograms have adequate tissue contrast and resolution to provide setup verification. Assuming only translational and rotational offset errors, it is also possible to determine the offsets directly from tomographic projections, avoiding the time-consuming image reconstruction operation. The offsets can be used to modify the leaf delivery pattern to match the beam to the patient's anatomy on each day of a course of treatment. If tomographic representations of the patient are generated, this information can also be used to perform dose reconstruction. In this way, the actual dose distribution delivered can be superimposed onto the tomographic representation of the patient obtained at the time of treatment. The results can be compared with the planned isodose on the planning CT. This comparison may be used as an accurate basis for adaptive radiotherapy whereby the optimized delivery is modified before subsequent fractions. The verification afforded tomotherapy allows more precise conformal therapy. It also enables conformal avoidance radiotherapy, the complement to conformal therapy, for cases in which the tumor volume is ill-defined, but the locations of sensitive structures are adequately determined. A clinical tomotherapy unit is under construction at the University of Wisconsin.

Dose Fractionation, Radiation↗
PubMed1999

Megavoltage CT on a tomotherapy system.

A megavoltage computed tomography (MVCT) system was developed on the University of Wisconsin tomotherapy benchtop. This system can operate either axially or helically, and collect transmission data without any bounds on delivered dose. Scan times as low as 12 s per slice are possible, and scans were run with linac output rates of 100 MU min(-1), although the system can be tuned to deliver arbitrarily low dose rates. Images were reconstructed with clinically reasonable doses ranging from 8 to 12 cGy. These images delineate contrasts below 2% and resolutions of 3.0 mm. Thus, the MVCT image quality of this system should be sufficient for verifying the patient's position and anatomy prior to radiotherapy. Additionally, synthetic data were used to test the potential for improved MVCT contrast using maximum-likelihood (ML) reconstruction. Specifically, the maximum-likelihood expectation-maximization (ML-EM) algorithm and a transmission ML algorithm were compared with filtered backprojection (FBP). It was found that for expected clinical MVCT doses enough imaging photons are used such that little benefit is conferred by the improved noise model of ML algorithms. For significantly lower doses, some quantitative improvement is achieved through ML reconstruction. Nonetheless, the image quality at those lower doses is not satisfactory for radiotherapy verification.

Algorithms↗
PubMed1999

Registration using tomographic projection files.

An algorithm has been developed and experimentally verified for tomographic registration--a patient positioning method using internal anatomy and standard external fiducial marks. This algorithm improves patient set-up and verification to an accuracy sufficient for tomotherapy. By implementation of this technique, the time-consuming reconstruction process is avoided. Instead, offsets in the x, y and z directions are determined directly from sinogram data by an algorithm that utilizes cross-correlations and Fourier transforms. To verify the efficiency and stability of the algorithm, data were collected on the University of Wisconsin's dedicated tomotherapy research workbench. The experiment indicates offset statistical errors of less than +/-0.8 mm for offsets up to 30 mm. With standard clinical techniques, initial patient offsets are expected to be less than 5 mm, so the 30 mm limitation is of no consequence. The angular resolution for the direction of patient translation is within the +/-2 degrees needed for tomotherapy.

Algorithms↗
PubMed1999

Delivery verification in sequential and helical tomotherapy.

Conformal and conformal avoidance radiation therapy are new therapeutic techniques that are generally characterized by high dose gradients. The success of this kind of treatment relies on quality assurance procedures in order to verify the delivery of the treatment. A delivery verification technique should consider quality assurance procedures for patient positioning and radiation delivery verification. A methodology for radiation delivery verification was developed and tested with our tomotherapy workbench. The procedure was investigated for two cases. The first treatment using a torus-shaped target was optimized for 72 beam directions and sequentially delivered as a single slice to a 33 cm diameter cylinder of homogeneous solid water. For the second treatment, a random pattern of energy fluence was helically delivered for two slices to a 9.0 cm diameter phantom containing inhomogeneities. The presented process provides the energy fluence (or a related quantity) delivered through the multileaf collimator (MLC) using the signal measured at the exit detector during the treatment delivery. As this information is created for every pulse of the accelerator, the energy fluence and state for each MLC leaf were verified on a pulse-by-pulse basis. The pulse-by-pulse results were averaged to obtain projection-by-projection information to allow for a comparison with the planned delivery. The errors between the planned and delivered energy fluences were concentrated between +/-2.0%, with none beyond +/-3.5%. In addition to accurately achieving radiation delivery verification, the process is fast, which could translate to radiation delivery verification in real time. This technique can also be extended to reconstruct the dose actually deposited in the patient or phantom (dose reconstruction).

Acrylates↗
PubMed1999

Image/patient registration from (partial) projection data by the Fourier phase matching method.

A technique for 2D or 3D image/patient registration, PFPM (projection based Fourier phase matching method), is proposed. This technique provides image/patient registration directly from sequential tomographic projection data. The method can also deal with image files by generating 2D Radon transforms slice by slice. The registration in projection space is done by calculating a Fourier invariant (FI) descriptor for each one-dimensional projection datum, and then registering the FI descriptor by the Fourier phase matching (FPM) method. The algorithm has been tested on both synthetic and experimental data. When dealing with translated, rotated and uniformly scaled 2D image registration, the performance of the PFPM method is comparable to that of the IFPM (image based Fourier phase matching) method in robustness, efficiency, insensitivity to the offset between images, and registration time. The advantages of the former are that subpixel resolution is feasible, and it is more insensitive to image noise due to the averaging effect of the projection acquisition. Furthermore, the PFPM method offers the ability to generalize to 3D image/patient registration and to register partial projection data. By applying patient registration directly from tomographic projection data, image reconstruction is not needed in the therapy set-up verification, thus reducing computational time and artefacts. In addition, real time registration is feasible. Registration from partial projection data meets the geometry and dose requirements in many application cases and makes dynamic set-up verification possible in tomotherapy.

Algorithms↗
PubMed1999

Development of a unique phantom to assess the geometric accuracy of magnetic resonance imaging for stereotactic localization.

OBJECTIVE: To test the spatial accuracy of coordinates generated from magnetic resonance imaging (MRI) scans, using the Brown-Roberts-Wells head frame and localizer system (Radionics, Inc., Burlington, MA). METHODS: An anthropomorphic head phantom, consisting of a two-dimensional lattice of acrylic spheres (4-mm diameter) spaced 10 mm apart and embedded in a brain tissue-mimicking gelatin-agar gel, was constructed. The intersphere distances for the target lattice positions in MRI and computed tomographic scan sets were compared. The data sets were fused, and differences in fiducial marker and intraphantom target positions were measured. RESULTS: Intersphere distances were identical for the MRI and computed tomographic scan sets (10 +/- 0.1 mm). Differences in fiducial marker positions [maximal lateral difference, 0.97 mm; mean absolute lateral difference, 0.69 +/- 0.22 mm; maximal anteroposterior (AP) difference, 1.99 mm; mean absolute AP difference, 1.29 +/- 0.67 mm] were correlated with differences in intraphantom target positions (maximal lateral difference, 0.83 mm; mean absolute lateral difference, 0.28 +/- 0.24 mm; maximal AP difference, -1.97 mm; mean absolute AP difference, 1.63 +/- 25 mm; maximal vertical difference, -0.73 mm; mean absolute vertical difference, 0.34 +/- 0.21 mm). This suggested that improper fiducial rod identification and the subsequent transformation to stereotactic coordinate space were the greatest sources of spatial uncertainty. CONCLUSION: With computed tomographic data as the standard, these differences resulted in maximal and minimal composite uncertainties of 2.06 and 1.17 mm, respectively. The measured uncertainties exceed recommended standards for radiosurgery but allow the possible use of MRI-based stereotactic treatment planning for certain intracranial lesions, if the errors are corrected using appropriate software. Clinicians must recognize that error magnitudes vary for different systems, and they should perform systematic, scheduled, institutional error analyses as part of their ongoing quality assurance processes. This phantom provides one tool for measuring such variances.

Brain Mapping↗
PubMed1999

Characterization of the output for helical delivery of intensity modulated slit beams.

The UW tomotherapy workbench utilizes a convolution/superposition based dose calculation and optimization program. It specifies the energy fluence that must be delivered from each leaf for each phantom projection angle. This requires that the spectrum of the radiation emitted from the one-dimensional MLC (multileaf collimator) attached to the linear accelerator be determined. The steps involved in that process are described. The spectrum along the central axis of the slit beam was determined, as well as the softening with off-axis position. Moreover, the magnitude of the energy-fluence output had to be quantified on a per MU (monitor unit) basis. This was done for a single leaf along the central axis of the beam. Factors, which modify that energy-fluence output, were investigated. The output increases with off-axis position due to the horns of the beam. The output for a leaf of interest will also increase if additional leaves are open due to the absence of the tongue-and-groove effect and penumbra blurring. The energy-fluence increase per leaf increase by 4.9% if an adjacent leaf is open. No other factors related to the state of additional leaves were found to significantly increase the energy-fluence output for an individual leaf.

Algorithms↗
PubMed1999

Multileaf collimator interleaf transmission.

Multileaf collimators (MLCs) have advanced past their original design purpose as a replacement for field shaping cerrobend blocks. Typically, MLCs incorporate an interlocking tongue-and-groove design between adjacent leaves to minimize leakage between leaves. They are beginning to be used to provide intensity modulation for conformal three-dimensional radiation therapy. It is possible that a critical target volume may receive an underdose due to the region of overlap if adjacent leaves are allowed to alternate between the open and closed positions, as they might if intensity modulation is employed. This work demonstrates the magnitude of that effect for a commercially available one-dimensional temporally modulated MLC. The magnitude of the transmission between leaves as a function of leaf separation was also studied, as well as the transmission as a function of leaf rotation away from the source. The results of this work were used for the design of a tomotherapy MLC. The radiation leakage considerations for a tomotherapy MLC are discussed.

Computer Simulation↗
PubMed1999

Monte Carlo investigation of electron beam output factors versus size of square cutout.

A major task in commissioning an electron accelerator is to measure relative output factors versus cutout size (i.e., cutout factors) for various electron beam energies and applicator sizes. We use the BEAM Monte Carlo code [Med Phys. 22, 503-524 (1995)] to stimulate clinical electron beams and to calculate the relative output factors for square cutouts. Calculations are performed for a Siemens MD2 linear accelerator with beam energies, 6, 9, 11, and 13 MeV. The calculated cutout factors for square cutouts in 10 X 10 cm2, 15 X 15 cm2, and 20 X 20 cm2 applicators at SSDs of 100 and 115 cm agree with the measurements made using a silicon diode within about 1% except for the smallest cutouts at SSD= 115 cm where they agree within 0.015. The details of each component of the dose, such as the dose from particles scattered off the jaws and the applicator, the dose from contaminant photons, the dose from direct electrons, etc., are also analyzed. The calculations show that inphantom side-scatter equilibrium is a major factor for the contribution from the direct component which usually dominates the output of a beam. It takes about 6 h of CPU time on a Pentium Pro 200 MHz computer to simulate an accelerator and additional 2 h to calculate the relative output factor for each cutout with a statistical uncertainty of 1%.

Monte Carlo Method↗
PubMed1999

A comparison of Monte Carlo and analytic first scatter dose spread arrays.

We compare first scattered point dose spread arrays generated by Monte Carlo and an analytic method. The analytic method models energy deposition using Klein-Nishina cross sections for Compton scatter and approximations for electron transport. Assumptions in the analytic method are shown to be valid within a region of the point dose spread kernel in which meaningful comparisons can be made. Differences between the models are less than 10% for the forward scatter directions for radii greater than the electron range associated with the first scattered Compton photon. Differences in the backscatter region are discussed and indicate that the analytic model is useful for identifying large errors that might be present in numerically generated first scatter point dose spread arrays. The analytic method is simple and useful for validating first scatter kernels.

Monte Carlo Method↗
PubMed1999

A simple model for examining issues in radiotherapy optimization.

Convolution/superposition software has been used to produce a library of photon pencil beam dose matrices. This library of pencil beams is designed to serve as a tool for both education and investigation in the field of radiotherapy optimization. The elegance of this pencil beam model stems from its cylindrical symmetry. Because of the symmetry, the dose distribution for a pencil beam from any arbitrary angle can be determined through a simple rotation of a pre-computed dose matrix. Rapid dose calculations can thus be performed while maintaining the accuracy of a convolution/superposition based dose computation. The pencil beam data sets have been made publicly available. It is hoped that the data sets will facilitate a comparison of a variety of optimization and delivery approaches. This paper will present a number of studies designed to demonstrate the usefulness of the pencil beam data sets. These studies include an examination of the extent to which a treatment plan can be improved through either an increase in the number of beam angles and/or a decrease in the collimator size. A few insights into the significance of heterogeneity corrections for treatment planning for intensity modulated radiotherapy will also be presented.

Computer Simulation↗