Tomotherapy: a new concept for the delivery of dynamic conformal radiotherapy

Cone-beam CT with megavoltage beams and an amorphous silicon electronic portal imaging device: potential for verification of radiotherapy of lung cancer

We investigate the potential of megavoltage (MV) cone-beam CT with an amorphous silicon electronic portal imaging device (EPID) as a tool for patient position verification and tumor/organ motion studies in radiation treatment of lung tumors. We acquire 25 to 200 projection images using a 22 x 29 cm EPID. The acquisition is automatic and requires 7 minutes for 100 projections; it can be synchronized with respiratory gating. From these images, volumetric reconstruction is accomplished with a filtered backprojection in the cone-beam geometry. Several important prereconstruction image corrections, such as detector sag, must be applied. Tests with a contrast phantom indicate that differences in electron density of 2% can be detected with 100 projections, 200 cGy total dose. The contrast-to-noise ratio improves as the number of projections is increased. With 50 projections (100 cGy), high contrast objects are visible, and as few as 25 projections yield images with discernible features. We identify a technique of acquiring projection images with conformal beam apertures, shaped by a multileaf collimator, to reduce the dose to surrounding normal tissue. Tests of this technique on an anthropomorphic phantom demonstrate that a gross tumor volume in the lung can be accurately localized in three dimensions with scans using 88 monitor units. As such, conformal megavoltage cone-beam CT can provide three-dimensional imaging of lung tumors and may be used, for example, in verifying respiratory gated treatments.

Methods for improving limited field-of-view radiotherapy reconstructions using imperfect a priori images

There are many benefits to having an online CT imaging system for radiotherapy, as it helps identify changes in the patient's position and anatomy between the time of planning and treatment. However, many current online CT systems suffer from a limited field-of-view (LFOV) in that collected data do not encompass the patient's complete cross section. Reconstruction of these data sets can quantitatively distort the image values and introduce artifacts. This work explores the use of planning CT data as a priori information for improving these reconstructions. Methods are presented to incorporate this data by aligning the LFOV with the planning images and then merging the data sets in sinogram space. One alignment option is explicit fusion, producing fusion-aligned reprojection (FAR) images. For cases where explicit fusion is not viable, FAR can be implemented using the implicit fusion of normal setup error, referred to as normal-error-aligned reprojection (NEAR). These methods are evaluated for multiday patient images showing both internal and skin-surface anatomical variation. The iterative use of NEAR and FAR is also investigated, as are applications of NEAR and FAR to dose calculations and the compensation of LFOV online MVCT images with kVCT planning images. Results indicate that NEAR and FAR can utilize planning CT data as imperfect a priori information to reduce artifacts and quantitatively improve images. These benefits can also increase the accuracy of dose calculations and be used for augmenting CT images (e.g., MVCT) acquired at different energies than the planning CT.

Reshapable physical modulator for intensity modulated radiation therapy

A new method of generating beam intensity modulation filters for intensity modulated radiation therapy (IMRT) is presented. The modulator was based on a reshapable material, which is not compressible but can be deformed under pressure. A two-dimensional (2D) piston array was used to repeatedly shape the attenuating material. The material is a mixture of tungsten powder and a silicon-based binder. The linear attenuation coefficient of the material was measured to be 0.409 cm(-1) for a 6 MV x-ray beam. The maximum thickness of the physical modulator is 10.2 cm, allowing a transmission of 1.5%. A 16 x 16 square piston array was used to generate a depth pattern in the deformable attenuating material. Each piston has a cross section of 6.37 x 6.37 mm2. The modulator was placed 65 cm from the radiation source of the linear accelerator in the position of the shielding tray. At this position, each piston projects to a 1.0 x 1.0 cm2 area at the isocenter, giving a treatment field of 16 x 16 cm2. The percent depth dose curve and output factor measurement show a slight beam hardening and a 1%-4% increase in scatter fraction when 2.2-4.4 cm uniform thickness filters are in the beam. The surface dose was decreased with the filter in the beam. Ion chamber and verification films were used to verify the entrance dose. The measured absolute and relative doses were compared with the calculated dose. The agreement of measurements and calculations is within 3%. In order to verify the spatial modulation of dose, 1-D dose profiles were obtained using dose calculations. Calculated and measured profiles were compared. The 20%-80% penumbra of the modulator was measured to be 5.5-10 mm. The results show that a physical modulator formed using a 16 x 16 piston array and a deformable attenuation material can provide intensity modulation for IMRT comparable with those provided by currently available commercial MLC techniques.

Intensity modulated arc deliveries approximated by a large number of fixed gantry position sliding window dynamic multileaf collimator fields

The intensity modulated arc has been proposed as an alternative to tomotherapy. Treatment planing systems more typically model the conventional step and shoot or sliding window dynamic multileaf collimator (DMLC) deliveries, and may not support intensity modulated arc therapy (IMAT). As well, another potential drawback to this technique is that increasing the number of intensity levels required to achieve certain dose distributions necessitates increasing the number of gantry passes, as may occur if the desired dose distribution is complex (e.g., concave or bifurcated), potentially increasing the overall treatment time. A technique is presented here for the delivery of tomotherapy like dose distributions in a single gantry pass by the use of a large number of fields modulated by a sliding window DMLC technique from fixed equally spaced gantry positions. This serves as a good approximation to either IMAT or tomotherapy deliveries. The planning of these fields is achieved using iterative filtered back projection. Measured results of deliveries of varying degrees of complexity on a homogeneous phantom are compared to desired distributions.

Re-optimization in adaptive radiotherapy

In routine clinical practice, radiotherapy treatment planning is performed based on the patient CT images obtained during the patient setup procedure. However, the actual delivered dose to the patient might be different from the planned dose because of various reasons such as patient motion. Under such situations, it is desirable to modify the original treatment plan in order to partially remedy the dose delivery errors in the subsequent dose delivery process. Such modification can be implemented by modifying the original treatment plan using re-optimization. In this work, issues such as the re-optimization dose prescription, optimization constraints in re-optimization, re-optimization in multiple fractionation schemes and re-optimization procedure with generalized dose-based objective functions were investigated and corresponding mathematical schemes proposed. The derived results were applied to a clinical case study in which it was shown that the proposed re-optimization method is able to remedy the cold spots in tumour while delivering low dose to normal structures. Thus the potential effectiveness of the method was demonstrated.

Independent monitor unit calculation for intensity modulated radiotherapy using the MIMiC multileaf collimator

A self-consistent monitor unit (MU) and isocenter point-dose calculation method has been developed that provides an independent verification of the MU for intensity modulated radiotherapy (IMRT) using the MIMiC (Nomos Corporation) multileaf collimator. The method takes into account two unique features of IMRT using the MIMiC: namely the gantry-dynamic arc delivery of intensity modulated photon beams and the slice-by-slice dose delivery for large tumor volumes. The method converts the nonuniform beam intensity planned at discrete gantry angles of 5 degrees or 10 degrees into conventional nonmodulated beam intensity apertures of elemental arc segments of 1 degree. This approach more closely simulates the actual gantry-dynamic arc delivery by MIMiC. Because each elemental arc segment is of uniform intensity, the MU calculation for an IMRT arc is made equivalent to a conventional arc with gantry-angle dependent beam apertures. The dose to the isocenter from each 1 degree elemental arc segment is calculated by using the Clarkson scatter summation technique based on measured tissue-maximum-ratio and output factors, independent of the dose calculation model used in the IMRT planning system. For treatments requiring multiple treatment slices, the MU for the arc at each treatment slice takes into account the MU, leakage and scatter doses from other slices. This is achieved by solving a set of coupled linear equations for the MUs of all involved treatment slices. All input dosimetry data for the independent MU/isocenter point-dose calculation are measured directly. Comparison of the MU and isocenter point dose calculated by the independent program to those calculated by the Corvus planning system and to direct measurements has shown good agreement with relative difference less than +/-3%. The program can be used as an independent initial MU verification for IMRT plans using the MIMiC multileaf collimators.

