A study of respiration-correlated cone-beam CT scans to correct target positioning errors in radiotherapy of thoracic cancer.

PURPOSE: There is increasingly widespread usage of cone-beam CT (CBCT) for guiding radiation treatment in advanced-stage lung tumors, but difficulties associated with daily CBCT in conventionally fractionated treatments include imaging dose to the patient, increased workload and longer treatment times. Respiration-correlated cone-beam CT (RC-CBCT) can improve localization accuracy in mobile lung tumors, but further increases the time and workload for conventionally fractionated treatments. This study investigates whether RC-CBCT-guided correction of systematic tumor deviations in standard fractionated lung tumor radiation treatments is more effective than 2D image-based correction of skeletal deviations alone. A second study goal compares respiration-correlated vs respiration-averaged images for determining tumor deviations. METHODS: Eleven stage II-IV nonsmall cell lung cancer patients are enrolled in an IRB-approved prospective off-line protocol using RC-CBCT guidance to correct for systematic errors in GTV position. Patients receive a respiration-correlated planning CT (RCCT) at simulation, daily kilovoltage RC-CBCT scans during the first week of treatment and weekly scans thereafter. Four types of correction methods are compared: (1) systematic error in gross tumor volume (GTV) position, (2) systematic error in skeletal anatomy, (3) daily skeletal corrections, and (4) weekly skeletal corrections. The comparison is in terms of weighted average of the residual GTV deviations measured from the RC-CBCT scans and representing the estimated residual deviation over the treatment course. In the second study goal, GTV deviations computed from matching RCCT and RC-CBCT are compared to deviations computed from matching respiration-averaged images consisting of a CBCT reconstructed using all projections and an average-intensity-projection CT computed from the RCCT. RESULTS: Of the eleven patients in the GTV-based systematic correction protocol, two required no correction, seven required a single correction, one required two corrections, and one required three corrections. Mean residual GTV deviation (3D distance) following GTV-based systematic correction (mean ± 1 standard deviation 4.8 ± 1.5 mm) is significantly lower than for systematic skeletal-based (6.5 ± 2.9 mm, p = 0.015), and weekly skeletal-based correction (7.2 ± 3.0 mm, p = 0.001), but is not significantly lower than daily skeletal-based correction (5.4 ± 2.6 mm, p = 0.34). In two cases, first-day CBCT images reveal tumor changes-one showing tumor growth, the other showing large tumor displacement-that are not readily observed in radiographs. Differences in computed GTV deviations between respiration-correlated and respiration-averaged images are 0.2 ± 1.8 mm in the superior-inferior direction and are of similar magnitude in the other directions. CONCLUSIONS: An off-line protocol to correct GTV-based systematic error in locally advanced lung tumor cases can be effective at reducing tumor deviations, although the findings need confirmation with larger patient statistics. In some cases, a single cone-beam CT can be useful for assessing tumor changes early in treatment, if more than a few days elapse between simulation and the start of treatment. Tumor deviations measured with respiration-averaged CT and CBCT images are consistent with those measured with respiration-correlated images; the respiration-averaged method is more easily implemented in the clinic.

Evaluation of deformable image registration and a motion model in CT images with limited features.

Deformable image registration (DIR) is increasingly used in radiotherapy applications and provides the basis for a previously described model of patient-specific respiratory motion. We examine the accuracy of a DIR algorithm and a motion model with respiration-correlated CT (RCCT) images of software phantom with known displacement fields, physical deformable abdominal phantom with implanted fiducials in the liver and small liver structures in patient images. The motion model is derived from a principal component analysis that relates volumetric deformations with the motion of the diaphragm or fiducials in the RCCT. Patient data analysis compares DIR with rigid registration as ground truth: the mean ± standard deviation 3D discrepancy of liver structure centroid positions is 2.0 ± 2.2 mm. DIR discrepancy in the software phantom is 3.8 ± 2.0 mm in lung and 3.7 ± 1.8 mm in abdomen; discrepancies near the chest wall are larger than indicated by image feature matching. Marker's 3D discrepancy in the physical phantom is 3.6 ± 2.8 mm. The results indicate that visible features in the images are important for guiding the DIR algorithm. Motion model accuracy is comparable to DIR, indicating that two principal components are sufficient to describe DIR-derived deformation in these datasets.

SU-D-BRA-03: Simultaneous MV-KV Imaging for Intra-Fractional Motion Management during Volume Modulated Arc Therapy (VMAT) Delivery on the Varian TrueBeam.

PURPOSE: To evaluate a MV-kV intra-fractional imaging technique for use during volume modulated arc therapy (VMAT) with the Varian TrueBeam. METHODS: MV-kV image pairs were acquired intra-fractionally during VMAT delivery. kV images (11 fps) were acquired throughout delivery using a standard pre-programmed imaging template. MV images (9.5 fps) were acquired simultaneously by deploying the EPID and passively collecting the resulting images using Varian proprietary software, iTools Capture. Localization accuracy was evaluated by imaging a Rando phantom implanted with 3 fiducials while moving the couch according to XML- programmed trajectories simulating typical prostate and respiratory motion. VMAT delivery was done using a single 360 degree arc in TrueBeam Developer mode. The effect on accuracy of total MU and gantry speed was studied. To improve image quality, MV frame averaging was performed and the MV and kV images were then registered to their corresponding DRRs using in-house registration software. From these 2D registrations, the 3D position at each MV-kV acquisition point was determined. RESULTS: Between 130 and 390 MV-kV pairs were acquired for each delivery. The mean difference between planned couch and measured fiducial 3D positions with prostate motion was less than 0.03 cm in each direction (SD 0.03 cm). Neither gantry speed nor MU significantly impacted accuracy. for respiratory motion, the mean difference between planned and measured position was less than 0.04 cm. Standard deviation averaged 0.06 cm but increased to 0.12 cm with large instantaneous motion and less MV dose per frame. MV frame averaging and inaccuracies in MV image gantry angle determination also affected accuracy, particularly with significant motion. CONCLUSIONS: With high quality MV imaging, MV-kV localization techniques can be highly accurate, even in the presence of significant motion. As clinical MV-kV methods become available, such techniques can provide an efficient and accurate method for monitoring intra-fractional motion. This work was partially supported through a research agreement with Varian Medical Systems, Palo Alto, CA.

TU-E-BRA-09: Evaluation of a Patient-Specific Respiratory Motion Model in Thoracic and Abdominal Phantom and Patient CT Images.

PURPOSE: We have previously described a model of patient-specific respiratory motion to predict organ deformations without assuming repeatable breath cycles. The model is derived from deformable image registration (DIR) between respiration-correlated images (RCCT), followed by a principal component analysis (PCA) which relates the first two principal components of 3D deformations to the position and direction of motion of the diaphragm or implanted fiducials. This study examines model accuracy in phantom and patient images. METHODS: We compare model and DIR accuracy using 3 types of image sets, each exhibiting different deformation patterns: (1) synthetic images in lung and abdomen from the 4D NURBS-based cardiac torso (NCAT) phantom with known deformations; (2) CT scans of physical deformable phantom with implanted markers in liver; and (3) liver structures in patient RCCT images using rigid registration in a small VOI as approximate ground truth. The model is calibrated by applying fast free-form DIR between a reference image set at end expiration and each of the other images at different motion states, defined by diaphragm or, in some patient cases, implanted fiducials as surrogate signals. Following PCA, the first two principal components are selected to yield a model-predicted displacement field for the given surrogate signal. RESULTS: Discrepancy between model prediction and ground truth (mean ± stand deviation) in 3D displacements is 3.3±2.0 mm in lung and 3.7±1.9 mm in abdomen in NCAT phantom, 3.8±2.7 mm in physical deformable phantom and 2.8±2.9 mm in patient data (N=7). Corresponding DIR discrepancies are 3.8±2.0 mm (NCAT lung), 3.7±1.8 mm (NCAT abdomen), 3.6±2.8 mm (physical phantom), and 2.0±2.2 mm (patient data). CONCLUSIONS: Motion model accuracy is found to be comparable to fast free-form in all three types of images, indicating that the assumption of two principal components is sufficient to describe the fast free-form DIR-derived deformations. NIH/NCI award R01 CA126993.

SU-E-J-119: Comparative Evaluation of Respiratory Motion-Corrected Cone-Beam CT Images Derived from Treatment-Day Vs. Simulation-Day Respiration-Correlated CT Scans.

