Dose-volume histograms

Estimating the dose variation in a volume of interest with explicit consideration of patient geometric variation

A method to measure the effects of internal organ motion and deformation and patient setup error on cumulative dose variation in a volume of interest is proposed. The method uses multiple CT scans and electronic portal images of a single patient to numerically simulate dose-volume effects over the entire course of the patient's external beam treatment. The results are expressed in the form of a novel dose-volume histogram, called an expected dose-volume histogram (EDVH).

A new three-dimensional dose distribution reduction scheme for tubular organs

In tubular structures, spatial aspects of the dose distribution may be important in determining the normal tissue response. Conventional dose-volume-histograms (DVHs) and dose-surface-histograms (DSHs) lack spatial information and may not be adequate to represent the three-dimensional (3D) dose data. A new 3D dose distribution data reduction scheme which preserves its longitudinal and circumferential character is presented. Dose distributions were generated at each axial level for esophagus or rectum in 123 patients with lung cancer or prostate cancer. Dose distribution histograms at each axial level were independently analyzed along the esophageal or rectal circumference to generate dose-circumference-histogram (DCH) sheets. Two types of plots were then generated from the DCH sheet. The first considered the percentage of the circumference at each axial level receiving various doses. The second considered the minimum dose delivered to any percentage of the circumference at each axial level. The DCH as a treatment planning tool can be easily implemented in a 3D planing system and is potentially useful for the study of the relationship between the complication risk and the longitudinal and circumferential dose distributions.

Analysis of biopsy outcome after three-dimensional conformal radiation therapy of prostate cancer using dose-distribution variables and tumor control probability models

PURPOSE

To investigate tumor control following three-dimensional conformal radiation therapy (3D-CRT) of prostate cancer and to identify dose-distribution variables that correlate with local control assessed through posttreatment prostate biopsies.

METHODS AND MATERIAL

Data from 132 patients, treated at Memorial Sloan-Kettering Cancer Center (MSKCC), who had a prostate biopsy 2.5 years or more after 3D-CRT for T1c-T3 prostate cancer with prescription doses of 64.8-81 Gy were analyzed. Variables derived from the dose distribution in the PTV included: minimum dose (Dmin), maximum dose (Dmax), mean dose (Dmean), dose to n% of the PTV (Dn), where n = 1%,...,99%. The concept of the equivalent uniform dose (EUD) was evaluated for different values of the surviving fraction at 2 Gy (SF(2)). Four tumor control probability (TCP) models (one phenomenologic model using a logistic function and three Poisson cell kill models) were investigated using two sets of input parameters, one for low and one for high T-stage tumors. Application of both sets to all patients was also investigated. In addition, several tumor-related prognostic variables were examined (including T-stage, Gleason score). Univariate and multivariate logistic regression analyses were performed. The ability of the logistic regression models (univariate and multivariate) to predict the biopsy result correctly was tested by performing cross-validation analyses and evaluating the results in terms of receiver operating characteristic (ROC) curves.

RESULTS

In univariate analysis, prescription dose (Dprescr), Dmax, Dmean, dose to n% of the PTV with n of 70% or less correlate with outcome (p < 0.01). The area under the ROC curve for Dmean is 0.64. In contrast, Dmin (p = 0.6), D98 (p = 0.2) or D95 (p = 0.1) are not significantly correlated with outcome. The results for EUD depend on the input parameter SF(2): EUD correlates significantly with outcome for SF(2) of 0.4 or more, but not for lower SF(2) values. Using either of the two input parameters sets, all TCP models correlate with outcome (p < 0.05; ROC areas 0.60-0.62). Using T-stage dependent input parameters, the correlation is improved (logistic function: p < 0.01, ROC area 0.67, Poisson models: p < 0.01, ROC areas 0.64-0.66). In comparison, the ROC area is 0.68 for the combination of Dmean and T-stage. After multivariate analysis, a model based on TCP, D20 and Gleason score is the best overall model (ROC area 0.73). However, an alternative model based on Dmean, Gleason score, and T-stage is competitive (ROC area 0.70).

CONCLUSION

Biopsy outcome after 3D-CRT of prostate cancer at MSKCC is not correlated with Dmin in the PTV and appears to be insensitive to cold spots in the dose distribution. This observation likely reflects the fact that much of the PTV, especially at the periphery, may not contain viable tumor cells and that the treatment margins were sufficiently large. Therefore, the predictive power of all variables which are sensitive to cold spots, like TCPs with Poisson models and EUD for low SF(2), is limited because the low dose region may not coincide with the tumor location. Instead, for MSKCC prostate cancer patients with their standardized CTV definition, substantial target motion and small dose inhomogeneities, Dmean (or any variable that downplays the effect of cold spots) is a very good predictor of biopsy outcome. While our findings may indicate a general problem in the application of current TCP models to clinical data, these conclusions should not be extrapolated to other disease sites without careful analysis.

Dose-volume analysis of different stereotactic radiotherapy mono-isocentric techniques

Several stereotactic irradiation techniques, using Linacs with the patient in lying and sitting position and a Gamma Knife Unit, were compared with regard to mono-isocentric three-dimensional dose distributions. Three types of target volumes, a sphere and two ellipsoids, were used for the comparisons. All three targets were centered on a real head, reconstructed from transversal CT scans. The ARTEMIS 3D Treatment Planning System, developed by the Tenon Hospital, Paris, was used for the dosimetry and the dose-volume histogram (DVH) calculation. For the comparative study, several quantitative parameters were used, derived from the dose-volume histogram calculation. Differential DVHs were plotted for each target volume and beam arrangement. Irradiation techniques were compared by deriving quantitative parameters from the DVHs such as mean and integral dose delivered to the target and normal tissue irradiated, as well as by the relative volume of the examined areas. All techniques used in this study produced very similar dose distributions. The small differences confirm the capability of the studied techniques to produce the same irradiation effects. By changing from the spherical target shape to a more elliptical shape, more of the normal tissue was irradiated with higher doses. For elliptical cases we therefore identified a need for more conformal stereotactic planning.

The probability of correct target dosage: dose-population histograms for deriving treatment margins in radiotherapy

PURPOSE

To provide an analytical description of the effect of random and systematic geometrical deviations on the target dose in radiotherapy and to derive margin rules.

METHODS AND MATERIALS

The cumulative dose distribution delivered to the clinical target volume (CTV) is expressed analytically. Geometrical deviations are separated into treatment execution (random) and treatment preparation (systematic) variations. The analysis relates each possible preparation (systematic) error to the dose distribution over the CTV and allows computation of the probability distribution of, for instance, the minimum dose delivered to the CTV.

RESULTS

