T Mackie

BEAM: a Monte Carlo code to simulate radiotherapy treatment units.

This paper describes BEAM, a general purpose Monte Carlo code to simulate the radiation beams from radiotherapy units including high-energy electron and photon beams, 60Co beams and orthovoltage units. The code handles a variety of elementary geometric entities which the user puts together as needed (jaws, applicators, stacked cones, mirrors, etc.), thus allowing simulation of a wide variety of accelerators. The code is not restricted to cylindrical symmetry. It incorporates a variety of powerful variance reduction techniques such as range rejection, bremsstrahlung splitting and forcing photon interactions. The code allows direct calculation of charge in the monitor ion chamber. It has the capability of keeping track of each particles history and using this information to score separate dose components (e.g., to determine the dose from electrons scattering off the applicator). The paper presents a variety of calculated results to demonstrate the codes capabilities. The calculated dose distributions in a water phantom irradiated by electron beams from the NRC 35 MeV research accelerator, a Varian Clinac 2100C, a Philips SL75-20, an AECL Therac 20 and a Scanditronix MM50 are all shown to be in good agreement with measurements at the 2 to 3% level. Eighteen electron spectra from four different commercial accelerators are presented and various aspects of the electron beams from a Clinac 2100C are discussed. Timing requirements and selection of parameters for the Monte Carlo calculations are discussed.

Image guidance for precise conformal radiotherapy

PURPOSE To review the state of the art in image-guided precision conformal radiotherapy and to describe how helical tomotherapy compares with the image-guided practices being developed for conventional radiotherapy. MATERIALS AND METHODS Image guidance is beginning to be the fundamental basis for radiotherapy planning, delivery, and verification. Radiotherapy planning requires more precision in the extension and localization of disease. When greater precision is not possible, conformal avoidance methodology may be indicated whereby the margin of disease extension is generous, except where sensitive normal tissues exist. Radiotherapy delivery requires better precision in the definition of treatment volume, on a daily basis if necessary. Helical tomotherapy has been designed to use CT imaging technology to plan, deliver, and verify that the delivery has been carried out as planned. The image-guided processes of helical tomotherapy that enable this goal are described. RESULTS Examples of the results of helical tomotherapy processes for image-guided intensity-modulated radiotherapy are presented. These processes include megavoltage CT acquisition, automated segmentation of CT images, dose reconstruction using the CT image set, deformable registration of CT images, and reoptimization. CONCLUSIONS Image-guided precision conformal radiotherapy can be used as a tool to treat the tumor yet spare critical structures. Helical tomotherapy has been designed from the ground up as an integrated image-guided intensity-modulated radiotherapy system and allows new verification processes based on megavoltage CT images to be implemented.

Fast free-form deformable registration via calculus of variations

In this paper, we present a fully automatic, fast and accurate deformable registration technique. This technique deals with free-form deformation. It minimizes an energy functional that combines both similarity and smoothness measures. By using calculus of variations, the minimization problem was represented as a set of nonlinear elliptic partial differential equations (PDEs). A Gauss-Seidel finite difference scheme is used to iteratively solve the PDE. The registration is refined by a multi-resolution approach. The whole process is fully automatic. It takes less than 3 min to register two three-dimensional (3D) image sets of size 256 x 256 x 61 using a single 933 MHz personal computer. Extensive experiments are presented. These experiments include simulations, phantom studies and clinical image studies. Experimental results show that our model and algorithm are suited for registration of temporal images of a deformable body. The registration of inspiration and expiration phases of the lung images shows that the method is able to deal with large deformations. When applied to the daily CT images of a prostate patient, the results show that registration based on iterative refinement of displacement field is appropriate to describe the local deformations in the prostate and the rectum. Similarity measures improved significantly after the registration. The target application of this paper is for radiotherapy treatment planning and evaluation that incorporates internal organ deformation throughout the course of radiation therapy. The registration method could also be equally applied in diagnostic radiology.

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 patients 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.

History of tomotherapy.

Tomotherapy is the delivery of intensity modulated radiation therapy using rotational delivery of a fan beam in the manner of a CT scanner. In helical tomotherapy the couch and gantry are in continuous motion akin to a helical CT scanner. Helical tomotherapy is inherently capable of acquiring CT images of the patient in treatment position and using this information for image guidance. This review documents technological advancements of the field concentrating on the conceptual beginnings through to its first clinical implementation. The history of helical tomotherapy is also a story of technology migration from academic research to a university-industrial partnership, and finally to commercialization and widespread clinical use.

Radiation characteristics of helical tomotherapy

