The first step in lung analysis by CT is the identification of the lung border. To deal with the increased number of sections per scan in thin-slice multidetector CT, it has been crucial to develop accurate and automated lung segmentation algorithms. In this study, an automated method for lung segmentation of thin-slice CT data is presented. The method exploits the advantages of a two-dimensional wavelet edge-highlighting step in lung border delineation. Lung volume segmentation is achieved with three-dimensional (3D) grey level thresholding, using a minimum error technique. 3D thresholding, combined with the wavelet pre-processing step, successfully deals with lung border segmentation challenges, such as anterior or posterior junction lines and juxtapleural nodules. Finally, to deal with mediastinum border under-segmentation, 3D morphological closing with a spherical structural element is applied. The performance of the proposed method is quantitatively assessed on a dataset originating from the Lung Imaging Database Consortium (LIDC) by comparing automatically derived borders with the manually traced ones. Segmentation performance, averaged over left and right lung volumes, for lung volume overlap is 0.983+/-0.008, whereas for shape differentiation in terms of mean distance it is 0.770+/-0.251 mm (root mean square distance is 0.520+/-0.008 mm; maximum distance is 3.327+/-1.637 mm). The effect of the wavelet pre-processing step was assessed by comparing the proposed method with the 3D thresholding technique (applied on original volume data). This yielded statistically significant differences for all segmentation metrics (p<0.01). Results demonstrate an accurate method that could be used as a first step in computer lung analysis by CT.

Identification and characterization of diffuse parenchyma lung disease (DPLD) patterns challenges computer-aided schemes in computed tomography (CT) lung analysis. In this study, an automated scheme for volumetric quantification of interstitial pneumonia (IP) patterns, a subset of DPLD, is presented, utilizing a multidetector CT (MDCT) dataset. Initially, lung-field segmentation is achieved by 3-D automated gray-level thresholding combined with an edge-highlighting wavelet preprocessing step, followed by a texture-based border refinement step. The vessel tree volume is identified and removed from lung field, resulting in lung parenchyma (LP) volume. Following, identification and characterization of IP patterns is formulated as a three-class pattern classification of LP into normal, ground glass, and reticular patterns, by means of k-nearest neighbor voxel classification, exploiting 3-D cooccurrence features. Performance of the proposed scheme in indentifying and characterizing ground glass and reticular patterns was evaluated by means of volume overlap (ground glass: 0.734 +/- 0.057, reticular: 0.815 +/- 0.037), true-positive fraction (ground glass: 0.638 +/- 0.055, reticular: 0.942 +/- 0.023) and false-positive fraction (ground glass: 0.361 +/- 0.027, reticular: 0.147 +/- 0.032) on five MDCT scans.

Accurate and automated lung field (LF) segmentation in high-resolution computed tomography (HRCT) is highly challenged by the presence of pathologies affecting lung borders, also affecting the performance of computer-aided diagnosis (CAD) schemes. In this work, a two-dimensional LF segmentation algorithm adapted to interstitial pneumonia (IP) patterns is presented. The algorithm employs k-means clustering followed by a filling operation to obtain an initial LF order estimate. The final LF border is obtained by an iterative support vector machine neighborhood labeling of border pixels based on gray level and wavelet coefficient statistics features. A second feature set based on gray level averaging and gradient features was also investigated to evaluate its effect on segmentation performance of the proposed method. The proposed method is evaluated on a dataset of 22 HRCT cases spanning a range of IP patterns such as ground glass, reticular, and honeycombing. The accuracy of the method is assessed using area overlap and shape differentiation metrics (d(mean), d(rms), and d(max)), by comparing automatically derived lung borders to manually traced ones, and further compared to a gray level thresholding-based (GLT-based) method. Accuracy of the methods evaluated is also compared to interobserver variability. The proposed method incorporating gray level and wavelet coefficient statistics demonstrated the highest segmentation accuracy, averaged over left and right LFs (overlap=0.954, d(mean)=1.080 mm, d(rms)=1.407 mm, and d(max)=4.944 mm), which is statistically significant (two-tailed student's t test for paired data, p<0.0083) with respect to all metrics considered as compared to the proposed method incorporating gray level averaging and gradient features (overlap=0.918, d(mean)=2.354 mm, d(rms)=3.711 mm, and d(max)=14.412 mm) and the GLT-based method (overlap=0.897, d(mean)=3.618 mm, d(rms)=5.007 mm, and d(max)=16.893 mm). The performance of the three segmentation methods, although decreased as IP pattern severity level (mild, moderate, and severe) was increased, did not demonstrate statistically significant difference (two-tailed student's t test for unpaired data, p>0.0167 for all metrics considered). Finally, the accuracy of the proposed method, based on gray level and wavelet coefficient statistics ranges within interobserver variability. The proposed segmentation method could be used as an initial stage of a CAD scheme for IP patterns.

No abstract

The automated segmentation of vessel tree structures is a crucial preprocessing stage in computer aided diagnosis (CAD) schemes of interstitial lung disease (ILD) patterns in multidetector computed tomography (MDCT). The accuracy of such preprocessing stages is expected to influence the accuracy of lung CAD schemes. Although algorithms aimed at improving the accuracy of lung fields segmentation in presence of ILD have been reported, the corresponding vessel tree segmentation stage is under-researched. Furthermore, previously reported vessel tree segmentation methods have only dealt with normal lung parenchyma. In this paper, an automated vessel tree segmentation scheme is proposed, adapted to the presence of pathologies affecting lung parenchyma. The first stage of the method accounts for a recently proposed method utilizing a 3-D multiscale vessel enhancement filter based on eigenvalue analysis of the Hessian matrix and on unsupervised segmentation. The second stage of the method is a texture-based voxel classification refinement to correct possible over-segmentation. The performance of the proposed scheme, and of the previously reported technique, in vessel tree segmentation was evaluated by means of area overlap (previously reported: 0.715 ± 0.082, proposed: 0.931 ± 0.027), true positive fraction (previously reported: 0.968 ± 0.019, proposed: 0.935 ± 0.036) and false positive fraction (previously reported: 0.400 ± 0.181, proposed: 0.074 ± 0.031) on a dataset of 210 axial slices originating from seven ILD affected patient scans (used for performance evaluation out of 15). The proposed method demonstrated a statistically significantly ( p < 0.05) higher performance as compared to the previously reported vessel tree segmentation technique. The impact of suboptimal vessel tree segmentation in a reticular pattern quantification system is also demonstrated.