PURPOSE: Respiration-induced motion artifacts in cone-beam CT (CBCT) can be corrected using a model of patient motion obtained from respiration-correlated CT (RCCT). This approach assumes that respiration-induced organ deformations at simulation, when RCCT scans are normally acquired, are still valid at treatment. The purpose of this study is to compare lung tumor image quality in motion-corrected CBCT images derived from treatment-day RCCT(tx) to simulation-day RCCT(sim) patient images. METHODS: In an IRB-approved study, lung cancer patients receive an RCCT at simulation, and an RCCT, gated CBCT and 1-minute CBCT at one treatment session. CBCT projections from the 1-minute scan are sorted according to breathing amplitude from an external monitor and reconstructed and warped to obtain a motion-corrected MC-CBCT at end expiration. Motion correction uses a model adapted from either RCCT(tx) or RCCT(sim), thus obtaining MC-CBCT(tx) and MC-CBCT(sim) images respectively. A gated CBCT, in which gantry rotation and projection acquisition occur within a gate at end expiration, serves as ground truth for comparison. Quality of MC-CBCT images is evaluated from tumor-to-background contrast ratio (TBCR) values measured by delineating the tumor and annular volume around it on the gated CBCT then transferring the contours and aligning them to each MC-CBCT. RESULTS: TBCR is found tobe lower in MC-CBCT(sim) images, relative to MC-CBCT(tx), in four out of five patients with mean 21% reduction in a range 9-39%. In the remaining case, where there was no change in TBCR, tumor motion observed in the RCCT was small (2mm). Tumor motion extent relative to diaphragm is observed to change between RCCT(tx) and RCCT(sim) scans. CONCLUSIONS: Preliminary results indicate that deformation patterns in lung do change between simulation and treatment. Such variations may reduce the validity of using simulation data for motion-corrected CBCT at treatment. The findings require confirmation with larger numbers of patients. NIH/NCI award R01 CA126993, research grant from Varian Medical Systems.

Phase and amplitude binning for 4D-CT imaging.

We compare the consistency and accuracy of two image binning approaches used in 4D-CT imaging. One approach, phase binning (PB), assigns each breathing cycle 2pi rad, within which the images are grouped. In amplitude binning (AB), the images are assigned bins according to the breathing signal's full amplitude. To quantitate both approaches we used a NEMA NU2-2001 IEC phantom oscillating in the axial direction and at random frequencies and amplitudes, approximately simulating a patient's breathing. 4D-CT images were obtained using a four-slice GE Lightspeed CT scanner operating in cine mode. We define consistency error as a measure of ability to correctly bin over repeated cycles in the same field of view. Average consistency error mue+/-sigmae in PB ranged from 18%+/-20% to 30%+/-35%, while in AB the error ranged from 11%+/-14% to 20%+/-24%. In PB nearly all bins contained sphere slices. AB was more accurate, revealing empty bins where no sphere slices existed. As a proof of principle, we present examples of two non-small cell lung carcinoma patients' 4D-CT lung images binned by both approaches. While AB can lead to gaps in the coronal images, depending on the patient's breathing pattern, PB exhibits no gaps but suffers visible artifacts due to misbinning, yielding images that cover a relatively large amplitude range. AB was more consistent, though often resulted in gaps when no data existed due to patients' breathing pattern. We conclude AB is more accurate than PB. This has important consequences to treatment planning and diagnosis.

Reduction of organ motion effects in IMRT and conformal 3D radiation delivery by using gating and tracking techniques.

Respiration-gated radiotherapy offers a significant potential for improvement in the irradiation of tumour sites affected by respiratory motion such as lung, breast and liver tumours. An increased conformality of irradiation fields leading to decreased complications rates of organs at risk (lung, heart) is expected. Four main strategies are used to reduce respiratory motion effects: integration of respiratory movements into treatment planning, breath-hold techniques, respiratory gating techniques, and tracking techniques. Measurements of respiratory movements can be performed either in a representative sample of the general population, or directly on the patient before irradiation. The measured amplitude could be applied to a geometrical margin or integrated into dosimetry. However, these strategies remain limited for very mobile tumours, in which this approach results in larger irradiated volumes. Reduction of breathing motion can be achieved by using either breath-hold techniques or respiration synchronized gating techniques. Breath-hold can be achieved with active techniques, in which a valve temporarily blocks airflow of the patient, or passive techniques, in which the patient voluntarily breath-holds. Synchronized gating techniques use external devices to predict the phase of the respiration cycle while the patient breaths freely. Another category is tumour tracking, which consists of two major aspects: real-time localization of, and real-time beam adaptation to, a constantly moving tumour. These techniques are presently being investigated in several medical centres worldwide. Although promising, the first results obtained in lung and liver cancer patients require confirmation. This paper describes the most frequently used gating and tracking techniques and the main published clinical reports.

Evaluation of an automated deformable image matching method for quantifying lung motion in respiration-correlated CT images.

We have evaluated an automated registration procedure for predicting tumor and lung deformation based on CT images of the thorax obtained at different respiration phases. The method uses a viscous fluid model of tissue deformation to map voxels from one CT dataset to another. To validate the deformable matching algorithm we used a respiration-correlated CT protocol to acquire images at different phases of the respiratory cycle for six patients with nonsmall cell lung carcinoma. The position and shape of the deformable gross tumor volumes (GTV) at the end-inhale (EI) phase predicted by the algorithm was compared to those drawn by four observers. To minimize interobserver differences, all observers used the contours drawn by a single observer at end-exhale (EE) phase as a guideline to outline GTV contours at EI. The differences between model-predicted and observer-drawn GTV surfaces at EI, as well as differences between structures delineated by observers at EI (interobserver variations) were evaluated using a contour comparison algorithm written for this purpose, which determined the distance between the two surfaces along different directions. The mean and 90% confidence interval for model-predicted versus observer-drawn GTV surface differences over all patients and all directions were 2.6 and 5.1 mm, respectively, whereas the mean and 90% confidence interval for interobserver differences were 2.1 and 3.7 mm. We have also evaluated the algorithm's ability to predict normal tissue deformations by examining the three-dimensional (3-D) vector displacement of 41 landmarks placed by each observer at bronchial and vascular branch points in the lung between the EE and EI image sets (mean and 90% confidence interval displacements of 11.7 and 25.1 mm, respectively). The mean and 90% confidence interval discrepancy between model-predicted and observer-determined landmark displacements over all patients were 2.9 and 7.3 mm, whereas interobserver discrepancies were 2.8 and 6.0 mm. Paired t tests indicate no significant statistical differences between model predicted and observer drawn structures. We conclude that the accuracy of the algorithm to map lung anatomy in CT images at different respiratory phases is comparable to the variability in manual delineation. This method has therefore the potential for predicting and quantifying respiration-induced tumor motion in the lung.

Developments in megavoltage cone beam CT with an amorphous silicon EPID: reduction of exposure and synchronization with respiratory gating.

We have studied the feasibility of a low-dose megavoltage cone beam computed tomography (MV CBCT) system for visualizing the gross tumor volume in respiratory gated radiation treatments of nonsmall-cell lung cancer. The system consists of a commercially available linear accelerator (LINAC), an amorphous silicon electronic portal imaging device, and a respiratory gating system. The gantry movement and beam delivery are controlled using dynamic beam delivery toolbox, a commercial software package for executing scripts to control the LINAC. A specially designed interface box synchronizes the LINAC, image acquisition electronics, and the respiratory gating system. Images are preprocessed to remove artifacts due to detector sag and LINAC output fluctuations. We report on the output, flatness, and symmetry of the images acquired using different imaging parameters. We also examine the quality of three-dimensional (3D) tomographic reconstruction with projection images of anthropomorphic thorax, contrast detail, and motion phantoms. The results show that, with the proper choice of imaging parameters, the flatness and symmetry are reasonably good with as low as three beam pulses per projection image. Resolution of 5% electron density differences is possible in a contrast detail phantom using 100 projections and 30 MU. Synchronization of image acquisition with simulated respiration also eliminated motion artifacts in a moving phantom, demonstrating the system's capability for imaging patients undergoing gated radiation therapy. The acquisition time is limited by the patient's respiration (only one image per breathing cycle) and is under 10 min for a scan of 100 projections. In conclusion, we have developed a MV CBCT system using commercially available components to produce 3D reconstructions, with sufficient contrast resolution for localizing a simulated lung tumor, using a dose comparable to portal imaging.

A method for determining the gantry angle for megavoltage cone beam imaging.

Accurate knowledge of gantry angle is essential in megavoltage cone beam imaging (MVCBI) with an electronic portal imager. We present a method for determining the gantry angle by detecting multileaf collimator (MLC) leaf positions in projection images. During image acquisition the gantry moves continuously and the MLC operates in dynamic arc mode. Our algorithm detects the leaf positions in the images and compares them with a stationary reference leaf. Comparison of the algorithm against angles determined from the locations of fiducial markers shows the accuracy (0.26 degrees rms error) to be sufficient for MVCBI.

Quantitation of respiratory motion during 4D-PET/CT acquisition.

We report on the variability of the respiratory motion during 4D-PET/CT acquisition. The respiratory motion for five lung cancer patients was monitored by tracking external markers placed on the abdomen. CT data were acquired over an entire respiratory cycle at each couch position. The x-ray tube status was recorded by the tracking system, for retrospective sorting of the CT data as a function of respiration phase. Each respiratory cycle was sampled in ten equal bins. 4D-PET data were acquired in gated mode, where each breathing cycle was divided into ten 500 ms bins. For both CT and PET acquisition, patients received audio prompting to regularize breathing. The 4D-CT and 4D-PET data were then correlated according to their respiratory phases. The respiratory periods, and average amplitude within each phase bin, acquired in both modality sessions were then analyzed. The average respiratory motion period during 4D-CT was within 18% from that in the 4D-PET sessions. This would reflect up to 1.8% fluctuation in the duration of each 4D-CT bin. This small uncertainty enabled good correlation between CT and PET data, on a phase-to-phase basis. Comparison of the average-amplitude within the respiration trace, between 4D-CT and 4D- PET, on a bin-by-bin basis show a maximum deviation of approximately 15%. This study has proved the feasibility of performing 4D-PET/CT acquisition. Respiratory motion was in most cases consistent between PET and CT sessions, thereby improving both the attenuation correction of PET images, and co-registration of PET and CT images. On the other hand, in two patients, there was an increased partial irregularity in their breathing motion, which would prevent accurately correlating the corresponding PET and CT images.

Four-dimensional (4D) PET/CT imaging of the thorax.

We have reported in our previous studies on the methodology, and feasibility of 4D-PET (Gated PET) acquisition, to reduce respiratory motion artifact in PET imaging of the thorax. In this study, we expand our investigation to address the problem of respiration motion in PET/CT imaging. The respiratory motion of four lung cancer patients were monitored by tracking external markers placed on the thorax. A 4D-CT acquisition was performed using a "step-and-shoot" technique, in which computed tomography (CT) projection data were acquired over a complete respiratory cycle at each couch position. The period of each CT acquisition segment was time stamped with an "x-ray ON" signal, which was recorded by the tracking system. 4D-CT data were then sorted into 10 groups, according to their corresponding phase of the breathing cycle. 4D-PET data were acquired in the gated mode, where each breathing cycle was divided into ten 0.5 s bins. For both CT and PET acquisitions, patients received audio prompting to regularize breathing. The 4D-CT and 4D-PET data were then correlated according to respiratory phase. The effect of 4D acquisition on improving the co-registration of PET and CT images, reducing motion smearing, and consequently increase the quantitation of the SUV, were investigated. Also, quantitation of the tumor motions in PET, and CT, were studied and compared. 4D-PET with matching phase 4D-CTAC showed an improved accuracy in PET-CT image co-registration of up to 41%, compared to measurements from 4D-PET with clinical-CTAC. Gating PET data in correlation with respiratory motion reduced motion-induced smearing, thereby decreasing the observed tumor volume, by as much as 43%. 4D-PET lesions volumes showed a maximum deviation of 19% between clinical CT and phase- matched 4D-CT attenuation corrected PET images. In CT, 4D acquisition resulted in increasing the tumor volume in two patients by up to 79%, and decreasing it in the other two by up to 35%. Consequently, these corrections have yielded an increase in the measured SUV by up to 16% over the clinical measured SUV, and 36% over SUV's measured in 4D-PET with clinical-CT Attenuation Correction (CTAC) SUV's. Quantitation of the maximum tumor motion amplitude, using 4D-PET and 4D-CT, showed up to 30% discrepancy between the two modalities. We have shown that 4D PET/CT is clinically a feasible method, to correct for respiratory motion artifacts in PET/CT imaging of the thorax. 4D PET/CT acquisition can reduce smearing, improve the accuracy in PET-CT co-registration, and increase the measured SUV. This should result in an improved tumor assessment for patients with lung malignancies.

Respiration-correlated spiral CT: a method of measuring respiratory-induced anatomic motion for radiation treatment planning.