Helical tomotherapy is a dedicated intensity modulated radiation therapy (IMRT) system with on-board imaging capability (MVCT) and therefore differs from conventional treatment units. Different design goals resulted in some distinctive radiation field characteristics. The most significant differences in the design are the lack of flattening filter, increased shielding of the collimators, treatment and imaging operation modes and narrow fan beam delivery. Radiation characteristics of the helical tomotherapy system, sensitivity studies of various incident electron beam parameters and radiation safety analyses are presented here. It was determined that the photon beam energy spectrum of helical tomotherapy is similar to that of more conventional radiation treatment units. The two operational modes of the system result in different nominal energies of the incident electron beam with approximately 6 MeV and 3.5 MeV in the treatment and imaging modes, respectively. The off-axis mean energy dependence is much lower than in conventional radiotherapy units with less than 5% variation across the field, which is the consequence of the absent flattening filter. For the same reason the transverse profile exhibits the characteristic conical shape resulting in a 2-fold increase of the beam intensity in the center. The radiation leakage outside the field was found to be negligible at less than 0.05% because of the increased shielding of the collimators. At this level the in-field scattering is a dominant source of the radiation outside the field and thus a narrow field treatment does not result in the increased leakage. The sensitivity studies showed increased sensitivity on the incident electron position because of the narrow fan beam delivery and high sensitivity on the incident electron energy, as common to other treatment systems. All in all, it was determined that helical tomotherapy is a system with some unique radiation characteristics, which have been to a large extent optimized for intensity modulated delivery.

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.

Stereotactic radiosurgery for glioblastoma multiforme: Report of a prospective study evaluating prognostic factors and analyzing long-term survival advantage

PURPOSE Prospective evaluation of the toxicity and efficacy of radiosurgery with external beam radiotherapy in the management of newly diagnosed glioblastoma. METHODS AND MATERIALS From 5/89 to 12/92, 31 out of 51 patients with glioblastoma multiforme underwent radiosurgery, in addition to 54 Gy in 1.8 Gy/fraction following biopsy (n = 12) or resection (n = 19). Eligibility required supratentorial glioblastoma, tumor not > 4 cm in > 1 axis, age > 18 years, and location > 1 cm from optic chiasm. Patient characteristics were: age 20-78 years (median = 57); 22 male, 9 female; Karnofsky score 20-90 (m = 70), and volume of 2.3-59.7 c.c. (m = 17.4). Eighteen patients were treated with 1 collimator, 5 with 2, 7 with 3, and 1 with 4; peripheral isodoses were 40-90% (m = 72.5) and minimum and maximum tumor dose ranges were 10-20 (m = 12) and 15-35 Gy (m = 18.75). Patients were followed clinically and radiographically every 8-12 weeks to analyze survival, quality of life, and toxicity. RESULTS With a follow-up of 12-171 weeks, 8 out of 31 (26%) patients are alive. Median survival is 42 weeks. Twelve and 24-month actuarial survival are 38 and 28%. Comparison of the 2-year survival with previous Radiation Therapy Oncology Group patients was carried out using a nonparametric recursive partitioning technique and the observed vs. expected values are 28 vs. 9.7% (p < 0.05). Extent of resection and performance status were associated with improved survival in a multivariate analysis. No significant acute toxicity was encountered. Four patients (13%) developed clinically significant necrosis verified by biopsy or positron emission tomography scan at 9-59 weeks after radiosurgery. CONCLUSION The improvement in median survival in broadly selected glioblastoma patients treated with radiosurgery is difficult to determine, but the 2-year survival may be superior. Future randomized trials of radiosurgery are recommended, and ad hoc use of this modality should be discouraged.

Accurate characterization of Monte Carlo calculated electron beams for radiotherapy

Monte Carlo studies of dose distributions in patients treated with radiotherapyelectron 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 scatteredbeam.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 Carlotreatment 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 Carlotreatment planning.

Application of the spirometer in respiratory gated radiotherapy

The signal from a spirometer is directly correlated with respiratory motion and is ideal for target respiratory motion tracking. However, its susceptibility to signal drift deters its application in radiotherapy. In this work, a few approaches are investigated to control spirometer signal drift for a Bernoulli-type spirometer. A method is presented for rapid daily calibration of the spirometer to obtain a flow sensitivity function. Daily calibration assures accurate airflow measurement and also reduces signal drift. Dynamic baseline adjustment further controls the signal drift. The accuracy of these techniques was studied and it was found that the spirometer is able to provide a long-term drift-free breathing signal. The tracking error is comprised of two components: calibration error and stochastic signal baseline variation error. The calibration error is very small (1% of 3 l) and therefore negligible. The stochastic baseline variation error can be as large as 20% of the normal breathing amplitude. In view of these uncertainties, the applications of spirometers in treatment techniques that rely on breathing monitoring are discussed. Spirometer-based monitoring is noted most suitable for deep inspiration breath-hold but less important for free breathing gating techniques.