Next Article in Journal
Study of Subjective Data Integrity for Image Quality Data Sets with Consumer Camera Content
Previous Article in Journal
Unveiling the Secrets of Escher’s Lithographs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lung Nodule Detection in CT Images Using Statistical and Shape-Based Features

1
Punjab University College of Information Technology, University of the Punjab, Lahore-54000, Pakistan
2
Department of Radiology, General Hospital, Lahore-54000, Pakistan
*
Author to whom correspondence should be addressed.
J. Imaging 2020, 6(2), 6; https://doi.org/10.3390/jimaging6020006
Submission received: 3 February 2020 / Revised: 20 February 2020 / Accepted: 21 February 2020 / Published: 24 February 2020

Abstract

:
The lung tumor is among the most detrimental kinds of malignancy. It has a high occurrence rate and a high death rate, as it is frequently diagnosed at the later stages. Computed Tomography (CT) scans are broadly used to distinguish the disease; computer aided systems are being created to analyze the ailment at prior stages productively. In this paper, we present a fully automatic framework for nodule detection from CT images of lungs. A histogram of the grayscale CT image is computed to automatically isolate the lung locale from the foundation. The results are refined using morphological operators. The internal structures are then extracted from the parenchyma. A threshold-based technique is proposed to separate the candidate nodules from other structures, e.g., bronchioles and blood vessels. Different statistical and shape-based features are extracted for these nodule candidates to form nodule feature vectors which are classified using support vector machines. The proposed method is evaluated on a large lungs CT dataset collected from the Lung Image Database Consortium (LIDC). The proposed method achieved excellent results compared to similar existing methods; it achieves a sensitivity rate of 93.75%, which demonstrates its effectiveness.

1. Introduction

Lung cancer has a high causality rate. According to a survey [1], more than 1.37 million people died from Lung cancer throughout the world only in 2008. The American Cancer Society has projected 1.74 million new cancer cases and 0.61 million cancer deaths in the year 2018 [2]. Two main reasons for the high mortality rate in lung cancer are the delay in early diagnosis and the poor prognosis [3]. The study reveals that 70% of lung cancers are diagnosed in too advanced stages, where the cancer prognosis is ineffective. Therefore, early diagnosis of cancer is momentous for increasing the patient’s chances of survival.
Computed Tomography (CT) scan images are utilized for cancer diagnosis; the radiologists examine them to detect and classify the nodules into malign and benign [4]. However, these methods require highly skilled radiologists who are not available, particularly, to the people of remote and poor regions. Moreover, there is a high risk of human error in manual examination, and thus Computer Aided Detection (CAD) systems are needed that can help the radiologists in the diagnosis and help decrease the rate of false reports. Digital image processing techniques can be used to detect the nodules, their type, size, and other features from CT scans.
Medical image processing has been extensively and increasingly applied to design expert support systems for the diagnosis of numerous diseases, e.g., arthritis detection [5,6], parasite detection [7,8,9], lung cancer detection [10,11,12], and rehabilitation [13,14,15]. A significant amount of research is being done in the early lung cancer diagnosis using CAD systems [16,17,18,19,20]. The need for automated systems arises because of the nature of the data—CT scans—used in lung cancer diagnosis. A lung CT scan contains usually more than 250 images per scan. Examining this extensive dataset for each patient is quite a challenging, time-consuming, and a tedious task for a radiologist. Moreover, the nature of nodules which decides the fate of a patient is also very complicated, as their shape and size varies from slice to slice. Sometimes they are attached to other pulmonary structures, such as vessels or bronchioles. The color in which they appear on CT scans may also differ. These factors add to the complexity of identifying them. However, the same elements, once recognized, help researchers in defining the course of their methodology.
A fundamental step in lung nodules detection is the accurate lung segmentation from the CT image. For this purpose, numerous techniques have been proposed for efficient lung segmentation. Some algorithms require a few seed pixels on the lungs region in the image and then utilize region-growing techniques to segment the lungs e.g., [21,22,23,24]. The algorithm in [21] extracts the chest out of complete CT images using a region-growing algorithm taking four seed points as input. The same technique was used to segment lungs from the chest. Cascio et al. [23] also use a region-growing algorithm for lung segmentation. They use the center pixel of the slice as seed. The method in [22] segments the nodules with the assistance of the radiologists. The technique proposed in [25] divides the chest CT scan into four classes: lung wall, parenchyma, bronchioles, and nodules. The active contour technique is used to segment lungs from the CT images.
Numerous methods in the literature use nodule intensity or color thresholding to detect nodules in CT images, e.g., [26,27,28,29,30,31,32]. Clifford et al. [26] preprocessed the CT image with bi-histogram equalization to improve its sharpness. The resultant image is thresholded to obtain connected components. The area and pixel values of these segmented regions are passed as features to fuzzy inference system for nodule detection. Another rule-based technique for lung segmentation is proposed in [27]. A nodule size is used in [28] as the standard to detect potential nodule regions. Local maxima are found for each subvolume with values and size larger than the standard nodule in 3D space. Messay et al. [29] use thresholding for initial segmentation, and then a rule-based analysis on the anatomical characteristics of bronchioles present between lungs. It detects and segments nodule candidates simultaneously, unlike the methods discussed previously. The algorithms in [30,31] are also thresholding-based methods for lung nodule classification.
Template matching has also been explored for nodule detection, e.g., [33,34,35]. In such methods, the template nodule images are searched in the target image in order to find the nodules. Ozekes et al. [33] developed a 3D template for the detection of nodule candidates. The template is convoluted to the Region Of Interests (ROIs) and extracted by applying an eight-directional search on the lungs region. In [36], thresholding and binarization of the CT image are used for lung segmentation. A multi-scale filtering is used for the detection of nodules.
Nodules have been shown to possess certain properties, such as shape, color, and intensity, which have been used to design a composite discriminative feature. The features are then classified to detect the nodule through different classification methods. The methods in [37,38,39,40,41] are a few such examples. The algorithm proposed in [37] extracts six features from each slice and forms feature vectors, which are classified using a support vector machine to identify the nodules. A two-dimensional (2D) multiscale filter is proposed in [39] for the detection of lung nodule candidates. The shape features for each segmented region are used to reduce the false positive rate through a classifier. The method proposed in [41] uses surface triangulation on different threshold values for the detection of nodule candidates and neural network classifiers are used to separate the nodules from the non-nodules. Deep learning techniques have also been explored lately for lung nodule detection, e.g., [42,43,44]. The research in [45,46,47] presents the recent developments in automatic lung nodule detection.
In this paper, we propose a CAD system that automatically extracts the lungs from chest CT scans and processes these segments to detect nodules. The major contributions of this paper are as follows:
  • In most existing CAD-based nodule detectors, the lungs are manually marked by the radiologist, which is a tedious and time-consuming task. In the proposed algorithm, the lungs are automatically segmented from the CT images without any user intervention;
  • Nodules can have different regular and irregular shapes and sizes. Some existing techniques use a few shape templates to detect the nodules, however the proposed algorithm is independent of the nodule shape and size;
  • The proposed system uses basic image processing techniques, e.g., histogram processing, morphological operators, connected components analysis, etc., which makes it implementable on simple computers, making it an efficient and cost effective solution;
  • In an experimental evaluation carried out on a standard LIDC dataset, the proposed system achieved high sensitivity and accuracy, outperforming the existing similar techniques.

