|
|
||||||||
Thoracic Imaging |
1 From the Departments of Medical Physics and Radiology (B.Z.), Radiology (L.H.S., M.S.G.), Epidemiology and Biostatistics (C.S.M.), and Medicine (N.A.R., M.G.K.), Memorial Sloan-Kettering Cancer Center, 1275 York Ave, New York, NY 10021. Received November 18, 2005; revision requested January 5, 2006; revision received January 26; accepted February 8; final version accepted March 28. Supported in part by grants from William H. Goodwin and Alice Goodwin, the Commonwealth Foundation for Cancer Research, the Experimental Therapeutics Center of Memorial Sloan-Kettering Cancer Center, and National Cancer Institute (R21 CA113653-01 and PO1 CA005826-39). Address correspondence to B.Z. (e-mail: zhaob{at}mskcc.org).
| ABSTRACT |
|---|
|
|
|---|
Materials and Methods: This HIPAA-compliant study was institutional review board approved; informed patient consent was waived. CT scans of 15 measurable nonsmall cell lung cancers (in five men and 10 women; mean age, 64 years; range, 3878 years) before and after gefitinib treatment were analyzed. A semiautomated three-dimensional lung cancer segmentation algorithm was developed and applied to each tumor at baseline and follow-up. The computer calculated the greatest diameter (unidimensional measurement), the product of the greatest diameter and the greatest perpendicular diameter (bidimensional measurement), and the volume of each tumor. Exact McNemar tests were used to analyze differences in the percentage change calculated with different measurement techniques.
Results: The computer accurately segmented 14 of the 15 tumors. One paramediastinal tumor required manual separation from the mediastinum. Eleven (73%) of the 15 patients had an absolute change in tumor volume of at least 20%, compared with one (7%) and four (27%) patients who had similar changes in unscaled unidimensional (P < .01) and bidimensional (P = .04) tumor measurements, respectively. Seven (47%) patients had an absolute change in tumor volume of at least 30%. In contrast, at unscaled analysis, no patients at unidimensional measurement (P = .02) and two (13%) patients at bidimensional measurement (P = .06) had a change of at least 30%.
Conclusion: Compared with the unidimensional and bidimensional techniques, semiautomated tumor segmentation enabled the identification of a larger number of patients with absolute changes in tumor volume of at least 20% and 30%.
© RSNA, 2006
| INTRODUCTION |
|---|
|
|
|---|
Current measurement techniques are limited by the inability to detect size changes in the axial direction, poor reproducibility due to complex disease morphology, and inter- and intraobserver variability (48). With the advent of thin-section CT, it is now possible to obtain image data sets with a spatial resolution adequate to measure tumor volumes. However, the radiologist faces challenges not only in producing accurate and reproducible measurements but also in performing measurements in a large number of tumors and carrying out efficient follow-up of the lesions to assess response to therapy.
Although isotropic imaging modalities are becoming more available in clinical settings and tumor volume measurements may enable more objective, accurate, and rapid assessment of therapy response than do traditional unidimensional and bidimensional techniques, further validation studies with large patient data sets of different cancer types are necessary. The development of computerized methods, both automated and semiautomated, for reliable quantification of tumor volume is essential to the successful validation of volumetric measurement techniques used to assess therapy responses and obtain routine diagnoses. Thus, the purpose of our study was to prospectively quantify tumor response or progression in patients with lung cancer by using thin-section CT and a semiautomated algorithm to calculate tumor volume and other parameter values.
| MATERIALS AND METHODS |
|---|
|
|
|---|
CT Imaging
Multidetector CT was performed with a LightSpeed 16 scanner (GE Medical Systems, Milwaukee, Wis); the parameters used are given in Table 1. The images were obtained without intravenous contrast material injection during a breath hold of 7 seconds. The entire thorax, from the lower neck and including the adrenal glands in their entirety, was imaged. The acquired thin-section images were downloaded directly from the CT acquisition workstation to our research picture archiving and communication system server by way of the hospital's fast networking system, where the Digital Imaging and Communications in Medicine images were stored with the patient identification information removed (9,10). The images were then transferred to an Ultra 10 workstation (Sun Microsystems, Santa Clara, Calif) for segmentation.
|
The three-dimensional algorithm was modified to adjust itself for segmentation of larger metastases in the lungs and used in two studies related to oncologic therapeutic trials (13,14). For our current study, this three-dimensional algorithm was modified for better measurement of the primary lung cancers seen on multiple thin-section CT images. Instead of using the shape constraint to automatically determine the size of the morphologic filter used to detach adjacent blood vessels, as proposed with the earlier versions of the algorithm, we predefined two sizes of the morphologic filter with the aim of removing surrounding blood vessels up to 5 mm in diameter. To maximally suppress the inconsistency in tumor segmentation that might occur on baseline and follow-up CT images owing to adjacent tissues with an attenuation similar to that of the tumor, tumor pairs were processed simultaneously. For the pleural lesions, the strategy proposed by Zhao et al (15) was adopted: The lungs were extracted by using a thresholding algorithm followed by a morphologic closing operator to repair the distortions on the lung pleural surface caused by the lesions. The lesion segmentation then started in the extracted lungs.
If a large lesion attached to the mediastinum or to the combined right diaphragm and liver had an attenuation similar to the attenuations of these structures, the boundaries lying between the lesion and either the mediastinum or the liver needed manual separation. Thus, a free-hand region of interest on each lesion section would be generated before the segmentation algorithm was applied. Last, nodules with ground-glass opacity were segmented by using a fixed threshold followed by the morphologic opening operation.
One author (B.Z.) used this algorithm to segment all 15 tumors depicted on the CT images obtained before and after gefitinib treatment. After segmentation, the greatest diameter, greatest perpendicular diameter, and volume of each tumor on both images were automatically calculated by the computer (14).
All segmentation results, including the segmented lesions, greatest diameters, and greatest perpendicular diameters, were visually inspected by a board-certified radiologist (L.H.S.) who has 14 years of experience in CT image interpretation. The evaluations were performed by overlapping the contours of each segmented lesion on the original transverse images, which were then viewed carefully in cine mode. To determine the consistency of the three-dimensional algorithm, the segmentation results of each lesion on the baseline and follow-up CT images were displayed side by side. In addition, the two greatest axial diameters of each lesion were placed on the segmented lesion for the assessment. No case required manual modification because of a suboptimal segmentation result.
Statistical Analyses
The unidimensional (RECIST), bidimensional (WHO), and three-dimensional volumetric tumor measurements were compared in two ways. The first way was a direct comparison of the percentage changes in each measurement; we called this the unadjusted/unscaled method. The second way was a comparison of the measurements after conversion of the unidimensional and bidimensional measurements to an equivalent volume by using the standard methods that James et al described with illustrations in the description of the RECIST (16); we called this the adjusted/scaled method. For example, a 30% decrease in the unidimensional measurement correlates geometrically to a 50% decrease in the bidimensional measurement and a 65% decrease in the volumetric measurement. These geometric correlations are determined with the assumption that tumors are spherical and are the basis for the relationships between the unidimensional RECIST, the bidimensional WHO response criteria, and the suggested volumetric response guidelines outlined in the RECIST. Differences in measurements determined at the unadjusted/unscaled comparison reflect clear disparities. In contrast, differences in measurements determined after they are converted into an equivalent volume reflect mathematically pure changes that are due to the fact that tumors are not perfect spheres and do not change in a spherical manner.
Exact McNemar tests were used to test for differences in the percentage change calculated with different techniques. Because of the distribution of data that suggested possible outliers, Wilcoxon sign-rank tests were used to compare the continuous measurements. These analyses were conducted by using Stata 8.0 for Windows, 2003 (Stata, College Station, Tex).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
A prerequisite for obtaining size measurementsespecially volume measurementsat serial thin-section CT is the automated or semiautomated (and reproducible) defining of the tumor's boundaries by means of separating the tumor from the surrounding normal tissues. Owing to the diversity of abnormalities, metastatic disease sites, and surrounding normal anatomic structures, the segmentation of lesions is very challenging; it requires the use of a priori knowledge and the development of task-specific strategies. In this study, we remodeled and applied a segmentation algorithm to measure the sizes of lung cancers in a clinical setting. The modified algorithm is based on an existing algorithm for segmenting small lung nodules depicted on CT images (11,12). The way the algorithm identifies nodule edges was unchanged: The edge is defined as the closed curve surrounding the lesion where the sharpest change in attenuation occurs around the boundary. Therefore, the modified algorithm had the same ability to identify a lesion's edge in regions where there were no abutting vascular structures.
If a lesion is attached to a surrounding soft tissue (eg, blood vessel) with similar attenuation, the soft tissue needs to be removed by using a morphologic filter. Although there is a default morphologic filter size in the modified algorithm, a larger filter can be chosen if the segmentation result is considered suboptimal. Larger morphologic filters were chosen manually in our study, and we anticipate further improvement in the algorithm with use of full automation. A limitation of our study was the lack of an independent reference standard for the tumor volume measurements.
Another factor that can affect volume measurements is the CT scanning technique used. Computer methods may perform differently in the identification of tumor contours on multiple CT images acquired with different parameters such as section thickness, pitch, reconstruction algorithm, peak voltage, and amperage. We specifically designed an imaging protocol for the evaluation of tumor volumes with use of our modified segmentation method. Investigating the effects of CT parameters on the performance of computer methods in the estimation of tumor volumes was beyond the scope of this study.
In our study, the unidimensional, bidimensional, and volumetric measurements of 15 measurable stage I or II nonsmall cell lung cancers depicted on baseline and follow-up CT images obtained 35 weeks apart were estimated by the computer. We selected cutoff values of 20% and 30% on an absolute scale because the RECIST criteria suggest that these are the unidimensional measurement values that should be used to identify tumor progression and response, respectively.
The differences in measurements seen at unadjusted/unscaled analysis in our study were greater than the differences seen at adjusted/scaled analysis. With either comparison method, the greatest changes were seen in the volumetric measurements. When the unadjusted scale was used, the changes seen differed significantly among the three measurement types. When the adjusted scale was used, the changes did not differ significantly among the three measurement types. These findings were not unexpected, because the unadjusted scale maximizes the differences between measurement methods by not correcting for the geometric scales. The differences seen at unadjusted/unscaled analysis reflect absolute differences in the measurement methods. Geometrically, a given percentage of unidimensional change in diameter corresponds to a larger change in volume. However, our unadjusted/unscaled analysis results show that the assumption that tumors are spherical is not always correct and that even when the known geometric relationship is taken into account, the various measurement methods may not yield equivalent results. Tumors do not necessarily have a spherical shape or change in a spherical manner.
Given the small sample size of 15 patients, differences between the volumetric estimations and the unidimensional and bidimensional measurements that were not significant at adjusted/scaled analysis were probably due to low statistical power. Further analysis of larger data sets is needed. Another limitation of our study is the fact that we do not yet have long-term follow-up data to assess whether patient benefitmeasured in terms of survival or time to progressioncorrelates better with one scale or the other. Our study results demonstrate that significantly greater changes in tumor size may be accurately measured at volumetric response assessment in the unadjusted scale and that a similar trend is seen in the adjusted scale.
Our initial study results suggest that with use of thin-section CT images and computer software, changes in tumor volume can be assessed as early as 3 weeks after initiation of gefitinib treatment, whereas a lower magnitude of changes in unidimensional and bidimensional measurements is seen during the same time period. Thus, volumetric changes in tumor size may have the potential to be an earlier or better marker of regression or progression. Such a marker, by enabling a lack of tumor response to be recognized earlier or more reliably than it is recognized with conventional unidimensional and bidimensional measurements, could help to limit the amount of ineffective chemotherapy patients receive. As more therapeutic agents and approaches become available, it will become increasingly important to reevaluate existing methods and investigate technologies for determining lung cancer response at imaging.
| ADVANCE IN KNOWLEDGE |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Abbreviations: RECIST = Response Evaluation Criteria in Solid Tumors WHO = World Health Organization
Authors stated no financial relationship to disclose.
Author contributions: Guarantors of integrity of entire study, B.Z., L.H.S.; study concepts/study design or data acquisition or data analysis/interpretation, all authors; manuscript drafting or manuscript revision for important intellectual content, all authors; manuscript final version approval, all authors; literature research, B.Z.; clinical studies, L.H.S., M.G., N.A.R., M.G.K.; experimental studies, B.Z., L.H.S.; statistical analysis, L.H.S., C.M.; and manuscript editing, B.Z., L.H.S., C.M., N.A.R., M.G.K.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Zhao, L. P. James, C. S. Moskowitz, P. Guo, M. S. Ginsberg, R. A. Lefkowitz, Y. Qin, G. J. Riely, M. G. Kris, and L. H. Schwartz Evaluating Variability in Tumor Measurements from Same-day Repeat CT Scans of Patients with Non-Small Cell Lung Cancer Radiology, July 1, 2009; 252(1): 263 - 272. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Gavrielides, L. M. Kinnard, K. J. Myers, and N. Petrick Noncalcified Lung Nodules: Volumetric Assessment with Thoracic CT Radiology, April 1, 2009; 251(1): 26 - 37. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Zhao, L. H. Schwartz, and S. M. Larson Imaging Surrogates of Tumor Response to Therapy: Anatomic and Functional Biomarkers J. Nucl. Med., February 1, 2009; 50(2): 239 - 249. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. M. Stadler Tumor Burden Endpoints and Phase II Clinical Trial Design ASCO Educational Book, January 1, 2008; 2008(1): 89 - 93. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. J. Riely, M. G. Kris, B. Zhao, T. Akhurst, D. T. Milton, E. Moore, L. Tyson, W. Pao, N. A. Rizvi, L. H. Schwartz, et al. Prospective Assessment of Discontinuation and Reinitiation of Erlotinib or Gefitinib in Patients with Acquired Resistance to Erlotinib or Gefitinib Followed by the Addition of Everolimus Clin. Cancer Res., September 1, 2007; 13(17): 5150 - 5155. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |