Thomas Guerrero

A framework for evaluation of deformable image registration spatial accuracy using large landmark point sets.

Expert landmark correspondences are widely reported for evaluating deformable image registration (DIR) spatial accuracy. In this report, we present a framework for objective evaluation of DIR spatial accuracy using large sets of expert-determined landmark point pairs. Large samples (>1100) of pulmonary landmark point pairs were manually generated for five cases. Estimates of inter- and intra-observer variation were determined from repeated registration. Comparative evaluation of DIR spatial accuracy was performed for two algorithms, a gradient-based optical flow algorithm and a landmark-based moving least-squares algorithm. The uncertainty of spatial error estimates was found to be inversely proportional to the square root of the number of landmark point pairs and directly proportional to the standard deviation of the spatial errors. Using the statistical properties of this data, we performed sample size calculations to estimate the average spatial accuracy of each algorithm with 95% confidence intervals within a 0.5 mm range. For the optical flow and moving least-squares algorithms, the required sample sizes were 1050 and 36, respectively. Comparative evaluation based on fewer than the required validation landmarks results in misrepresentation of the relative spatial accuracy. This study demonstrates that landmark pairs can be used to assess DIR spatial accuracy within a narrow uncertainty range.

Implementation and evaluation of various demons deformable image registration algorithms on a GPU

Online adaptive radiation therapy (ART) promises the ability to deliver an optimal treatment in response to daily patient anatomic variation. A major technical barrier for the clinical implementation of online ART is the requirement of rapid image segmentation. Deformable image registration (DIR) has been used as an automated segmentation method to transfer tumor/organ contours from the planning image to daily images. However, the current computational time of DIR is insufficient for online ART. In this work, this issue is addressed by using computer graphics processing units (GPUs). A gray-scale-based DIR algorithm called demons and five of its variants were implemented on GPUs using the compute unified device architecture (CUDA) programming environment. The spatial accuracy of these algorithms was evaluated over five sets of pulmonary 4D CT images with an average size of 256 x 256 x 100 and more than 1100 expert-determined landmark point pairs each. For all the testing scenarios presented in this paper, the GPU-based DIR computation required around 7 to 11 s to yield an average 3D error ranging from 1.5 to 1.8 mm. It is interesting to find out that the original passive force demons algorithms outperform subsequently proposed variants based on the combination of accuracy, efficiency and ease of implementation.

Dynamic ventilation imaging from four-dimensional computed tomography

A novel method for dynamic ventilation imaging of the full respiratory cycle from four-dimensional computed tomography (4D CT) acquired without added contrast is presented. Three cases with 4D CT images obtained with respiratory gated acquisition for radiotherapy treatment planning were selected. Each of the 4D CT data sets was acquired during resting tidal breathing. A deformable image registration algorithm mapped each (voxel) corresponding tissue element across the 4D CT data set. From local average CT values, the change in fraction of air per voxel (i.e. local ventilation) was calculated. A 4D ventilation image set was calculated using pairs formed with the maximum expiration image volume, first the exhalation then the inhalation phases representing a complete breath cycle. A preliminary validation using manually determined lung volumes was performed. The calculated total ventilation was compared to the change in contoured lung volumes between the CT pairs (measured volume). A linear regression resulted in a slope of 1.01 and a correlation coefficient of 0.984 for the ventilation images. The spatial distribution of ventilation was found to be case specific and a 30% difference in mass-specific ventilation between the lower and upper lung halves was found. These images may be useful in radiotherapy planning.

Four-dimensional deformable image registration using trajectory modeling

A four-dimensional deformable image registration (4D DIR) algorithm, referred to as 4D local trajectory modeling (4DLTM), is presented and applied to thoracic 4D computed tomography (4DCT) image sets. The theoretical framework on which this algorithm is built exploits the incremental continuity present in 4DCT component images to calculate a dense set of parameterized voxel trajectories through space as functions of time. The spatial accuracy of the 4DLTM algorithm is compared with an alternative registration approach in which component phase to phase (CPP) DIR is utilized to determine the full displacement between maximum inhale and exhale images. A publically available DIR reference database (http://www.dir-lab.com) is utilized for the spatial accuracy assessment. The database consists of ten 4DCT image sets and corresponding manually identified landmark points between the maximum phases. A subset of points are propagated through the expiratory 4DCT component images. Cubic polynomials were found to provide sufficient flexibility and spatial accuracy for describing the point trajectories through the expiratory phases. The resulting average spatial error between the maximum phases was 1.25 mm for the 4DLTM and 1.44 mm for the CPP. The 4DLTM method captures the long-range motion between 4DCT extremes with high spatial accuracy.

Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method

The purpose of this work was to develop and validate an automated method for intrathoracic tumour motion estimation from breath-hold computed tomography (BH CT) imaging using the three-dimensional optical flow method (3D OFM). A modified 3D OFM algorithm provided 3D displacement vectors for each voxel which were used to map tumour voxels on expiration BH CT onto inspiration BH CT images. A thoracic phantom and simulated expiration/inspiration BH CT pairs were used for validation. The 3D OFM was applied to the measured inspiration and expiration BH CT images from one lung cancer and one oesophageal cancer patient. The resulting displacements were plotted in histogram format and analysed to provide insight regarding the tumour motion. The phantom tumour displacement was measured as 1.20 and 2.40 cm with full-width at tenth maximum (FWTM) for the distribution of displacement estimates of 0.008 and 0.006 cm, respectively. The maximum error of any single voxels motion estimate was 1.1 mm along the z-dimension or approximately one-third of the z-dimension voxel size. The simulated BH CT pairs revealed an rms error of less than 0.25 mm. The displacement of the oesophageal tumours was nonuniform and up to 1.4 cm, this was a new finding. A lung tumour maximum displacement of 2.4 cm was found in the case evaluated. In conclusion, 3D OFM provided an accurate estimation of intrathoracic tumour motion, with estimated errors less than the voxel dimension in a simulated motion phantom study. Surprisingly, oesophageal tumour motion was large and nonuniform, with greatest motion occurring at the gastro-oesophageal junction.

Ventilation from four-dimensional computed tomography: density versus Jacobian methods

Two calculation methods to produce ventilation images from four-dimensional computed tomography (4DCT) acquired without added contrast have been reported. We reported a method to obtain ventilation images using deformable image registration (DIR) and the underlying CT density information. A second method performs the ventilation image calculation from the DIR result alone, using the Jacobian determinant of the deformation field to estimate the local volume changes resulting from ventilation. For each of these two approaches, there are variations on their implementation. In this study, two implementations of the Jacobian-based methodology are evaluated, as well as a single density change-based model for calculating the physiologic specific ventilation from 4DCT. In clinical practice, (99m)Tc-labeled aerosol single photon emission computed tomography (SPECT) is the standard method used to obtain ventilation images in patients. In this study, the distributions of ventilation obtained from the CT-based ventilation image calculation methods are compared with those obtained from the clinical standard SPECT ventilation imaging. Seven patients with 4DCT imaging and standard (99m)Tc-labeled aerosol SPECT/CT ventilation imaging obtained on the same day as part of a prospective validation study were selected. The results of this work demonstrate the equivalence of the Jacobian-based methodologies for quantifying the specific ventilation on a voxel scale. Additionally, we found that both Jacobian- and density-change-based methods correlate well with global measurements of the resting tidal volume. Finally, correlation with the clinical SPECT was assessed using the Dice similarity coefficient, which showed statistically higher (p-value < 10(-4)) correlation between density-change-based specific ventilation and the clinical reference than did either Jacobian-based implementation.

Promising early local control of malignant pleural mesothelioma following postoperative intensity modulated radiotherapy (IMRT) to the chest

PURPOSEMalignant pleural mesothelioma often recurs locally in spite of aggressive resection by extrapleural pneumonectomy and conventional radiotherapy. This may be due to failure to recognize the extent of clinical target volume (CTV) or suboptimal dose delivery to a target that abuts the heart, esophagus, liver, lung, kidney, and spinal cord. We report how these geometric/dosimetric constraints were overcome by exploiting intensity-modulated radiotherapy in the first cohort patient. MATERIALS AND METHODSTwenty-eight patients who had undergone extrapleural pneumonectomy were treated with intensity-modulated radiotherapy. The CTV included the surgically violated inner chest wall, insertion of diaphragm, pleural reflections, and deep margin of the incision. CTV delineation was facilitated by intraoperative radio-opaque marking. Motion was assessed. CTV doses were 45–50 Gy with boosts taken to 60 Gy. RESULTSDespite the large, irregular CTV (median, 4151 cc; range, 2667–7286 cc), an average of 97% of the CTV was covered to the target dose (range, 92%–100%). Respiratory motion was minimal because of immobility of the prosthetic diaphragm. Normal tissue dose constraints were met. The commonest effects were nausea/vomiting (89%) and dyspnea (80%). Esophagitis was absent (59% of patients) or mild (34% grade 1/2). At median follow-up of 9 months (range, 5–27 months), local control within the contoured target was 100%. One-year survival, disease-specific survival, and disease-free survival are 65%, 91%, and 88%, respectively. CONCLUSIONSIntensity-modulated radiotherapy after extrapleural pneumonectomy is tolerable and seems effective, at least at this early point. As local control improves, systemic metastases become more common, and it may be appropriate to add novel agents to further improve the therapeutic ratio.

Four‐dimensional cone beam CT with adaptive gantry rotation and adaptive data sampling

We have developed a new four-dimensional cone beam CT (4D-CBCT) on a Varian image-guided radiation therapy system, which has radiation therapy treatment and cone beam CT imaging capabilities. We adapted the speed of gantry rotation time of the CBCT to the average breath cycle of the patient to maintain the same level of image quality and adjusted the data sampling frequency to keep a similar level of radiation exposure to the patient. Our design utilized the real-time positioning and monitoring system to record the respiratory signal of the patient during the acquisition of the CBCT data. We used the full-fan bowtie filter during data acquisition, acquired the projection data over 200 deg of gantry rotation, and reconstructed the images with a half-scan cone beam reconstruction. The scan time for a 200-deg gantry rotation per patient ranged from 3.3 to 6.6 min for the average breath cycle of 3-6 s. The radiation dose of the 4D-CBCT was about 1-2 times the radiation dose of the 4D-CT on a multislice CT scanner. We evaluated the 4D-CBCT in scanning, data processing and image quality with phantom studies. We demonstrated the clinical applicability of the 4D-CBCT and compared the 4D-CBCT and the 4D-CT scans in four patient studies. The contrast-to-noise ratio of the 4D-CT was 2.8-3.5 times of the contrast-to-noise ratio of the 4D-CBCT in the four patient studies.

Four-dimensional computed tomography-based treatment planning for intensity-modulated radiation therapy and proton therapy for distal esophageal cancer.

PURPOSE To compare three-dimensional (3D) and four-dimensional (4D) computed tomography (CT)-based treatment plans for proton therapy or intensity-modulated radiation therapy (IMRT) for esophageal cancer in terms of doses to the lung, heart, and spinal cord and variations in target coverage and normal tissue sparing. METHODS AND MATERIALS The IMRT and proton plans for 15 patients with distal esophageal cancer were designed from the 3D average CT scans and then recalculated on 10 4D CT data sets. Dosimetric data were compared for tumor coverage and normal tissue sparing. RESULTS Compared with IMRT, median lung volumes exposed to 5, 10, and 20 Gy and mean lung dose were reduced by 35.6%, 20.5%, 5.8%, and 5.1 Gy for a two-beam proton plan and by 17.4%, 8.4%, 5%, and 2.9 Gy for a three-beam proton plan. The greater lung sparing in the two-beam proton plan was achieved at the expense of less conformity to the target (conformity index [CI], 1.99) and greater irradiation of the heart (heart-V40, 41.8%) compared with the IMRT plan(CI, 1.55, heart-V40, 35.7%) or the three-beam proton plan (CI, 1.46, heart-V40, 27.7%). Target coverage differed by more than 2% between the 3D and 4D plans for patients with substantial diaphragm motion in the three-beam proton and IMRT plans. The difference in spinal cord maximum dose between 3D and 4D plans could exceed 5 Gy for the proton plans partly owing to variations in stomach gas filling. CONCLUSIONS Proton therapy provided significantly better sparing of lung than did IMRT. Diaphragm motion and stomach gas-filling must be considered in evaluating target coverage and cord doses.

A Phase I Clinical Trial of Thoracic Radiotherapy and Concurrent Celecoxib for Patients with Unfavorable Performance Status Inoperable/Unresectable Non–Small Cell Lung Cancer

Objectives: Preclinical observations that selective cyclooxygenase-2 inhibitors enhance in vitro cell radiosensitivity and in vivo tumor radioresponse led to clinical trials testing therapeutic efficacy of these agents. Our study was designed to determine whether the COX-2 inhibitor celecoxib could be safely administered in doses within those approved by the Food and Drug Administration when used concurrently with thoracic radiotherapy in patients with poor prognosis non–small cell lung cancer (NSCLC). Patients and Methods: The trial consisted of three cohorts of patients: (a) locally advanced NSCLC with obstructive pneumonia, hemoptysis, and/or minimal metastatic disease treated with 45 Gy in 15 fractions; (b) medically inoperable early-stage NSCLC treated with definitive radiation of 66 Gy in 33 fractions; and (c) patients who received induction chemotherapy but who were not eligible for concurrent chemoradiotherapy trials. These patients received 63 Gy in 35 fractions. Celecoxib was administered p.o. on a daily basis 5 days before and throughout the course of radiotherapy. Celecoxib doses were escalated from 200, 400, 600, to 800 mg/d given in two equally divided doses. Two to eight patients of each cohort were assigned to each dose level of celecoxib. Results: Forty-seven patients were enrolled in this protocol (19 in cohort I, 22 in cohort II, and 6 in cohort III). The main toxicities were grades 1 and 2 nausea and esophagitis, and they were independent of the dose of celecoxib or radiotherapy schedule. Only two patients in group II developed grade 3 pneumonitis 1 month after treatment, one on 200 mg, and the other on 400 mg celecoxib. Celecoxib-related toxicity developed in 3 of 47 patients: an uncontrolled hypertension in one patient on 800 mg celecoxib and hemorrhagic episodes in 2 patients (shoulder hematoma in one and hemoptysis in the other) on 200 mg celecoxib who were on warfarin for other medical reasons. Of 37 patients evaluable for tumor response, 14 had complete response, 13 partial responses, and 10 stable or progressive disease. The actuarial local progression-free survival was 66.0% at 1 year and 42.2% at 2 years following initiation of radiotherapy. Conclusions: These results show that celecoxib can be safely administered concurrently with thoracic radiotherapy when given up to the highest Food and Drug Administration–approved dose of 800 mg/d, which we used. A maximal tolerated dose was not reached in this study. The treatment resulted in actuarial local progression-free survival of 66.0% at 1 year and 42.2% at 2 years, an encouraging outcome that warrants further assessment in a phase II/III trial.