2. Materials and Methods

The strategy proposed for lung nodule detection comprises two major phases. In the first phase, lung are segmented from the Digital Imaging and Communications in Medicine (DICOM) CT scans, and in the second phase nodules are detected from the lungs. In most existing techniques, the lungs segmentation is performed in semi-automatic ways, where the radiologist assists the system by either specifying the region of interest or by drawing a few scribbles on the target object. Subsequently, based on the user input, the lungs are extracted from the rest of the image. In the proposed framework, the lung segmentation is totally automatic and no assistance from the user is needed. The proposed segmentation technique constructs a histogram of the given image and analyzes it to automatically select a threshold. Based on this threshold the outer region in the image is identified and dropped, and the rest of the region is further processed to extract the lungs. This process involves morphological operations and connected components analysis.
In the second phase, the nodules are detected from the segmented lungs. It is achieved by separating the inner structures of nodules, bronchi, and blood vessels from the parenchyma region. A model is trained using the statistical and shape-based features of the nodules, and a support vector machine-based classifier is used to separate the nodules from the bronchi and the blood vessels. A block diagram of the proposed algorithm is shown in Figure 1.

2.1. Lung Segmentation

The medical images are not in the conventional image formats such as PNG or JPEG, as they are conducted under a specified constrained environment, which has a direct impact on the image attained. Digital Imaging and Communications in Medicine (DICOM) is a standard for storing and transmitting digital images, enabling the integration of different medical imaging devices, e.g., scanners, servers, workstations, and printers. A DICOM image may also contain the information of the patient, date, and many other data that are not required in lung segmentation. To conveniently perform the image processing tasks, we convert the DICOM images to the loss-less Portable Network Graphics (PNG) format. When a DICOM image is converted to PNG format, the personal information of the patient and all tags which come with a DICOM format get removed so as to preserve the individual’s privacy.
The lung segmentation is considered a fundamental activity in nodule CAD systems, as the performance of the later stages in such an analysis largely depends on the segmentation accuracy. In this section, we propose a lung segmentation algorithm that utilizes the histogram and morphological image processing techniques. The converted PNG image has four components: (i) a black background; (ii) a dark gray circular region; (iii) a brighter region; and (iv) the lungs in a dark gray shade, as shown in Figure 2. Our region of interest is the lungs, thus here we remove the first two components.
To reduce the image to our ROI, we perform thresholding. Thresholding is mainly dependent on the value of the threshold, and this value is usually user specified. In our case, we compute the threshold for each slice from the histogram of the input image. Figure 3b shows the histogram of a sample lung image shown in Figure 3a. It can be noted that the histogram has four prominent peaks. One very high peak is at 0, which corresponds to the black background in the image. The second peak around the gray level 60 is formed from the dark gray circular region covering the bright region. These two peaks correspond to the regions (i) and (ii) as discussed above. A high peak at 255 corresponds to the white region mainly around the lungs, and the third peak around gray level of 210 in this example is formed by the intensities of the lungs and patches inside the bright region around the lungs. Thus, by dropping the pixels that fall in the first two peaks, we can remove the background regions (i) and (ii). The valley where the second peak ends serves as a separator between regions i–ii and regions iii–iv. This valley can be estimated by computing the second minima of the histogram. The value of the second minima is used as a threshold to remove the regions (i) and (ii) from the CT scan. Figure 3c shows the result achieved after thresholding the image (Figure 3a) with the estimated threshold.
Let I be an input lung image of size M × N after pre-processing. A histogram H of I is computed using the step size κ , and the 2nd local minima λ of H is computed and used in thresholding the image I to remove the background from the image.
I ^ ( x , y ) = I ( x , y ) if I ( x , y ) > λ , 0 otherwise
where I ^ is the thresholded image. Conventionally, the white color is used to represent the foreground and the black color is used for the background; we then complement the resultant image. Figure 3d shows the complemented image I ¯ obtained from I ^ . The next step is to separate the lung region from its surrounding bright region. The histogram-based thresholding applied to remove the background regions is not effective in this case, as the bright region covering the lungs contains patches of significantly different intensities (Figure 3c). For this purpose, we use the Otsu method [48] to separate the lungs from the surrounding region. Let τ be the threshold obtained from [48], which is used to obtain the binary mask B of lungs:
B ( x , y ) = 1 if I ¯ ( x , y ) > τ , 0 otherwise
The resultant binary mask is shown in Figure 3e. It can be observed that the mask still contains a few unwanted objects. If we see this image as components of pixels, it can be concluded that our region of interest is the two separate components enclosed in the largest component. To this end, we compute the connected components [49] of the mask to achieve our goal of getting a binary map of the lungs.
We simply determine the largest component that corresponds to the region enclosing the lungs, as evident from Figure 3e. Hence after determining this component, we shred it off to obtain the lungs mask M. Moreover, a morphological dilation operation is also applied to the extracted map to include the false negatives inside the lungs. The resultant map is shown in Figure 3f. Using this map, we segment the lungs out of the original image I, shown in Figure 3g. Salt and pepper noise can be noted in the mask (Figure 3f) and in the final segmented lungs (Figure 3g). We apply a median filter of size 3 × 3 to remove this noise. The final segmented lungs are shown in Figure 3h.

2.2. Nodules Detection

In this phase, first the inner structures of the lungs, i.e., nodules, bronchi, and blood vessels, are separated from the parenchyma region. The inner structures in the lungs appear as bright spots (Figure 3h), which can be easily separated through thresholding as the intensity levels of the parenchyma region and inner structures are distinguishably different. The Otsu method [48] is used on the segmented lungs to separate the inner structure vessels, bronchi, and nodules (if there are any) from the rest of the region.
The nodules differ from other structures present in lungs in many aspects. One key difference is their shape, and we exploit this property to isolate the nodules from the non-nodules structures, i.e., vessels and bronchi. The nodules are spherical in shape, whereas the vessels and the bronchioles are cylindrical, as shown in Figure 4. We use the size invariant round/near-round shape detection algorithm proposed in [50] to identify the circular shapes in the detected set of structures. It returns the centers of the potential nodule locations, which are used as seed points in a region-growing algorithm [51] to extract the nodule candidates. In contrast to existing nodule template-based techniques, the proposed strategy enables us to extract nodules of any shape, making our methodology independent of any nodule template.
Let { c 1 , c 2 , c 3 , , c n } be the centers of the shapes extracted using algorithm [50]. The centers are passed to the region-growing algorithm as seeds, which returns us the corresponding n nodule candidate regions { A 1 , A 2 , A 3 , , A n } . We compute different features of nodules and construct a feature vector to discriminate nodules from the other inner structures. We exploit different statistical properties, shape-based features, and across-slice characteristics of the candidate regions to design a discriminative feature vector. In particular we use the following statistical and shape properties of candidates regions:
  • Mean ( μ i ) represents the average value of the region A i :
    μ i = 1 n x A i x ;
  • Median ( m e i ) is the mid-point of A i when arranged in non-decreasing order;
  • Mode ( m o i ) is the most repetitive element of the data in A i ;
  • Variance ( σ 2 ) represents to what extent the data varies from the mean value. For region A i , σ i 2 is:
    σ i 2 = 1 A i x A i ( x μ i ) 2
    where A i and μ i represent the size and mean of the region A i , respectively;
  • Standard deviation σ is the square root of variance:
    σ i = σ i 2 ;
  • Consistency feature: one more important feature is based on the shape of the lesion and its appearance in the colocated slices of the CT scan. That is, if a nodule exists in one slice, it must also appear in the preceding slices or in the succeeding slices of a CT scan. On the other hand, the vessels and bronchi transform further into new shapes, so if they are detected in one slice there are high chances that they will not be present in the exact location in the next slice of the series. This property is an important feature of nodules. Therefore, the center points detected in slice S j are traced in a window of size 2 k + 1 in adjacent slices. That is, the center points detected in slice S i will be compared with the center points identified in k previous and k next slices,
    S j k , , S j 2 , S j 1 , S j , S j + 1 , , S j + k Slice Search Window
    and we assign a center point 1 if it exists in any of those 2 k slices, and 0 if it does not exist in any of them. This consistency feature for candidate A i in the current slice S j is represented as t i . In this paper, we use a window of size 3, ( k = 1 ), to determine the the value of feature t i for each location A i in a given slice.
Figure 5 shows an example of computing the consistency feature for a sample slide S j . The preceding slide is S j 1 and the next slide is S j + 1 . The three slices are shown in column (a), the lungs detected using the proposed algorithm are shown in (b), and the inner structures identified by the proposed technique are shown in (c). The centers of the round/near-round regions detected by our method are highlighted using small red circles in (c). These centers are passed to the region-growing algorithm to mark the shape of the objects. The spherical objects detected in each slice are shown in (d) separately. One potential nodule is detected in slice S j 1 , four candidates are detected in slice S j , and three candidates are detected in slice S j + 1 , as is shown in column (d) of Figure 5. It can be noted that if only one candidate of slice S j is found in the previous and also in the next slices of S j , it is present in the last figure (from left-to-right) of each row. This object of slice S j is taken as a nodule candidate, the other three regions which could not be traced either in S j 1 or in S j + 1 are marked as vessels or bronchi and are dropped from further processing.
All the features described above are computed for each nodule candidate region A i in a slice S j and then combined to obtain a feature vector F i ,
F i = μ i , m e i , m o i , σ i 2 , σ i , t i .
Based on feature F, the selected regions are classified as nodules and non-nodules using a Support Vector Machine (SVM). Training is done by using features of a training dataset, and then these same features are calculated for the testing dataset and passed to the model for classification. We used the LIBLINEAR SVM library [52] in our implementation for SVM classification.

3. Experiments and Results

To assess the performance of the proposed algorithm, we performed a large set of experiments on a standard lung CT dataset. The performance is objectively computed and the results are also compared with 10 existing similar techniques, including [11,21,29,30,31,36,53,54,55,56]. In the region-growing algorithm, the maximum intensity distance was set to 0.18. In all experiments, κ was set to 5 and the window of size 3 was used in the computation of the nodule consistency feature vector. Figure 6 presents a few more results of nodule detection achieved by our method on images from the test dataset.

3.1. Evaluation Dataset

The dataset of lungs CT scans was collected from the Lung Image Database Consortium (LIDC) database [57]. This collection contained 70 cases of Lung scan acquired using different CT scanners. Four radiologists tagged these scans and the tagging was done in two phases. In the first phase, each radiologist tagged the scans independently, and in next phase, results from all radiologists were compiled together and then given to each radiologist for a second review. In this phase, all radiologist were able to review their previous annotations as well as annotations done by the other radiologists. Each of these 70 cases is a series of 250–350 images.

3.2. Performance Evaluation

In this section, we evaluate the performance of the proposed lung nodule detection algorithm using different statistical metrics. There are four possible outcomes of the proposed algorithm when run on a test image: true positive (TP), true negative (TN), false positive (FP), and false negative (FN). True positives (TP) means that nodules exist in the image and they are detected correctly, and true negative (TN) means that there are no nodules in the images and this is correctly identified. False positive (FP) occurs when no nodule exists in the image but is incorrectly detected by the algorithm, and false negative (FN) occurs when a nodule is missed by the algorithm.
To objectively quantify the performance of the proposed algorithm, we chose various statistical measurement parameters: sensitivity, specificity, precision, accuracy, and F score. Sensitivity, also called Recall, is a measurement of true positive rate, i.e., a nodule is tagged as the nodule:
S e n s i t i v i t y or R e c a l l = T P T P + F N .
Specificity measures the true negative rate, i.e., non-nodule is tagged as non-nodule. It is computed as:
S p e c i f i c i t y = T N T N + F P .
Precision demonstrates how much is the algorithm precise in detecting true positive results:
P r e c i s i o n = T P T P + F P .
Accuracy is the measurement of how well the binary classifier correctly identifies or excludes a condition:
A c c u r a c y = T P + T N T P + T N + F P + F N .
When the positive and negative classes in the binary classification are highly unbalanced, the accuracy and precision metrics can be delusive [58], and F measure is considered to be more reliable in such situations. The F measure is a weighted harmonic average of the precision and recall. Therefore, this score takes both false positives and false negatives into account and provides the overall accuracy of the model. The F β score is used in our analysis. F β measures the effectiveness of retrieval with respect to a user that attaches β times as much importance to recall as precision [59,60]:
F β = ( 1 + β 2 ) · precision · recall β 2 · precision + recall
We chose the widely used F 0.5 in our evaluation.
To further investigate the performance of the proposed method, we also performed a Matthews Correlation Coefficient (MCC) [61] test. It measures the quality of the binary classification and is considered to be more truthful and informative than other parametric statistical measures [62]. The value of MCC varies between −1 and +1, where the maximum value +1 represents a perfect prediction, −1 indicates total disagreement between prediction and ground truth, and 0 represents that it is no better than random prediction. It is computed as:
MCC = ( T P × T N ) ( F P × F N ) ( T P + F P ) ( T P + F N ) ( T N + F P ) ( T N + F N )
This dataset was divided into two subsets, one used for training and other used for testing. We used different divisions—40:60, 50:50, 60:40, 70:30—for training and testing, respectively, to analyze the performance of the proposed algorithm. The results are reported in Table 1. The statistics show that the best results were gained at the 70:30 percent division. The same division trend was observed in [21], therefore in our experiments we used 70% of the dataset for training and 30% for testing. In our experiments, we performed a 4-fold cross-validation. The statistics presented in Table 1 show that our algorithm achieves an accuracy of 0.92, F score of 1.0976, and MCC value 0.8385. The results reveal that the proposed method is reliable for lung nodule detection.

3.3. Performance Comparison

We also compared the results of our method with ten existing similar nodule detection algorithms. The list of compared methods includes [11,21,29,30,31,36,53,54,55,56]. Although it is difficult to compare their performance because it depends on the image datasets and detection parameters, it is still important to attempt making a relative comparison. The Sensitivity measure is widely used to report such results, as can be found in [11,21,29,30,31,36,53,54,55,56], therefore we also used a sensitivity metric to present the comparison. The results of the comparison are presented in Table 2. To further analyze the performance of the proposed and compared methods, we also report the average FPI (False Positive per Image) and average FPE (False Positive per Exam) metric values. The FPI is the ratio of falsely accepted negative samples to the total number of images, that is,
F P I = False Positives ( FN ) Total number of tested images
Similarly, the false positive per exam (FPE) is the ratio of false positives to the total number of cases evaluated in the experiment.
The statistics presented in Table 2 reveal that our method achieved convincing results. In the sensitivity measure, our method achieved 93.75% sensitivity, outperforming all compared methods. In terms of FPI and FPE, our method achieved 0.13 and 0.22 scores, respectively, and the best results were achieved by Stelmo [21]. However, the results were computed only for 29 scans compared to our dataset of 75 scans. Moreover, our method has a better sensitivity rate than [21]. From the results of the objective performance evaluation, one can conclude that the proposed method is effective and accurate for lung nodule detection. Moreover, in contrast to most compared methods which are semi-automatic, our method is fully automatic. All the thresholds and other parameters used in our method are automatically estimated, and no external assistance is needed at any stage of the algorithm. These characteristics make the proposed algorithm ideal for lung nodule detection.

3.4. Computational Complexity Analysis