We describe a method for generating CT images at multiple respiratory phases with a single spiral CT scan, referred to as respiratory-correlated spiral CT (RCCT). RCCT relies on a respiration wave form supplied by an external patient monitor. During acquisition this wave form is recorded along with the initiation time of the CT scan, so as to "time stamp" each reconstructed slice with the phase of the respiratory cycle. By selecting the appropriate slices, a full CT image set is generated at several phases, typically 7-11 per cycle. The CT parameters are chosen to optimize the temporal resolution while minimizing the spatial gap between slices at successive respiratory cycles. Using a pitch of 0.5, a gantry rotation period of 1.5 s, and a 180 degrees reconstruction algorithm results in approximately 5 mm slice spacing at a given phase for typical respiration periods, and a respiratory motion within each slice that is acceptably small, particularly near end expiration or end inspiration where gated radiotherapy is to occur. We have performed validation measurements on a phantom with a moving sphere designed to simulate respiration-induced tumor motion. RCCT scans of the phantom at respiratory periods of 4, 5, and 6 s show good agreement of the sphere's motion with that observed under fluoroscopic imaging. The positional deviations in the sphere's centroid between RCCT and fluoroscopy are 1.1+/-0.9 mm in the transaxial direction (average over all scans at all phases +/-1 s.d.) and 1.2+/-1.0 mm in the longitudinal direction. Reconstructed volumes match those expected on the basis of stationary-phantom scans to within 5% in all cases. The surface distortions of the reconstructed sphere, as quantified by deviations from a mathematical reference sphere, are similar to those from a stationary phantom scan and are correlated with the speed of the phantom. A RCCT scan of the phantom undergoing irregular motion, demonstrates that successful reconstruction can be achieved even with irregular respiration. Limitations from x-ray tube heating in our current CT unit restrict the length of the scan region to 9 cm for the RCCT settings used, though this will not be a limitation for a multislice scanner. RCCT offers an alternative to the current method of respiration-triggered axial scans. Multiple phases of respiration are imaged with RCCT in approximately the same scanning time required to image a single phase with a triggered axial scan. RCCT scans can be used in connection with respiratory-gated treatment to identify the patient-specific phase of minimum tumor motion, determine residual tumor motion within the gate interval, and compare treatment plans at different phases.

Effect of respiratory gating on reducing lung motion artifacts in PET imaging of lung cancer.

Positron emission tomography (PET) has shown an increase in both sensitivity and specificity over computed tomography (CT) in lung cancer. However, motion artifacts in the 18F fluorodioxydoglucose (FDG) PET images caused by respiration persists to be an important factor in degrading PET image quality and quantification. Motion artifacts lead to two major effects: First, it affects the accuracy of quantitation, producing a reduction of the measured standard uptake value (SUV). Second, the apparent lesion volume is overestimated. Both impact upon the usage of PET images for radiation treatment planning. The first affects the visibility, or contrast, of the lesion. The second results in an increase in the planning target volume, and consequently a greater radiation dose to the normal tissues. One way to compensate for this effect is by applying a multiple-frame capture technique. The PET data are then acquired in synchronization with the respiratory motion. Reduction in smearing due to gating was investigated in both phantoms and patient studies. Phantom studies showed a dependence of the reduction in smearing on the lesion size, the motion amplitude, and the number of bins used for data acquisition. These studies also showed an improvement in the target-to-background ratio, and a more accurate measurement of the SUV. When applied to one patient, respiratory gating showed a 28% reduction in the total lesion volume, and a 56.5% increase in the SUV. This study was conducted as a proof of principle that a gating technique can effectively reduce motion artifacts in PET image acquisition.

Evaluation of respiratory movement during gated radiotherapy using film and electronic portal imaging.

PURPOSE: To evaluate the effectiveness of a commercial system(1) in reducing respiration-induced treatment uncertainty by gating the radiation delivery. METHODS AND MATERIALS: The gating system considered here measures respiration from the position of a reflective marker on the patient's chest. Respiration-triggered planning CT scans were obtained for 8 patients (4 lung, 4 liver) at the intended phase of respiration (6 at end expiration and 2 at end inspiration). In addition, fluoroscopic movies were recorded simultaneously with the respiratory waveform. During the treatment sessions, gated localization films were used to measure the position of the diaphragm relative to the vertebral bodies, which was compared to the reference digitally reconstructed radiograph derived from the respiration-triggered planning CT. Variability was quantified by the standard deviation about the mean position. We also assessed the interfraction variability of soft tissue structures during gated treatment in 2 patients using an amorphous silicon electronic portal imaging device. RESULTS: The gated localization films revealed an interfraction patient-averaged diaphragm variability of 2.8 +/- 1.0 mm (error bars indicate standard deviation in the patient population). The fluoroscopic data yielded a patient-averaged intrafraction diaphragm variability of 2.6 +/- 1.7 mm. With no gating, this intrafraction excursion became 6.9 +/- 2.1 mm. In gated localization films, the patient-averaged mean displacement of the diaphragm from the planning position was 0.0 +/- 3.9 mm. However, in 4 of the 8 patients, the mean (over localization films) displacement was >4 mm, indicating a systematic displacement in treatment position from the planned one. The position of soft tissue features observed in portal images during gated treatments over several fractions showed a mean variability between 2.6 and 5.7 mm. The intrafraction variability, however, was between 0.6 and 1.4 mm, indicating that most of the variability was due to patient setup errors rather than to respiratory motion. CONCLUSIONS: The gating system evaluated here reduces the intra- and interfraction variability of anatomy due to respiratory motion. However, systematic displacements were observed in some cases between the location of an anatomic feature at simulation and its location during treatment. Frequent monitoring is advisable with film or portal imaging.

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.

Fluoroscopic evaluation of diaphragmatic motion reduction with a respiratory gated radiotherapy system.

We report on initial patient studies to evaluate the performance of a commercial respiratory gating radiotherapy system. The system uses a breathing monitor, consisting of a video camera and passive infrared reflective markers placed on the patient's thorax, to synchronize radiation from a linear accelerator with the patient's breathing cycle. Six patients receiving treatment for lung cancer participated in a study of system characteristics during treatment simulation with fluoroscopy. Breathing synchronized fluoroscopy was performed initially without instruction, followed by fluoroscopy with recorded verbal instruction (i.e., when to inhale and exhale) with the tempo matched to the patient's normal breathing period. Patients tended to inhale more consistently when given instruction, as assessed by an external marker movement. This resulted in smaller variation in expiration and inspiration marker positions relative to total excursion, thereby permitting more precise gating tolerances at those parts of the breathing cycle. Breathing instruction also reduced the fraction of session times having irregular breathing as measured by the system software, thereby potentially increasing the accelerator duty factor and decreasing treatment times. Fluoroscopy studies showed external monitor movement to correlate well with that of the diaphragm in four patients, whereas time delays of up to 0.7 s in diaphragm movement were observed in two patients with impaired lung function. From fluoroscopic observations, average patient diaphragm excursion was reduced from 1.4 cm (range 0.7-2.1 cm) without gating and without breathing instruction, to 0.3 cm (range 0.2-0.5 cm) with instruction and with gating tolerances set for treatment at expiration for 25% of the breathing cycle. Patients expressed no difficulty with following instruction for the duration of a session. We conclude that the external monitor accurately predicts internal respiratory motion in most cases; however, it may be important to check with fluoroscopy for possible time delays in patients with impaired lung function. Furthermore, we observe that verbal instruction can improve breathing regularity, thus improving the performance of gated treatments with this system.

An iterative EPID calibration procedure for dosimetric verification that considers the EPID scattering factor.

There has been an increasing interest in the application of electronic portal imaging devices (EPIDs) to dosimetric verification, particularly for intensity modulated radiotherapy. Although not water equivalent, the phantom scatter factor of an EPID, Spe, is generally assumed to be that of a full phantom, Sp, a slab phantom, Sps, or a mini phantom. This assumption may introduce errors in absolute dosimetry using EPIDs. A calibration procedure that iteratively updates Spe and the calibration curve (pixel value to dose rate) is presented. The EPID (Varian Portal Vision) is irradiated using a 20 x 20 cm2 field with different beam intensities. The initial guess of dose rates in the EPID is calculated from ionization chamber measurements in air, multiplied by Sp or Sps. The calibration curve is obtained by fitting EPID readings from pixels near the beam central axis and dose rates in EPID to a quadratic equation. The Spe is obtained from EPID measurements in 10 X 10 cm2 and 20 x 20 cm2 field and from the calibration curve, and is in turn used to adjust the dose rate measurements and hence the calibration curve. The above procedure is repeated until it converges. The final calibration curve is used to convert portal dose to dose in the slab phantom, using the calibrated Spe, or assuming Spe = Sp or Spe=Sps . The converted doses are then compared with the dose measured using an ionization chamber. We also apply this procedure to off-axis points and study its dependence on the energy spectrum. The hypothesis testing results (on the 95% significance level) indicate that systematic errors are introduced when assuming Spe = Sp or Spe=Sps and the dose calculated using Spe is more consistent with ionization chamber measurements. Differences between Spe and Sps are as large as 2% for large field sizes. The measured relative dose profile at dmax using the EPID agrees well with the measured profile at dmax of the isocentric plane using film in a polystyrene phantom with full buildup and full backup, for open and wedged fields, and for a broad range of field sizes of interest. The dependence of the EPID response on the energy spectrum is removed once the calibration is performed under the same conditions as the actual measurements.

Technical aspects of the deep inspiration breath-hold technique in the treatment of thoracic cancer.

PURPOSE: The goal of this paper is to describe our initial experience with the deep inspiration breath-hold (DIBH) technique in conformal treatment of non-small-cell lung cancer with particular emphasis on the technical aspects required for implementation. METHODS AND MATERIALS: In the DIBH technique, the patient is verbally coached through a modified slow vital capacity maneuver and brought to a reproducible deep inspiration breath-hold level. The goal is to immobilize the tumor and to expand normal lung out of the high-dose region. A physicist or therapist monitors and records patient breathing during simulation, verification, and treatment using a spirometer with a custom computer interface. Examination of internal anatomy during fluoroscopy over multiple breath holds establishes the reproducibility of the DIBH maneuver for each patient. A reference free-breathing CT scan and DIBH planning scan are obtained. To provide an estimate of tumor motion during normal tidal breathing, additional scan sets are obtained at end inspiration and end expiration. These are also used to set the spirometer action levels for treatment. Patient lung inflation is independently verified over the course of treatment by comparing the distance from the isocenter to the diaphragm measured from the DIBH digitally reconstructed radiographs to the distance measured on the portal films. Patient breathing traces obtained during treatment were examined retrospectively to assess the reproducibility of the technique. RESULTS: Data from the first 7 patients, encompassing over 250 treatments, were analyzed. The inferred displacement of the centroid of gross tumor volume from its position in the planning scan, as calculated from the spirometer records in over 350 breath holds was 0.02 +/- 0.14 cm (mean and standard deviation). These data are consistent with the displacements of the diaphragm (-0.1 +/- 0.4 cm; range, from -1.2 to 1.1 cm) relative to the isocenter, as measured on the (92) portal films. The latter measurements include the patient setup error. The patient averaged displacement of the tumor during free breathing, determined from the tumor displacement between end inspiration and end expiration, was 0.8 +/- 0.5 cm in both the superior-inferior and anterior-posterior directions and 0.1 cm (+/- 0.1 cm) medial-laterally. CONCLUSION: Treatment of patients with the DIBH technique is feasible in a clinical setting. With this technique, consistent lung inflation levels are achieved in patients, as judged by both spirometry and verification films. Breathing-induced tumor motion is significantly reduced using DIBH compared to free breathing, enabling better target coverage.

The deep inspiration breath-hold technique in the treatment of inoperable non-small-cell lung cancer.

PURPOSE: Conventional radiotherapeutic techniques are associated with lung toxicity that limits the treatment dose. Motion of the tumor during treatment requires the use of large safety margins that affect the feasibility of treatment. To address the control of tumor motion and decrease the volume of normal lung irradiated, we investigated the use of three-dimensional conformal radiation therapy (3D-CRT) in conjunction with the deep inspiration breath-hold (DIBH) technique. METHODS AND MATERIALS: In the DIBH technique, the patient is initially maintained at quiet tidal breathing, followed by a deep inspiration, a deep expiration, a second deep inspiration, and breath-hold. At this point the patient is at approximately 100% vital capacity, and simulation, verification, and treatment take place during this phase of breath-holding. RESULTS: Seven patients have received a total of 164 treatment sessions and have tolerated the technique well. The estimated normal tissue complication probabilities decreased in all patients at their prescribed dose when compared to free breathing. The dose to which patients could be treated with DIBH increased on average from 69.4 Gy to 87.9 Gy, without increasing the risk of toxicity. CONCLUSIONS: The DIBH technique provides an advantage to conventional free-breathing treatment by decreasing lung density, reducing normal safety margins, and enabling more accurate treatment. These improvements contribute to the effective exclusion of normal lung tissue from the high-dose region and permit the use of higher treatment doses without increased risks of toxicity.