A Segmentation Method for Lung Parenchyma Image Sequences Based on Superpixels and a Self-Generating Neural Forest.

BACKGROUND: Lung parenchyma segmentation is often performed as an important pre-processing step in the computer-aided diagnosis of lung nodules based on CT image sequences. However, existing lung parenchyma image segmentation methods cannot fully segment all lung parenchyma images and have a slow processing speed, particularly for images in the top and bottom of the lung and the images that contain lung nodules.
METHOD: Our proposed method first uses the position of the lung parenchyma image features to obtain lung parenchyma ROI image sequences. A gradient and sequential linear iterative clustering algorithm (GSLIC) for sequence image segmentation is then proposed to segment the ROI image sequences and obtain superpixel samples. The SGNF, which is optimized by a genetic algorithm (GA), is then utilized for superpixel clustering. Finally, the grey and geometric features of the superpixel samples are used to identify and segment all of the lung parenchyma image sequences.
RESULTS: Our proposed method achieves higher segmentation precision and greater accuracy in less time. It has an average processing time of 42.21 seconds for each dataset and an average volume pixel overlap ratio of 92.22 ± 4.02% for four types of lung parenchyma image sequences.

Introduction

Lung cancer is one of the most common causes of cancer-related death worldwide [1]. Computed tomography (CT) [2] scanning technology has good density resolution for lesions in the human body and is currently the most effective and direct imaging method for the early diagnosis of lung cancer. However, as the accuracy requirements for clinical imaging of lesions increase, the CT scanning thickness decreases, and a large number of CT image sequences need to be produced [3]. The massive amount of image data will inevitably increase the challenge of CT image processing, leading to a slow processing speed and decreased efficiency. In addition, because each pulmonary CT image presents a different morphological structure from the top to the bottom of the lung in CT image sequences, the general segmentation algorithm is not effective. Therefore, determining how to segment the lung parenchyma image sequences quickly without reducing accuracy is of great significance for the subsequent segmentation of pulmonary nodules and benign and malignant diagnoses.

Lung segmentation can be an important component of computer-aided diagnosis (CAD) systems [4]. Geng H’s group used an iterative gray threshold to select seed points automatically and then extract each lung parenchyma image with the region growing method, which is sensitive to background noise [5]. Liming D and colleague present a new form of lung parenchyma segmentation. The optimal threshold value method and the boundary tracking method are used to segment the lung region and can effectively eliminate the influence of background noise but may lose some of the lung parenchyma [6]. Mansoor A and coworkers segmented the lung parenchyma in two steps [7] by using the fuzzy connectedness (FC) image segmentation algorithm to perform the initial lung parenchyma extraction and then texture-based local descriptors to segment abnormal imaging patterns using a near-optimal keypoint analysis. However, this method is not effective for processing irregular images. Wavelet transform has been applied by Shojaii R [8] to decompose an image into several regions, and the regions with low pixel intensities are kept and grown to segment the honeycomb regions. This method can effectively segment irregular lung parenchyma images. Yan-hua and coworkers used several methods, including the optimal iterative threshold, three-dimensional connectivity labeling, and three-dimensional region growing methods, for the initial segmentation of the lung parenchyma and used the morphological method to repair the lung parenchyma [9]. Luo X [10] and others used an improved active contour model, which can obtain better segmentation result with the help of artificial segmentation but is very time consuming. There are also some scholars who use superpixel to segment medical images. Yu N and Weinstein S P [11] proposed a novel automatic segmentation framework for tumor on breast DCE-MRI images by using graph-cuts and superpixel classification, which can achieve a classification accuracy of 96%.

Superpixel was also used for bacteria cell segmentation by Song Y’s group [12]. Features of superpixels are extracted and trained by supervised deep learning method with an accuracy of 99% and a sensitivity of 100% for four types of different bacteria.

In general, lung segmentation methods are based on threshold, region, and mathematical morphology. However, lung CT images are sequential, and the existing methods in the study of lung parenchyma segmentation algorithms are generally for single-image segmentation of CT images and ignore the before-to-after image correlation. A few scholars have studied sequential image segmentation, but this often involves a long processing time, low efficiency and extensibility.

In this paper, we use the position particularity of the lung parenchyma in lung CT images, fully consider the strong correlation between adjacent slices of CT image sequences, and put forward a segmentation method for lung parenchyma image sequences based on superpixels and a self-generating neural forest. The experimental results show that our proposed method can significantly increase the speed of segmentation for four types of lung parenchyma images, which guarantees accuracy and integrity.

Materials and Methods

2.1 Materials

2.1.1 Ethics statement

This study was approved by the institutional review board (IRB) of the Coal Center Hospital in Shanxi. The study was conducted in accordance with the hospital’s ethics requirements. Informed consent was obtained from all patients for being included in the study.

2.1.2 Datasets

The CT image datasets used in this study were obtained from a hospital in Shanxi Province, China. All data can be accessed at https://figshare.com/s/254e3467efd57a442334. We used a Discovery ST16 PET-CT scanner from the General Electric Company of America (150 mA, 140 kV, with a slice thickness of 3.75 mm). In the experiment, we select lung CT sequence image datasets from 80 people with a total of 4812 CT images, and the size of each image was 512 × 512. Based on the physician’s prior knowledge and the morphological perspective of lung CT image sequences, the 80 datasets were divided into four categories: without nodules, benign nodules, malignant SPN (solitary pulmonary nodules) and pleural nodules. Each category had 20 datasets and approximately 1200 CT images.

2.2 Proposed Method

We propose a segmentation method for lung parenchyma image sequences based on superpixels and a self-generating neural forest that mainly involves a gradient and sequential linear iterative clustering algorithm (GSLIC) to obtain superpixels, clustering of superpixels with a self-generating neural forest (SGNF), and lung parenchyma image sequences segmentation. A block diagram of the lung parenchyma image sequences segmentation is shown in .

2.2.1 Gradient and sequential linear iterative clustering (GSLIC)

It has been difficult to obtain better segmentation in a relatively short time in lung parenchyma image sequence segmentation. To solve this problem, we used multiple CT images and drew on the prior knowledge of the physicians to propose a superpixel segmentation algorithm for image sequence. Our method is based on a gradient and sequential linear iterative clustering (GSLIC) algorithm and includes lung ROI image sequence extraction and ROI sequence superpixel segmentation.

1. Lung ROI sequence extraction: Because the position of the lung parenchyma region in CT images is relatively fixed, we increased the running speed by first extracting the lung ROI sequences using a statistical method for ROI extraction that is adopted in this paper. More than 4800 CT lung images of 80 individuals were analyzed, and we determined that the rectangles in the upper-left (100, 60) and lower-right corners (400,420) could include all lung parenchyma regions. Therefore, we can obtain the entire lung ROI sequence based on the two points in the CT image sequence. In , the original lung CT image (a) was used to extract its ROI image (b). Extracting the lung ROI sequences can reduce the processing time and simultaneously eliminate some noise.

Lung ROI extraction result (b) of the original lung CT image (a).

2. Superpixel segmentation on ROI image sequences: The concept of superpixels was first put forward by Ren [13] in 2003. A superpixel is a collection of pixels with similar characteristics, such as color, brightness, and texture. An image can be composed of a certain number of superpixels that contain multiple combination characteristics of the pixels and can preserve the edge information of the original image. Compared with a single pixel, a superpixel contains rich characteristic information and can greatly reduce image post-processing complexity and significantly increase the speed of image segmentation.

Traditional superpixel segmentation uses a process of simple, linear iterative clustering. After this method was improved by Lucchi, Hammoudi, and Wang J [14-16], it was applied to a single image segmentation. In this paper, an algorithm based on gradient and sequential linear iterative clustering (GSLIC) is proposed to segment image sequences. Each pixel in the sequence of lung CT images can be represented by a six-dimensional feature vector ([ ). The similarity between the pixels can be measured by the Euclidean distance between them. A pixel’s feature vector is made up of its color vector [ in color space and its space coordinate vector [, where and are the pixel coordinates, and is the serial number of the image. The GSLIC procedure is shown in .

In an original CT image with pixels that needs to be divided into superpixels, each superpixel contains approximately pixels, and therefore, the average length of each superpixel is about (1).

We first take an initial clustering center every pixels and then select cluster centers using a 3 * 3 nuclear window to the lowest gradient position. When selecting the initial clustering center, a method of gradient descent is adopted to sample pixels at a regular grid so that the edge points are not selected as the cluster centers. By using a nuclear window, a pixel’s gradient G(x, y) can be defined as (2). The coordinates in each grid can be chose as the cluster center. Following that, each clustering center can be search for neighboring similar pixels around the search space for * based on the similarity of between the pixels.

In our method, for the same CT image with serial number , the similarity of between the pixels ([, , , , , ) to the clustering center ([, , , , , ) can be calculated by their color feature distance and space feature distance . The calculation formulas of , and are as follows in (3), (4) and (5).

In (5), is a parameter to adjust the weight of and . The larger the value, the bigger the weight of to calculate will be, which is generally between 1 and 20. The result of superpixels segmentation on ROI image is shown in

Superpixel segmentation result (b) of the original lung ROI image (a).

The proposed method can segment an image into a series of superpixels, each of which can also be expressed by feature vector ([ ). In the process of obtaining superpixels, effective clustering centers can be chose by gradient descent, and the blocks can be obtained through the clustering algorithm. In addition, taking the correlation between the sequences of CT images into account, the coordinate information of the clustering centers in the previous image is directly transmitted to the next image, which can significantly improve the image’s superpixel segmentation speed.

2.2.2 Superpixel clustering based on SGNF and GA

In this paper, a method using a self-generated neural forest (SGNF) algorithm optimized by a genetic algorithm (GA) was proposed to cluster the superpixels. The genetic algorithm is used to select the optimal clustering centers, which are used to generate the neural trees that form the neural forest. Our method effectively overcomes the instability of the primary SGNN and improves the efficiency and accuracy of clustering.

1. Self-generating neural tree (SGNT): Self-generated neural networks (SGNNs) were developed in 1992 [17] using a competitive learning mechanism for samples learning. A self-generating neural tree (SGNT) is generated by an SGNN using unsupervised learning.

An SGNT includes neurons, weights, and connections. In this paper, we use an ordered pair < {n}, {l} > to express an SGNT, where {n} is the set of neurons, and {l} is the set of connections. Each neuron also can be represented as ordered pair < w, {n>, where is the weight of the neuron, and {n} is the set of child neurons of the neuron. Each leaf neuron corresponds to a sample, and each root neuron is a cluster center. All leaf neurons of the root neuron belong to the same cluster, and the weight of every neuron is the average attribute of all the leaf neurons it covers.

Therefore, the structure of an SGNN shows simplicity and a good self-organizing capability and learning speed. It is beneficial to learn clustering that has high performance. shows the structure of an SGNT with five samples. lists a clustering sample set where W, j = a, b, …, e are the sample attributes. is the SGNT generated by following SGNT generating rules [18-20].

The structure of the SGNT (b) generated by five samples (a).

2. Self-generating neural forest (SGNF) optimized by GA: Despite its good capacity for clustering, an SGNT is influenced by the input order of the samples [21]. To solve this problem, we propose an adaptive clustering algorithm that is optimized by GA and in which the SGNT is generalized to a self-generating neural forest (SGNF), and the GA is applied to select optimal superpixel seeds as the initial input into the SGNF. The process of clustering with an SGNF optimized by a GA is shown in .

Each SGNT in an SGNF corresponds to a cluster, and all the leaf neurons in an SGNT belong to the same cluster. In section 2.2.1, we show that each superpixel can be expressed by a feature vector ([ ). For given sample {} where = 1, 2, …, L, the distance between sample and clustering center can be calculated as (6): where is the sequence number of the element in the feature vector, and w is the weight of the first k attribute. When processing the superpixels segmented by GSLIC, the attributes of superpixels, such as color and coordinate feature values, can be used to generate an SGNF. The result of superpixel clustering with GA-SGNF is shown in .

Superpixel clustering result (b) of the superpixel segmentation image (a).

Different seed samples generate a SGNF with different structures. Therefore, the choice of K seed samples can be seen as an optimization problem. To obtain preferable clustering results, the GA is used to search for K seed samples to optimize the clustering results as described in the next section.

3. Genetic Algorithm: The genetic algorithm (GA) is a method based on the probability search technology of population optimization [22-24]. The GA has a good ability for global searching and can search for the optimal solution quickly. In the previous section, each superpixel is expressed by a feature vector that can also correspond to the process of chromosome encoding in the GA. In addition, the sample capacity and clustering numbers are small, making the algorithm converge rapidly. Therefore, the optimal seed points will be found by the GA to obtain the best clustering results in a short time. The process of choosing optimal seed points using the GA is shown in .

We first define a chromosome structure C = (, , ), where the ci, = 1, 2, …, K, including initial K superpixels, and K is user-specified. For each superpixel with gene code string , , ), where in , = 1, 2, …, L, represents the number of superpixel attributes in section 2.2.1.

For the given superpixels {,}, we finally obtain K classes of superpixels with cluster center R = (, , ). The number of superpixels for cluster r is n, and x are all superpixels in cluster r. We define a fitness function φ(C) that can be obtained by the between-class variance δ2 to evaluate the goodness of a chromosome. φ(C) and δ2 can be calculated using (7) and (8).

The higher the value of φ(C), the better the chromosome quality is assumed to be. The chromosome with the maximum value of φ(C) is considered the optimal one in the population, and the K superpixel seeds are chose to generate the SNGF.

The selection process copies individual strings with high fitness function values into the next population based on the ‘‘roulette wheel” selection approach. The main purpose of crossover is to exchange genetic information of the selected chromosomes. Mutation is the process of a random alteration in the genetic structure of a chromosome, which can introduce genetic diversity into the population. In our method, the probabilities of crossover and mutation are τ and η, respectively, and the termination criteria are as follows:

The biggest fitness function value is obtained and the algorithm converges.

The fixed number of generations is reached.

2.2.3 Feature extraction and lung parenchyma segmentation

1. Feature extraction and lung identification: After clustering the superpixel samples using the optimized SGNF algorithm, four superpixel sample sets that include the left and the right lung parenchyma images, pleural tissue and extrathoracic area are obtained. We still need to identify and segment the lung parenchyma from the image sequences. Because the average greyscale value of each superpixel sample set is equal to the average value of all superpixels, the sample set with the highest value should be the pleural tissue, and the value for the lung parenchyma is close to that of the extrathoracic area. As the distribution of superpixel coordinates of the lung parenchyma is relatively concentrated, the sample set with the highest value of coordinate variance should be the extrathoracic area. And the left two sample sets should be the left and the right lung parenchyma. In this paper we will mainly extract two features of superpixel samples: the average grayscale value and the coordinate variance. An overview of the detection method is shown in .

As mentioned previously, each superpixel can be expressed by feature vector ([ ), , , ). For a sample set {X}, the formulas to determine the average grayscale value, centroid coordinate and coordinate variance are as follows:

The sample set with the smaller coordinate variance should be the left and the right lung parenchyma images, which are assumed to be S1 and S2. We still need to traverse all superpixel samples in S1 and S2. As the superpixels in the same image will have the same attribute value z, all superpixel samples can be connected according to the value of attribute z. By sequentially outputting all images, the sequential coarse lung parenchyma images were fully segmented.

2. Removing the trachea/bronchus and refining the lung contour: After coarse segmentation of lung CT image sequences, there are still the trachea/bronchus at the top of the lung image. To ensure the integrity of the lung parenchyma segmentation, we adopt an improved region growing method [25] to remove them. The description of improved region growing is given as follows.

Step 1: Binarization for the coarse lung image sequences

Step 2: Extract the minimum bounding rectangle of the lung.

Step 3: Select seed points by using LRS algorithm.

Step 4: Refine the lung contours with erosion and dilation.

Step 5: Acquire the final lung mask sequences.

We first use adaptive threshold method for image binarization and extract the minimum bounding rectangle of the lung. And then LRS algorithm is employed to select left and right seed points. In LRS, scan the minimum bounding rectangle image along the left and right sides simultaneously until there are more than consecutive points on direction with the pixel value 255, and record the middle (third) one’s ordinate value as the seeds. Next we adopt the improved region growing method based on these seeds to discard disconnected trachea, bronchus and other noise and extract lung out. Finally, dilation and erosion are used to smooth the contour and eliminate some vessels, small nodules as well as bones. Thus we acquire the final lung mask sequences with which to segment lung parenchyma image sequences accurately. The process of removing the trachea/bronchus and refining the lung contour is shown in

The process of trachea/bronchus removing and lung contour refining.

(a) Binarization of the coarse lung image; (b) Extraction of minimum bounding rectangle; (c) Select seed points with the LRS algorithm; (d) Final lung mask.

Results

To verify our method’s validity and universality on lung parenchyma segmentation for these four types of lung CT image sequences, we compare the results of our method with those of some existing algorithms, such as the active contour model (ACM) [26], the watershed (Watershed) [27], region growing (RG) [28] and the level set (Level Set) [29], and with manual segmentation by two experts. All our implementations were programmed in the Microsoft Visual Studio 10.0 environment and executed on a personal computer equipped with a 3.40GHz Intel Core i7-3770 processor with 8 GB RAM. The software packages we used for medical image processing and 3D visualization are ITK 4.4.2 and VTK 6.1.0.

3.1 Qualitative evaluation

A solitary pulmonary nodule (SPN) is one of the most common types of pulmonary nodules. In this paper, for a series of lung CT image sequences with solitary pulmonary nodules, we use our algorithm and the ACM, watershed, RG and level set algorithms for lung segmentation. We must set up some of the necessary parameters to ensure the accuracy and effectiveness of the segmentation method. The values of these parameters are shown in .

Because of the large number of image sequences, we select five lung CT images from the top to the bottom of the lung in a dataset with SPN and then use one image out of every twelve to demonstrate the process and the results of lung image segmentation.

The process and results of our method are shown in . Column (b) is the ROI extraction result of the original lung CT image sequences column (a); column (c) is the result using the SGLIC algorithm for superpixel segmentation; column (d) is the result of using the SGNF algorithm for clustering, which is optimized by the genetic algorithm; column (e) and (f) are the coarse and the final lung parenchyma mask; and columns (g) and (h) are the final segmentation results of our method and the artificial segmentation.

The segmentation results of our proposed method.

Column (a) and (b) are five original lung CT images and ROI images from the top to the bottom; (c) and (d) shows the results of GSLIC and SGNF; (e) and (f) are the coarse and the final lung parenchyma mask; (g) and (h) present the final results of the proposed method and manual segmentation.

When using the RG, watershed and active contour model algorithms for sequence image segmentation, we set up the left and right lung seed points and select five lung CT images based on the result of the experiment. In this paper, the coordinates of these from top to bottom are (235, 272) and (206, 276); (212, 277) and (331, 260); (187, 249) and (341, 251); (194, 245) and (343, 265); and (195, 304) and (333, 327), which will be used to segment the left and the right lung parenchyma images. The process and results of using the watershed and RG segmentation algorithms are shown in and . We have observed that when the watershed algorithm is used to segment the images, different values will have different results (, Column (b)—(d)). Compared with the segmentation results of the level values of 0.05, 0.1 and 0.15, the best level value of the best segmentation results is 0.15, and the coarse and final segmentation results are shown in , Column(e) and (f). In addition, when using the RG algorithm, the of the best segmentation result is between 100 and 120 (, Column (f)).

The segmentation results of using the watershed method.

Column (a) is the five original lung CT images from top to bottom; (b)-(d) show the results of the level values of 0.05, 0.1 and 0.15, respectively; and (e) and (f) present the final results using the watershed algorithm and manual segmentation.

The segmentation results using the RG (region growing) method.

Column (a) is the five original lung CT images from top to bottom; (b)-(d) show the mask of the left, right and whole lung; (e) are the coarse lung parenchyma mask; (f) and (g) present the final results of using RG algorithm and manual segmentation.

In , for the CT images with SPN, we give the comparison of the artificial segmentation results of the 5 methods. The best segmentation result is obtained by using the ACM to segment, the we select is 5 pixels, and the value is 2.0 (, Column (d)). When we use the level set algorithm for segmentation, the best segmentation result is obtained when the is 100 and the is 500 (, Column (g)). shows the front and back of the lung in the 3D reconstruction of the lung parenchyma image sequences segmentation results with our proposed method using VTK.

The comparison of the final segmentation results lung parenchyma scans with malignant SPN.

Column (a) is the five original lung CT images from top to bottom, and (b)-(g) show manual segmentation; our proposed method; and the ACM, watershed, RG and level set methods, respectively.

The results of the lung parenchyma image 3D reconstruction.

(a) and (b) show the front and back of the whole lung.

For the other three types of lung sequence images, we also compare the segmentation results of these five segmentation methods without nodules (), with benign nodules (), and with pleural nodules ().

The comparison of the final segmentation results of lung parenchyma scans without nodules.

Column (a) is the five original lung CT images from top to bottom, and (b)-(g) show manual human segmentation; our proposed method; and the ACM, watershed, RG and level set methods, respectively.

The comparison of the final segmentation results of lung parenchyma scans with benign nodules.

Column (a) is the five original lung CT images from top to bottom, and (b)-(g) show manual human segmentation; our proposed method; and the ACM, watershed, RG and level set methods, respectively.

The comparison of the final segmentation results of lung parenchyma scans with pleural nodules.

Column (a) is the five original lung CT images from top to bottom, and (b)-(g) show manual human segmentation; our proposed method; and the ACM, watershed, RG and level set methods, respectively.

When the experimental results are compared with the results of manual segmentation, it is shown that the method in this paper has the best segmentation results, particularly for the segmentation of irregular lung images, such as the top and the bottom of the lung and the pleural nodules in CT image sequences, and the advantages are more obvious.

For the lung sequences images without nodules (), the segmentation results of our method (, Column (c))and the region growing method (, Column (f)) have good segmentation results, but the segmentation results of the watershed (, Column (e)) and level set algorithms (, Column (g)) are relatively poor.

For the lung image sequences with benign nodules (, Row 2), the ACM and watershed algorithm will lose some of the lung parenchyma (, Column (d) and (e)); the RG and level set algorithms will lose some of the pulmonary nodules (, Column (f) and (g)); and the method in this paper has the best segmentation results (, Column (c)).

For the lung sequence images with malignant SPN (, Row 2), the algorithm in this paper can ensure the integrity of the segmentation of the lung parenchyma (, Column (c)), while ACM will lose some of the lung parenchyma (, Column (d)) and the watershed, RG and level set algorithms (, Column (e)—(g)) will lose some of the pulmonary nodules.

For the lung sequence images with pleural nodules (, Row 4), the ACM, watershed, RG and level set algorithms will miss the retraction part of the pulmonary and pleural nodules (, Column (d)—(g)), while the method in this paper can guarantee much of the retraction part of the pulmonary and pleural nodules (, Column (c)) and is the closest to the manual segmentation results (, Column (b)).

In addition, for all CT images from the top to the bottom of the lung, the method in this paper is the only one that can ensure the integrity of the segmentation. Therefore, our proposed method has a better segmentation result and a higher generality on segmentation of lung parenchyma images.

3.2 Quantitative comparisons

Quantitative evaluation has significant importance in objectively assessing the effectiveness of an algorithm. The probabilistic rand index (PRI) [30], variation of information (VoI) [31] and Jaccard similarity coefficient [32] (Kim et al., 2005) are used to objectively assess the performance of the proposed algorithm.

Assuming that the original lung image S contains M pixels, the referential and actual segmentation results are expressed as and , respectively, and the following conditions theoretically should be have met (12): where K and N are the number of segmented regions in the referential and the actual segmentation results, respectively.

The probabilistic rand index (PRI) is a parameter to evaluate the consistency of attribute symbiosis between the actual segmentation results and the reference.

For a pixel pair (x, x) in the original lung image S marked (s, s) with the same attributes in the referential segmentation result , which should be the same in , the value of PRI [30] can be calculated as (13): where is a discriminant function that is used to determine whether the pixel pair has the same label. The value of PRI is in the range of [0, 1], and the larger the value, the better the result.

The variation of information (VoI) [31] is a measure of information content that depicts how much one segmentation reflects the information of the other segmentation. It is the conditional entropy among the distributions of the segments labels. Therefore, the VOI value can be calculated as: where H(S) and H(S) represent the entropy, and I(S, S) represents the mutual information. H(S) and I(S, S) also can be calculated as (15) (16) and (17).

The VoI values lie in [0, ∞). The 0 indicates that the two segmentations match perfectly. The smaller the value of VoI, the less information changes and the better the results will be.

The Jaccard similarity coefficient (Jaccard) is a measure to compare the similarity between the sample sets, which can indicate the coincidence degree of two images. The value of Jaccard can be calculated as (18).

The Jaccard values lie in [0, 1], and a higher Jaccard similarity coefficient indicates a better segmentation result.

shows the average scores of the PRI, VoI and Jaccard measures for the five algorithms on four types of lung CT image sequences. It is clear from that the proposed method outperforms the other state-of-the-art algorithms in terms of PRI, VoI and Jaccard.

For the lung sequences images of without nodules, five methods have better segmentation results. The average PRI and Jaccard values of the RG algorithm are close to the proposed method but the average VoI values are much greater. The watershed algorithm has a minimum PRI value of 0.9474, while the level set algorithm has a maximum VoI value and the lowest Jaccard value. Moreover, for lung sequence images with benign nodules, malignant SPN and pleural nodules, the RG algorithm’s performance drops rapidly with the highest VoI value, while the watershed and level set algorithms are relatively stable. Our method has better results in term of PRI, VoI and Jaccard values.

In general, our proposed method has the best segmentation performance in terms of the parameter comparison of the five types of image segmentation. From lungs without nodules to lungs with pleural nodules, our method’s performance declines slightly. The main cause perhaps is that the lung parenchyma images contain lung nodules, and the types of lung nodules are becoming more and more complicated. It is clear from that the RG algorithm is the most sensitive to lung nodules. Unless otherwise stipulated in the image sequences without nodules, the RG algorithm is very close to the proposed method, but once lung nodules are included, there is a drastic decline in the indicators.

Moreover, the performance of each of the five algorithms in terms of PRI, VoI and Jaccard values on each image with malignant SPN (a) and without nodules (b) is graphically represented in , and . It is clear from that the proposed method performs better than the other methods in terms of the PRI value, the RG algorithm is close to the proposed method and better than the other algorithms. shows that our proposed method can always keep the global minimum value of VOI, while the level set value nearly approaches ours. Both the Watershed and RG algorithm have higher VoI values. It is obvious that our Jaccard value is higher than that of RG and is far better than for the other methods.

The PRI values of the five methods on the segmentation results of lung parenchyma scans with malignant SPN (a) and without nodules (b).

The VOI values of the five methods on the segmentation results of lung parenchyma scans with malignant SPN (a) and without nodules (b).

The Jaccard values of the five methods on the segmentation results of lung parenchyma scans with malignant SPN (a) and without nodules (b).

In consequence, based on the comparison and analysis of the three measures PRI, VoI and Jaccard, our method is the most close to the artificial segmentation results, which can further reflect our method’s high preformation and wide generality in the segmentation of lung parenchyma images.

For the four types of lung parenchyma images, we also analyzed the time performance of the five methods, as shown in , and the average processing times for the five methods are shown in . The “Average Dataset size” row signifies the average number of lung CT images in a dataset. It is clear from that the average processing time for each dataset using our method is 42.21 seconds; i.e., it will take 0.71 seconds to process a single slice, which is far better than in other four methods. Therefore, the proposed method has obvious advantages over the other methods in terms of segmentation speed of lung CT images.

Conclusion

Our work indicates that our proposed method can segment various types of lung parenchyma image sequences effectively. This method is more accurate and universally applicable than any of the traditional methods. Based on the segmentation of the four different types of sequences of lung CT images, which included 4812 images from 80 datasets, we compare the results of our method and those of the existing algorithms with manual segmentation. The experimental results show that our method can achieve accurate segmentation of the lung parenchyma and in particular accurate segment the lung CT images, which have complex morphological structures such as the top and bottom of the lung and contain pulmonary nodules. Our method can achieve an average volume pixel overlap ratio of 92.22 ± 4.02% for the four types of lung parenchyma image sequences. Moreover, our method is less time consuming, with an average processing time of 42.21 seconds for each dataset, meaning it takes approximately 0.71 seconds to process a single slice. Therefore, in the segmentation of lung parenchyma image sequences, taking the high correlation between adjacent slices of CT image sequences into consideration can significantly improve the speed of segmentation, while superpixels can guarantee the quality and the post-processing of image segmentation, and the SGNF optimized by the GA can be more effective at maintaining the integrity of the lung parenchyma segmentation.

Table 1

GSLIC algorithm.

Algorithm 1 GSLIC Algorithm
1: Initialize cluster centers Ck = [Lk, ak, bk, xk, yk, zk]T by sampling pixels at regular grid steps.
2: Choose cluster centers in an n * n neighborhood, to the lowest gradient position.
3: Repeat
4: for each cluster center Ck in zk do
5: Assign the best matching pixels from an 2S * 2S square neighborhood around the cluster center according to the distance measure.
6: end for
7: Compute new cluster centers and residual error E{distance between previous centers and recomputed centers}.
8: until E < = threshold. Save cluster centers Ck and transmit Ck to the next image zk+1.
9: until all image sequences are segmented.
10: Enforce connectivity.

Table 2

SGNF optimized by GA.

Algorithm 2 SGNF Optimized by GA
1: Input the initial sample set {X}, and the K initial seeds are selected randomly.
2: Obtain the optimal K superpixel seeds to generate the initial SNGF with K neural tree using genetic algorithm.
3: Repeat
4: For each superpixel sample in {X} do
5: Generate neuron nj from {X} and traverse the SGNF to find a neuron nwin with the shortest distance to nj.
6: Connect nj to the SGNT.
If nwin is a leaf neuron in the current SGNT
Create a new neuron nj+1 and copy weight wwin to wj+1;
Connect nj and nj+1 to nwin as a child.
Else Connect nj to nwin as a child, update node weight.
7: End for.
8: Until all samples have been input into the SGNF.

Table 3

Genetic Algorithm.

Algorithm 3 Genetic Algorithm
1: Chromosome encoding.
2: Initialize the maximum iterations T and generate an initial population G (0) randomly.
3: Compute the fitness value φ(C) of each sample in initial population G (0).
4: If Δφ(C) < μ or iterations > T
5: Return to step 10.
6: Else perform genetic operators: selection, crossover and mutation.
7: Generate the next generation of population G (t+1).
8: Compute the fitness value φ(C) of each sample in initial population G (t+1).
9: Return to step 4.
10: Output the optimal solution Copt.

Table 4

Feature extraction and lung parenchyma segmentation.

Algorithm 4 Feature Extraction Algorithm
1: Assume that the final four types of sample sets are {Si}, i = 1, 2, 3, 4, then
∪i=14S_i=X,∀(1≤i,j≤4,i≠j)(S_i∩S_j)=∅
2: Calculate the average grayscale values φ(Si) for four sample sets using (9)
3: max(φ(Si)) → pleural tissue(S4)
4: Compute the centroid coordinates (xi,yi) and the coordinate variances ξ2(Si) of the three left sample sets {Si}using (10) and (11)
5: max(ξ2(Si)) → extrathoracic area(S3)
6: two left sample sets → the left(S1) and the right(S2) lung parenchyma
7: Traverse superpixel samples in {S1, S2}
8: For Xi,Xj ∈ S1 ∪ S2, If aj6 = ai6, Xi,Xj → N (The same CT image N whose slice number is ai6), Connect Xi and Xj to Nk.
9: Until all superpixel samples have been input.
10: Sequential output all Nk corresponding to each lung parenchyma image.

Table 5

Parameter values setting.

Methods	Parameters	Values
Our method	N, n, S, δ, E in GSLIC	1000, 3, 16, 10, 0.0001
	k, p, T, τ, η, μ in SGNF	3, 6, 800, 0.7, 0.001, 0.0001
RG	Threshold	100–120
Watershed	Level	0.05, 0.1 and 0.15
ACM	image spacing, expansion coefficient	5, 2.0
Level set	time threshold, stop time	100,500

Table 6

Average values of PRI, VOI and Jaccard for the five algorithms on four types of image sequences.

Types	Measures	ACM	Watershed	RG	Level set	GSLIC-SGNF
Without-	PRI	0.9581	0.9486	0.9611	0.9474	0.9624
nodules	VoI	1.779	2.283	1.754	1.4535	1.3815
	Jaccard	0.9375	0.9252	0.9385	0.9267	0.941
Benign-	PRI	0.8624	0.8964	0.9285	0.8695	0.9465
nodules	VoI	1.9532	2.369	2.7471	1.8382	1.6515
	Jaccard	0.8845	0.8972	0.9238	0.8758	0.9342
Malignant-	PRI	0.8948	0.9025	0.9094	0.8989	0.9169
SPN	VoI	1.9875	2.6475	2.8305	1.8975	1.8195
	Jaccard	0.8708	0.8636	0.902	0.8361	0.9257
	PRI	0.8556	0.8629	0.8505	0.8636	0.8926
Pleural-	VoI	3.825	3.323	4.152	2.8136	2.2052
nodules	Jaccard	0.8122	0.8535	0.7934	0.8242	0.8821

Table 7

Average processing time(s) for the five algorithms on CT image sequences.

Types	Average Dataset size	ACM	Watershed	RG	Level set	GSLIC-SGNF
Without nodules	512*512*60	254.57	169.26	83.35	117.56	43.22
Benign nodules	512*512*59	250.42	164.45	81.22	114.28	41.69
Malignant SPN	512*512*60	235.51	162.88	78.92	111.36	41.18
Pleural nodules	512*512*62	246.84	168.23	83.25	118.23	42.76