The proposed algorithm is implemented in Matlab and is made freely available for peer and public use on the project web-page (http://www.di.unito.it/~farid/Research/hls.html). We performed an execution time analysis of the proposed and the compared methods. To this end, the proposed algorithm was executed on the test dataset and average execution time was computed. The experiment was executed on an Intel® Core™ i5 processor with 4GB RAM and a 64-bit operating system. Our method takes approximately 12 s to detect nodules from a slice, which is quite efficient considering that a nodule feature is computed temporally by locating the nodule candidates in the current slice and tracing them in adjacent images. By contrast, the Stelmo [23] and Froz [53] methods take on average 90 s and 9.2 s, respectively, to examine a CT image. An efficient implementation of the proposed algorithm can further reduce its running time.

4. Summary and Conclusions

In this paper, we presented an automated system for the detection of lung nodules from CT images. The functionality of our system can be divided into two phases: first, the lung segmentation from the chest CT scan, and second, the nodule detection. The lung segmentation in the proposed algorithm is performed automatically, a novel histogram-based threshold estimation technique is proposed in this regard. This method shows that the classification of structures is based on their dimensionalities. A large number of ROIs extracted from lungs after the lung segmentation phase is a challenge for accurate classification. This problem is addressed by using the shape feature and the property that a nodule exists in consecutive slices. The testing stage of the SVM classifier resulted in 0.13 false positives per slice. The proposed algorithm achieved excellent results with a sensitivity of 0.9375, accuracy of 0.92, and a Matthews correlation coefficient of 0.8385. These results and the comparison with the existing CAD systems reveal the effectiveness of the proposed method.
The high incident rate of lung cancer and the late diagnosis show that this automated system can be conducive in the early scan stages. An exam of lung CT scan consists of a long series of images, and this system can analyze these images fast and reduce the risk of human error. This system can be used as the first step of a diagnosis; the marked cases can be passed for medical analysis for further studies and confirmation. We need to mention here that public hospitals lack in number of specialists and the number of patients visiting the hospital is enormous as compared to the number of doctors available. The proposed method can be used to lower this burden. The method can be implemented for preliminary scans and a radiologist can validate these results. Thus this system will be an assistant to the radiologists. This system is economical to design as it requires regular computers for deployment, which are usually already available in hospitals and clinics or can be easily procured.

Author Contributions

Conceptualization, N.K.; methodology, N.K., M.S.F., S.B., and M.H.K.; software, N.K. and M.S.F.; validation, N.K., S.B., and M.S.F.; writing—original draft preparation, N.K. and M.S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors acknowledge the National Cancer Institute and the Foundation for the National Institutes of Health, and their critical role in the creation of the free publicly available LIDC/IDRI database used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ferlay, J.; Shin, H.; Bray, F.; Forman, D.; Mathers, C.; Parkin, D.M. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int. J. Cancer 2010, 127, 2893–2917. [Google Scholar] [CrossRef]
  2. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2018. CA-Cancer J. Clin. 2018, 68, 7–30. [Google Scholar] [CrossRef] [PubMed]
  3. Capocaccia, R.; Gatta, G.; Dal Maso, L. Life expectancy of colon, breast, and testicular cancer patients: An analysis of US-SEER population-based data. Ann. Oncol. 2015, 26, 1263–1268. [Google Scholar] [CrossRef] [PubMed]
  4. Ning, J.; Zhao, H.; Lan, L.; Sun, P.; Feng, Y. A Computer-Aided Detection System for the Detection of Lung Nodules Based on 3D-ResNet. Appl. Sci. 2019, 9, 5544. [Google Scholar] [CrossRef] [Green Version]
  5. Saleem, M.; Farid, M.S.; Saleem, S.; Khan, M.H. X-ray image analysis for automated knee osteoarthritis detection. Signal Image Video Process. 2020. [Google Scholar] [CrossRef]
  6. Thomson, J.; O’Neill, T.; Felson, D.; Cootes, T. Automated shape and texture analysis for detection of osteoarthritis from radiographs of the knee. In International Conference on Medical Image Computing and Computer-Assisted Intervention; Springer: Cham, Switzerland, 2015; pp. 127–134. [Google Scholar]
  7. Das, D.; Mukherjee, R.; Chakarborty, C. Computational microscopic imaging for malaria parasite detection: A systematic review. J. Microsc. 2015, 260, 1–19. [Google Scholar] [CrossRef]
  8. Fatima, T.; Farid, M.S. Automatic detection of Plasmodium parasites from microscopic blood images. J. Parasit. Dis. 2019. [Google Scholar] [CrossRef]
  9. Mehanian, C.; Jaiswal, M.; Delahunt, C.; Thompson, C.; Horning, M.; Hu, L.; Ostbye, T.; McGuire, S.; Mehanian, M.; Champlin, C.; et al. Computer-Automated Malaria Diagnosis and Quantitation Using Convolutional Neural Networks. In Proceedings of the IEEE International Conference on Computer Vision Workshops (ICCVW), Venice, Italy, 22–29 October 2017. [Google Scholar]
  10. Nasrullah, N.; Sang, J.; Alam, M.S.; Mateen, M.; Cai, B.; Hu, H. Automated Lung Nodule Detection and Classification Using Deep Learning Combined with Multiple Strategies. Sensors 2019, 19, 3722. [Google Scholar] [CrossRef] [Green Version]
  11. Saien, S.; Moghaddam, H.A.; Fathian, M. A unified methodology based on sparse field level sets and boosting algorithms for false positives reduction in lung nodules detection. Int. J. Comput. Assist. Radiol. Surg. 2018, 13, 397–409. [Google Scholar] [CrossRef]
  12. Pena, D.M.; Luo, S.; Abdelgader, A.M.S. Auto Diagnostics of Lung Nodules Using Minimal Characteristics Extraction Technique. Diagnostics 2016, 6, 13. [Google Scholar] [CrossRef] [Green Version]
  13. Gu, Y.; Pandit, S.; Saraee, E.; Nordahl, T.; Ellis, T.; Betke, M. Home-Based Physical Therapy with an Interactive Computer Vision System. In Proceedings of the IEEE International Conference on Computer Vision Workshops (ICCVW), Seoul, Korea, 27 October–2 November 2019. [Google Scholar]
  14. Khan, M.H.; Helsper, J.; Farid, M.S.; Grzegorzek, M. A computer vision-based system for monitoring Vojta therapy. Int. J. Med. Inform. 2018, 113, 85–95. [Google Scholar] [CrossRef] [PubMed]
  15. Khan, M.H.; Schneider, M.; Farid, M.S.; Grzegorzek, M. Detection of Infantile Movement Disorders in Video Data Using Deformable Part-Based Model. Sensors 2018, 18, 3202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Narayanan, B.N.; Hardie, R.C.; Kebede, T. M Performance analysis of a computer-aided detection system for lung nodules in CT at different slice thicknesses. J. Med. Imaging 2018, 5, 5–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Nishio, M.; Nishizawa, M.; Sugiyama, O.; Kojima, R.; Yakami, M.; Kuroda, T.; Togashi, K. Computer-aided diagnosis of lung nodule using gradient tree boosting and Bayesian optimization. PLoS ONE 2018, 13, e0195875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zhang, G.; Yang, Z.; Gong, L.; Jiang, S.; Wang, L.; Cao, X.; Wei, L.; Zhang, H.; Liu, Z. An Appraisal of Nodule Diagnosis for Lung Cancer in CT Images. J. Med. Syst. 2019, 43, 181. [Google Scholar] [CrossRef]
  19. El-Regaily, S.A.; Salem, M.A.; Abdel Aziz, M.H.; Roushdy, M.I. Survey of Computer Aided Detection Systems for Lung Cancer in Computed Tomography. Curr. Med. Imaging Rev. 2018, 14, 3–18. [Google Scholar] [CrossRef]
  20. Rajan, J.R.; Chelvan, A.C.; Duela, J.S. Multi-Class Neural Networks to Predict Lung Cancer. J. Med. Syst. 2019, 43, 211. [Google Scholar] [CrossRef]
  21. Netto, S.M.B.; Silva, A.C.; Nunes, R.A.; Gattass, M. Automatic segmentation of lung nodules with growing neural gas and support vector machine. Comput. Biol. Med. 2012, 42, 1110–1121. [Google Scholar] [CrossRef]
  22. Lee, S.L.A.; Kouzani, A.Z.; Hu, E.J. Automated identification of lung nodules. In Proceedings of the 10th Workshop on Multimedia Signal Processing (MMSP), Cairns, Australia, 8–10 October 2008; pp. 497–502. [Google Scholar]
  23. Cascio, D.; Magro, R.; Fauci, F.; Iacomi, M.; Raso, G. Automatic detection of lung nodules in CT datasets based on stable 3D mass–spring models. Comput. Biol. Med. 2012, 42, 1098–1109. [Google Scholar] [CrossRef] [Green Version]
  24. Akram, S.; Javed, M.Y.; Akram, M.U.; Qamar, U.; Hassan, A. Pulmonary Nodules Detection and Classification Using Hybrid Features from Computerized Tomographic Images. J. Med. Imaging Health Inform. 2016, 6, 252–259. [Google Scholar] [CrossRef]
  25. Keshani, M.; Azimifar, Z.; Tajeripour, F.; Boostani, R. Lung nodule segmentation and recognition using SVM classifier and active contour modeling: A complete intelligent system. Comput. Biol. Med. 2013, 43, 287–300. [Google Scholar] [CrossRef] [PubMed]
  26. Samuel, C.C.; Saravanan, V.; Devi, M.R.V. Lung Nodule Diagnosis from CT Images Using Fuzzy Logic. In Proceedings of the ICCIMA, Sivakasi, India, 13–15 December 2007; Volume 3, pp. 159–163. [Google Scholar]
  27. Özekes, S. Rule-Based Lung Region Segmentation and Nodule Detection via Genetic Algorithm Trained Template Matching. İstanbul Ticaret Üniversitesi Fen Bilim. Derg. 2007, 6, 17–30. [Google Scholar]
  28. Pu, J.; Zheng, B.; Leader, J.; Wang, X.; Gur, D. An automated CT based lung nodule detection scheme using geometric analysis of signed distance field. Med. Phys. 2008, 35, 3453–3461. [Google Scholar] [CrossRef] [PubMed]
  29. Messay, T.; Hardie, R.C.; Rogers, S.K. A new computationally efficient CAD system for pulmonary nodule detection in CT imagery. Med. Image Anal. 2010, 14, 390–406. [Google Scholar] [CrossRef]
  30. Tan, M.; Deklerck, R.; Jansen, B.; Bister, M.; Cornelis, J. A novel computer-aided lung nodule detection system for CT images. Med. Phys. 2011, 38, 5630–5645. [Google Scholar] [CrossRef]
  31. Golosio, B.; Masala, G.L.; Piccioli, A.; Oliva, P.; Carpinelli, M.; Cataldo, R.; Cerello, P.; De Carlo, F.; Falaschi, F.; Fantacci, M.E.; et al. A novel multithreshold method for nodule detection in lung CT. Med. Phys. 2009, 36, 3607–3618. [Google Scholar] [CrossRef]
  32. Shi, Z.; Ma, J.; Zhao, M.; Liu, Y.; Feng, Y.; Zhang, M.; He, L.; Suzuki, K. Many is better than one: An integration of multiple simple strategies for accurate lung segmentation in CT Images. BioMed Res. Int. 2016, 2016, 1480423. [Google Scholar] [CrossRef] [Green Version]
  33. Ozekes, S.; Osman, O.; Ucan, O.N. Nodule detection in a lung region that’s segmented with using genetic cellular neural networks and 3D template matching with fuzzy rule based thresholding. Korean J. Radiol. 2008, 9, 1–9. [Google Scholar] [CrossRef] [Green Version]
  34. Jo, H.H.; Hong, H.; Goo, J.M. Pulmonary nodule registration in serial CT scans using global rib matching and nodule template matching. Comput. Biol. Med. 2014, 45, 87–97. [Google Scholar] [CrossRef]
  35. El-Baz, A.; Elnakib, A.; Abou El-Ghar, M.; Gimel’farb, G.; Falk, R.; Farag, A. Automatic detection of 2D and 3D lung nodules in chest spiral CT scans. Int. J. Biomed. Imaging 2013, 2013, 517632. [Google Scholar] [CrossRef]
  36. Dolejsi, M.; Kybic, J.; Polovincak, M.; Tuma, S. The Lung TIME: Annotated lung nodule dataset and nodule detection framework. In Proceedings of the Medical Imaging 2009: Computer-Aided Diagnosis, Lake Buena Vista (Orlando Area), FL, USA, 7–12 February 2009; p. 72601U. [Google Scholar]
  37. Opfer, R.; Wiemker, R. Performance analysis for computer-aided lung nodule detection on LIDC data. In Proceedings of the Medical Imaging 2007: Image Perception, Observer Performance, and Technology Assessment, San Diego, CA, USA, 17–22 February 2007; Volume 6515, p. 65151C. [Google Scholar]
  38. Narayanan, B.N.; Hardie, R.C.; Kebede, T.M.; Sprague, M.J. Optimized Feature Selection-Based Clustering Approach for Computer-Aided Detection of Lung Nodules in Different Modalities. Pattern Anal. Appl. 2019, 22, 559–571. [Google Scholar] [CrossRef]
  39. Pei, X.; Guo, H.; Dai, J. Computerized Detection of Lung Nodules in CT Images by Use of Multiscale Filters and Geometrical Constraint Region Growing. In Proceedings of the 4th International Conference on Bioinformatics and Biomedical Engineering, Chengdu, China, 18–20 June 2010; pp. 1–4. [Google Scholar]
  40. Soliman, A.; Khalifa, F.; Elnakib, A.; Abou El-Ghar, M.; Dunlap, N.; Wang, B.; Gimel’farb, G.; Keynton, R.; El-Baz, A. Accurate Lungs Segmentation on CT Chest Images by Adaptive Appearance-Guided Shape Modeling. IEEE Trans. Med. Imaging 2017, 36, 263–276. [Google Scholar] [CrossRef] [PubMed]
  41. da Silva Sousa, J.R.F.; Silva, A.C.; de Paiva, A.C.; Nunes, R.A. Methodology for automatic detection of lung nodules in computerized tomography images. Comput. Methods Programs Biomed. 2010, 98, 1–14. [Google Scholar] [CrossRef] [PubMed]
  42. Negahdar, M.; Beymer, D.; Syeda-Mahmood, T. Automated volumetric lung segmentation of thoracic CT images using fully convolutional neural network. In Proceedings of the Medical Imaging 2018: Computer-Aided Diagnosis, Houston, TX, USA, 10–15 February 2018; Petrick, N., Mori, K., Eds.; Volume 10575, pp. 356–361. [Google Scholar]
  43. Gruetzemacher, R.; Gupta, A.; Paradice, D. 3D deep learning for detecting pulmonary nodules in CT scans. J. Am. Med. Inf. Assoc. 2018, 25, 1301–1310. [Google Scholar] [CrossRef] [PubMed]
  44. Jiang, H.; Ma, H.; Qian, W.; Gao, M.; Li, Y. An Automatic Detection System of Lung Nodule Based on Multigroup Patch-Based Deep Learning Network. IEEE J. Biomed. Health Inform. 2018, 22, 1227–1237. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, G.; Jiang, S.; Yang, Z.; Gong, L.; Ma, X.; Zhou, Z.; Bao, C.; Liu, Q. Automatic nodule detection for lung cancer in CT images: A review. Comput. Biol. Med. 2018, 103, 287–300. [Google Scholar] [CrossRef]
  46. Shaukat, F.; Raja, G.; Frangi, A.F. Computer-aided detection of lung nodules: A review. J. Med. Imaging 2019, 6, 020901. [Google Scholar] [CrossRef]
  47. Pehrson, L.M.; Nielsen, M.B.; Ammitzbøl Lauridsen, C. Automatic Pulmonary Nodule Detection Applying Deep Learning or Machine Learning Algorithms to the LIDC-IDRI Database: A Systematic Review. Diagnostics 2019, 9, 29. [Google Scholar] [CrossRef] [Green Version]
  48. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef] [Green Version]
  49. Samet, H.; Tamminen, M. Efficient component labeling of images of arbitrary dimension represented by linear bintrees. IEEE Trans. Pattern Anal. Mach. Intell. 1988, 10, 579–586. [Google Scholar] [CrossRef] [Green Version]
  50. Atherton, T.; Kerbyson, D. Size invariant circle detection. Image Vis. Comput. 1999, 17, 795–803. [Google Scholar] [CrossRef]
  51. Adams, R.; Bischof, L. Seeded region growing. IEEE Trans. Pattern Anal. Mach. Intell. 1994, 16, 641–647. [Google Scholar] [CrossRef] [Green Version]
  52. Fan, R.E.; Chang, K.W.; Hsieh, C.J.; Wang, X.R.; Lin, C.J. LIBLINEAR: A library for large linear classification. J. Mach. Learn. Res. 2008, 9, 1871–1874. [Google Scholar]
  53. Froz, B.R.; de Carvalho Filho, A.O.; Silva, A.C.; de Paiva, A.C.; Nunes, R.A.; Gattass, M. Lung nodule classification using artificial crawlers, directional texture and support vector machine. Expert Syst. Appl. 2017, 69, 176–188. [Google Scholar] [CrossRef]
  54. Teramoto, A.; Fujita, H. Fast lung nodule detection in chest CT images using cylindrical nodule-enhancement filter. Int. J. Comput. Assist. Radiol. Surg. 2013, 8, 193–205. [Google Scholar] [CrossRef]
  55. Wu, P.; Xia, K.; Yu, H. Correlation coefficient based supervised locally linear embedding for pulmonary nodule recognition. Comput. Methods Programs Biomed. 2016, 136, 97–106. [Google Scholar] [CrossRef] [Green Version]
  56. Bergtholdt, M.; Wiemker, R.; Klinder, T. Pulmonary nodule detection using a cascaded SVM classifier. In Proceedings of the Medical Imaging 2016: Computer-Aided Diagnosis, San Diego, CA, USA, 27 February–3 March 2016; Volume 9785, pp. 268–278. [Google Scholar]
  57. Armato, S.G., III; McLennan, G.; McNitt-Gray, M.F.; Meyer, C.R.; Yankelevitz, D.; Aberle, D.R.; Henschke, C.I.; Hoffman, E.A.; Kazerooni, E.A.; MacMahon, H.; et al. Lung Image Database Consortium: Developing a Resource for the Medical Imaging Research Community. Radiology 2004, 232, 739–748. [Google Scholar]
  58. Farid, M.S.; Lucenteforte, M.; Grangetto, M. DOST: A distributed object segmentation tool. Multimed. Tools Appl. 2018, 77, 20839–20862. [Google Scholar] [CrossRef] [Green Version]
  59. Sasaki, Y. The truth of the F-measure. Teach Tutor Mater. 2007, 1, 1–5. [Google Scholar]
  60. Chinchor, N. The Statistical Significance of the MUC-4 Results. In Proceedings of the 4th Conference on Message Understanding, MUC4 ’92; Association for Computational Linguistics: Stroudsburg, PA, USA, 1992; pp. 30–50. [Google Scholar]
  61. Matthews, B. Comparison of the predicted and observed secondary structure of T4 phage lysozyme. Biochim. Biophys. Acta BBA Protein Struct. 1975, 405, 442–451. [Google Scholar] [CrossRef]
  62. Chicco, D.; Jurman, G. Machine learning can predict survival of patients with heart failure from serum creatinine and ejection fraction alone. BMC Med. Inform. Decis. Mak. 2020, 20, 16. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Block diagram of the proposed algorithm.