A review of intensity modulated radiation therapy: incorporating a report on the seventh education workshop of the ACPSEM--ACT/NSW branch. Australasian College of Physical Scientists and Engineers in Medicine

Intensity modulated radiation therapy (IMRT) is an evolving treatment technique that has become a clinical treatment option in several radiotherapy centres around the world. In August 2001 the ACT/NSW branch of the ACPSEM held its seventh education workshop, the subject was IMRT. This review considers the current use of IMRT and reports on the proceedings of the workshop. The workshop provided some of the theory behind IMRT, discussion of the practical issues associated with IMRT, and also involved presentations from Australian centres that had clinically implemented IMRT. The main topics of discussion were patient selection, plan assessment, multi-disciplinary approach, quality assurance and delivery of IMRT. Key points that were emphasised were the need for a balanced multi-disciplinary approach to IMRT, in both the establishment and maintenance of an IMRT program; the importance of the accuracy of the final dose distribution as compared to the minor in-field fluctuations of individual beams; and that IMRT is an emerging treatment technique, undergoing continuing development and refinement.

Fast and accurate leaf verification for dynamic multileaf collimation using an electronic portal imaging device

A prerequisite for accurate dose delivery of IMRT profiles produced with dynamic multileaf collimation (DMLC) is highly accurate leaf positioning. In our institution, leaf verification for DMLC was initially done with film and ionization chamber. To overcome the limitations of these methods, a fast, accurate and two-dimensional method for daily leaf verification, using our CCD-camera based electronic portal imaging device (EPID), has been developed. This method is based on a flat field produced with a 0.5 cm wide sliding gap for each leaf pair. Deviations in gap widths are detected as deviations in gray scale value profiles derived from the EPID images, and not by directly assessing leaf positions in the images. Dedicated software was developed to reduce the noise level in the low signal images produced with the narrow gaps. The accuracy of this quality assurance procedure was tested by introducing known leaf position errors. It was shown that errors in leaf gap as small as 0.01-0.02 cm could be detected, which is certainly adequate to guarantee accurate dose delivery of DMLC treatments, even for strongly modulated beam profiles. Using this method, it was demonstrated that both short and long term reproducibility in leaf positioning were within 0.01 cm (1sigma) for all gantry angles, and that the effect of gravity was negligible.

Exploring the limits of spatial resolution in radiation dose delivery

Flexibility and complexity in patient treatment due to advances in radiotherapy techniques necessitates a simple method for evaluating spatial resolution capabilities of the dose delivery device. Our purpose in this investigation is to evaluate a model that describes the ability of a radiation therapy device to deliver a desired dose distribution. The model is based on linear systems theory and is analogous to methods used to describe resolution degradation in imaging systems. A qualitative analysis of spatial resolution degradation using the model is presented in the spatial and spatial frequency domains. The ability of the model to predict the effects of geometric dose conformity to treatment volumes is evaluated by varying multileaf collimator leaf width and magnitude of dose spreading. Dose distributions for three clinical treatment shapes, circular shapes of varying diameter and one intensity modulated shape are used in the evaluation. We show that the model accurately predicts the dependence of dose conformity on these parameters. The spatial resolution capabilities of different radiation therapy devices can be quantified using the model, providing a simple method for comparing different treatment machine characteristics. Also, as different treatment sites have different resolution requirements this model may be used to tailor machine characteristics to the specific site.

Helical tomotherapy: an innovative technology and approach to radiation therapy

Helical tomotherapy represents both a novel radiation treatment device and an innovative means of delivering radiotherapy. The helical tomotherapy unit itself is essentially a hybrid between a linear accelerator and a helical CT scanner for the purpose of delivering intensity-modulated radiation therapy (IMRT). The imaging capacity conferred by the CT component allows targeted regions to be visualized prior to, during, and immediately after each treatment. The megavoltage CT (MVCT) images supplant the port-films used in conventional radiotherapy, providing unprecedented anatomical detail. Image-guidance through MVCT will allow the development and refinement of the concept of "adaptive radiotherapy", the reconstruction of the actual daily delivered dose (as opposed to planned dose) accompanied by prescription and delivery adjustments when appropriate. In addition to this unique feature, helical tomotherapy appears capable of further improvements over 3-dimensional conformal radiation therapy and non-helical IMRT in the specific avoidance of critical normal structures, i.e "conformal avoidance", the counterpart of conformal radiation therapy. Based on radiobiological principles that exploit the physical advantages of helical tomotherapy, several dosimetric and clinical investigations are underway.

Quality assurance of intensity-modulated radiotherapy

Intensity-modulated radiation therapy (IMRT) requires the use of inverse treatment planning and nonuniform fluence beams delivered by a series of complex radiation portals. The quality assurance procedures for conventional three-dimensional conformal radiation therapy (3D-CRT) have been developed and are in worldwide clinical use, but the more complex nature of IMRT limits the application of much of the quality assurance (QA) procedures developed for IMRT. Although consensus has not yet been reached regarding which procedures will eventually become recommended by official organizations, the field is rapidly coming to agreement on a basic set of procedures. This manuscript describes some of the novel techniques recently developed for IMRT QA, both for the validation of the calculated dose distribution and for assuring that the dose distribution reaches its intended targets.

Three-dimensional conformal radiotherapy for lung cancer: promises and pitfalls

Lung cancer represents a major source of morbidity and mortality. Despite recent advances, long-term survival remains elusive in most patients with locally advanced cancer. A substantial proportion of these patients experience a relapse at the original site of disease within the thorax, making radiotherapy an important component of treatment. Of several approaches investigated to improve the therapeutic ratio in radiotherapy, three-dimensional conformal radiotherapy holds the most promise, primarily because it allows higher doses to be delivered to the target by improved shaping of radiation portals and conformal avoidance of normal structures. The rationale and evolution of this technology and its potential pitfalls are presented in this review.

A comparison of physically and radiobiologically based optimization for IMRT

Many optimization techniques for intensity modulated radiotherapy have now been developed. The majority of these techniques including all the commercial systems that are available are based on physical dose methods of assessment. Some techniques have also been based on radiobiological models. None of the radiobiological optimization techniques however have assessed the clinically realistic situation of considering both tumor and normal cells within the target volume. This study considers a ratio-based fluence optimizing technique to compare a dose-based optimization method described previously and two biologically based models. The biologically based methods use the values of equivalent uniform dose calculated for the tumor cells and integral biological effective dose for normal cells. The first biologically based method includes only tumor cells in the target volume while the second considers both tumor and normal cells in the target volume. All three methods achieve good conformation to the target volume. The biologically based optimization without the normal tissue in the target volume shows a high dose region in the center of the target volume while this is reduced when the normal tissues are also considered in the target volume. This effect occurs because the normal tissues in the target volume require the optimization to reduce the dose and therefore limit the maximum dose to that volume.

Evaluation of dose calculation algorithms for dynamic arc treatments of head and neck tumors

BACKGROUND AND PURPOSE

To investigate if the Pencil Beam (PB) algorithm takes the disturbance of the dose distribution due to tissue inhomogeneities sufficiently into account in dynamic field shaping rotation therapy (called the dynamic arc treatment modality) for fractionated stereotactic radiation therapy of head and neck tumors.

MATERIAL AND METHODS

A treatment plan using the dynamic arc treatment modality of an oropharynx lesion on a humanoid phantom was evaluated. The same plan was calculated with three different calculation algorithms: the Clarkson and the PB algorithm (both available on the planning system of the Novalis system used for dynamic arc treatments), and the Collapsed Cone Convolution Superposition (CC) algorithm (used by the Pinnacle planning system). The three resulting plans are compared using isodose distributions and cumulative dose volume histograms (CDVHs). An intercomparison of the results of the three algorithms was performed to investigate how accurately each of them takes the influence of tissue inhomogeneities into account such as bony structures and air cavities often appearing in the head and neck region. Additionally, the resulting plans were compared with absolute and relative dosimetric measurements of the treatment plan on the humanoid phantom with thermoluminescent detectors and radiographic film, respectively.

RESULTS

All calculated dose distributions show a good agreement with the measured distribution except in the planning target volume (PTV) in and at the border of the air cavity. All three algorithms overestimate the dose in the PTV at the boundary with the low-density tissue, with 12, 10 and 7% for the Clarkson, the PB and the CC algorithm, respectively. The correspondence between the calculated dose distributions is reflected in the graphs of the CDVHs. They show similar curves for the PTV and the structures except for the left parotic gland and the myelum.

CONCLUSIONS

The PB algorithm of the Novalis system calculates a treatment plan for the dynamic arc treatment modality adequately for fractionated stereotactic radiation therapy of head and neck tumors, except in the PTV in and at the border of the air cavity where the actual dose is overestimated. Care needs to be taken in clinical cases where it is critical to irradiate the air-tissue boundary to a sufficient dose.

Dosimetric validation for multileaf collimator-based intensity-modulated radiotherapy: a review

The creation of intricate dose distributions produced by intensity-modulated radiotherapy (IMRT) depends on complex planning systems and specialized mechanical devices. The many possible sources of inaccuracy and the complexity of the dose maps themselves require that a substantial effort be made to ensure that calculated and delivered dose distributions agree. This review provides an overview of the current status of the validation of dose predictions of IMRT planning systems by comparisons with measurements. Emphasis is placed on multileaf collimator- (MLC) based IMRT. Discrepancies between calculations and measurements may be due to any of 3 causes: errors and uncertainties in the dose calculation algorithm, in measurements, or in beam delivery by the accelerator/MLC combination. Some of the factors affecting dosimetry include: the technique employed for modulating the fluence, the dose calculation algorithm and other aspects of the planning system, mechanical limitations of the MLC hardware, dosimetric characteristics of the MLC, such as MLC leakage and rounded leaf ends, the choice of dosimeter, and the measurement geometry and technique. The advantages and drawbacks of various dosimeters including film, ion chambers, thermoluminescent dosimetry, and electronic portal imaging devices are discussed. The steps involved in validating dosimetrically a planning system are outlined, including the various fields that need to be measured, the phantoms that may be used, and measurement techniques. The achievable accuracy of dosimetry for IMRT is discussed.

Monte Carlo simulation for MLC-based intensity-modulated radiotherapy

This article describes photon beam Monte Carlo simulation for multi leaf collimator (MLC)-based intensity-modulated radiotherapy (IMRT). We present the general aspects of the Monte Carlo method for the non-Monte Carloist with an emphasis given to patient-specific radiotherapy application. Patient-specific application of the Monte Carlo method can be used for IMRT dose verification, inverse planning, and forward planning in conventional conformal radiotherapy. Because it is difficult to measure IMRT dose distributions in heterogeneous phantoms that approximate a patient, Monte Carlo methods can be used to verify IMRT dose distributions that are calculated using conventional methods. Furthermore, using Monte Carlo as the dose calculation method for inverse planning results in better-optimized treatment plans. We describe both aspects and present our recent results to illustrate the discussion. Finally, we present current issues related to clinical implementation of Monte Carlo dose calculation. Monte Carlo is the most recent, and most accurate, method of radiotherapy dose calculation. It is currently in the process of being implemented by various treatment planning vendors and will be available for clinical use in the immediate future.

Clinical implementation of intensity-modulated arc therapy

PURPOSE

Intensity-modulated arc therapy (IMAT) is a method for delivering intensity-modulated radiation therapy (IMRT) using rotational beams. During delivery, the field shape, formed by a multileaf collimator (MLC), changes constantly. The objectives of this study were to (1) clinically implement the IMAT technique, and (2) evaluate the dosimetry in comparison with conventional three-dimensional (3D) conformal techniques.

METHODS AND MATERIALS

Forward planning with a commercial system (RenderPlan 3D, Precision Therapy International, Inc., Norcross, GA) was used for IMAT planning. Arcs were approximated as multiple shaped fields spaced every 5-10 degrees around the patient. The number and ranges of the arcs were chosen manually. Multiple coplanar, superimposing arcs or noncoplanar arcs with or without a wedge were allowed. For comparison, conventional 3D conformal treatment plans were generated with the same commercial forward planning system as for IMAT. Intensity-modulated treatment plans were also created with a commercial inverse planning system (CORVUS, Nomos Corporation). A leaf-sequencing program was developed to generate the dynamic MLC prescriptions. IMAT treatment delivery was accomplished by programming the linear accelerator (linac) to deliver an arc and the MLC to step through a sequence of fields. Both gantry rotation and leaf motion were enslaved to the delivered MUs. Dosimetric accuracy of the entire process was verified with phantoms before IMAT was used clinically. For each IMAT treatment, a dry run was performed to assess the geometric and dosimetric accuracy. Both the central axis dose and dose distributions were measured and compared with predictions by the planning system.

RESULTS

By the end of May 2001, 50 patients had completed their treatments with the IMAT technique. Two to five arcs were needed to achieve highly conformal dose distributions. The IMAT plans provided better dose uniformity in the target and lower doses to normal structures than 3D conformal plans. The results varied when the comparison was made with fixed gantry IMRT. In general, IMAT plans provided more uniform dose distributions in the target, whereas the inverse-planned fixed gantry treatments had greater flexibility in controlling dose to the critical structures. Because the field sizes and shapes used in the IMAT were similar to those used in conventional treatments, the dosimetric uncertainty was very small. Of the first 32 patients treated, the average difference between the measured and predicted doses was -0.54 +/- 1.72% at isocenter. The 80%-95% isodose contours measured with film dosimetry matched those predicted by the planning system to within 2 mm. The planning time for IMAT was slightly longer than for generating conventional 3D conformal plans. However, because of the need to create phantom plans for the dry run, the overall planning time was doubled. The average time a patient spent on the table for IMAT treatment was similar to conventional treatments.

CONCLUSION

Initial results demonstrated the feasibility and accuracy of IMAT for achieving highly conformal dose distributions for different sites. If treatment plans can be optimized for IMAT cone beam delivery, we expect IMAT to achieve dose distributions that rival both slice-based and fixed-field IMRT techniques. The efficient delivery with existing linac and MLC makes IMAT a practical choice.

IMRT in the treatment of head and neck cancer: is the present already the future?

Disease outcome in locally advanced head and neck cancer patients is far from satisfactory. The main causes of failure remain linked to locoregional recurrences, which are due to incomplete eradication of clonogenic cells. Conventional radiation therapy or 3-dimensional conformal radiation therapy are currently carried out at their extreme possibilities due to their intrinsic limitation--namely the impossibility to generate concave dose distributions without compromising tumor irradiation. Approximately a third of patients treated with radiotherapy and most head and neck cancer cases present concave shapes of the target volumes. With the advent of intensity modulated radiation therapy--clinically available for only few years--head and neck patients can now benefit from strategies based on highly conformal techniques. It is possible to exploit efficiently dose-escalation protocols to increase probabilities to eradicate clonogens, to reduce overall treatment time, to control repopulation problems and to keep as low as reasonably necessary the irradiation of healthy tissues minimizing acute and late complications. Today, both planning and clinical studies demonstrate these advantages but larger controlled trials are necessary to assess the true potentialities of techniques based on intensity modulation for head and neck cancers. In a speculative view, proton therapy, possibly with intensity modulation, or light ion therapy should be considered for selected cases or for reirradiation due to their higher biological efficacy and their degree of dose-conformation to target volumes.

A performance comparison of flat-panel imager-based MV and kV cone-beam CT

The use of cone-beam computed tomography (CBCT) has been proposed for guiding the delivery of radiation therapy, and investigators have examined the use of both kilovoltage (kV) and megavoltage (MV) x-ray beams in the development of such CBCT systems. In this paper, the inherent contrast and signal-to-noise ratio (SNR) performance for a variety of existing and hypothetical detectors for CBCT are investigated analytically as a function of imaging dose and object size. Theoretical predictions are compared to the results of experimental investigations employing largearea flat-panel imagers (FPIs) at kV and MV energies. Measurements were performed on two different FPI-based CBCT systems: a bench-top prototype incorporating an FPI and kV x-ray source (100 kVp x rays), and a system incorporating an FPI mounted on the gantry of a medical linear accelerator (6 MV x rays). The SNR in volume reconstructions was measured as a function of dose and found to agree reasonably with theoretical predictions. These results confirm the theoretically predicted advantages of employing kV energy x rays in imaging soft-tissue structures found in the human body. While MV CBCT may provide a valuable means of correcting 3D setup errors and may offer an advantage in terms of simplicity of mechanical integration with a linear accelerator (e.g., implementation in place of a portal imager), kV CBCT offers significant performance advantages in terms of image contrast and SNR per unit dose for visualization of soft-tissue structures. The relatively poor SNR performance at MV energies is primarily a result of the low x-ray quantum efficiencies (approximately a few percent or less) that are currently achieved with FPIs at high energies. Furthermore, kV CBCT with an FPI offers the potential of combined volumetric and radiographic/fluoroscopic imaging using the same device.

A deterministic iterative least-squares algorithm for beam weight optimization in conformal radiotherapy

Currently, inverse treatment planning in conformal radiotherapy is, in part, a trial-and-error process due to the interplay of many competing criteria for obtaining a clinically acceptable dose distribution. A new method is developed for beam weight optimization that incorporates clinically relevant nonlinear and linear constraints. The process is driven by a nonlinear, quasi-quadratic objective function and the solution space is defined by a set of linear constraints. At each step of iteration, the optimization problem is linearized by a self-consistent approximation that is local to the existing dose distribution. The dose distribution is then improved by solving a series of constrained least-squares problems using an established method until all prescribed constraints are satisfied. This differs from the current approaches in that it does not rely on the search for the global minimum of a specific objective function. Essentially, our proposed objective function can be construed as a functional that comprises a class of dose-based quadratic objective functions. Empirical adjustment for appropriate model parameters in the construction of objective function is minimized, since these parameters are in effect adaptively adjusted during optimization. The method is robust in solving difficult clinical cases using either aperture or pencil beam based planning techniques for intensity-modulated radiation therapy.

Intensity-modulated arc therapy simplified

PURPOSE

We present a treatment planning strategy for intensity-modulated radiation therapy using gantry arcs with dynamic multileaf collimator, previously termed intensity-modulated arc therapy (IMAT).

METHODS AND MATERIALS

The planning strategy is an extension of the photon bar arc and asymmetric arc techniques and is classified into three levels of complexity, with increasing number of gantry arcs. This principle allows us to generalize the analysis of the number of arcs required for intensity modulation for a given treatment site. Using a phantom, we illustrate how the current technique is more flexible than the photon bar arc technique. We then compare plans from our strategy with conventional three-dimensional conformal treatment plans for three sites: prostate (prostate plus seminal vesicles), posterior pharyngeal wall, and chest wall.

RESULTS

Our strategy generates superior IMAT treatment plans compared to conventional three-dimensional conformal plans. The IMAT plans spare critical organs well, and the trade-off for simplicity is that the dose uniformity in the target volume may not rival that of true inverse treatment plans.

CONCLUSIONS

The analyses presented in this paper give a better understanding of IMAT plans. Our strategy is easier to understand and more efficient in generating plans than inverse planning systems; our plans are also simpler to modify, and quality assurance is more intuitive.

Tomographic motion detection and correction directly in sinogram space

Patient motion, especially respiratory motion, results in various artefacts such as blurring and streaks in tomographic images. The interplay of the movement of the beam aperture and variations of organ anatomy during delivery can create 'hot' and 'cold' spots throughout the field in intensity-modulated radiation therapy (IMRT). Detection and correction of patient motion is extremely important in tomographic imaging and IMRT. Tomographic projection data (sinogram) encode not only the patient anatomy information, but also the intra-scanning motion information. In this paper, we developed an algorithm to detect and correct the in-plane respiratory motion directly in sinogram space. The respiratory motion is modelled as time-varying scaling along the x and y directions. Its effects on the sinogram are discussed. Based on the traces of some nodal points in the sinogram, the intra-scanning motion is determined. The motion correction is also implemented in sinogram space. The motion-corrected sinogram is used for reconstruction by the filtered back-projection (FBP) method. Computer simulations validate the motion detection and correction algorithm. The reconstructed images from the motion-corrected sinogram eliminate the majority of the artefacts. The method could be applied to projection data used in CT and ECT, as well as in tomotherapy delivery modification and dose reconstruction.

Assessment of the acceptability of the Elekta multileaf collimator (MLC) within the Corvus planning system for static and dynamic delivery of intensity modulated beams (IMBs)

The sliding window technique used for static and dynamic segmentation of intensity modulated beams is evaluated. Dynamic delivery is preferred since the resulting distributions correspond better with the calculated distributions, the treatment beam is used more efficiently and the delivery is less sensitive to small variations in the accuracy of the multileaf collimator (MLC).

Quasi-independent monitor unit calculation for intensity modulated sequential tomotherapy

The number of linac monitor units (MU) from intensity modulated sequential tomotherapy (IMST) is substantially larger than the MU delivered in conventional radiation therapy, and the relation between MU and dose is obscure due to complicated variation of the beam intensities. The purpose of this work was to develop a practical method of verifying the MU and dose from IMST so that the MU of each arced beam could be double-checked for accuracy. MU calculations for 41 arced beams from 14 IMST patients were performed using the variables of vane open fraction time, field size, target depth, output factor, TMR, and derived intensity distribution. Discrepancy between planned and checked MU was quantified as 100 (MU(cal)-MU(plan))/MU(plan) percent. All 41 discrepancies were clustered between -5% to +4%, illustrated in a Gaussian-shaped histogram centered at -1.0+/-3.5% standard deviation indicating the present MU calculations are in agreement with the planned expectations. To confirm the correctness of the present calculated MUs of the IMST plans, eight of the calculated IMST plans are performed dose verifications using their hybrid plans, which are created by transporting patient's IMST plan beams onto a spherical polystyrene Phantom for dose distribution within the Phantom. The dose was measured with a 0.07 cc ionization chamber inserted in the spherical Phantom during the hybrid plan irradiation. Average discrepancy between planned and measured doses was found to be 0.6+/-3.4% with single standard deviation uncertainty. The spread of the discrepancies of present calculated MUs relative to their planned ones are attributed to uncertainties of effective field size, effective planned dose corresponding to each arc, and inaccuracy of quantification of scattered dose from adjacent arced beams. Overall, the present calculation of MUs is consistent with what derived from treatment plans. Since the MUs are verified by actual dose measurements, therefore the present MU calculation technique is considered adequate for double-checking planned IMST MUs.

Potential of tomotherapy for total scalp treatment

PURPOSE

Total scalp radiotherapy is required in a variety of clinical situations. We compared conventional lateral photon-electron (LPE) technique with tomotherapy (intensity-modulated radiotherapy [IMRT]).

METHODS AND MATERIALS

A patient with Merkel cell carcinoma was treated at our institution using conventional treatment techniques. Treatment plans were conducted using conventional three-dimensional treatment planning and IMRT. The clinical target volume included the entire scalp tissue volumes to the surface of underlying cranial bone, as well as superficial and deep neck nodes in the bilateral neck. To provide a consistent comparison between the IMRT and three-dimensional conventional treatment plans, the dose distributions were normalized such that 90% of the target volumes received the prescription dose.

RESULTS

Examination of our results revealed an acceptable dose-volume histogram and adequate coverage of the clinical target volume using the conventional LPE technique. The IMRT plan provided a more homogeneous dose to the target volume; however, critical structure doses were uniformly higher than for the conventional treatment plan.

CONCLUSIONS

The IMRT plan resulted in a substantial dose to the lens, brain, and orbit, making it clinically unacceptable compared with the LPE technique. Overall, the LPE technique was superior.

Monte Carlo study of a highly efficient gas ionization detector for megavoltage imaging and image-guided radiotherapy

The imaging characteristics of an arc-shaped xenon gas ionization chamber for the purpose of megavoltage CT imaging were investigated. The detector consists of several hundred 320 microm thick gas cavities separated by thin tungsten plates of the same thickness. Dose response, efficiency and resolution parameters were calculated using Monte Carlo simulations. The calculations were compared to measurements taken in a 4 MV photon beam, assuming that the measured signal in the chambers corresponds to the therein absorbed dose. The measured response profiles for narrow and broad incident photon beams could be well reproduced with the Monte Carlo calculations. They show, that the quantum efficiency is 29.2% and the detective quantum efficiency at zero frequency DQE(0) is 20.4% for the detector arc placed in focus with the photon source. For a detector placed out of focus, these numbers even increase. The efficiency of this kind of radiation detector for megavoltage radiation therefore surpasses the reported efficiency of existing detector technologies. The resolution of the detector is quantified with calculated and measured line spread functions. The corresponding modulation transfer functions were determined for different thicknesses of the tungsten plates. They show that the resolution is only slightly dependent on the plate thickness but is predominantly determined by the cell size of the detector. The optimal plate thickness is determined by a tradeoff between quantum efficiency, total signal generation and resolution. Thicker plates are more efficient but the total signal and the resolution decrease with plate thickness. In conclusion, a gas ionization chamber of the described type is a highly efficient megavoltage radiation detector, allowing to obtain CT images with very little dose for a sufficient image quality for anatomy verification. This kind of detector might serve as a model for a future generation of highly efficient radiation detectors.

An investigation of the reproducibility and usefulness of automatic couch motion in complex radiation therapy techniques

Radiation therapy techniques that incorporate multiple couch motions are becoming more common, and they often involve an increasing level of complexity along with a need for automatic motion. The reproducibility of automatic couch motion is thus a growing concern. In this work we carried out various tests to assess the automatic motion of a commercial treatment couch, including tests to evaluate the digital readout reproducibility, as well as an independent verification of the reproducibility of the couch positions on repeated motions, using phantoms as well as a volunteer subject. It was shown that the couch motion is highly reproducible, with no discomfort to the patient, and can greatly improve treatment times as well as reduce errors in couch positioning.

Analysis of the effects of the delivery technique on an IMRT plan: comparison for multiple static field, dynamic and NOMOS MIMiC collimation

The process of delivering an IMRT treatment may involve various beam-modifying techniques such as multileaf collimators (MLCs), the NOMOS MIMiC, blocks, wedges, etc. In the case of the MLC, the spatial/temporal variation of the position of the leaves and diaphragms in the beam allows the delivery of modulated beam profiles either by the multiple-static-field (MSF) method or by the dynamic multileaf collimator (DMLC) method. The constraints associated with the IMRT delivery technique are usually neglected in the process of obtaining the 'optimal' inverse treatment plan. Consequently, dose optimization may be significantly reduced when the 'optimal' beam profiles are converted to leaf/diaphragm positions via a leaf-sequencing interpreter. The paper presented here assesses the effects on the optimum treatment plan of the following leaf-sequencing algorithms: MSF, DMLC and NOMOS MIMiC. The results obtained suggest that the delivery of an 'optimum' plan produces an overdosage of the PTV region due to various factors such as leaf/diaphragm transmission effects, head-scatter and phantom-scatter contributions. The overdosage observed for a cohort of ten patients was 2.5, 3.7 and 5.7%, respectively, for the DMLC, MSF and NOMOS MIMiC, after normalizing the delivered fluence to account for IMRT effects (using the method of Convery et al (Convery D J, Cogrove V P and Webb S 2000 Proc. 13th Int. Conf. on Computers in Radiotherapy (Heidelberg, 2000)) such as to obtain 70 Gy at the isocentre. The IMRT techniques DMLC, MIMiC and MSF were compared for the organs at risk: rectum, bladder, and left and right femoral heads.

Intensity-modulated radiotherapy: current status and issues of interest

PURPOSE

To develop and disseminate a report aimed primarily at practicing radiation oncology physicians and medical physicists that describes the current state-of-the-art of intensity-modulated radiotherapy (IMRT). Those areas needing further research and development are identified by category and recommendations are given, which should also be of interest to IMRT equipment manufacturers and research funding agencies.

METHODS AND MATERIALS

The National Cancer Institute formed a Collaborative Working Group of experts in IMRT to develop consensus guidelines and recommendations for implementation of IMRT and for further research through a critical analysis of the published data supplemented by clinical experience. A glossary of the words and phrases currently used in IMRT is given in the. Recommendations for new terminology are given where clarification is needed.

RESULTS

IMRT, an advanced form of external beam irradiation, is a type of three-dimensional conformal radiotherapy (3D-CRT). It represents one of the most important technical advances in RT since the advent of the medical linear accelerator. 3D-CRT/IMRT is not just an add-on to the current radiation oncology process; it represents a radical change in practice, particularly for the radiation oncologist. For example, 3D-CRT/IMRT requires the use of 3D treatment planning capabilities, such as defining target volumes and organs at risk in three dimensions by drawing contours on cross-sectional images (i.e., CT, MRI) on a slice-by-slice basis as opposed to drawing beam portals on a simulator radiograph. In addition, IMRT requires that the physician clearly and quantitatively define the treatment objectives. Currently, most IMRT approaches will increase the time and effort required by physicians, medical physicists, dosimetrists, and radiation therapists, because IMRT planning and delivery systems are not yet robust enough to provide totally automated solutions for all disease sites. Considerable research is needed to model the clinical outcomes to allow truly automated solutions. Current IMRT delivery systems are essentially first-generation systems, and no single method stands out as the ultimate technique. The instrumentation and methods used for IMRT quality assurance procedures and testing are not yet well established. In addition, many fundamental questions regarding IMRT are still unanswered. For example, the radiobiologic consequences of altered time-dose fractionation are not completely understood. Also, because there may be a much greater ability to trade off dose heterogeneity in the target vs. avoidance of normal critical structures with IMRT compared with traditional RT techniques, conventional radiation oncology planning principles are challenged. All in all, this new process of planning and treatment delivery has significant potential for improving the therapeutic ratio and reducing toxicity. Also, although inefficient currently, it is expected that IMRT, when fully developed, will improve the overall efficiency with which external beam RT can be planned and delivered, and thus will potentially lower costs.

CONCLUSION

Recommendations in the areas pertinent to IMRT, including dose-calculation algorithms, acceptance testing, commissioning and quality assurance, facility planning and radiation safety, and target volume and dose specification, are presented. Several of the areas in which future research and development are needed are also indicated. These broad recommendations are intended to be both technical and advisory in nature, but the ultimate responsibility for clinical decisions pertaining to the implementation and use of IMRT rests with the radiation oncologist and radiation oncology physicist. This is an evolving field, and modifications of these recommendations are expected as new technology and data become available.

On the verification of the incident energy fluence in tomotherapy IMRT

For any radiotherapy verification technique, it is desirable that issues with the accelerator, multileaf collimator and patient position be detected. In previous works, an effective method for this level of verification was presented. This paper identifies second-order issues affecting the part of the process in which the incident energy fluence is verified. These problems will affect any rotational intensity-modulated radiotherapy delivery that divides each rotation or arc into projections: however the solutions offered in this paper are specific to the method previously developed. The issues affecting the energy fluence verification method include leaf bouncing. delivery implementation and leaf latency. All three matters were found to introduce small errors in the verified energy fluence values for a small fraction of leaf states. The overall effect on the deposited dose over the course of a rotational delivery involving thousands of beam pulses per rotation is negligible. Regardless, effective correction strategies are presented; these are utilized in order to characterize both the delivered energy fluence and deposited dose as accurately as possible.

Abutment dosimetry for serial tomotherapy

The dose distributions at the abutment region for serial tomotherapy are reviewed. While tomotherapy provides unparalleled dose distributions, precise couch motion and good patient immobilization are required because the dose in the abutment region changes by 25% for each millimeter of misalignment. The process of delivering intensity-modulated radiation therapy using sequentially delivered modulated arcs yields hot spots below and cold spots above the machine isocenter when arc angles of less than 360 degrees are used. The magnitude of the hot and cold spots increases significantly as the arc angle is reduced 180 degrees such as when limited by couch clearance restrictions. Placement of isocenter also significantly affects the dose heterogeneity in the abutment region, with the hot and cold spots increasing nearly linearly with off-axis distance in the vertical direction. Reduction of the magnitude of the abutment region dose heterogeneities is possible if helical delivery is provided by moving the couch during arc delivery. The dose heterogeneity can also be reduced by creating 2 treatment plans, each with slightly different abutment region positions, or by using multiple couch angles.

Comparative study between IMRT with NOMOS BEAK and linac-based radiosurgery in the treatment of intracranial lesions

A comparative study was undertaken to examine intracranial irradiation using intensity-modulation radiation therapy (IMRT) and linear accelerator-based radiosurgery. The IMRT was examined using the Peacock system with a BEAK attachment. A clinical case involving a metastatic brain lesion, treated with 3 radiosurgery isocenters, was planned for IMRT. The radiosurgery was planned using the Leibinger planning system. The IMRT was planned using the CORVUS planning system. The CORVUS planning system uses an inverse planning algorithm, a recent development in radiotherapy. Isodose distributions and dose volume histograms were generated and compared. Analysis of the dosimetry shows that the dose conformity and homogeneity within the target using the RTOG guidelines are superior for IMRT. The advantages of IMRT using inverse planning system include the ease of planning and execution of treatment, especially for cases that involve concave targets that require multiple isocenters using radiosurgery.

NOMOS Peacock IMRT utilizing the Beak post collimation device

Intensity-modulated radiation therapy (IMRT), an exciting recent development in the field of radiation therapy, is widely anticipated by many to make possible significant improvements in the quality of radiation treatments delivered to patients. The NOMOS Peacock method of delivery, often referred to as serial tomotherapy because of its "slice-wise" treatment of a tumor, has been used since 1994 to treat some 8000+ patients worldwide. This slice-wise method of treatment is known to produce extremely conformal dose distributions due to its ability to specifically match the dose distribution on each slice to the shape of the target volume on that same slice. Based on the belief of this institution, and the NOMOS Corporation, that an increase in the number of treatment slices into which the target is segmented would lead directly to an improvement in three-dimensional (3D) dose conformality, a joint effort was undertaken to develop a new MIMIC collimator treatment mode. Inherent to the original design of the NOMOS MIMIC binary multileaf collimator were 2 treatment modes: a 2-cm mode with a slice thickness of approximately 1.7 cm and a 1-cm mode with a slice thickness of approximately 0.85 cm. As a result of this collaborative effort, a new MIMIC treatment mode has been developed. The method employs a slit collimator, post-collimation device known as the BEAK, enabling the treatment mode referred to as Beak Mode. The device imposes a distal redefinition of the slice thickness, or length, by effectively blocking the full retraction of the MIMIC vanes. The end result is a newly available slice thickness of approximately 4 mm, which is shown in this work to yield significant improvements in dose conformality for 2 representative patients. The comparative analysis of these 2 patient plans includes, in addition to a comparison of isodose distributions, an evaluation of dose-volume histogram (DVH) information, and a comparison of indices of conformality (CI) and homogeneity (HI).

Where goest the Peacock?

Since the treatment of the first patient in 1994, the Peacock system has maintained its presence as the dominant method of intensity-modulated radiation therapy (IMRT) delivery. Currently in use at nearly 80 institutions, over 8000 patients have been treated using the system. Peacock treatments have been delivered to sites throughout the body, including CNS, head & neck, prostate, liver, kidney, lung, mediastinum, and extremities. IMRT, however, is a young and rapidly evolving treatment methodology. As institutions have explored new ways of improving radiation therapy with intensity-modulated techniques, the requirements for the Peacock system have also expanded. More sophisticated planning algorithms have been implemented to satisfy these new requirements, as well as better tools for treatment verification and quality assurance. In addition, new delivery techniques are being examined to improve the ability of IMRT to increase target volume doses while limiting organ-at-risk doses. One such technique, using helical tomotherapy (Peacock is an example of sequential tomotherapy), is currently being evaluated at one institution. Both techniques use narrow, modulated delivery beams. However, helical tomotherapy requires continuous movement of the couch during radiation, similar to helical CT. This work reviews the development of tomotherapy with the Peacock system. It then looks at current IMRT treatment techniques using tomotherapy, and how the field has broadened since the first treatments were delivered. Finally, it looks at the future of tomotherapy techniques, and how these techniques will adapt to the changing requirements for radiation therapy.

Non-coplanar inverse planning IMRT using the MIMiC system: clinical significance in choice of 2-cm/1-cm mode and single couch vs. multiple couch angles

The IMRT delivery process using MIMiC is extremely effective. In the hands of an experienced practitioner, MIMiC can provide excellent conformity to tumor volume with minimal damage to adjacent critical structures. However, the process is not completely automatic, even though it is an inverse planning process. The choice of machine treatment mode, 1- or 2-cm mode, and the choice of single couch angle vs. multiple, will yield quite different end results, with all prescription parameters being equal. The choice of mode and number of couch angles that yield excellent results of good coverage of the PTV with minimal hot spots and appropriate degree of normal tissue sparing still depends on the clinical experience of the planner. The size, shape, and total length of the PTV and those of the organs of interest dictate the choice. The authors demonstrated some clinical examples with comparative studies of the results of different choice of these parameters. In general, smaller beamlets with multiple couch angles give better PTV conformity, and critical organ avoidance. The down side of these choices are: (a) relatively more scattered hot spots within PTV, (b) longer time required to deliver the treatment, and (c) increased complexity of the process that could increase the chance of human error in treatment delivery.

The effect of source-axis distance on integral dose: implications for IMRT

The source-axis distance (SAD) is a treatment machine design parameter that affects integral dose, dose rate and patient clearance. The aim of this work was to investigate the effect of source-axis distance on integral dose for conformal arc therapy. This work is part of a larger project to determine the ideal characteristics of a dedicated IMRT machine. The sensitivity of SAD to beam energy, PTV size, body size and PTV position were determined for conformal arc therapy. For the calculations performed here it was assumed that dose equals terma. The integral dose ratio (IDR) was used to quantify the calculation results. It was found that the IDR increases as both SAD and photon energy increase, though the dependence of IDR on SAD decreases as energy increases. The PTV size was found to have a negligible effect on the relationship between the SAD and IDR, however the body size does affect the relationship between the SAD and IDR. The position of the PTV within the body also affects the IDR. From dosimetric considerations alone, the larger the SAD, the better the possible dose distribution. The IDR for a very large SAD is increased by approximately 5% when compared with the IDR for 100 cm SAD. Similarly, the IDR for 100 cm SAD is approximately 5% higher than the IDR at 50 cm SAD.

Monte Carlo dose verification for intensity-modulated arc therapy

Intensity-modulated arc therapy (IMAT), a technique which combines beam rotation and dynamic multileaf collimation, has been implemented in our clinic. Dosimetric errors can be created by the inability of the planning system to accurately account for the effects of tissue inhomogeneities and physical characteristics of the multileaf collimator (MLC). The objective of this study is to explore the use of Monte Carlo (MC) simulation for IMAT dose verification. The BEAM/DOSXYZ Monte Carlo system was implemented to perform dose verification for the IMAT treatment. The implementation includes the simulation of the linac head/MLC (Elekta SL20), the conversion of patient CT images and beam arrangement for 3D dose calculation, the calculation of gantry rotation and leaf motion by a series of static beams and the development of software to automate the entire MC process. The MC calculations were verified by measurements for conventional beam settings. The agreement was within 2%. The IMAT dose distributions generated by a commercial forward planning system (RenderPlan. Elekta) were compared with those calculated by the MC package. For the cases studied, discrepancies of over 10% were found between the MC and the RenderPlan dose calculations. These discrepancies were due in part to the inaccurate dose calculation of the RenderPlan system. The computation time for the IMAT MC calculation was in the range of 20-80 min on 15 Pentium-Ill computers. The MC method was also useful in verifying the beam apertures used in the IMAT treatments.

On few-view tomographic reconstruction with megavoltage photon beams

Currently portal imaging devices are used to obtain information on patient localization during radiation therapy treatments. Such obtained information is two dimensional in nature, limited to the plane of the captured image. It has been proposed that megavoltage computed tomography images be reconstructed to overcome this limitation. This study explores the feasibility of reconstructing tomographic images from fan-beam projection data acquired with a commercial portal imaging device on a standard radiotherapy linear accelerator. Several CT reconstruction algorithms are examined as to their performance and suitability for applications in radiation therapy verification. The results show that it is possible, using some of the iterative reconstruction techniques, to obtain an image useful for patient localization from only several (< or =10) projection views.

An intercomparison of IMRT delivery techniques: a case study for breast treatment

Intensity-modulated radiotherapy beams can be delivered using a multileaf collimator by one of two methods: either by superposition of a series of multiple-static fields, or by moving the collimators while the beam is on to produce 'dynamically' modulated beams. The leaf trajectories in this dynamic mode are given by a series of linear steps between control points defining each collimator position at known intervals throughout an exposure. The complexity of the resulting modulation is limited in the first case by the number of fields superposed and in the second case by the number of control points defined. Results are presented for an experimental study that investigates the effect of changing both the number of fields for the multiple-static technique, and the number of control points for a dynamic 'close-in' technique. All deliveries studied are clinical intensity-modulated breast fields. The effect of using a universal wedge in conjunction with the multileaf collimator is also studied, together with a comparison of the relative efficiency, time taken and the absolute dosimetric accuracy of the various delivery options. It is shown that all delivery techniques produce equivalent dose distributions when using 15 control points, with 10 control points being sufficient to produce an adequate breast compensator distribution. Except for the case of a four-control-point dynamic delivery, the universal wedge makes no significant difference to the dose distribution. However, it makes the delivery less efficient. The close-in interpreter consistently produces deliveries that are more efficient than the more conventional sliding-window technique and faster than the multiple-static-field technique. Finally the close-in technique is compared to the more 'standard' leaf-sweep technique and shown to be equivalent.

Directional dependence in film dosimetry: radiographic and radiochromic film

The trend towards conformal, intensity modulated radiotherapy treatments has established the need for a true integrating dosimeter. In traditional radiotherapy, radiographic film dosimetry is commonly used. The accuracy and reproducibility of film optical density as an indicator of dose is influenced by several variables, including the chemical processing conditions. As a result radiochromic film, with all the advantages of radiographic film but without the need for chemical processing, has increased in popularity, although the low-dose sensitivity of radiochromic film does remain a disadvantage for some experiments. Several studies have investigated the reproducibility of radiochromic film results, but none have specifically addressed the well-known directional dependence seen with traditional radiographic film. In this study, the directional dependence of radiographic (Kodak X-omat V) and radiochromic (Gafchromic) films were measured. It was found that both films over responded when exposed parallel to the central axis of the beam as opposed to perpendicular exposure. An attempt is made to explain the reason for the responses of both films in terms of spectral effects and the air gap between the phantom segments. Although radiographic film exposed parallel rather than perpendicular to the central axis of the beam exhibits a measured difference in film response at depth, this over response does not occur when the extent of the film is restricted to a small region at the centre of the phantom (in this case an air gap is not introduced across the phantom). This suggests that it is the air gap rather than the orientation of the film that is the cause of the over response. Furthermore, when film occupies a slice through the entire phantom an over response occurs for both radiographic and radiochromic film, indicating that spectral effects are not the cause.

A feasible method for clinical delivery verification and dose reconstruction in tomotherapy

Delivery verification is the process in which the energy fluence delivered during a treatment is verified. This verified energy fluence can be used in conjunction with an image in the treatment position to reconstruct the full three-dimensional dose deposited. A method for delivery verification that utilizes a measured database of detector signal is described in this work. This database is a function of two parameters, radiological path-length and detector-to-phantom distance, both of which are computed from a CT image taken at the time of delivery. Such a database was generated and used to perform delivery verification and dose reconstruction. Two experiments were conducted: a simulated prostate delivery on an inhomogeneous abdominal phantom, and a nasopharyngeal delivery on a dog cadaver. For both cases, it was found that the verified fluence and dose results using the database approach agreed very well with those using previously developed and proven techniques. Delivery verification with a measured database and CT image at the time of treatment is an accurate procedure for tomotherapy. The database eliminates the need for any patient-specific, pre- or post-treatment measurements. Moreover, such an approach creates an opportunity for accurate, real-time delivery verification and dose reconstruction given fast image reconstruction and dose computation tools.

Organ/patient geometric variation in external beam radiotherapy and its effects

Treatment variation in positioning of the organ/patient with respect to the radiation beams causes a temporal dose variation in critical normal tissues adjacent to the treatment target. This temporal variation induces uncertainties in understanding the normal tissue dose response, thereby limiting reliable treatment evaluation and optimization. The aim of this study is to model and analyze the temporal variation of organ dose distribution, and its effect on the biological effective dose. The study mainly focuses on the temporal dose variation caused by intertreatment organ motion/ deformation and daily setup error. Sensitivity of the biological effective dose to organ/patient geometric variation, dose distribution, and treatment fractionation will be investigated. Significant deviation of the biological effective dose could be expected in a critical normal structure, even if the cumulative dose deviation in this structure is negligible. Patients with similar geometric variation characteristics can experience significantly different biological effective dose, and the differences are sensitive to the dose distribution and the total number of treatment fractions.

On the accuracy and effectiveness of dose reconstruction for tomotherapy

Dose reconstruction is a process that re-creates the treatment-time dose deposited in a patient provided there is knowledge of the delivered energy fluence and the patient's anatomy at the time of treatment. A method for reconstructing dose is presented. The process starts with delivery verification, in which the incident energy fluence from a treatment is computed using the exit detector signal and a transfer matrix to convert the detector signal to energy fluence. With the verified energy fluence and a CT image of the patient in the treatment position, the treatment-time dose distribution is computed using any model-based algorithm such as convolution/superposition or Monte Carlo. The accuracy of dose reconstruction and the ability of the process to reveal delivery errors are presented. Regarding accuracy, a reconstructed dose distribution was compared with a measured film distribution for a simulated breast treatment carried out on a thorax phantom. It was found that the reconstructed dose distribution agreed well with the dose distribution measured using film: the majority of the voxels were within the low and high dose-gradient tolerances of 3% and 3 mm respectively. Concerning delivery errors, it was found that errors associated with the accelerator, the multileaf collimator and patient positioning might be detected in the verified energy fluence and are readily apparent in the reconstructed dose. For the cases in which errors appear in the reconstructed dose, the possibility for adaptive radiotherapy is discussed.

Intensity modulated conformal therapy for intracranial lesions

The Peacock planning and delivery system was used to create treatment plans and deliver these plans to patients. The system involves an arc therapy delivery of small (2 cm long) slices of radiation combined with indexing of the couch to achieve target coverage. Two clinical examples are shown to demonstrate the system's capability and evaluate the resources required to produce and deliver the plans. One plan is an optic sheath meningioma and the other is a craniopharyngioma that surrounded the optic chiasm. The optic sheath meningioma was treated to 50 Gy in 25 fractions. The treatment involved delivery of two arcs. The total time to set up the patient and deliver the treatment was less than 15 min. Planning and plan validation after computed tomography required approximately 3 days. The patient had 100% restoration of her field of vision and is stable 3 years post therapy. The second patient is a 9-year-old who had a craniopharyngioma which surrounded the optic chiasm. The tumor was treated to 50.4 Gy in 28 fractions and the dose to the optic chiasm was limited to 45 Gy. The treatment required three arcs and total treatment time was less than 20 min. The patient is stable 15 months post therapy. The system is able to create and deliver radiation patterns that are unique. These plans can be created and delivered in times that rival conventional forward planning conformal radiotherapy systems that cannot produce or conveniently deliver such plans.

A method to incorporate leakage and head scatter corrections into a tomotherapy inverse treatment planning algorithm

A detailed tomotherapy inverse treatment planning method is described which incorporates leakage and head scatter corrections during each iteration of the optimization process, allowing these effects to be directly accounted for in the optimized dose distribution. It is shown that the conventional inverse planning method for optimizing incident intensity can be extended to include a 'concurrent' leaf sequencing operation from which the leakage and head scatter corrections are determined. The method is demonstrated using the steepest-descent optimization technique with constant step size and a least-squared error objective. The method was implemented using the MATLAB scientific programming environment and its feasibility demonstrated for 2D test cases simulating treatment delivery using a single coplanar rotation. The results indicate that this modification does not significantly affect convergence of the intensity optimization method when exposure times of individual leaves are stratified to a large number of levels (>100) during leaf sequencing. In general, the addition of aperture dependent corrections, especially 'head scatter', reduces incident fluence in local regions of the modulated fan beam, resulting in increased exposure times for individual collimator leaves. These local variations can result in 5% or greater local variation in the optimized dose distribution compared to the uncorrected case. The overall efficiency of the modified intensity optimization algorithm is comparable to that of the original unmodified case.

Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware

Algebraic reconstruction methods, such as the algebraic reconstruction technique (ART) and the related simultaneous ART (SART). reconstruct a two-dimensional (2-D) or three-dimensional (3-D) object from its X-ray projections. The algebraic methods have, in certain scenarios, many advantages over the more popular Filtered Backprojection approaches and have also recently been shown to perform well for 3-D cone-beam reconstruction. However, so far the slow speed of these iterative methods have prohibited their routine use in clinical applications. In this paper, we address this shortcoming and investigate the utility of widely available 2-D texture mapping graphics hardware for the purpose of accelerating the 3-D algebraic reconstruction. We find that this hardware allows 3-D cone-beam reconstructions to be obtained at almost interactive speeds, with speed-ups of over 50 with respect to implementations that only use general-purpose CPUs. However, we also find that the reconstruction quality is rather sensitive to the resolution of the framebuffer, and to address this critical issue we propose a scheme that extends the precision of a given framebuffer by 4 bits, using the color channels. With this extension, a 12-bit framebuffer delivers useful reconstructions for 0.5% tissue contrast, while an 8-bit framebuffer requires 4%. Since graphics hardware generates an entire image for each volume projection, it is most appropriately used with an algebraic reconstruction method that performs volume correction at that granularity as well, such as SART or SIRT. We chose SART for its faster convergence properties.

Monitor unit settings for intensity modulated beams delivered using a step-and-shoot approach

Two linear accelerators have been commissioned for delivering IMRT treatments using a step-and-shoot approach. To assess beam startup stability for 6 and 18 MV x-ray beams, dose delivered per monitor unit (MU), beam flatness, and beam symmetry were measured as a function of the total number of MU delivered at a clinical dose rate of 400 MU per minute. Relative to a 100 MU exposure, the dose delivered per MU by both linear accelerators was found to be within +/-2% for exposures larger than 4 MU. Beam flatness and symmetry also met accepted quality assurance standards for a minimum exposure of 4 MU. We have found that the performance of the two machines under study is well suited to the delivery of step-and-shoot IMRT. A system of dose calculation has also been commissioned for applying head scatter corrections to fields as small as 1x1 cm2. The accuracy and precision of the relative output calculations in water was validated for small fields and fields offset from the axis of collimator rotation. For both 6 and 18 MV x-ray beams, the dose per MU calculated in a water phantom agrees with measured data to within 1% on average, with a maximum deviation of 2.5%. The largest output factor discrepancies were seen when the actual radiation field size deviated from the set field size. The measured output in water can vary by as much 16% for 1x1 cm2 fields, when the measured field size deviates from the set field size by 2 mm. For a 1 mm deviation, this discrepancy was reduced to 8%. Steps should be taken to ensure collimator precision is tightly controlled when using such small fields. If this is not possible, very small fields should not contribute to a significant portion of the treatment, or uncertainties in the collimator position may effect the accuracy of the dose delivered.