The probability distributions of the cumulative dose over a population of patients are called dose-population histograms in short. Large execution (random) variations lead to CTV underdosage for a large number of patients, while the same level of preparation (systematic) errors leads to a much larger underdosage for some of the patients. A single point on the histogram gives a simple "margin recipe." For example, to ensure a minimum dose to the CTV of 95% for 90% of the patients, a margin between CTV and planning target volume (PTV) is required of 2.5 times the total standard deviation (SD) of preparation (systematic) errors (Sigma) plus 1.64 times the total SD of execution (random) errors (sigma') combined with the penumbra width, minus 1.64 times the SD describing the penumbra width (sigma(p)). For a sigma(p) of 3.2 mm, this recipe can be simplified to 2.5 Sigma + 0.7 sigma'. Because this margin excludes rotational errors and shape deviations, it must be considered as a lower limit for safe radiotherapy.

CONCLUSION

Dose-population histograms provide insight into the effects of geometrical deviations on a population of patients. Using a dose-probability based approach, simple algorithms for choosing margins were derived.

Residual plaque burden, delivered dose, and tissue composition predict 6-month outcome after balloon angioplasty and beta-radiation therapy

BACKGROUND

Inhomogeneity of dose distribution and anatomic aspects of the atherosclerotic plaque may influence the outcome of irradiated lesions after balloon angioplasty (BA). We evaluated the influence of delivered dose and morphological characteristics of coronary stenoses treated with beta-radiation after BA.

METHODS AND RESULTS

Eighteen consecutive patients treated according to the Beta Energy Restenosis Trial 1.5 were included in the study. The site of angioplasty was irradiated with the use of a beta-emitting (90)Sr/(90)Y source. With the side branches used as anatomic landmarks, the irradiated area was identified and volumetric assessment was performed by 3D intracoronary ultrasound imaging after treatment and at 6 months. The type of tissue, the presence of dissection, and the vessel volumes were assessed every 2 mm within the irradiated area. The minimal dose absorbed by 90% of the adventitial volume (D(v90)Adv) was calculated in each 2-mm segment. Diffuse calcified subsegments and those containing side branches were excluded. Two hundred six coronary subsegments were studied. Of those, 55 were defined as soft, 129 as hard, and 22 as normal/intimal thickening. Plaque volume showed less increase in hard segments as compared with soft and normal/intimal thickening segments (P<0.0001). D(v90)Adv was associated with plaque volume at follow-up after a polynomial equation with linear and nonlinear components (r = 0.71; P = 0.0001). The multivariate regression analysis identified the independent predictors of the plaque volume at follow-up: plaque volume after treatment, D(v90)Adv, and type of plaque.

CONCLUSIONS

Residual plaque burden, delivered dose, and tiss composition play a fundamental role in the volumetric outcome at 6-month follow-up after beta-radiation therapy and BA.

A semi-empirical treatment planning model for optimization of multiprobe cryosurgery

A model is presented for treatment planning of multiprobe cryosurgery. In this model a thermal simulation algorithm is used to generate temperature distribution from cryoprobes, visualize isotherms in the anatomical region of interest (ROI) and provide tools to assist estimation of the amount of freezing damage to the target and surrounding normal structures. Calculations may be performed for any given freezing time for the selected set of operation parameters. The thermal simulation is based on solving the transient heat conduction equation using finite element methods for a multiprobe geometry. As an example, a semi-empirical optimization of 2D placement of six cryoprobes and their thermal protocol for the first freeze cycle is presented. The effectiveness of the optimized treatment protocol was estimated by generating temperature-volume histograms and calculating the objective function for the anatomy of interest. Two phantom experiments were performed to verify isotherm locations predicted by calculations. A comparison of the predicted 0 degrees C isotherm with the actual iceball boundary imaged by x-ray CT demonstrated a spatial agreement within +/-2 mm.

The role of intravascular ultrasound imaging in vascular brachytherapy

Intracoronary brachytherapy has recently emerged as a new therapy to prevent restenosis. Initial experimental work was achieved in animal models and the results were assessed by histomorphometry. Initial clinical trials used angiography to guide dosimetry and to assess efficacy. Intravascular ultrasound (IVUS) permits tomographic examination of the vessel wall, elucidating the true morphology of the lumen and transmural components, which cannot be investigated on the lumenogram obtained by angiography. This paper reviews the use of IVUS in the clinical studies of brachytherapy conducted to date. IVUS allows clinicians to make a thorough assessment of the remodeling of the vessel and appears to have a major role to play in facilitating understanding of the underlying mechanisms of action in this emerging field. The authors propose that state-of-the-art IVUS techniques should be employed to further knowledge of the mechanisms of action of brachytherapy in atherosclerotic human coronary arteries.

Optimum beam configurations in tomographic intensity modulated radiation therapy

We review and extend the theory of tomographic dose reconstruction for intensity modulated radiotherapy (IMRT). We derive the basis for a saturation with beam number of dose conformation, and provide an analysis which ranks particular beam orientations in terms of the contribution to the delivered dose. Preferred beam directions are found which effectively reduce the number of beams necessary to achieve a given level of dose conformation. The analysis is a new application of the tomographic projection-slice theorem to the problem of beam orientation determination. The effects of the beam front filter and the positivity constraint arising from the tomographic approach are analysed, and modifications of the beam front filter for small beam numbers are suggested. The theory is applied to simple geometric shapes in two dimensions. A Gaussian ellipse, where analytical results are obtained, and simple hard-edged convex prescribed dose shapes are examined to illustrate beam selection based on the beam overlap metric. More complex concave prescribed dose shapes which contain a sensitive organ are also analysed and for low beam numbers are found to have preferred beam directions.

Applying the equivalent uniform dose formulation based on the linear-quadratic model to inhomogeneous tumor dose distributions: Caution for analyzing and reporting

We apply the concept of equivalent uniform dose (EUD) to our data set of model distributions and intensity modulated radiotherapy (IMRT) treatment plans as a method for analyzing large dose inhomogeneities within the tumor volume. For large dose nonuniformities, we find that the linerar-quadratic based EUD model is sensitive to the linear-quadratic model parameters, alpha and beta, making it necessary to consider EUD as a function of these parameters. This complicates the analysis for inhomogeneous dose distributions. EUD provides a biological estimate that requires interpretation and cannot be used as a single parameter for judging an inhomogeneous plan. We present heuristic examples to demonstrate the dose volume effect associated with EUD and the correlation to statistical parameters used for describing dose distributions. From these examples and patient plans, we discuss the risk of incorrectly applying EUD to IMRT patient plans.

Approaches to local salvage of soft tissue sarcoma after primary site failure

Improvements in pretreatment evaluation and the wider application of specialized multidisciplinary care have substantially reduced the risk of local recurrence for patients with soft tissue sarcomas arising at any site, and the recurrences that are still seen are often those that are most difficult to manage effectively. The management strategy for an isolated local recurrence of soft tissue sarcoma is usually similar to that used for primary disease. With appropriate pretreatment evaluation and salvage therapy that includes a multidisciplinary approach, most patients with local recurrence can expect local disease control and a good functional outcome. The development of effective management of a local recurrence is often a complex problem. The necessary decisions are influenced by the tumor location, disease extent, and previous local therapy. The need for specialized care is stressed. Patient evaluation, management strategies, and expected outcome for various clinical scenarios are discussed.

Quantitative pulmonary single photon emission computed tomography for radiotherapy applications

Pulmonary imaging using single photon emission computed tomography (SPECT) is the focus of current radiotherapy research, including dose-response analysis and three-dimensional (3D) radiation treatment planning. Improvement in the quantitative capability of SPECT may help establish its potential role in this application as well as others requiring accurate knowledge of pulmonary blood flow. The purposes of this study were to quantitatively evaluate SPECT filtered backprojection (FBP) and ordered subset-expectation maximization (OS-EM) reconstruction implementations for measuring absolute activity concentration in lung phantom experiments, and to incorporate quantitative SPECT techniques in 3D-RTP for lung cancer. Quantitative FBP (nonuniform iterative Chang attenuation compensation, scatter correction, and 3D postreconstruction Metz filtering) and OS-EM implementations were compared with a "clinical" implementation of FBP (uniform multiplicative Chang attenuation compensation and post-reconstruction von Hann filtering), for their ability to improve quantification of inactive and active spherical defects in the lungs of an anthropomorphic torso phantom. Activity concentration estimates were found to depend on many factors, such as region of interest size, scatter subtraction constant (k), postreconstruction deconvolution filtering and, in the case of OS-EM, total number of iterations. In general, reconstruction implementations incorporating compensation for nonuniform attenuation and scatter provided reduced bias relative to the clinical implementation. Potential applications to lung radiotherapy, including dose-functional histograms and treatment planning are also discussed. SPECT has the potential to provide accurate estimates of lung activity distributions that, together with improved image quality, may be useful for the study and prediction of therapeutic response.

An improved technique for comparing Gamma Knife dose-volume distributions in stereotactic radiosurgery

A function derived from the geometry of brachytherapy dose distributions is applied to stereotactic radiosurgery and an algorithm for the production of a novel dose-volume histogram, the Anderson inverse-square shifted dose-volume histogram (DVH), is proposed. The expected form of the function to be plotted is checked by calculating its value for single focus exposures, and its application to clinical examples of Gamma Knife treatments described. The technique is shown to provide a valuable tool for assessing the adequacy of radiosurgical plans and comparing and reporting dose distributions.

A quality assurance phantom for three-dimensional radiation treatment planning

PURPOSE

Three-dimensional (3D) radiation treatment planning is facilitated through the use of computerized radiation treatment planning systems (RTPSs) and CT simulators (CT-sims). Quality assurance (QA) of these systems is necessary for ensuring that they fulfill their potential. However, comprehensive tools for the systematic QA of these systems have not been developed. We present a phantom that facilitates the evaluation of a large number of nondosimetric functions. These include CT image acquisition and transfer, graphical displays of 3D radiation beams, multiplanar CT image reconstructions, digitally reconstructed radiographs, the representation and manipulation of contoured patient anatomy, dose volume histograms, and the conversion of CT numbers to relative electron densities.

METHODS AND MATERIALS

A phantom was constructed which contains materials and geometries that are appropriate for the routine QA of the features described above. The anatomy of the phantom is used as a standard against which the performance of the 3D-RTPS or CT-sim is evaluated. The phantom was used to evaluate three different 3D-RTPSs and a CT-sim at four institutions.

RESULTS

Using this phantom, clinically significant errors and limitations in commercially available 3D treatment planning software were discovered. No errors were discovered in the beam display or image reconstructions in the systems examined. Problems were found in the anatomy display, automatic tools, and the CT number to relative electron density conversion data used in some of the systems.

CONCLUSION

This phantom is a unique tool designed explicitly for the QA of 3D treatment planning software. Errors and limitations discovered through its use indicate that the QA of commercial treatment planning software is necessary, and that this phantom is an effective device for this task.

The delta-TCP concept: a clinically useful measure of tumor control probability

PURPOSE

The aim of this article is to provide a quantitative tool to evaluate the influence of the different dose regions in a non-uniformly irradiated tumour upon the probability of controlling that tumor.

METHODS AND MATERIALS

First, a method to generate a distribution of the probability of controlling the cells in a voxel (VCP) is explored and found not to be useful. Second, we introduce the concept of delta-TCP, which represents the gain or loss in the overall TCP as a result of each particular bin in a DVH not receiving the prescribed dose (the same concept is applicable to dose cubes or to a fraction of the bin). The delta-TCP method presented here is based on the Poisson TCP model, but any other model could also be used. Third, using this tool, with parameters appropriate to Stage C prostate tumors, the consequences of "cold" and "hot" dose regions have been explored.

RESULTS

We show that TCP is affected by the minimum dose, even if it is delivered to a very small volume (20% dose deficit to 5% of the volume makes the TCP decrease by 18%), and that a hot region may be "wasted" unless the boost is to the bulk of the volume. An example of the application of the delta-TCP concept to a prostate radiotherapy plan is also given.

CONCLUSION

The delta-TCP distribution adds more objective information to the original DVH by enabling the clinician or planner to directly evaluate the effects of a non-uniform dose distribution on local control.

Comparison of brachytherapy strategies based on dose-volume histograms derived from quantitative intravascular ultrasound

PURPOSE

We present in this paper the comparison, by simulation, of different treatment strategies based either on beta- or gamma-sources, both with and without a centering device. Ionizing radiation to prevent restenosis is an emerging modality in interventional cardiology. Numerous clinical studies are presently being performed or planned, but there is variability in dose prescription, and both gamma- and beta-emitters are used, leading to a wide range of possible dose distributions over the arterial vessel wall. This paper discusses the potential merits of dose-volume histograms (DVH) based on three-dimensional (3-D) reconstruction of electrocardiogram (ECG)-gated intravascular ultrasound (IVUS) to compare brachytherapy treatment strategies.

MATERIALS AND METHODS

DVH describe the cumulative distribution of dose over three specific volumes: (1) at the level of the luminal surface, a volume was defined with a thickness of 0.1 mm from the automatically detected contour of the highly echogenic blood-vessel interface; (2) at the level of the IVUS echogenic media-adventitia interface (external elastic lamina [EEL]), an adventitial volume was computed considering a 0.5-mm thickness from EEL; and (3) the volume encompassed between the luminal surface and the EEL (plaque + media). The IVUS data used were recorded in 23 of 31 patients during the Beta Energy Restenosis Trial (BERT) conducted in our institution.

RESULTS

On average, the minimal dose in 90% of the adventitial volume was 37 +/- 16% of the prescribed dose; the minimal dose in 90% of the plaque + media volume was 58 +/- 24% and of the luminal surface volume was 67 +/- 31%. The minimal dose in the 10% most exposed luminal surface volume was 296 +/- 42%. Simulations of the use of a gamma-emitter and/or a radioactive source train centered in the lumen are reported, with a comparison of the homogeneity of the dose distribution.

CONCLUSIONS

It is possible to derive DVH from IVUS, to evaluate the dose delivered to different parts of the coronary wall. This process should improve our understanding of the mechanisms of action of brachytherapy.

Treatment plan evaluation using dose-volume histogram (DVH) and spatial dose-volume histogram (zDVH)

OBJECTIVE

The dose-volume histogram (DVH) has been accepted as a tool for treatment-plan evaluation. However, DVH lacks spatial information. A new concept, the z-dependent dose-volume histogram (zDVH), is presented as a supplement to the DVH in three-dimensional (3D) treatment planning to provide the spatial variation, as well as the size and magnitude of the different dose regions within a region of interest.

MATERIALS AND METHODS

Three-dimensional dose calculations were carried out with various plans for three disease sites: lung, breast, and prostate. DVHs were calculated for the entire volume. A zDVH is defined as a differential dose-volume histogram with respect to a computed tomographic (CT) slice position. In this study, zDVHs were calculated for each CT slice in the treatment field. DVHs and zDVHs were compared.

RESULTS

In the irradiation of lung, DVH calculation indicated that the treatment plan satisfied the dose-volume constraint placed on the lung and zDVH of the lung revealed that a sizable fraction of the lung centered about the central axis (CAX) received a significant dose, a situation that warranted a modification of the treatment plan due to the removal of one lung. In the irradiation of breast with tangential fields, the DVH showed that about 7% of the breast volume received at least 110% of the prescribed dose (PD) and about 11% of the breast received less than 98% PD. However, the zDVHs of the breast volume in each of seven planes showed the existence of high-dose regions of 34% and 15%, respectively, of the volume in the two caudal-most planes and cold spots of about 40% in the two cephalic planes. In the treatment planning of prostate, DVHs showed that about 15% of the bladder and 40% of the rectum received 102% PD, whereas about 30% of the bladder and 50% of the rectum received the full dose. Taking into account the hollow structure of both the bladder and the rectum, the dose-surface histograms (DSH) showed larger hot-spot volume, about 37% of the bladder wall and 43% of the rectal wall. The zDVHs of the bladder revealed that the hot-spot region was superior to the central axis. The zDVHs of the rectum showed that the high-dose region was an 8-cm segment mostly superior to the central axis. The serial array-like of the rectum warrants a closer attention with regard to the complication probability of the organ.

CONCLUSIONS

Although DVH provides an averaged dose-volume information, zDVH provides differential dose-volume information with respect to the CT slice position. zDVH is a 2D analog of a 3D DVH and, in some situations, more superior. It provides additional information on plan evaluation that otherwise could not be appreciated. The zDVH may be used along with DVH for plan evaluation and for the correlation of radiation outcome.

Incorporation of functional status into dose-volume analysis

The dose-volume histogram (DVH) has gained wide acceptance as a mechanism for reducing the voluminous data of a three-dimensional dose distribution into a two-dimensional graph. These graphs are often converted to a single figure of merit. This data reduction technique is used both for clinical treatment plan evaluation and as part of proposed systems for estimating control and complication probabilities. It has long been recognized that a major shortcoming of the DVH as an analysis tool is that all spatial information is discarded. A subtler problem, which is addressed in this work, is that the DVH also implies homogeneity of biological consequence of irradiation in what may be a functionally heterogeneous volume of tissue. An extension to the DVH, the functional dose-volume histogram, or dose-function histogram (DFH), is proposed, that explicitly includes quantitative three-dimensional functional information. The concept is illustrated by the use of SPECT imaging to assess the functional status of irradiated lung.

Quality control of dose volume histogram computation characteristics of 3D treatment planning systems

Detailed quality control (QC) protocols are a necessity for modern radiotherapy departments. The established QC protocols for treatment planning systems (TPS) do not include recommendations on the advanced features of three-dimensional (3D) treatment planning, like the dose volume histograms (DVH). In this study, a test protocol for DVH characteristics was developed. The protocol assesses the consistency of the DVH computation to the dose distribution calculated by the same TPS by comparing DVH parameters with values obtained by the isodose distributions. The computation parameters (such as the dimension of the computation grid) that are applied to the TPS during the tests are not fixed but set by the user as if the test represents a typical clinical case. Six commercial TPS were examined with this protocol within the frame of the EC project Dynarad (Biomed I). The results of the intercomparison prove the consistency of the DVH results to the isodose values for most of the examined TPS. However, special attention should be paid when working with cases of adverse conditions such as high dose gradient regions. In these cases, higher errors are derived, especially when an insufficient number of dose calculation points are used for the DVH computation.

Guidance of intracoronary radiation therapy based on dose-volume histograms derived from quantitative intravascular ultrasound

Application of ionizing radiation to prevent restenosis in atherosclerotic vessels treated by balloon angioplasty is a new treatment under investigation in interventional cardiology and radiology. There is variability in dose prescription, and both gamma- and beta-emitters are used, leading to a wide range of dose distribution over the arterial vessel wall. We present a new modality of dosimetry based on a method that three-dimensional (3-D) image reconstruction of electrocardiogram (ECG)-gated intravascular ultrasound (IVUS) images. Dose volume histograms (DVH) are used to describe the cumulative distribution of dose over two specific volumes: i) at the level of the luminal surface, defined with a thickness of 0.1 mm from the automatically detected contour of the highly echogenic blood-vessel interface, and ii) the adventitia volume is computed considering a 0.5-mm thickness from the echogenic media-adventitia interface. DVH provide a tool for reporting the actual delivered dose at the site believed to be the target: the adventitia, and to detect excessive radiation which could lead to vascular complications. Simulation of a gamma-emitter or of a radioactive source train in the center of the lumen are possible. The data obtained from the first ten patients included in the beta-irradiation trial (BERT 1.5) conducted in our institution are presented, supporting the use of DVH based on quantitative IVUS measurements for optimal dose prescription in vascular interventional radiation therapy.

Developing a dose-volume histogram computation program for brachytherapy

A dose-volume histogram (DVH) computation program was developed for brachytherapy treatment planning in an attempt to benefit from the DVH's ability to present graphically information on 3D dose distributions. The program is incorporated into a planning system that utilizes a pair of orthogonal radiographs to localize the radiation sources. DVHs are calculated for the volume of tissue enclosed by an isodose surface (e.g. half the value of the reference isodose). The calculation algorithm is based on a non-uniform random sampling that gives a denser point distribution at the centre of the implants. Our program was tested and proved to be fast enough for clinical use and sufficiently accurate (i.e. computation time of 20 s and less than 2% relative error for one point source, for 100,000 calculation points). The accuracy improves when a larger calculation point number is used, but the computation time also increases proportionally. The DVH is presented in the form of a simple graph or table, or as Anderson's 'natural' DVH graph. The cumulative DVH tables can be used to extract a series of indexes characterizing the homogeneity and the dose levels of the distribution in the treatment volume and the surrounding tissues. If a reference plan is available, the DVH results can be assessed relative to the reference plan's DVH.

[Conformational radiotherapy with multi-leaf collimators: one year experience at the Leon-Berard Centre]

Taking advantage of the renewal of a linear accelerator, the Radiation Therapy Department of the Centre Léon Bérard implemented, in collaboration with Philips Systèmes Médicaux, a conformal therapy set-up procedure using CT-scan for 3D treatment planning and a multileaf collimator that allows achievement of numerous irregular-shaped beams via the multileaf preparation system. The various elements of this equipment make possible well defined and structured procedures for treatment planning with different steps and essential tools used by this technique. We describe the means used and indicate future improvements that will lead to automation in order to provide good quality assurance, better security and substantial time saving. During the first year, 115 patients were treated with this new technique. They presented with central nervous system tumors (32 patients), lung cancer (29 patients), prostate cancer (20 patients), paranasal sinus tumors (14 patients) and tumors located in other sites (13 patients with soft sarcoma, hepato-bilary tumor, etc).

