Adaptive Multiple Feature Method (AMFM) for the early detection of parenchymal pathology in a smoking population

Renuka Uppaluri¹, Geoffrey McLennan², Paul Enright³,
James R. Standen⁴, Pamela Boyer-Pfersdorf³, and Eric A. Hoffman^1,5

¹Department of Radiology University of Iowa, Iowa City, IA 52242

²Department of Internal Medicine University of Iowa, Iowa City, IA 52242

³Department of Internal Medicine University of Arizona, Tucson, AZ 85724

⁴Department of Radiology University of Arizona, Tucson, AZ 85724

⁵Department of Biomedical Engineering University of Iowa, Iowa City, IA 52242

Table of Contents

• ABSTRACT
• INTRODUCTION
• METHODS
• EXPERIMENTAL METHODS
  • Subjects and Image Acquisition
  • Measurements and Validation
• RESULTS
  • Global analysis
  • Regional analysis (observer-defined regions)
  • Regional analysis (whole lung slices)
• DISCUSSION
• CONCLUSION
• ACKNOWLEDGEMENTS
• REFERENCES

Abstract:

Tissue characterization, texture, smoking, pulmonary, Bayesian classifier

INTRODUCTION

Smoking is one of the main causes of debilitating diseases such as emphysema. Emphysema is a form of chronic obstructive pulmonary disorder (COPD). The onset of emphysema starts early but it presents itself clinically at a stage when it is incurable and care is largely supportive. In recent years, lung volume reduction surgery has become a viable intervention for end-stage emphysema subjects. However, patient selection criteria are poorly defined and, more importantly, sensitive measures for outcome evaluation are lacking. Novel anti-elastase and other anti-inflammatory modulators are being developed to halt the progression of early emphysema. Therefore, the motivation for this study is to identify a means of detecting changes caused by smoking before clinical symptoms appear. This can then be used either to advise the patient to adopt lifestyle changes such as smoking cessation or to determine the outcome of therapeutic interventions. It should be noted, once established, emphysema may progress even with smoking cessation.

In this study, three groups of subjects are included - healthy non-smokers with normal lung function, healthy smokers with normal lung function, and smokers with abnormal lung function or COPD. The AMFM will be used to demonstrate the differences between these three groups.

METHODS

The optimal feature selection is performed by using the divergence measure along with correlation analysis [4] and the classification is performed by a Bayesian approach. The texture measures used are derived from various classes. The first order measures based on the grey level histogram, run-length measures, the co-occurrence matrix measures, and the geometric fractal dimension computations are presented in [5,6,1]. In addition to these 17, 10 new measures, developed by us, are included in this study. They are the percentile measures, also based on the grey level histogram, and the stochastic fractal dimension measures [7]. All measures used and their abbreviations are provided in Table 1.

EXPERIMENTAL METHODS

Subjects and Image Acquisition

The subjects included in this study were: 41 healthy non-smokers (NONSMK), 24 smokers (SMK), and 19 smokers with abnormal lung function (COPD).

A Siemens DRH Somatom scanner was used to acquire CT scans of all subjects in supine position. Two slices (1 mm slice thickness) were acquired at full inspiration 5 cm above and 5 cm below the carina. The images were reconstructed to 512512 pixels and the grey levels were rescaled to lie between 0-2048.

Measurements and Validation

1.

Global analysis:

All three subject groups were used in this experiment. The left and right lungs from each image were isolated and used individually as the region of interest for texture computation and classification. Half the lungs were put into a training set and the remaining were put into a test set. The goal was to discriminate a NONSMK lung from SMK lung, a NONSMK lung from a COPD lung, and a SMK lung from a COPD lung. The ability to discriminate simultaneously between a NONSMK lung, SMK lung, and COPD lung was also tested.

On the test set, classification accuracy of each subject group and the overall accuracy of all subject groups was computed for each discrimination task. The classification accuracy is the percent of samples belonging to a group correctly classified as such. The overall accuracy is the percent of correctly classified samples of all groups together.

2.

Regional analysis (observer-defined regions):

In this experiment, NONSMK and COPD subject groups were used. A trained observer identified areas of four tissue types, namely, normal lung from NONSMK subjects (NN), emphysema-like lung from COPD subjects (EC), ``normal" lung from COPD subjects (NC), and mild emphysema-like lung from COPD subjects (MC). In each defined region, pixel blocks of 3131 overlapping by 1515 were placed. Half the blocks were put into a training set and the remaining were put into a test set. The goal was to discriminate NN samples from each of EC, NC, and MC samples, individually. The ability to discriminate between NN from EC, NC, MC simultaneously was also tested.

As before, classification accuracies and overall accuracies were computed.

3.

Regional analysis (whole lung slices):

All three subject groups were included in this experiment. On each CT slice, 3131 pixel blocks overlapping by 1515 were translated through the lung regions and the pattern classification obtained for each block was assigned to the 1515 block at its center. The classifier assigned either a NN, EC, NC, or MC pattern using the training acquired in the previous experiment. Each block was color-coded based on its pattern and the images were displayed. The goal of this experiment was to determine what percent of a NONSMK, SMK, and COPD lung is comprised of each of the four patterns.

RESULTS

Global analysis

Regional analysis (observer-defined regions)

Regional analysis (whole lung slices)

  

Figure 1: Examples of lungs with their color-coded computer output. (a) An healthy non-smoker's lung (NONSMK) (b) Computer output of NONSMK lung (c) An smoker's lung who has normal lung function (SMK) (d) Computer output of SMK lung (e) Lung of a smoker with abnormal lung function (COPD) (f) Computer output of the COPD lung. NN: Normal lung region from healthy non-smokers, EC: Emphysematous lung region from smoking COPD subjects, NC: Normal appearing lung region from smoking COPD subjects, MC: Mild emphysematous lung region from smoking COPD subjects.

DISCUSSION

The AMFM has demonstrated to be successful in differentiating normal subjects from subjects with severe emphysema [1], and more recently with interstitial lung diseases such as idiopathic pulmonary fibrosis, and sarcoidosis [2]. Since, in this study, NONSMK (non-smokers with normal lung function) can be differentiated from SMK (smokers with normal lung function) with 80.0% accuracy, this further demonstrates that the AMFM is sensitive enough to pick up subtle differences. Our regional analysis experiments were used to localize and characterize the differences between the groups. It was first demonstrated that visually normal appearing lung regions from NONSMK and COPD subjects (smokers with abnormal lung function) were differentiated with 89.5% accuracy. An explanation for this could be that normal appearing lung in COPD subjects appears ``normal" to the observer only because it is juxtaposed with severely diseased lung. The mild emphysematous regions and emphysema regions were also discriminated from normal lung of NONSMK with accuracies higher than 94.0%.

In analyzing whole lung slices, it was found that the majority of a SMK lung was composed of mild emphysema and ``normal from COPD", i.e., lung which visually appears normal but is different from normal lung of non-smokers. It is also noted that 22.9% of a NONSMK lung is composed of mild emphysema. We attribute this to the presence of dark regions juxtaposed to ribs which may be occurring due to beam hardening. Since such areas were not included in the training set as being normal, the computer assigned labels of mild emphysema or emphysema. The computer can perform only as well as the training it receives. Therefore, representative training examples along with good data acquisition using standardized protocols are of utmost importance. Futhermore, it is critical that motion be eliminated during scanning and that lung volume be standardized.

This study also demonstrates the applicability of the AMFM across scanners. While our previous studies were performed on Imatron C-150 scanner with 3 mm thick slices, this study uses 1 mm thick slices from Siemens Somatom DRH scanner. It should be noted that, for best results, the AMFM should be re-trained for each scanner and for each new set of scanning and reconstruction parameter.

CONCLUSION

ACKNOWLEDGEMENTS

References

1

R. Uppaluri, T. Mitsa, M. Sonka, E. A. Hoffman, and G. McLennan, ``Quantification of pulmonary emphysema from lung computed tomography images,'' Am. J. Respir. Crit. Care Med. 156(1), pp. 248-254, 1997.

2

R. Uppaluri, M. Sonka, E. A. Hoffman, T. Mitsa, C. Dayton, G. W. Hunninghake, and G. McLennan, ``Quantitative assessment in interstitial lung disease using the multiple feature (MF) method (abstract),'' Am. J. Respir. Crit. Care Med. 155(4), p. A322, 1997.

3

R. Uppaluri, P. G. Hartley, E. A. Hoffman, M. Sonka, and G. McLennan, ``Validation of AMFM: A quantitative, regional assessment of pulmonary parenchymal pathology from HRCT (abstract),'' Radiology 205((P)), p. 481, 1997.

4

H. C. Andrews, Introduction to mathematical techniques in pattern recognition, Wiley, New York, 1972.

5

R. Uppaluri, T. Mitsa, and J. R. Galvin, ``Fractal analysis of high-resolution CT images as a tool for quantification of lung diseases,'' in Medical Imaging 1995: Physiology and Function from Multidimensional Images, Vol. 2433, San Diego, CA, pp. 133-142, SPIE, (Bellingham, WA), 1995.

6

R. Uppaluri, T. Mitsa, E. A. Hoffman, G. McLennan, and M. Sonka, ``Texture analysis of pulmonary parenchyma in normal and emphysematous lung,'' in Medical Imaging 1996: Physiology and Function from Multidimensional Images, Vol. 2709, Newport Beach, CA, pp. 456-467, SPIE, (Bellingham, WA), 1996.

7

R. Uppaluri, Automated analysis of pulmonary parenchyma from CT images.
PhD thesis, University of Iowa, Iowa City, Iowa, USA, 1997.