Figure 1. Block diagram of the proposed algorithm.
Jimaging 06 00006 g001
Figure 2. A CT image of lung and proposed division into four regions.
Figure 2. A CT image of lung and proposed division into four regions.
Jimaging 06 00006 g002
Figure 3. Automatic histogram based on initial segmentation: (a) Input CT slice; (b) histogram of the image in (a); (c) the results after thresholding the image with the threshold estimated from the histogram in (b); (d) the complemented image of (c); (e) binary segmentation map; (f) map after connected components-based refinement; (g) detected lungs; (h) lungs after noise removal.
Figure 3. Automatic histogram based on initial segmentation: (a) Input CT slice; (b) histogram of the image in (a); (c) the results after thresholding the image with the threshold estimated from the histogram in (b); (d) the complemented image of (c); (e) binary segmentation map; (f) map after connected components-based refinement; (g) detected lungs; (h) lungs after noise removal.
Jimaging 06 00006 g003
Figure 4. The region circled in red is a nodule and the region circled in green is a non-nodule.
Figure 4. The region circled in red is a nodule and the region circled in green is a non-nodule.
Jimaging 06 00006 g004
Figure 5. Spherical shape detection (reader is requested to magnify the images to see the details). The 2nd row shows the processing of the slice under consideration, whereas the 1st and the 3rd rows show the processing of previous and next slices, respectively. (a) Segmented lungs, (b) the extracted internal structures, (c) the detected circle in red boundaries, and (d) the final detected nodules.
Figure 5. Spherical shape detection (reader is requested to magnify the images to see the details). The 2nd row shows the processing of the slice under consideration, whereas the 1st and the 3rd rows show the processing of previous and next slices, respectively. (a) Segmented lungs, (b) the extracted internal structures, (c) the detected circle in red boundaries, and (d) the final detected nodules.
Jimaging 06 00006 g005
Figure 6. Nodule detection results: (a,d) input CT slice, (b,e) segmented lungs, (c,f) detected nodule.
Figure 6. Nodule detection results: (a,d) input CT slice, (b,e) segmented lungs, (c,f) detected nodule.
Jimaging 06 00006 g006
Table 1. Performance evaluation of the proposed algorithm on different training and testing dataset divisions. The division shows the ratio of training to test dataset size. The best results are marked in bold.
Table 1. Performance evaluation of the proposed algorithm on different training and testing dataset divisions. The division shows the ratio of training to test dataset size. The best results are marked in bold.
DivisionSensitivitySpecificityPrecisionAccuracyF ScoreMCC
40:600.86110.88240.83780.87361.02760.3491
50:500.86210.88640.83330.87671.02740.3908
60:400.92000.88890.85190.90161.08660.5838
70:300.93750.91180.83330.92001.09760.8385
Table 2. Performance evaluation of the proposed algorithm and the compared methods. Size represents the number of scans in the dataset. FPI: False Positive per Image; FPE: False Positive per Exam.
Table 2. Performance evaluation of the proposed algorithm and the compared methods. Size represents the number of scans in the dataset. FPI: False Positive per Image; FPE: False Positive per Exam.
MethodYearDatabaseSizeSensitivityFPIFPE
Dolejsi [36]2009TIME-LIDC-ANODE3889.6012.03-
Golosio [31]2009LIDC48471.00-4
Messay [29]2010LIDC8482.66-3
Tan [30]2011LIDC39987.50-4
Stelmo [21]2012LIDC2985.930.0010.14
Teramoto [54]2013LIDC8480.00-4.2
Bergtholdt [56]2016LIDC-IDRI24385.90-2.5
Wu [55]2017LIDC-IDRI6079.23--
Froz [53]2017LIDC-IDRI83391.86--
Saien [11]2018LIDC/LIDC-IDRI7083.980.02-
Ours2019LIDC7593.750.130.22

Share and Cite

MDPI and ACS Style

Khehrah, N.; Farid, M.S.; Bilal, S.; Khan, M.H. Lung Nodule Detection in CT Images Using Statistical and Shape-Based Features. J. Imaging 2020, 6, 6. https://doi.org/10.3390/jimaging6020006

AMA Style

Khehrah N, Farid MS, Bilal S, Khan MH. Lung Nodule Detection in CT Images Using Statistical and Shape-Based Features. Journal of Imaging. 2020; 6(2):6. https://doi.org/10.3390/jimaging6020006

Chicago/Turabian Style

Khehrah, Noor, Muhammad Shahid Farid, Saira Bilal, and Muhammad Hassan Khan. 2020. "Lung Nodule Detection in CT Images Using Statistical and Shape-Based Features" Journal of Imaging 6, no. 2: 6. https://doi.org/10.3390/jimaging6020006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop