Sasa Mutic

A technique for the quantitative evaluation of dose distributions

The commissioning of a three-dimensional treatment planning system requires comparisons of measured and calculated dose distributions. Techniques have been developed to facilitate quantitative comparisons, including superimposed isodoses, dose-difference, and distance-to-agreement (DTA) distributions. The criterion for acceptable calculation performance is generally defined as a tolerance of the dose and DTA in regions of low and high dose gradients, respectively. The dose difference and DTA distributions complement each other in their useful regions. A composite distribution has recently been developed that presents the dose difference in regions that fail both dose-difference and DTA comparison criteria. Although the composite distribution identifies locations where the calculation fails the preselected criteria, no numerical quality measure is provided for display or analysis. A technique is developed to unify dose distribution comparisons using the acceptance criteria. The measure of acceptability is the multidimensional distance between the measurement and calculation points in both the dose and the physical distance, scaled as a fraction of the acceptance criteria. In a space composed of dose and spatial coordinates, the acceptance criteria form an ellipsoid surface, the major axis scales of which are determined by individual acceptance criteria and the center of which is located at the measurement point in question. When the calculated dose distribution surface passes through the ellipsoid, the calculation passes the acceptance test for the measurement point. The minimum radial distance between the measurement point and the calculation points (expressed as a surface in the dose-distance space) is termed the gamma index. Regions where gamma > 1 correspond to locations where the calculation does not meet the acceptance criteria. The determination of gamma throughout the measured dose distribution provides a presentation that quantitatively indicates the calculation accuracy. Examples of a 6 MV beam penumbra are used to illustrate the gamma index.

A method for the reconstruction of four-dimensional synchronized CT scans acquired during free breathing

Breathing motion is a significant source of error in radiotherapy treatment planning for the thorax and upper abdomen. Accounting for breathing motion has a profound effect on the size of conformal radiation portals employed in these sites. Breathing motion also causes artifacts and distortions in treatment planning computed tomography (CT) scans acquired during free breathing and also causes a breakdown of the assumption of the superposition of radiation portals in intensity-modulated radiation therapy, possibly leading to significant dose delivery errors. Proposed voluntary and involuntary breath-hold techniques have the potential for reducing or eliminating the effects of breathing motion, however, they are limited in practice, by the fact that many lung cancer patients cannot tolerate holding their breath. We present an alternative solution to accounting for breathing motion in radiotherapy treatment planning, where multislice CT scans are collected simultaneously with digital spirometry over many free breathing cycles to create a four-dimensional (4-D) image set, where tidal lung volume is the additional dimension. An analysis of this 4-D data leads to methods for digital-spirometry, based elimination or accounting of breathing motion artifacts in radiotherapy treatment planning for free breathing patients. The 4-D image set is generated by sorting free-breathing multislice CT scans according to user-defined tidal-volume bins. A multislice CT scanner is operated in the ciné mode, acquiring 15 scans per couch position, while the patient undergoes simultaneous digital-spirometry measurements. The spirometry is used to retrospectively sort the CT scans by their correlated tidal lung volume within the patients normal breathing cycle. This method has been prototyped using data from three lung cancer patients. The actual tidal lung volumes agreed with the specified bin volumes within standard deviations ranging between 22 and 33 cm3. An analysis of sagittal and coronal images demonstrated relatively small (<1 cm) motion artifacts along the diaphragm, even for tidal volumes where the rate of breathing motion is greatest. While still under development, this technology has the potential for revolutionizing the radiotherapy treatment planning for the thorax and upper abdomen.

The ViewRay System: Magnetic Resonance- Guided and Controlled Radiotherapy

A description of the first commercially available magnetic resonance imaging (MRI)-guided radiation therapy (RT) system is provided. The system consists of a split 0.35-T MR scanner straddling 3 (60)Co heads mounted on a ring gantry, each head equipped with independent doubly focused multileaf collimators. The MR and RT systems share a common isocenter, enabling simultaneous and continuous MRI during RT delivery. An on-couch adaptive RT treatment-planning system and integrated MRI-guided RT control system allow for rapid adaptive planning and beam delivery control based on the visualization of soft tissues. Treatment of patients with this system commenced at Washington University in January 2014.

Quality assurance for computed-tomography simulators and the computed-tomography-simulation process : report of the AAPM Radiation Therapy Committee Task Group No. 66

This document presents recommendations of the American Association of Physicists in Medicine (AAPM) for quality assurance of computed-tomography- (CT) simulators and CT-simulation process. This report was prepared by Task Group No. 66 of the AAPM Radiation Therapy Committee. It was approved by the Radiation Therapy Committee and by the AAPM Science Council.

Clinical Outcomes of Definitive Intensity-Modulated Radiation Therapy With Fluorodeoxyglucose-Positron Emission Tomography Simulation in Patients With Locally Advanced Cervical Cancer

PURPOSE This study aimed to evaluate the toxicity and clinical outcomes for cervical cancer patients treated definitively with intensity-modulated radiation therapy (IMRT) compared with non-IMRT treatment. METHODS AND MATERIALS This prospective cohort study included 452 patients with newly diagnosed cervical cancer treated with curative intent (135 IMRT and 317 non-IMRT). Treatment involved external irradiation and brachytherapy, and 85% of patients received concurrent chemotherapy. All IMRT patients underwent an F-18 fluorodeoxyglucose positron emission tomography (FDG-PET/CT) simulation. A 3-month post-therapy PET was obtained to evaluate treatment response. Toxicity was scored by the Common Terminology Criteria for Adverse Events Version 3.0. RESULTS The IMRT and non-IMRT groups had similar stage distribution and histology. For all patients, the post-therapy FDG-PET response correlated with overall recurrence risk (p < 0.0001) and cause-specific survival (p < 0.0001). Post-treatment FDG-PET findings were not significantly different between the IMRT and non-IMRT patients (p = 0.9774). The mean follow-up for all patients alive at the time of last follow-up was 52 months (72 months non-IMRT, 22 months IMRT). At last follow-up, 178 patients (39 IMRT, 139 non-IMRT) had developed a recurrence. The difference in recurrence-free survival between the two groups did not reach statistical significance (p = 0.0738), although the IMRT group showed better overall and cause-specific survivals (p < 0.0001). Of the patients, 62 patients (8 IMRT and 54 non-IMRT) developed Grade 3 or greater bowel or bladder complications, and by cumulative hazard function analysis the risk was significantly less for patients treated with IMRT (p = 0.0351). CONCLUSION Cervical cancer patients treated with FDG-PET/CT-guided IMRT have improved survival and less treatment-related toxicity compared with patients treated with non-IMRT radiotherapy.

EXPERIENCE-BASED QUALITY CONTROL OF CLINICAL INTENSITY-MODULATED RADIOTHERAPY PLANNING

PURPOSE To incorporate a quality control tool, according to previous planning experience and patient-specific anatomic information, into the intensity-modulated radiotherapy (IMRT) plan generation process and to determine whether the tool improved treatment plan quality. METHODS AND MATERIALS A retrospective study of 42 IMRT plans demonstrated a correlation between the fraction of organs at risk (OARs) overlapping the planning target volume and the mean dose. This yielded a model, predicted dose = prescription dose (0.2 + 0.8 [1 - exp(-3 overlapping planning target volume/volume of OAR)]), that predicted the achievable mean doses according to the planning target volume overlap/volume of OAR and the prescription dose. The model was incorporated into the planning process by way of a user-executable script that reported the predicted dose for any OAR. The script was introduced to clinicians engaged in IMRT planning and deployed thereafter. The scripts effect was evaluated by tracking δ = (mean dose-predicted dose)/predicted dose, the fraction by which the mean dose exceeded the model. RESULTS All OARs under investigation (rectum and bladder in prostate cancer; parotid glands, esophagus, and larynx in head-and-neck cancer) exhibited both smaller δ and reduced variability after script implementation. These effects were substantial for the parotid glands, for which the previous δ = 0.28 ± 0.24 was reduced to δ = 0.13 ± 0.10. The clinical relevance was most evident in the subset of cases in which the parotid glands were potentially salvageable (predicted dose <30 Gy). Before script implementation, an average of 30.1 Gy was delivered to the salvageable cases, with an average predicted dose of 20.3 Gy. After implementation, an average of 18.7 Gy was delivered to salvageable cases, with an average predicted dose of 17.2 Gy. In the prostate cases, the rectum model excess was reduced from δ = 0.28 ± 0.20 to δ = 0.07 ± 0.15. On surveying dosimetrists at the end of the study, most reported that the script both improved their IMRT planning (8 of 10) and increased their efficiency (6 of 10). CONCLUSIONS This tool proved successful in increasing normal tissue sparing and reducing interclinician variability, providing effective quality control of the IMRT plan development process.

Evaluation of polymer gels and MRI as a 3-D dosimeter for intensity-modulated radiation therapy.

BANG gel (MGS Research, Inc., Guilford, CT) has been evaluated for measuring intensity-modulated radiation therapy (IMRT) dose distributions. Treatment plans with target doses of 1500 cGy were generated by the Peacock IMRT system (NOMOS Corp., Sewickley, PA) using test target volumes. The gels were enclosed in 13 cm outer diameter cylindrical glass vessels. Dose calibration was conducted using seven smaller (4 cm diameter) cylindrical glass vessels irradiated to 0-1800 cGy in 300 cGy increments. Three-dimensional maps of the proton relaxation rate R2 were obtained using a 1.5 T magnetic resonance imaging (MRI) system (Siemens Medical Systems, Erlangen, Germany) and correlated with dose. A Hahn spin echo sequence was used with TR = 3 s, TE = 20 and 100 ms, NEX = 1, using 1 x 1 x 3 mm3 voxels. The MRI measurements were repeated weekly to identify the gel-aging characteristics. Ionization chamber, thermoluminescent dosimetry (TLD), and film dosimetry measurements of the IMRT dose distributions were obtained to compare against the gel results. The other dosimeters were used in a phantom with the same external cross-section as the gel phantom. The irradiated R2 values of the large vessels did not precisely track the smaller vessels, so the ionization chamber measurements were used to normalize the gel dose distributions. The point-to-point standard deviation of the gel dose measurements was 7.0 cGy. When compared with the ionization chamber measurements averaged over the chamber volume, 1% agreement was obtained. Comparisons against radiographic film dose distribution measurements and the treatment planning dose distribution calculation were used to determine the spatial localization accuracy of the gel and MRI. Spatial localization was better than 2 mm, and the dose was accurately determined by the gel both within and outside the target. The TLD chips were placed throughout the phantom to determine gel measurement precision in high- and low-dose regions. A multidimensional dose comparison tool that simultaneously examines the dose-difference and distance-to-agreement was used to evaluate the gel in both low-and high-dose gradient regions. When 3% and 3 mm criteria were used for the comparisons, more than 90% of the TLD measurements agreed with the gel, with the worst of 309 TLD chip measurements disagreeing by 40% of the criteria. All four MRI measurement session gel-measured dose distributions were compared to evaluate the time behavior of the gel. The low-dose regions were evaluated by comparison with TLD measurements at selected points, while high-dose regions were evaluated by directly comparing measured dose distributions. Tests using the multidimensional comparison tool showed detectable degradation beyond one week postirradiation, but all low-dose measurements passed relative to the test criteria and the dose distributions showed few regions that failed.

Quantitative dosimetric verification of an IMRT planning and delivery system

BACKGROUND AND PURPOSE The accuracy of dose calculation and delivery of a commercial serial tomotherapy treatment planning and delivery system (Peacock. NOMOS Corporation) was experimentally determined. MATERIALS AND METHODS External beam fluence distributions were optimized and delivered to test treatment plan target volumes, including three with cylindrical targets with diameters ranging from 2.0 to 6.2 cm and lengths of 0.9 through 4.8 cm, one using three cylindrical targets and two using C-shaped targets surrounding a critical structure, each with different dose distribution optimization criteria. Computer overlays of film-measured and calculated planar dose distributions were used to assess the dose calculation and delivery spatial accuracy. A 0.125 cm3 ionization chamber was used to conduct absolute point dosimetry verification. Thermoluminescent dosimetry chips, a small-volume ionization chamber and radiochromic film were used as independent checks of the ion chamber measurements. RESULTS Spatial localization accuracy was found to be better than +/-2.0 mm in the transverse axes (with one exception of 3.0 mm) and +/-1.5 mm in the longitudinal axis. Dosimetric verification using single slice delivery versions of the plans showed that the relative dose distribution was accurate to +/-2% within and outside the target volumes (in high dose and low dose gradient regions) with a mean and standard deviation for all points of -0.05% and 1.1%, respectively. The absolute dose per monitor unit was found to vary by +/-3.5% of the mean value due to the lack of consideration for leakage radiation and the limited scattered radiation integration in the dose calculation algorithm. To deliver the prescribed dose, adjustment of the monitor units by the measured ratio would be required. CONCLUSIONS The treatment planning and delivery system offered suitably accurate spatial registration and dose delivery of serial tomotherapy generated dose distributions. The quantitative dose comparisons were made as far as possible from abutment regions and examination of the dosimetry of these regions will also be important. Because of the variability in the dose per monitor unit and the complex nature of the calculation and delivery of serial tomotherapy, patient-specific quality assurance procedures will include a measurement of the delivered target dose.

Quantitation of the reconstruction quality of a four-dimensional computed tomography process for lung cancer patients

We have developed a four-dimensional computed tomography (4D CT) technique for mapping breathing motion in radiotherapy treatment planning. A multislice CT scanner (1.5 mm slices) operated in ciné mode was used to acquire 12 contiguous slices in each couch position for 15 consecutive scans (0.5 s rotation, 0.25 s between scans) while the patient underwent simultaneous quantitative spirometry measurements to provide a sorting metric. The spirometry-sorted scans were used to reconstruct a 4D data set. A critical factor for 4D CT is quantifying the reconstructed data set quality which we measure by correlating the metric used relative to internal-object motion. For this study, the internal air content within the lung was used as a surrogate for internal motion measurements. Thresholding and image morphological operations were applied to delineate the air-containing tissues (lungs, trachea) from each CT slice. The Hounsfield values were converted to the internal air content (V). The relationship between the air content and spirometer-measured tidal volume (v) was found to be quite linear throughout the lungs and was used to estimate the overall accuracy and precision of tidal volume-sorted 4D CT. Inspection of the CT-scan air content as a function of tidal volume showed excellent correlations (typically r>0.99) throughout the lung volume. Because of the discovered linear relationship, the ratio of internal air content to tidal volume was indicative of the fraction of air change in each couch position. Theoretically, due to air density differences within the lung and in room, the sum of these ratios would equal 1.11. For 12 patients, the mean value was 1.08 +/- 0.06, indicating the high quality of spirometry-based image sorting. The residual of a first-order fit between v and V was used to estimate the process precision. For all patients, the precision was better than 8%, with a mean value of 5.1% +/- 1.9%. This quantitative analysis highlights the value of using spirometry as the metric in sorting CT scans. The 4D reconstruction provides the CT data required to measure the three-dimensional trajectory of tumor and lung tissue during free breathing.

PET-GUIDED IMRT FOR CERVICAL CARCINOMA WITH POSITIVE PARA-AORTIC LYMPH NODES—A DOSE-ESCALATION TREATMENT PLANNING STUDY

PURPOSE To evaluate a treatment planning method for dose escalation to the para-aortic lymph nodes (PALNs) based on positron emission tomography (PET) with intensity-modulated radiotherapy (IMRT) for cervical cancer patients with PALN involvement. One goal of this process was not to modify the traditional treatment of the pelvic region. METHODS AND MATERIALS PET images for 4 cervical cancer patents with PALN involvement were registered with their corresponding CT scans. Positive PALNs were identified on PET images, and the surrounding critical structures were delineated on CT images. The treatment machine central axis (CAX) was placed at the level of the L4-L5 vertebral body interspace. There were two distinct treatment regions: the para-aortic bed superior to the CAX and the whole pelvis region inferior to the CAX. IMRT was used for treatment planning of PALN bed irradiation. The positive PALNs identified on PET images were defined as the gross target volume, and the para-aortic bed was defined as the clinical target volume. The radiation doses were escalated from the conventional 45 Gy to 59.4 Gy for the gross target volume and 50.4 Gy for the clinical target volume in 33 fractions. The pelvis area was treated with conventional treatment methods, AP-PA beams to 50.4 Gy in 28 fractions with a brachytherapy implant boost. The placement of the CAX allowed the two treatment regions to be abutted using the treatment machines independent jaws. RESULTS Dose escalation to positive PALNs, as identified on PET images, and the PALN bed is feasible with IMRT. Treatment plans for 4 patients revealed that escalated prescription doses could be delivered to target volumes while maintaining acceptable doses to the surrounding critical structures. Strategic placement of the treatment isocenter allows the IMRT region (PALN bed) and whole pelvis fields to be treated with a relatively uniform dose distribution in the abutment region. CONCLUSION This study indicates that PET-guided IMRT could be used in a clinical trial in an attempt to escalate doses delivered to patients with cervical cancer who have positive PALNs.