Pulmonary Lobe Segmentation in CT Images Based on Lung Anatomy Knowledge
In computed tomography (CT) images, pulmonary lobe segmentation is an arduous task due to its complex structures. To remedy the problem, we introduce a new framework based on lung anatomy knowledge for lung lobe segmentation. Firstly, the priori knowledge of lung anatomy is used to identify the fissure region of interest. Then, an oriented derivative of stick filter is applied to isolate plate-like structures from clutters for lobar fissure verification. Finally, a surface fitting model is employed to complete the incomplete fissure surface for lung lobe segmentation. Compared with manually segmented fissure references, the designed approach obtained a high median F1-score of 0.8865 in the left lung and obtained a high median F1-score of 0.9200 in the right lung. The average percentages of the segmented lung lobes in the lung lobe ground truth are 0.960, 0.989, 0.973, 0.920, and 0.985 for the left upper, left lower, right upper, right middle, and right lower lobes, respectively. The perfect performance of the proposed scheme is tested by visual inspection and quantitative evaluation.
In medical practice, recognition of lung lobes is useful for clinical diagnosis and lung disease assessment . Hence, accurate segmentation of lung lobes is urgently needed. Anatomically, the left lung is composed of two lobes, whereas the right lung is composed of three lobes . Clinically, lung lobe segmentation is extremely important in surgical treatment of lung diseases [3, 4]. However, manual segmentation of lung lobes is impracticable in CT images due to the large number of images. Therefore, automatic segmentation methods for lung lobes are necessary.
A schematic diagram is shown in Figure 1; the human lungs are composed of five lobes. The right lung is subdivided into the right upper lobe (RUL), the right middle lobe (RML), and the right lower lobe (RLL), separated by the right oblique fissure and the right horizontal fissure, whereas the left lung is subdivided into the left upper lobe (LUL) and the left lower lobe (LLL), separated by the left oblique fissure. Lung lobes are served by separate bronchial trees and their corresponding pulmonary arteries and are anatomically independent. Unfortunately, lung lobe segmentation is an arduous task due to complex lung anatomy.
Pulmonary lobe segmentation algorithms can be divided into three categories. The first approaches consist of methods that search for the fissure surface based on shape features and gray-level information. The second approaches consist of methods that use this information from airways, vessels, and fissures to estimate lobar boundaries. The third approaches consist of methods that construct pulmonary atlases to guide lung lobe segmentation. Although many methods have been presented, no standard method has yet emerged .
Considering the fact that fissures are physical boundaries between different lobes, Pu et al.  combined a threshold with a marching cube algorithm to identify fissures and utilized the direction difference for fissure decomposition; then, an implicit surface fitting model was applied for pulmonary lobe segmentation. Shamonin et al.  designed a fissure enhancement filter to distinguish between fissures and noise; then, an iterative B-spline fitting model was applied to segment pulmonary lobes. However, their algorithms [6, 7] may not be able to accurately segment pulmonary lobes in some cases due to their fissure decomposition error and confounding structures. For the same purpose, Ross et al.  employed a particle system and ridge surface sampling to identify lobar fissures; then, a thin plate spline was used to segment pulmonary lobes. The approach fails if the particle signal is extremely weak. Using a different strategy, Lo et al.  used a fissure confidence function to verify pulmonary fissures; then, a surface evolution algorithm was applied for pulmonary lobe segmentation. Lassen-Schmidt et al.  firstly selected a number of representative lobar boundary points; then, a plane was constructed by the given points to segment pulmonary lobes, but their algorithms [9, 10] need manual interaction. Recently, Wei and Hu  used a texture analysis method based on the neural network and a stack projection to verify lobar fissures; then, a fissure surface model was used for pulmonary lobe segmentation. Qi et al.  used a line enhancement filter followed by a uniform cost search to verify pulmonary fissures; then, an implicit surface model was used to segment pulmonary lobes, but their methods [11, 12] may not be able to accurately verify pulmonary fissures. The reason is that the orientation of fissures can vary too much in some cases .
Motivated by the assumption that the information of airways and blood vessels may contribute to segment pulmonary lobes, Ukil and Reinhardt  used airways and vascular trees to construct a region of interest (ROI) encompassing fissures and applied a ridgeness measure to enhance fissures; then, a fissure decomposition technique and a minimum curvature biharmonic spline interpolation algorithm were employed to segment pulmonary lobes. In a similar way, Doel et al.  applied a similar approach to identify fissure ROI and developed a fissureness filter to highlight the fissure representation; then, a multilevel B-spline method was used for lung lobe segmentation. Their algorithms [13, 14] fail if pulmonary veins cross lobar fissures. Using a different strategy, Lassen et al.  constructed a cost image based on fissures, vessels, and bronchi for the watershed transformation to segment pulmonary lobes. Its improved version was presented by Bragman et al.  using a groupwise fissure prior approach for pulmonary lobe segmentation. Later, Bauer et al.  presented a subject-specific prior method for lung lobe separation. Recently, Giuliani et al.  proposed a lung lobe segmentation method using alpha-expansion. However, their methods [15–18] heavily depend on airway segmentation.
In addition, a number of studies constructed pulmonary atlases to guide lung lobe segmentation. Following this strategy, Zhang et al.  combined the ridgeness measure and an atlas to verify pulmonary fissures; then, a fissure surface was constructed for pulmonary lobe segmentation. Later, van Rikxoort et al.  improved the method and presented a multiatlas approach to segment pulmonary lobes. Recently, Schmidt-Richberg et al.  combined the level sets with atlases to segment pulmonary lobes, but their methods [19–21] have a longer computation time.
The proposed scheme aims at automatically segmenting pulmonary lobes. In order to identify fissure ROI, many studies [13, 14, 22, 23] used bronchial trees and pulmonary vessels to approximately construct lobar shapes. Their algorithms [13, 14, 22, 23] fail in some cases because of pulmonary veins may cross fissures. Different from the algorithms [13, 14, 22, 23] which have been previously proposed, we present a novel and effective framework to automatically identify pulmonary fissure ROI. For subsequent segmentation, we improved and extended our previous works [24–26] to accurately verify pulmonary fissures. Finally, a surface fitting model [6, 11] is used to complete the incomplete fissure surface for lung lobe segmentation. Experimental results show that the proposed method has a good performance in pulmonary fissure segmentation and lung lobe segmentation.
2. Materials and Methods
2.1. Data and References
Fifteen CT images were used to evaluate the validation of our scheme and acquired from the LObe and Lung Analysis 2011 (LOLA11) challenge , including various pathologic abnormalities. To better evaluate the segmented right fissures and left fissures, we divided CT images into left lungs and right lungs.
A manual segmentation of lung lobes in 9 coronal slices was made by two human observers. Unfortunately, the ground truth is not publicly available. To make full use of the publicly available data for experimental evaluation, we manually segmented 15 lung lobes and regarded them as the lung lobe ground truth. In addition, lobar fissure segmentation is different from lung lobe segmentation . Therefore, two internationally renowned experts verified the LOLA11 reference  and removed parts of the invisible reference and added parts of the visible reference that were lost in the LOLA11 reference. In this study, we regarded the fissure reference  as the fissure ground truth for experimental evaluation.
2.2. Overview of the Proposed Method
Our scheme (developed in MATLAB and C++) uses a new framework for lung lobe segmentation. The designed algorithm employs three stages: (1) fissure ROI identification, (2) lobar fissure verification, and (3) lung lobe segmentation. Figure 2 depicts the flowchart of this computerized scheme. The following context presents these three stages.
2.3. Fissure ROI Identification
In this paper, we present a new computerized scheme for fissure ROI identification. Firstly, an effective method based on the graphic processing unit (GPU) is used to extract the centerline of airway trees , and a graph search approach is employed to obtain the subtrees in airway trees. Then, a shape prior and graph cut method  is applied to separate pulmonary arteries and veins. Furthermore, the priori knowledge of airways and pulmonary arteries is fused to segment pulmonary arteries. Finally, an alpha-shape algorithm [31, 32] is applied to identify fissure ROI.
2.3.1. Airway Segmentation
To accurately segment airways, a useful and valuable method  based on GPU is applied for tube extraction. Firstly, Gaussian smoothing with the parameter = 0.5 is used to blur the dataset I. Subsequently, a gradient vector field is designed:
Here, controls the normalization. To prevent different tubular structures diffusing together, a gradient vector flow (GVF) method  is used:where and denote the initial gradient vector field and the weighting constant. Xu and Prince  proposed an iterative approach to solve the Euler equation:
The purpose of the GVF model is to make the following tube detection filter (TDF) [34, 35] invariant to contrast of airways . Subsequently, the TDF is applied to detect tubular structures. Firstly, three eigenvectors of the Hessian matrix corresponding to the tube are estimated: one of three eigenvectors is related to the direction of the tube, and the other two eigenvectors and are associated with the cross-sectional plane of the tube. The main idea is to create a circle with a radius in the cross-sectional plane, and N is the number of spaced points i on the circle. Therefore, the direction from the center to each point i is
The average dot product is
When is increased, the average dot product is calculated again. In this paper, ranges from 0.5 to 25. Lastly, the airway centerline is extracted by traversing a ridge in the TDF results. All values in TDF results above a threshold are extracted and treated as candidate starting points, starting from the candidate point with the positive tangent direction and negative . For a current centerline point and the corresponding tangent direction , all local neighbors are considered:
The neighbor point with the highest TDF response is chosen as the next centerline point . Finally, the airway centerline is acquired and marked with S.
Different subtrees in the airway tree belong to different lobes . Based on this theory, a graph search method is used to search for major bifurcations for bronchial tree separation in appropriate orientations in Figure 3. In the first step, a centerline extraction method  is used to transfer the airway tree segmentation to its skeleton representation. Then, a directed graph is modeled from the airway tree with the trachea as the root. To extract the bifurcations, a simple but effective criterion [37, 38] is used:where p and N26 denote the pixel in the centerline of the airway tree S and its 26-neighborhood area, respectively. If there are 4 or more points within a 3 3 3 cube (p and its 26 neighborhood), we consider the pixel p as a bifurcation. Finally, a graph search method is used to remove major bifurcations labeled with yellow dots in Figure 3(b) and calculate for each subtree in Figure 3(c) along with the particular separation direction in Figure 3(b), in which EndRMB and LMB denote the ends in the right main bronchus and the left main bronchus, respectively.
2.3.2. Pulmonary Artery and Vein Separation
Airways are tubular structures and less adhered with each other. Furthermore, the GVF model has the advantage that it is feature preserving. Therefore, a method based on the GVF model is used to segment the airway tree, but the approach is not applicable for pulmonary artery and vein separation. As shown in Figure 4, pulmonary arteries may adhere with pulmonary veins, and pulmonary veins may cross lobar fissures (marked with yellow arrows). To remedy the problem, a shape prior and graph cut method  is applied to separate and segment multiple interwoven vessel trees.
Pulmonary artery and vein separation is a challenging task due to its complex structure. To solve the problem, the shape prior and graph cut method  is applied to segment the interwoven three-dimensional (3D) tubular trees. The first step is to identify all tubular objects by TDF. Three eigenvectors of the Hessian matrix corresponding to tubular objects are estimated: one of the three eigenvectors is related to the direction of the tube, and the other two eigenvectors and are associated with the cross-sectional plane of the tubular objects. The main idea is to create a circle with a radius in the cross-sectional plane, and N is the number of spaced points i on the circle. Therefore, the direction from the center to each point i is
To segment the interwoven 3D tubular trees, a medialness method  is applied to separate pulmonary arteries and veins. The initial boundariness contribution and averaging boundariness contribution arewhere N represents the number of spaced points on the circle, is the circle radius, is the tube direction, and the gradient vector field is with . The variance can be calculated as follows:
Thus, the response of the final offset medialness is
To suppress the response away from the tube center, Bauer et al.  proposed a method based on an adaptive threshold:
The multiscale medialness response is
Here, and . After acquiring the multiscale medialness response, we set a threshold and take the medialness response larger than the threshold as seed points. Finally, a tree construction algorithm is used to separate pulmonary arteries and veins. The main idea is to extract the centerline of the multiscale medialness response by traversing a ridge map.
Pulmonary arteries may adhere with pulmonary veins in Figure 5(a). To separate the multiple interwoven vessel trees, we use the medialness method to extract the centerline of vessel trees, and the results are shown in Figure 5(b). As observed in the regions marked with yellow arrows in Figure 5, pulmonary arteries and veins have been separated.
2.3.3. Fissure ROI Identification
Based on the reality that pulmonary arteries are close to bronchial trees [39–41], we combine the bronchial trees (BT) with vessel subtrees (VST) to remove mounts of pulmonary veins. We dilate BT and use a connected-component analysis method (CCAM) to select 5 biggest objects for pulmonary artery segmentation. To give a mathematical expression,where D denotes the morphological dilation.
To illustrate the validation of the proposed scheme, we integrate with bronchial trees, pulmonary arteries, and lobar fissures to describe the distribution map from different angles of view. Bronchial trees and lobar fissures are labeled with yellow color and green color in Figure 6, respectively. It can be seen that pulmonary arteries are close to bronchial trees [40, 41], and pulmonary arteries cannot cross lobar fissures [15, 39].
As a result, we consider different bronchial trees and the corresponding pulmonary arteries belonging to different pulmonary lobes as a whole in Figure 7(a). Then, an alpha-shape algorithm [31, 32] is used to construct lobe shapes. The alpha-shape algorithm was presented by Edelsbrunner and Mücke  that constructs an entire family of shapes for a given set of points. We adopt for the implementation of alpha-shape algorithm to construct the geometric shape of different bronchial trees and their corresponding pulmonary arteries. Finally, surface fitting processing and region filling processing are used to approximately estimate the lobe distribution map for fissure ROI identification in Figure 7(b).
Fissure ROI identification is useful in lung lobe segmentation. As shown in Figure 8(a), the estimated right upper, middle, and lower lobe distribution maps are marked with , , and , respectively. We just regard and as an entirety and as another part; then, the alpha-shape algorithm is applied to identify the right oblique fissure ROI, and its corresponding results are shown in Figure 8(b).
2.4. Lobar Fissure Verification
To effectively verify pulmonary fissures, an oriented derivative of stick (ODoS) filter  is applied to enhance pulmonary fissures. Then, an orientation partition scheme is used to separate pulmonary fissures and adhering clutters. Next, a connected-component analysis approach is adopted to remove noise for plate-like structure extraction. Finally, the identified fissure ROI is employed to isolate lobar fissures from adhering clutters.
2.4.1. Retrospection of the DoS Filter and ODoS Filter
To efficiently detect the weak and thin pulmonary fissures, Xiao et al.  presented a derivative of stick (DoS) filter for pulmonary fissure enhancement. The filtering kernel consists of three parallel sticks (LS, MS, and RS), and θ is the stick orientation. Two 2D line strength measures and arewhere , , and denote the intensity of the jth pixel along the left, middle, and right stick, represents a positive coefficient, and is an expected value operator. To suppress noise and other structures, Xiao et al.  integrated the line strength measures and in multidirections:
The DoS filtering response can be formulated as
represents the cascading operator. To integrate the DoS filtering response in multidirections, Xiao et al. defined a shape-tuned fissureness function:
Here, , , and denote the DoS filtering response across the axial, sagittal, and coronal sections. However, the original DoS filter cannot efficiently suppress adhering clutters. To solve the problem, Peng and Xiao  improved and extended it by constructing its orientation response. The orientation response is
A vector representation is
Therefore, the response of the ODoS filter can be described as follows:
2.4.2. Plate-Like Structure Extraction
To suppress noise and other structures, Peng and Xiao  firstly used a threshold T to modulate the ODoS filtering response:where the threshold T we set is equal to 1. Based on the fact that pulmonary fissures appear as plate-like structures and have similar directions, Peng and Xiao  presented an orientation partition scheme to separate lobar fissures and noise in Figure 8. They divided the enhanced 3D CT image into n overlapped subregions, and the subregions are marked with and Rn. To give a mathematical expression,
In 3D partitioned images , clutters such as airways and vessels usually appear as small objects. As a result, they used a connected-component analysis method to select the candidate fissures, and the results are displayed in Figure 9.
2.4.3. Oblique Fissure and Horizontal Fissure Verification
From the above processing, some plate-like structures are easily acquired, including lobar fissures and accessory fissures. As is shown in Figure 10, accessory fissures have bad influence on fissure segmentation, but few methods have been proposed by using the priori knowledge of pulmonary structures to remove them. In this study, the identified fissure ROI is used to limit the search area and remove those unrelated accessory fissures for right fissure verification. Similarly, we can use the previously identified fissure ROI to verify right oblique fissures and left oblique fissures.
After acquiring the right fissure (RF) and the right oblique fissure (ROF), we use a simple but effective criterion to verify the right horizontal fissure (RHF):
As shown in Figure 11, it can be seen that the proposed method has a good performance in lobar fissure verification.
2.5. Lung Lobe Segmentation
Lung lobe segmentation is an arduous task due to incomplete, deformed, and disrupted fissures. To remedy the problem, a surface fitting model [6, 11, 12] is applied to construct a complete surface. Subsequently, the complete surface is employed to segment lung lobes.
2.5.1. Surface Fitting Model
Once oblique fissures and horizontal fissures are verified, we modeled a fissure surface in the form ofwhere x, y, and z are the coordinate points in the 3D space. Imaging a thin and flexible plate for a given set of points by elastic bands is a nontrivial task. Each of these elastic bands draws the plate towards the given set of points, stretching only in the Z-direction. The given set of points was to fall on a simple planar surface [6, 11]:
A matrix expression is
At every point of the surface, the surface fitting model attempts to force the partial derivatives of the surface in neighboring cells to be equal. A mathematical expression in the second set of linear equations is of the form
To construct a smooth fissure surface, a Laplacian function is treated as the regularization term. It can be described aswhere α is a scale factor that controls the effect of regularization. Then, we solve for the vector W by performing the following minimization:
Therefore, the surface fitting model is employed to construct a complete fissure surface, as shown in Figure 12.
2.5.2. Cutting Model
To segment pulmonary lobes, we present a cutting model to subdivide human lungs into five lobes. As shown in Figure 13, the right complete oblique fissure surface is used to divide the right lung into the special part and the lower lobe, and then the right complete horizontal fissure surface is employed to separate the special part into the upper and middle lobes. In the same way, we use the left complete fissure surface to divide the left lung into left upper and lower lobes. As a result, the human lung is subdivided into five lobes.
3. Experimental Results
The proposed scheme is tested on 15 clinical CT images. All experiments are performed in 64 bit Windows 10 operating system, running on a computer with an Intel (R) Xeon (R) E5-1607 v2 3.00GHZ CPU and 20 GB RAM. The runtime of the proposed method for a 400 × 512 × 512-size 3D image is approximately 1580 s. In addition, the 3D visualization results are rendered by MeVisLab software .
3.1. Evaluation of Pulmonary Fissure Segmentation
3.1.1. Evaluation Criteria
The manual reference and its adjacent 3 mm region are labeled with , and the segmented fissure and its adjacent 3 mm region are labeled with . The voxels of the segmented fissure inside are considered true positive (), and the rest are false positive (). Similarly, the voxels of the manual reference inside are considered , and the rest are false negative (). As a result, the false discovery rate (), false negative rate (), and F1-score () are given as follows:
3.1.2. Visual Inspection
As shown in Figure 14, the segmentation results of the proposed method are shown in Figure 14(a). For comparison, the ODoS  and DoS  filtering segmentations are, respectively, displayed in Figures 14(b) and 14(c). The segmented results, the manual reference, and their overlapped regions are, respectively, rendered in green, yellow, and purple. In the areas labeled with red ellipses in Figure 14, the weak objects are treated as the undersegmentation error, which may lead to a high FNR index. As observed, the proposed method has a good performance in weak object detection.
3.1.3. Quantitative Evaluation
The proposed method is validated on 15 CT examinations in Figure 15. The box plots of indices (F1, FDR, and FNR) corresponding to the proposed method (p), ODoS (od), and DoS (d) filtering scheme are drawn next to each other. In the left lung, the median values of box plots are 0.8865, 0.8828, 0.8374, 0.0755, 0.0489, 0.1971, 0.0978, 0.1368, and 0.1478. In the right lung, the median values of box plots are 0.9200, 0.8997, 0.8271, 0.0381, 0.0344, 0.1544, 0.0933, 0.1221, and 0.1906. Both visual inspection and quantitative evaluation exhibited that our scheme can outperform the compared methods [25, 26].
3.2. Evaluation of Lung Lobe Segmentation
3.2.1. Evaluation Criteria
The volumes of the lung lobe ground truth and the segmented lung lobes are, respectively, marked with and . The volumetric overlap results between and are marked with . We regard the percentage as the evaluation criteria.
3.2.2. Visual Inspection
To provide a global impression, four representative segmentations are chosen from the LOLA11 dataset. As shown in Figure 16, the upper, middle, and lower lobes are rendered in green, red, and yellow, respectively. The first, second, and third rows represent the visual lobar boundaries, the proposed method, and the constrained interpolation profile method . As observed in the areas labeled with red arrows, the proposed method has a good performance in lung lobe segmentation.
3.2.3. Quantitative Evaluation
Table 1 presents the evaluation results of lung lobes in 15 patients. The average percentages of the segmented LUL, LLL, RUL, RML, and RLL in the lung lobe ground truth are 0.960, 0.989, 0.973, 0.920, and 0.985, respectively. Lung lobe segmentation is useful in clinical diagnosis, especially in rapid diagnosis of new coronary pneumonia , in which we divide all patients into two parts: normal and abnormal groups. In abnormal groups, parts of bronchus and fissures are missed due to various pathological abnormalities. As shown in Figure 17(a), the disrupted bronchus has no power to guide and segment pulmonary arteries, and it may cause parts of clutters not to be removed. Another limitation of our scheme is drawn on incomplete fissures. As observed in Figure 17(b), incomplete fissures are used to construct complete lobe boundaries by using the surface fitting model; it may lead to a wrong lobar segmentation. Furthermore, the results can be seen in Table 1, and the abnormal cases have low volumetric overlap values with bold type.
In this paper, a unique scheme based on lung anatomy knowledge is presented for lung lobe segmentation. The described scheme has many benefits. First, the priori knowledge of airways and pulmonary arteries is used to identify fissure ROI. Many methods [13, 14, 22, 23] use only pulmonary airways and vessels to identify fissure ROI. Second, the identified fissure ROI and our previous works [25, 26] are combined to accurately verify oblique fissures and horizontal fissures. While parts of methods use the orientation information to verify oblique fissures and horizontal fissures, this will influence the results where parts of pulmonary fissures cannot be detected because the orientation can vary too much between patients . Third, the proposed algorithm can effectively eliminate accessory fissures, which is a well-known challenge about pulmonary fissure detection [24–26, 45]. Fourth, the surface fitting model  is applied to complete the incomplete surface of pulmonary fissures. The approach ensures the local continuity of the estimated lobar boundaries. Last, the computerized scheme has a good performance in lobar fissure detection and lung lobe segmentation.
Compared with manually segmented fissure references, our scheme achieved a median F1-score of 0.8865 in the left lung and a median F1-score of 0.8920 in the right lung. The perfect performance of our computerized scheme can be exhibited from the lowest FNR index, which indicates the completeness of lobar fissure detection has been mended by our scheme. The major contribution is that the described scheme applied the priori knowledge of airways and pulmonary arteries to accurately identify fissure ROI and used a lower threshold to suppress clutters. While Xiao et al.  applied a sophisticated postprocessing pipeline and a volume sorting scheme for lobar fissure detection, a part of fissures was considered as clutters and removed. Using a different strategy, Peng et al.  used an appropriate threshold to suppress clutters and applied the volume sorting scheme to segment lobar fissures. However, the DoS and ODoS methods [25, 26] often faced a tradeoff between weak object detection and noise suppression.
Our scheme is validated on 15 CT examinations. The average percentages of the segmented lung lobes in the lung lobe ground truth are 0.960, 0.989, 0.973, 0.920, and 0.985 for the left upper, left lower, right upper, right middle, and right lower lobes, respectively. The computerized scheme has a good performance in lung lobe segmentation. This is attributed to an ingeniously designed fusion of the surface fitting model, identified fissure ROI, and the oblique and horizontal fissure verification.
In spite of the aforementioned benefits, the described scheme has many limitations. First, our scheme is heavily depending on airway segmentation, which is an arduous task and sensitive to the image quality. Second, the lung anatomy knowledge is employed to segment pulmonary fissures, which is time consuming. Third, our scheme is derived from the ODoS filter and inherits some drawbacks. For example, some deformed fissures cannot be detected by the computerized scheme due to their orientation distribution. Fourth, poor segmented airways and pulmonary arteries may cause parts of fissures to be undetected in some cases. Finally, the fitting accuracy of the estimated lobe boundaries is largely determined by verified fissures; a better segmented method of fissure detection and lung lobe segmentation is needed to pursue.
In this paper, we present a new lung lobe segmentation scheme. Motivated by the fact that pulmonary veins may cross pulmonary fissures, the priori knowledge of airways and pulmonary arteries is used to identify pulmonary fissure ROI. The second contribution of our scheme emphasizes on pulmonary fissure segmentation. Using an ingeniously designed framework, adhering clutters are reasonably removed. Finally, the surface fitting model, fissure line strength, and orientation are fused to complete the incomplete fissure surface for lobar boundary estimation. The proposed method is validated on 15 CT examinations. Experimental results showed that the described method has a good performance in pulmonary fissure detection and lung lobe segmentation.
The LOLA11 dataset can be downloaded from https://lola11.grand-challenge.org/.
Conflicts of Interest
The authors declare no conflicts of interest.
Yuanyuan Peng and Lan Peng conceptualized the study, developed the methodology, and contributed to data curation and software. Hualan Zhong, Zheng Xu, and Hongbin Tu validated the study. Xiong Li performed formal analysis and supervised the study. Yuanyuan Peng wrote the original draft. Yuanyuan Peng and Hongbin Tu reviewed and edited the article. Hualan Zhong and Zheng Xu administered the project. Yuanyuan Peng, Hongbin Tu, and Xiong Li were responsible for funding acquisition. All authors read and agreed to the published version of the manuscript.
This research was funded by the Jiangxi Provincial Natural Science Foundation (nos. 20202BAB212004, 20204BCJL23035, 20192ACB21004, and 20181BAB202017), the Educational Science Research Project of China Institute of Communications Education (no. JTYB20-33), the Scientific and Technological Research Project of Education Department in Jiangxi Province (no. GJJ190356), and the Science and Technology Project of Changsha City (no. kq2001014).
K. Cederlund, U. Tylén, L. Jorfeldt, and P. Aspelin, “Classification of emphysema in candidates for lung volume reduction surgery,” Chest, vol. 122, no. 2, pp. 590–596, 2002.View at: Publisher Site | Google Scholar
Q. Li, L. Chen, X. Li, S. Xia, and Y. Kang, “An improved random forests approach for interactive lobar segmentation on emphysema detection,” Granular Computing, vol. 12, pp. 1–10, 2019.View at: Google Scholar
H. Nakamura, Y. Taniguchi, and T. Haruki, “Securing and treatment of interlobar pulmonary artery in severe insufficient interlobar fissure of the lung,” The Japanese Journal of Thoracic Surgery, vol. 72, no. 10, pp. 829–833, 2019.View at: Google Scholar
A. R. Nene, K. S. Gajendra, and M. V. R. Sarma, “Lung lobes and fissures: a morphological study,” Anatomy (International Journal of Experimental and Clinical Anatomy), vol. 5, no. 1, pp. 30–38, 2011.View at: Publisher Site | Google Scholar
T. Doel, D. J. Gavaghan, and V. Grau, “Review of automatic pulmonary lobe segmentation methods from CT,” Computerized Medical Imaging and Graphics, vol. 40, pp. 14–29, 2015.View at: Publisher Site | Google Scholar
J. Pu, B. Zheng, J. K. Leader et al., “Pulmonary lobe segmentation in CT examinations using implicit surface fitting,” IEEE Transactions on Medical Imaging, vol. 28, no. 12, pp. 1986–1996, 2009.View at: Publisher Site | Google Scholar
D. P. Shamonin, M. Staring, M. E. Bakker et al., “Automatic lung lobe segmentation of copd patients using iterative B-spline fitting,” International Society for Optics and Photonics, vol. 20, 2012.View at: Google Scholar
J. C. Ross, G. L. Kindlmann, Y. Okajima et al., “Pulmonary lobe segmentation based on ridge surface sampling and shape model fitting,” Medical Physics, vol. 40, no. 12, p. 121903, 2013.View at: Publisher Site | Google Scholar
P. Lo, E. M. van Rikxoort, F. Abtin et al., “Automated segmentation of pulmonary lobes in chest CT scans using evolving surfaces,” International Society for Optics and Photonics, vol. 20, 2013.View at: Google Scholar
B. C. Lassen-Schmidt, J.-M. Kuhnigk, O. Konrad, B. van Ginneken, and E. M. van Rikxoort, “Fast interactive segmentation of the pulmonary lobes from thoracic computed tomography data,” Physics in Medicine & Biology, vol. 62, no. 16, pp. 6649–6665, 2017.View at: Publisher Site | Google Scholar
Q. Wei and Y. Hu, “A hybrid approach to segmentation of diseased lung lobes,” IEEE Journal of Biomedical and Health Informatics, vol. 18, no. 5, pp. 1696–1706, 2014.View at: Publisher Site | Google Scholar
S. Qi, H. J. W. van Triest, Y. Yue, M. Xu, and Y. Kang, “Automatic pulmonary fissure detection and lobe segmentation in CT chest images,” Biomedical Engineering Online, vol. 13, no. 1, p. 59, 2014.View at: Publisher Site | Google Scholar
S. Ukil and J. M. Reinhardt, “Anatomy-guided lung lobe segmentation in X-ray CT images,” IEEE Transactions on Medical Imaging, vol. 28, no. 2, pp. 202–214, 2009.View at: Publisher Site | Google Scholar
T. Doel, T. N. Matin, F. V. Gleeson, D. J. Gavaghan, and V. Grau, “Pulmonary Lobe Segmentation from CT Images Using Fissureness, Airways, Vessels and Multilevel B-Splines,” 2012.View at: Google Scholar
B. Lassen, E. M. van Rikxoort, M. Schmidt, S. Kerkstra, B. van Ginneken, and J.-M. Kuhnigk, “Automatic segmentation of the pulmonary lobes from chest CT scans based on fissures, vessels, and bronchi,” IEEE Transactions on Medical Imaging, vol. 32, no. 2, pp. 210–222, 2013.View at: Publisher Site | Google Scholar
F. J. S. Bragman, J. R. McClelland, J. Jacob, J. R. Hurst, and D. J. Hawkes, “Pulmonary lobe segmentation with probabilistic segmentation of the fissures and a groupwise fissure prior,” IEEE Transactions on Medical Imaging, vol. 36, no. 8, pp. 1650–1663, 2017.View at: Publisher Site | Google Scholar
C. Bauer, M. Eberlein, and R. R. Beichel, “Pulmonary lobe separation in expiration chest CT scans based on subject-specific priors derived from inspiration scans,” Journal of Medical Imaging, vol. 5, no. 1, 2018.View at: Publisher Site | Google Scholar
N. Giuliani, C. Payer, M. Pienn, H. Olschewski, and M. Urschler, “Pulmonary lobe segmentation in CT images using alpha-expansion,” Processing of VISIGRAPP, vol. 23, pp. 387–394, 2018.View at: Google Scholar
L. Zhang, E. A. Hoffman, and J. M. Reinhardt, “Atlas-driven lung lobe segmentation in volumetric X-ray CT images,” IEEE Transactions on Medical Imaging, vol. 25, no. 1, pp. 1–16, 2006.View at: Google Scholar
E. M. van Rikxoort, M. Prokop, B. de Hoop, M. A. Viergever, J. Pluim, and B. van Ginneken, “Automatic segmentation of pulmonary lobes robust against incomplete fissures,” IEEE Transactions on Medical Imaging, vol. 29, no. 6, pp. 1286–1296, 2010.View at: Publisher Site | Google Scholar
A. Schmidt-Richberg, J. Ehrhardt, M. Wilms, R. Werner, and H. Handels, “Pulmonary Lobe Segmentation with Level Sets,” International Society for Optics and Photonics, vol. 20, 2012.View at: Google Scholar
X. Zhou, T. Hayashi, T. Hara et al., “Automatic segmentation and recognition of anatomical lung structures from high-resolution chest CT images,” Computerized Medical Imaging and Graphics, vol. 30, no. 5, pp. 299–313, 2006.View at: Publisher Site | Google Scholar
Y. Iwao, T. Gotoh, S. Kagei, T. Iwasawa, and M. D. S. G. Tsuzuki, “Integrated lung field segmentation of injured region with anatomical structure analysis by failure-recovery algorithm from chest CT images,” Biomedical Signal Processing and Control, vol. 12, pp. 28–38, 2014.View at: Publisher Site | Google Scholar
C. Xiao, M. Staring, J. Wang, D. P. Shamonin, and B. C. Stoel, “A Derivative of Stick Filter for Pulmonary Fissure Detection in CT Images,” 2013.View at: Google Scholar
C. Xiao, B. C. Stoel, M. E. Bakker, Y. Peng, J. Stolk, and M. Staring, “Pulmonary fissure detection in CT images using a derivative of stick filter,” IEEE Transactions on Medical Imaging, vol. 35, no. 6, pp. 1488–1500, 2016.View at: Publisher Site | Google Scholar
Y. Peng and C. Xiao, “An oriented derivative of stick filter and post-processing segmentation algorithms for pulmonary fissure detection in CT images,” Biomedical Signal Processing and Control, vol. 43, no. 5, pp. 278–288, 2018.View at: Publisher Site | Google Scholar
E. M. van Rikxoort and B. van Ginneken, “Automatic Segmentation of the Lungs and Lobes from Thoracic CT Scans,” 2011.View at: Google Scholar
J. Pu, J. K. Leader, B. Zheng et al., “A computational geometry approach to automated pulmonary fissure segmentation in CT examinations,” IEEE Transactions on Medical Imaging, vol. 28, no. 5, pp. 710–719, 2009.View at: Publisher Site | Google Scholar
E. Smistad, A. C. Elster, and F. Lindseth, “GPU accelerated segmentation and centerline extraction of tubular structures from medical images,” International Journal of Computer Assisted Radiology and Surgery, vol. 9, no. 4, pp. 561–575, 2014.View at: Publisher Site | Google Scholar
C. Bauer, T. Pock, E. Sorantin, H. Bischof, and R. Beichel, “Segmentation of interwoven 3d tubular tree structures utilizing shape priors and graph cuts,” Medical Image Analysis, vol. 14, no. 2, pp. 172–184, 2010.View at: Publisher Site | Google Scholar
H. Edelsbrunner and E. P. Mücke, “Three-dimensional alpha shapes,” ACM Transactions on Graphics, vol. 13, no. 1, pp. 43–72, 1994.View at: Publisher Site | Google Scholar
Y. Ge, Z. Lin, H. Tang, P. Zhong, and B. Cao, “Measurement of particle size of loose accumulation based on alpha shapes (AS) and hill climbing-region growing (HC-RG) algorithms,” Sensors, vol. 20, no. 3, p. 883, 2020.View at: Publisher Site | Google Scholar
C. Xu and J. L. Prince, “Snakes, shapes, and gradient vector flow,” IEEE Transactions on Image Processing:A Publication of the IEEE Signal Processing Society, vol. 7, no. 3, pp. 359–369, 1998.View at: Publisher Site | Google Scholar
T. Pock, R. Beichel, and H. Bischof, A Novel Robust Tube Detection Filter for 3D Centerline Extraction, Springer, Berlin, Germany, 2005.
K. Krissian, G. Malandain, N. Ayache, R. Vaillant, and Y. Trousset, “Model-based detection of tubular structures in 3D images,” Computer Vision and Image Understanding, vol. 80, no. 2, pp. 130–171, 2000.View at: Publisher Site | Google Scholar
J. Tschirren, G. McLennan, K. Palagyi, E. A. Hoffman, and M. Sonka, “Matching and anatomical labeling of human airway tree,” IEEE Transactions on Medical Imaging, vol. 24, no. 12, pp. 1540–1547, 2005.View at: Publisher Site | Google Scholar
T. C. Lee, R. L. Kashyap, and C. N. Chu, “Building skeleton models via 3-D medial surface Axis thinning algorithms,” CVGIP: Graphical Models and Image Processing, vol. 56, no. 6, pp. 462–478, 1994.View at: Publisher Site | Google Scholar
M. Kerschnitzki, P. Kollmannsberger, M. Burghammer et al., “Architecture of the osteocyte network correlates with bone material quality,” Journal of Bone and Mineral Research, vol. 28, no. 8, pp. 1837–1845, 2013.View at: Publisher Site | Google Scholar
C. Zhang, M. Sun, Y. Wei et al., “Automatic segmentation of arterial tree from 3D computed tomographic pulmonary angiography (CTPA) scans,” Computer Assisted Surgery, vol. 24, no. 2, pp. 79–86, 2019.View at: Publisher Site | Google Scholar
P. Nardelli, D. Jimenez-Carretero, D. Bermejo-Pelaez et al., “Pulmonary artery-vein classification in CT images using deep learning,” IEEE Transactions on Medical Imaging, vol. 37, no. 11, pp. 2428–2440, 2018.View at: Publisher Site | Google Scholar
D. Jimenez-Carretero, D. Bermejo-Peláez, P. Nardelli et al., “A graph-cut approach for pulmonary artery-vein segmentation in noncontrast CT images,” Medical Image Analysis, vol. 52, 2019.View at: Publisher Site | Google Scholar
F. Ritter, T. Boskamp, A. Homeyer et al., “Medical image analysis,” IEEE Pulse, vol. 2, no. 6, pp. 60–70, 2011.View at: Publisher Site | Google Scholar
S. Kazuhiko, M. Yoshitake, M. Takahiro, and S. Satoru, “Simulation of the transient current of radiation detector materials using the constrained profile interpolation method,” Nuclear Instruments and Methods in Physics Research Section A, vol. 971, p. 164128, 2020.View at: Google Scholar
E. H. Maria, C. Michael, B. Adam et al., “Review of chest CT manifestations of COVID-19 infection,” European Journal of Radiology Open, vol. 7, p. 100239, 2020.View at: Google Scholar
S. Anitha and T. G. Babu, “An efficient method for the detection of oblique fissures from computed tomography images of lungs,” Journal of Medical Systems, vol. 43, no. 8, p. 252, 2019.View at: Publisher Site | Google Scholar
L. Zhang, T. Gu, J. Zhao et al., “An improved moving least squares method for curve and surface fitting,” Mathematical Problems in Engineering, vol. 23, 2013.View at: Google Scholar