[Evaluation tests of computer systems concerning tri-dimensional dose calculations]

The development of irradiation techniques in radiotherapy shows a clear tendency towards the systematic use of three-dimensional (3D) information. Great efforts are being made to set up 3D conformal radiotherapy. Consequently, in the aim of greater coherence and accuracy, "the dosimetric tool" must also meet the requirements of 3D radiotherapy, as it plays a role in the treatment chain. To know if the treatment planning system is a "3D", "2D" or even "1D" system, one should not be satisfied with reading the technical documentation and the program algorithm description nor entirely trust the constructor's assertions. It is essential to clearly and precisely evaluate the possibilities of the treatment planning system. Even if it is proved not to satisfy perfectly all the tests which would qualify it as a real 3D calculation system, the study of the test results helps to give clear explanations of the dosimetric results. Two series of test cases are proposed. The first series allows us to understand in which conditions the treatment planning system takes into account the scatter influence in a volume. The second series is designed to inform us about the capability of the dose calculation algorithm when the medium encloses non-homogeneities. These test cases do not constitute an exhaustive "check-list" able to tackle completely the question of 3D calculation. They are submitted as examples and should be considered as an evaluation methodology for the software implanted in the treatment planning system.

The use of medical images in planning and delivery of radiation therapy

The authors provide a survey of how images are used in radiation therapy to improve the precision of radiation therapy plans, and delivery of radiation treatment. In contrast to diagnostic radiology, where the focus is on interpretation of the images to decide if disease is present, radiation therapy quantifies the extent of the region to be treated, and relates it to the proposed treatment using a quantitative modeling system called a radiation treatment planning (RTP) system. This necessitates several requirements of image display and manipulation in radiation therapy that are not usually important in diagnosis. The images must have uniform spatial fidelity: i.e., the pixel size must be known and consistent throughout individual images, and between spatially related sets. The exact spatial relation of images in a set must be known. Radiation oncologists draw on images to define target volumes; dosimetrists use RTP systems to superimpose quantitative models of radiation beams and radiation dose distributions on the images and on the sets of organ and target contours derived from them. While this mainly uses transverse cross-sectional images, projected images are also important, both those produced by the radiation treatment simulator and the treatment machines, and so-called "digital reconstructed radiographs," computed from spatially related sets of cross-sectional images. These requirements are not typically met by software produced for radiologists but are addressed by RTP systems. This review briefly summarizes ongoing work on software development in this area at the University of Washington Department of Radiation Oncology.

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.

Dose-volume histograms for bladder and rectum

PURPOSE

A careful examination of the foundation upon which the concept of the Dose-Volume Histogram (DVH) is built, and the implications of this set of parameters on the clinical application and interpretation of the DVH concept has not been conducted since the introduction of DVHs as a tool for the quantitative evaluation of treatment plans. The purpose of the work presented herein is to illustrate problems with current methods of implementing and interpreting DVHs when applied to hollow anatomic structures such as the bladder and rectum.

METHODS AND MATERIALS

A typical treatment plan for external beam irradiation of a patient with prostate cancer was chosen to provide a data set from which DVH curves for both the bladder and rectum were calculated. The two organs share the property of being shells with contents that are of no clinical importance. DVHs for both organs were computed using a solid model and using a shell model. Typical treatment plans for prostate cancer were used to generate DVH curves for both models. The Normal Tissue Complication Probability (NTCP) for these organs is discussed in this context.

RESULTS

For an eight-field conformal treatment plan of the prostate, a bladder DVH curve generated using the shell model is higher than the corresponding curve generated using the solid model. The shell model also has a higher NTCP. A six-field conformal treatment plan also results in a higher DVH curve for the shell model. A treatment plan consisting of bilateral 120-degree arcs, results in a higher DVH curve for the shell model, as well as a higher NTCP.

CONCLUSION

The DVH concept currently used in evaluation of treatment plans is problematic because current practices of defining exactly what constitutes "bladder" and "rectum." Commonly used methods of tracing the bladder and rectum imply use of a solid structure model for DVHs. In reality, these organs are shells and the critical structure associated with NTCP is obviously and indisputably the shell, as opposed to its contents. Treatment planning algorithms for DVH computation should thus be modified to utilize the shell model for these organs.

Functional dose-volume histograms for functionally heterogeneous normal organs

Functional dose-volume histograms are proposed as an extension of the conventional dose-volume histograms, for quantitative assessment of three-dimensional radiation dose coverage of functionally heterogeneous normal organs. Examples are given to illustrate possible applications of this approach to the treatment of a brain tumour or a lung tumour, in which cases the distribution of the normal organ function can be obtained from functional dose-volume modalities. It is shown that a significant difference exists between the functional dose-volume histograms and the conventional dose-volume histograms when the normal organ function is non-uniformly distributed within the organ. Utilization of functional dose-volume histograms as the input for the calculation of normal tissue complication probabilities is discussed for different normal tissue structures.

A numerical simulation of organ motion and daily setup uncertainties: implications for radiation therapy

PURPOSE

In radiotherapy planning, the clinical target volume (CTV) is typically enlarged to create a planning target volume (PTV) that accounts for uncertainties due to internal organ and patient motion as well as setup error. Margin size clearly determines the volume of normal tissue irradiated, yet in practice it is often given a set value in accordance with a clinical precedent from which variations are rare. The (CTV/PTV) formalism does not account for critical structure dose. We present a numerical simulation to assess (CTV) coverage and critical organ dose as a function of treatment margins in the presence of organ motion and physical setup errors. An application of the model to the treatment of prostate cancer is presented, but the method is applicable to any site where normal tissue tolerance is a dose-limiting factor.

METHODS AND MATERIALS

A Monte Carlo approach was used to simulate the cumulative effect of variation in overall tumor position, for individual treatment fractions, relative to a fixed distribution of dose. Distributions of potential dose-volume histograms (DVHs), for both tumor and normal tissues, are determined that fully quantify the stochastic nature of radiotherapy delivery. We introduce the concept of Probability of Prescription Dose (PoPD) isosurfaces as a tool for treatment plan optimization. Outcomes resulting from current treatment planning methods are compared with proposed techniques for treatment optimization. The standard planning technique of relatively large uniform margins applied to the CTV, in the beam's eye view (BEV), was compared with three other treatment strategies: (a) reduced uniform margins, (b) nonuniform margins adjusted to maximize normal tissue sparing, and (c) a reduced margin plan in which nonuniform fluence profiles were introduced to compensate for potential areas of reduced dose.

RESULTS

Results based on 100 simulated full course treatments indicate that a 10 mm CTV to PTV margin, combined with an additional 5 mm dosimetric margin, provides adequate CTV coverage in the presence of known treatment uncertainties. Nonuniform margins can be employed to reduce dose delivered to normal tissues while preserving CTV coverage. Nonuniform fluence profiles can also be used to further reduce dose delivered to normal tissues, though this strategy does result in higher dose levels delivered to a small volume of the CTV and normal tissues.

CONCLUSIONS

Monte Carlo-based treatment simulation is an effective means of assessing the impact of organ motion and daily setup error on dose delivery via external beam radiation therapy. Probability of Prescription Dose (PoPD) isosurfaces are a useful tool for the determination of nonuniform beam margins that reduce dose delivered to critical organs while preserving CTV dose coverage. Nonuniform fluence profiles can further alter critical organ dose with potential therapeutic benefits. Clinical consequences of this latter approach can only be assessed via clinical trials.

Reporting and analyzing dose distributions: a concept of equivalent uniform dose

Modern treatment planning systems for three-dimensional treatment planning provide three-dimensionally accurate dose distributions for each individual patient. These data open up new possibilities for more precise reporting and analysis of doses actually delivered to irradiated organs and volumes of interest. A new method of summarizing and reporting inhomogeneous dose distributions is reported here. The concept of equivalent uniform dose (EUD) assumes that any two dose distributions are equivalent if they cause the same radiobiological effect. In this paper the EUD concept for tumors is presented, for which the probability of local control is assumed to be determined by the expected number of surviving clonogens, according to Poisson statistics. The EUD can be calculated directly from the dose calculation points or, from the corresponding dose-volume distributions (histograms). The fraction of clonogens surviving a dose of 2 Gy (SF2) is chosen to be the primary operational parameter characterizing radiosensitivity of clonogens. The application of the EUD concept is demonstrated on a clinical dataset. The causes of flattening of the observed dose-response curves become apparent since the EUD concept reveals the finer structure of the analyzed group of patients in respect to the irradiated volumes and doses actually received. Extensions of the basic EUD concept to include nonuniform density of clonogens, dose per fraction effects, repopulation of clonogens, and inhomogeneity of patient population are discussed and compared with the basic formula.

Real-time 3D dose calculation and display: a tool for plan optimization

PURPOSE

Both human and computer optimization of treatment plans have advantages; humans are much better at global pattern recognition, and computers are much better at detailed calculations. A major impediment to human optimization of treatment plans by manipulation of beam parameters is the long time required for feedback to the operator on the effectiveness of a change in beam parameters. Our goal was to create a real-time dose calculation and display system that provides the planner with immediate (fraction of a second) feedback with displays of three-dimensional (3D) isodose surfaces, digitally reconstructed radiographs (DRRs), dose-volume histograms, and/or a figure of merit (FOM) (i.e., a single value plan score function). This will allow the experienced treatment planner to optimize a plan by adjusting beam parameters based on a direct indication of plan effectiveness, the FOM value, and to use 3D display of target, critical organs, DRRs, and isodose contours to guide changes aimed at improving the FOM value.

METHODS AND MATERIALS

We use computer platforms that contain easily utilized parallel processors and very tight coupling between calculation and display. We ported code running on a network of two workstations and an array of transputers to a single multiprocessor workstation. Our current high-performance graphics workstation contains four 150-MHz processors that can be readily used in a shared-memory multithreaded calculation.

RESULTS

When a 10 x 10-cm beam is moved, using an 8-mm dose grid, the full 3D dose matrix is recalculated using a Bentley-Milan-type dose calculation algorithm, and the 3D dose surface display is then updated, all in < 0.1s. A 64 x 64-pixel DRR calculation can be performed in < 0.1 s. Other features, such as automated aperture calculation, are still required to make real-time feedback practical for clinical use.

CONCLUSION

We demonstrate that real-time plan optimization using general purpose multiprocessor workstations is a practical goal. Parallel processing technology provides this capability for 3D planning systems, and when combined with objective plan ranking algorithms should prove effective for optimizing 3D conformal radiation therapy. Compared to our earlier transputer work, multiprocessor workstations are more easily programmed, making software development costs more reasonable compared with uniprocessor development costs. How the dose calculation is partitioned into parallel tasks on a multiprocessor work station can make a significant difference in performance. Shared-memory multiprocessor workstations are our first choice for future work, because they require minimum programming effort and continue to be driven to higher performance by competition in the workstation arena.

The "critical volume tolerance method" for estimating the limits of dose escalation during three-dimensional conformal radiotherapy for prostate cancer

PURPOSE

To describe the "Critical Volume Tolerance" (CVT) method for defining normal tissue tolerance during 3D-based dose escalation studies for prostate cancer.

METHODS AND MATERIALS

The CVT method predicts the tolerance to radiation for "in series"-type functional units based on the assumption that tolerance depends on a critical threshold "low-volume high-dose region." The data used for describing this model were generated from 3D analysis of randomly selected patients with prostate cancer. Commonly used coplanar four-and six-field conformal (SFC) techniques were chosen as the comparison techniques. For purposes of comparison, rectal tolerance was assumed to be reached following whole pelvic irradiation using a four-field box technique to 50 Gy, followed by a conedown boost to 70 Gy using bilateral 9 x 9 cm 120 degree arcs as popularized by investigators from Stanford University (SUH).

RESULTS

Based on the average dose volume histograms for the patients studied, the maximum safe increase in dose for the SFC technique compared to the SUH technique, would be 10% if 30% of the rectal volume was the critical dose limiting volume (CVT = 30%), 5% if the CVT = 10%, or greater than 20% if the CVT = 40%. Commonly used four-field conformal techniques would not be expected to allow significant escalation of the dose without increasing the risk of complications.

CONCLUSIONS

The CVT method is relatively simple, and data generated based on it can be used to support normal tissue complication probability equations. The CVT method can be verified or modified as partial tolerance data become available. Based on the CVT model, sophisticated treatment techniques should allow a modest increase in the total dose of radiation delivered to the prostate without an increase in late complications.

Spatial information on dose distribution using multisectional dose-volume histograms

Dose-volume histograms are useful tools to summarize the information on the dose profiles resulting within a target volume. However, the spatial relationships of the hot and the cold spots are blunted in the dose-volume histograms. This study tries to circumvent this problem using multisectional dose-volume histograms and highlight the utility of these in the optimization of a radiation therapy plan.

Conventional vs. conformal radiotherapy for prostate cancer: preliminary results of dosimetry and acute toxicity

PURPOSE

To compare conformal radiotherapy using three dimensional treatment planning (3D-CRT) to conventional radiotherapy (Conven-RT) for patients with Stages T2-T4 adenocarcinoma of the prostate.

METHODS AND MATERIALS

A Phase III randomized study was activated in May 1993, to compare treatment toxicity and patient outcome after 78 Gy in 39 fractions using 3D-CRT to that after 70 Gy in 35 fractions using Conven-RT. The first 46 Gy were administered using the same nonconformal field arrangement (four field) in both arms. The boost was given nonconformally using four fields in the Conven-RT arm and conformally using six fields in the 3D-CRT arm. The dose was specific to the isocenter. The first 60 patients, 29 in the 3D-CRT arm and 31 in the Conven-RT arm, are the subject of this preliminary analysis.

RESULTS

The two treatment arms were first compared in terms of dosimetry by dose-volume histogram analysis. Using a subgroup of patients in the 3D-CRT arm (n=15), both Conven-RT and 3D-CRT plans were generated and the dose-volume histogram data compared. The mean volumes treated to doses above 60 Gy for the bladder and rectum were 28 and 36% for the 3D-CRT plans, and 43 and 38% for the Conven-RT plans, respectively (p < 0.05 for the bladder volumes). The mean clinical target volume (prostate and seminal vesicles) treated to 95% of the prescribed dose was 97.5% for the 3D-CRT arm, and 95.6% for the Conven-RT arm (p < 0.05). There were no significant differences in the acute reactions between the two arms, with the majority experiencing Grade 2 or less toxicity (92%). Moreover, no relationship was seen between acute toxicity and the volume of bladder and rectum receiving in excess of 60 Gy for those in the 3D-CRT arm. There was also no difference between the groups in terms of early biochemical response. Prostate-specific antigen levels at 3 and 6 months after completion of radiotherapy were similar in the two treatment arms. There was only one biochemical failure in the study population at the time of the analysis.

CONCLUSIONS

Comparison of the Conven-RT and 3D-RT treatment plans revealed that significantly less bladder was in the high dose volume in the 3D-CRT plans, while the volume of rectum receiving doses over 60 Gy was equivalent. There were no differences between the two treatment arms in terms of acute toxicity or early biochemical response. Longer follow-up is needed to determine the impact of 3D-CRT on long-term patient outcome and late reactions.

Head and neck cancer

This synthesis of the literature on radiotherapy for head and neck cancer is based on 424 scientific articles, including 3 meta-analyses, 38 randomized studies, 45 prospective studies, and 246 retrospective studies. These studies involve 79174 patients. The literature review shows that radiotherapy, either alone or in combination with surgery, plays an essential role in treating head and neck cancers. When tumors are localized, many tumor patients can be cured by radiotherapy alone and thereby maintain full organ function (1, 2). Current technical advancements in radiotherapy offer the potential for better local tumor control with lower morbidity (3). This, however, will require more sophisticated dose planning resources. To further improve treatment results for advanced tumors, other fractionation schedules, mainly hyperfractionation, should be introduced (5). This mainly increases the demands on staff resources for radiotherapy. The combination of radiotherapy and chemotherapy should be subjected to further controlled studies involving a sufficiently large number of patients (4, 5). Interstitial treatment (in the hands of experienced radiotherapists) yields good results for selected cancers. The method should be more generally accessible in Sweden. Intraoperative radiotherapy should be targeted for further study and development.

Three dimensional conformal radiation therapy in pediatric parameningeal rhabdomyosarcomas

PURPOSE

We evaluated the utility of three dimensional (3D) treatment planning in the management of children with parameningeal head and neck rhabdomyosarcomas.

METHODS AND MATERIALS

Five children with parameningeal rhabdomyosarcoma were referred for treatment at our radiation oncology center from May 1990 through January 1993. Each patient was evaluated, staged, and treated according to the Intergroup Rhabdomyosarcoma Study. Patients were immobilized and underwent a computed tomography scan with contrast in the treatment position. Tumor and normal tissues were identified with assistance from a diagnostic radiologist and defined in each slice. The patients were then planned and treated with the assistance of a 3D treatment planning system. A second plan was then devised by another physician without the benefit of the 3D volumetric display. The target volumes designed with the 3D system and the two-dimensional (2D) method were then compared. The dosimetric coverage to tumor, tumor plus margin, and normal tissues was also compared with the two methods of treatment planning.

RESULTS

The apparent size of the gross tumor volume was underestimated with the conventional 2D planning method relative to the 3D method. When margin was added around the gross tumor to account for microscopic extension of disease in the 2D method, the expected area of coverage improved relative to the 3D method. In each circumstance, the minimum dose that covered the gross tumor was substantially less with the 2D method than with the 3D method. The inadequate dosimetric coverage was especially pronounced when the necessary margin to account for subclinical disease was added. In each case, the 2D plans would have delivered substantial dose to adjacent normal tissues and organs, resulting in a higher incidence of significant complications.

CONCLUSIONS

3D conformal radiation therapy has a demonstrated advantage in the treatment of sarcomas of the head and neck. The improved dosimetric coverage of the tumor and its margin for subclinical extensions may result in improvement in local control of these tumors. In addition, lowering of radiation dose to adjacent critical structures may help lower the incidence of adverse late effects in children.

Comparison of treatment techniques for conformal radiotherapy of the prostate using dose-volume histograms and normal tissue complication probabilities

The aim of this study was to evaluate the relative merits of the coplanar field arrangements most frequently used for conformal radiotherapy of the prostate using dose-volume histograms and normal tissue complication probabilities (NTCPs). Twelve patients with early prostate cancer underwent a planning CT scan of the pelvis. Isocentric plans for each patient were devised using three, four, six and eight conformal fields and beam-weights optimised using fast simulated annealing to give a dose homogeneity across the planning target volume of +/- 5% or better while minimising irradiation of the relevant organs at risk. The plans were then evaluated using dose-volume histograms of the organs at risk (bladder, rectum and both femoral heads) and the Lyman model of normal tissue complication probability for the rectum. Analysis of dose-volume histogram data averaged over the 12 patients indicates an advantage for six (p = 0.002) and eight (p = 0.0001) fields with respect to the percentage volume of the femoral heads receiving > 50% of the prescribed dose compared with three fields. There was a similar advantage for six (p = 0.0007) and eight (p = 0.0001) fields compared with four fields. Ranking of the treatment techniques indicates that the four-field technique is the worst with respect to femoral head irradiation but the best with respect to reducing rectal irradiation. A higher dose can be prescribed to the isocentre with the four-field technique for a 5% rectal NTCP. The six-field technique led to sparing of the bladder when the different treatment techniques were ranked using bladder dose-volume histogram data. We conclude that none of the techniques studied consistently proved to be superior when applied to this sample of patients with prostate cancer with respect to sparing all the organs at risk. The absolute differences between techniques are small and would be very difficult to detect with respect to clinically relevant endpoints.

Dosimetric considerations in treating mediastinal disease with mantle fields: characterization of the dose under mantle blocks

PURPOSE

While the rationale for using mantle fields is well understood and the prescription of these fields is straightforward, the underlying complexity of the dose distributions that result is not generally appreciated. This is especially true in the choice of lung block design, which affects the dose to both the target volume as well as to the normal lung tissue. The key to the design of optimal lung blocks is the physician's perception of the complex relationship between the geometric and dosimetric aspects of heavily modified fields, as well as how the physical and anatomical properties of the target volume and the shape of the patient's lungs relate to the images visualized on simulator films.

METHODS AND MATERIALS

Depth doses and cross-beam profiles of blocks ranging in width from 1 cm to 10 cm were taken using an automated beam scanning system. These data were then converted to "shadow fields." The results were compared to open fields of the same size using standard methodology.

RESULTS

Shadow fields behave quite similarly to small, open fields in terms of x-ray-light field congruence, flatness, symmetry, and penumbra. There is a 2-3 mm rim between the edge of the block and the point at which it becomes nominally effective. The dose at the center of a block, which gives the normalization of the shadow fields, is given by a block transmission factor (BTF), which produces results in excellent agreement with measurements over a wide variety of block sizes and tissue depths.

CONCLUSION

The radiation dose under shielding blocks can be considerably higher than expected, and care must be exercised when drawing blocks close to critical structures. The effects of blocks can be described in terms of normalized shadow fields, which behave similar to narrow, open fields, but with a divergence characteristic of their position relative to the radiation source. The normalization value for these fields, which gives the relative dose under the block, can be obtained from a straightforward analytical expression, the BTF.

Dose-volume histogram analysis of techniques for irradiating pituitary adenomas

PURPOSE

Three-dimensional treatment planning was performed to evaluate three standard coplanar irradiation techniques (two-field parallel-opposed, three-field, and 110 degrees bilateral arcs), the 330 degrees single rotational arc, and a four noncoplanar arc technique for the treatment of pituitary adenomas. We sought to identify the optimal technique for minimizing the dose delivered to the normal tissues around the pituitary gland.

METHODS AND MATERIALS

Contours of the pituitary tumor and normal tissues were traced onto computed axial tomography (CT) scans and reconstructed in three dimensions using a three-dimensional planning system. A total dose of 45 Gy was delivered to the pituitary lesion with the five techniques using 6 MV and 18 MV photons, and dose-volume histograms were generated.

RESULTS

The 18 MV photons delivered a lower dose to the temporal lobe than did the 6 MV photons in the two-field technique, but this advantage was not evident for the other techniques. The three-field technique improved dose distribution throughout the temporal lobes with low doses being delivered to the frontal lobe. The bilateral arc and the 330 degrees arc techniques were superior to stationary two- and three-fields techniques for sparing the temporal lobes. The four noncoplanar arc technique delivered less doses to the temporal and frontal lobes than did the other techniques. However, the lens dose (3.6 Gy/25 fractions) was higher compared to the other techniques.

CONCLUSION

Analysis of the dose-volume histograms shows the various dosimetric advantages and disadvantages of the five techniques. Based upon individual considerations, including the patient's age and medical history, one can decide the optimal technique for treatment.

External beam radiotherapy of localized prostatic adenocarcinoma. Evaluation of conformal therapy, field number and target margins

The purpose of the present study was to identify factors of importance in the planning of external beam radiotherapy of prostatic adenocarcinoma. Seven patients with urogenital cancers were planned for external radiotherapy of the prostate. Four different techniques were used, viz. a 4-field box technique and four-, five- or six-field conformal therapy set-ups combined with three different margins (1-3 cm). The evaluations were based on the doses delivered to the rectum and the urinary bladder. A normal tissue complication probability (NTCP) was calculated for each plan using Lyman's dose volume reduction method. The most important factors that resulted in a decrease of the dose delivered to the rectum and the bladder were the use of conformal therapy and smaller margins. Conformal therapy seemed more important for the dose distribution in the urinary bladder. Five- and six-field set-ups were not significantly better than those with four fields. NTCP calculations were in accordance with the evaluation of the dose volume histograms. To conclude, four-field conformal therapy utilizing reduced margins improves the dose distribution to the rectum and the urinary bladder in the radiotherapy of prostatic adenocarcinoma.

Non-uniform dwell times in line source high dose rate brachytherapy: physical and radiobiological considerations

The ability to vary source dwell times in high dose rate (HDR) brachytherapy allows for the use of non-uniform dwell times along a line source. This may have advantages in the radical treatment of tumours depending on individual tumour geometry. This study investigates the potential improvements in local tumour control relative to adjacent normal tissue isoeffects when intratumour source dwell times are increased along the central portion of a line source (technique A) in radiotherapy schedules which include a relatively small component of HDR brachytherapy. Such a technique is predicted to increase the local control for tumours of diameters ranging between 2 cm and 4 cm by up to 11% compared with a technique in which there are uniform dwell times along the line source (technique B). There is no difference in the local control rates for the two techniques when used to treat smaller tumours. Normal tissue doses are also modified by the technique used. Technique A produces higher normal tissue doses at points perpendicular to the centre of the line source and lower doses at points nearer the ends of the line source if the prescription point is not in the central plane of the line source. Alternatively, if the dose is prescribed at a point in the central plane of the line source, the dose at all the normal tissue points are lower when technique A is used.

Integrated software tools for the evaluation of radiotherapy treatment plans

PURPOSE

This article announces the availability of a convenient and useful software environment for the evaluation of three-dimensional (3D) radiotherapy treatment plans.

MATERIALS AND METHODS

Using standards such as American National Standards for Information Systems C and the X Window System allowed us to bring the computation and display of dose-volume histograms, dose statistics, tumor control probabilities, normal tissue complication probabilities, and a figure of merit together under one user interface. These plan evaluation tools are not stand alone, but must interact with a 3D radiation therapy planning system to obtain the required dose matrices and patient anatomical contours. Installation of the software involves a programmer who writes a software bridge between the radiation therapy planning system and the tools, thereby providing access to local data files. This design strategy confines portability issues to one area of the software.

RESULTS

Access to the other tools is through the Graphical Plan Evaluation Tool (GPET). GPET coordinates the use of each of the tools and provides graphical facilities for display of their results. Importantly, GPET assures that the displayed results of each tool have been computed with the same input specifications for all treatment plans being compared. For added convenience, the user can rearrange the resultant data to be reviewed in various ways on the video screen. The software design also allows incorporation of customized algorithms and input data for computing tumor control probability and normal tissue complication probabilities, since those currently available are controversial.

CONCLUSION

The Graphical Plan Evaluation Tool unifies the simultaneous computation for several analytical tools and graphical display of their results. Within the constraints of the X Window System environment, this assemblage of software tools provides a portable, flexible, and convenient method for the quantitative evaluation of several radiotherapy treatment plans.

Assurance of high quality linac-based stereotactic radiosurgery

PURPOSE

Stereotactic radiosurgery is generally a single, high-dose radiation treatment for the brain requiring targeting accuracy on the order of a millimeter. From the initial implementation of radiosurgery, therefore, quality assurance is an ongoing process of paramount importance. In this paper, we outline the basic elements of a quality assurance program for our linear accelerator that has been in use at Washington University Medical Center over the past 2 years.

METHODS AND MATERIALS

Various devices and procedures have been developed to verify the accuracy and safety of the stereotactic radiosurgery regimen. Specifically, we present methods for assessing the attainment of spatially correct patient images, the reliability of the computerized treatment planning system, achieving physical safety for the patient, as well as the proper operation of the radiation treatment device.

RESULTS

Our procedures have allowed us to assure quality patient treatments and, additionally, has permitted monitoring our performance for continual improvement. For example, a plot of targeting accuracy with the number of patients shows an asymptotic approach to a value within 0.6 mm of that ideally expected.

CONCLUSION

To maintain high-quality patient care, one must review critical aspects of the treatment regimen on a periodic basis. Providing for the appropriate level of staff training, periodic reviews of procedures and maintenance of forms are also very important.