Graph-Based Airway Tree Reconstruction From Chest CT Scans: Evaluation of Different Features on Five Cohorts

Accurate Airway Tree Segmentation in CT Scans via Anatomy-Aware Multi-Class Segmentation and Topology-Guided Iterative Learning

Intrathoracic airway segmentation in computed tomography is a prerequisite for various respiratory disease analyses such as chronic obstructive pulmonary disease, asthma and lung cancer. Due to the low imaging contrast and noises execrated at peripheral branches, the topological-complexity and the intra-class imbalance of airway tree, it remains challenging for deep learning-based methods to segment the complete airway tree (on extracting deeper branches). Unlike other organs with simpler shapes or topology, the airway's complex tree structure imposes an unbearable burden to generate the "ground truth" label (up to 7 or 3 hours of manual or semi-automatic annotation per case). Most of the existing airway datasets are incompletely labeled/annotated, thus limiting the completeness of computer-segmented airway. In this paper, we propose a new anatomy-aware multi-class airway segmentation method enhanced by topology-guided iterative self-learning. Based on the natural airway anatomy, we formulate a simple yet highly effective anatomy-aware multi-class segmentation task to intuitively handle the severe intra-class imbalance of the airway. To solve the incomplete labeling issue, we propose a tailored iterative self-learning scheme to segment toward the complete airway tree. For generating pseudo-labels to achieve higher sensitivity (while retaining similar specificity), we introduce a novel breakage attention map and design a topology-guided pseudo-label refinement method by iteratively connecting breaking branches commonly existed from initial pseudo-labels. Extensive experiments have been conducted on four datasets including two public challenges. The proposed method achieves the top performance in both EXACT'09 challenge using average score and ATM'22 challenge on weighted average score. In a public BAS dataset and a private lung cancer dataset, our method significantly improves previous leading approaches by extracting at least (absolute) 6.1% more detected tree length and 5.2% more tree branches, while maintaining comparable precision.

A Precise Pulmonary Airway Tree Segmentation Method Using Quasi-Spherical Region Constraint and Tracheal Wall Gap Sealing

Accurate segmentation of the pulmonary airway tree is crucial for diagnosing lung diseases. To tackle the issues of low segmentation accuracy and frequent leaks in existing methods, this paper proposes a precise segmentation method using quasi-spherical region-constrained wavefront propagation with tracheal wall gap sealing. Based on the characteristic that the surface formed by seed points approximates the airway cross-section, the width of the unsegmented airway is calculated, determining the initial quasi-spherical constraint region. Using the wavefront propagation method, seed points are continuously propagated and segmented along the tracheal wall within the quasi-spherical constraint region, thus overcoming the need to determine complex segmentation directions. To seal tracheal wall gaps, a morphological closing operation is utilized to extract the characteristics of small holes and locate low-brightness tracheal wall gaps. By filling the CT values at these gaps, the method seals the tracheal wall gaps. Extensive experiments on the EXACT09 dataset demonstrate that our algorithm ranks third in segmentation completeness. Moreover, its performance in preventing airway leaks is significantly better than the top-two algorithms, effectively preventing large-scale leak-induced spread.

Tree Recovery by Dynamic Programming

Tree-like structures are common, naturally occurring objects that are of interest to many fields of study, such as plant science and biomedicine. Analysis of these structures is typically based on skeletons extracted from captured data, which often contain spurious cycles that need to be removed. We propose a dynamic programming algorithm for solving the NP-hard tree recovery problem formulated by (Estrada et al. 2015), which seeks a least-cost partitioning of the graph nodes that yields a directed tree. Our algorithm finds the optimal solution by iteratively contracting the graph via node-merging until the problem can be trivially solved. By carefully designing the merging sequence, our algorithm can efficiently recover optimal trees for many real-world data where (Estrada et al. 2015) only produces sub-optimal solutions. We also propose an approximate variant of dynamic programming using beam search, which can process graphs containing thousands of cycles with significantly improved optimality and efficiency compared with (Estrada et al. 2015).

Hybrid Airway Segmentation Using Multi-Scale Tubular Structure Filters and Texture Analysis on 3D Chest CT Scans

Airway diseases are frequently related to morphological changes that may influence lung physiology. Accurate airway region segmentation may be useful for quantitative evaluation of disease prognosis and therapy efficacy. The information can also be applied to understand the fundamental mechanisms of various lung diseases. We present a hybrid method to automatically segment the airway regions on 3D volume chest computed tomography (CT) scans. This method uses multi-scale filtering and support vector machine (SVM) classification. The proposed scheme is comprised of two hybrid steps. First, a tubular structure-based multi-scale filter is applied to find the initial candidate airway regions. Second, for identifying candidate airway regions using the fuzzy connectedness technique, the small and disconnected branches of airway regions are detected using SVM classification trained to differentiate between airway and non-airway regions through texture analysis of user-defined landmark points. For development and evaluation of the method, two datasets were incorporated: (1) 55 lung-CT volumes from the Korean Obstructive Lung Disease (KOLD) Cohort Study and (2) 20 cases from the publicly open database (EXACT'09). The average tree-length detection rates of EXACT'09 and KOLD were 56.9 ± 11.0 and 70.5 ± 8.98, respectively. Comparison of the results for the EXACT'09 data set between the presented method and other methods revealed that our approach was a high performer. The method limitations were higher false-positive rates than those of the other methods and risk of leakage. In future studies, application of a convolutional neural network will help overcome these shortcomings.

Segmentation of distal airways using structural analysis

Segmentation of airways in Computed Tomography (CT) scans is a must for accurate support of diagnosis and intervention of many pulmonary disorders. In particular, lung cancer diagnosis would benefit from segmentations reaching most distal airways. We present a method that combines descriptors of bronchi local appearance and graph global structural analysis to fine-tune thresholds on the descriptors adapted for each bronchial level. We have compared our method to the top performers of the EXACT09 challenge and to a commercial software for biopsy planning evaluated in an own-collected data-base of high resolution CT scans acquired under different breathing conditions. Results on EXACT09 data show that our method provides a high leakage reduction with minimum loss in airway detection. Results on our data-base show the reliability across varying breathing conditions and a competitive performance for biopsy planning compared to a commercial solution.

Chest wall strapping increases expiratory airflow and detectable airway segments in computer tomographic scans of normal and obstructed lungs

Chest wall strapping (CWS) induces breathing at low lung volumes but also increases parenchymal elastic recoil. In this study, we tested the hypothesis that CWS dilates airways via airway-parenchymal interdependence. In 11 subjects (6 healthy and 5 with mild to moderate COPD), pulmonary function tests and lung volumes were obtained in control (baseline) and the CWS state. Control and CWS-CT scans were obtained at 50% of control (baseline) total lung-capacity (TLC). CT lung volumes were analyzed by CT volumetry. If control and CWS-CT volumetry did not differ by more than 25%, airway dimensions were analyzed via automated airway segmentation. CWS-TLC was reduced on average to 71% of control-TLC in normal subjects and 79% of control-TLC in subjects with COPD. CWS increased expiratory airflow at 50% of control-TLC by 41% (3.50 ± 1.6 vs. 4.93 ± 1.9 l/s, P = 0.04) in normals and 316% in COPD(0.25 ± 0.05 vs 0.79 ± 0.39 l/s, P = 0.04). In 10 subjects (5 normals and 5 COPD), control and CWS-CT scans at 50% control-TLC did not differ more than 25% on CT volumetry and were included in the airway structure analysis. CWS increased the mean number of detectable airways with a diameter of ≤2 mm by 32.5% (65 ± 10 vs. 86 ± 124, P = 0.01) in normal subjects and by 79% (59 ± 19 vs. 104 ± 16, P = 0.01) in subjects with COPD. There was no difference in the number of detectable airways with diameters 2-4 mm and >4 mm in normal or in COPD subjects. In conclusion, CWS enhances the detection of small airways via automated CT airway segmentation and increases expiratory airflow in normal subjects as well as in subjects with mild to moderate COPD. NEW & NOTEWORTHY In normal and COPD subjects, chest wall strapping(CWS) increased the number of detectable small airways using automated CT airway segmentation. The concept of dysanapsis expresses the physiological variation in the geometry of the tracheobronchial tree and lung parenchyma based on development. We propose a dynamic concept to dysanapsis in which CWS leads to breathing at lower lung volumes with a corresponding increase in the size of small airways, a potentially novel, nonpharmacological treatment for COPD.

Higher BMI is associated with higher expiratory airflow normalised for lung volume (FEF25-75/FVC) in COPD

INTRODUCTION

The obesity paradox in chronic obstructive pulmonary disease (COPD), whereby patients with higher body mass index (BMI) fare better, is poorly understood. Higher BMIs are associated with lower lung volumes and greater lung elastic recoil, a key determinant of expiratory airflow. The forced expiratory flow (25-75) (FEF25-75)/forced vital capacity (FVC) ratio reflects effort-independent expiratory airflow in the context of lung volume and could be modulated by BMI.

METHODS

We analysed data from the COPDGene study, an observational study of 10 192 subjects, with at least a 10 pack-year smoking history. Data were limited to subjects with BMI 20-40 kg/m2 (n=9222). Subjects were stratified according to forced expiratory volume in 1 s (FEV1) (%predicted)-quintiles. In regression analyses and Cox proportional hazard models, we analysed the association between BMI, the FEF25-75/FVC ratio, the imaging phenotype, COPD exacerbations, hospitalisations and death.

RESULTS

There was no correlation between BMI and FEV1(%predicted). However, a higher BMI is correlated with a higher FEF25-75/FVC ratio. In CT scans, a higher BMI was associated with less emphysema and less air trapping. In risk-adjusted models, the quintile with the highest FEF25-75/FVC ratio was associated with a 46% lower risk of COPD exacerbations (OR 0.54, p<0.001) and a 40% lower risk of death (HR 0.60, p=0.02), compared with the lowest quintile. BMI was not independently associated with these outcomes.

CONCLUSIONS

A higher BMI is associated with lower lung volumes and higher expiratory airflows when normalised for lung volume, as quantified by the FEF25-75/FVC ratio. A higher FEF25-75/FVC ratio is associated with a lower risk of COPD exacerbations and death and might quantify functional aspects of the paradoxical effect of higher BMIs on COPD.

Airway tree reconstruction in expiration chest CT scans facilitated by information transfer from corresponding inspiration scans

PURPOSE

Analysis and comparison of airways imaged in pairs of inspiration and expiration lung CT scans provides important information for quantitative assessment of lung diseases like chronic obstructive pulmonary disease. However, airway tree reconstruction in expiration CT scans is a challenging problem. Typically, only a low number of airway branches are found in expiration scans, compared to inspiration scans. To detect more airways in expiration CT scans, the authors introduce a novel airway reconstruction approach and assess its performance.

METHODS

The method requires a pair of inspiration and expiration CT scans and utilizes information from the inspiration scan to facilitate reconstructing the airway tree in the expiration lung CT scan. First, an initial airway tree (high confidence) and airway candidates (limited confidence) are reconstructed in the expiration scan by utilizing a 3D graph-based reconstruction method. Then, the 3D airway tree is reconstructed in the inspiration scan. Second, correspondences between expiration and inspiration tree structures are identified by utilizing a novel hierarchical tree matching algorithm that utilizes a local CT image-based similarity criterion. Third, the tree information from the inspiration airway tree is used to select expiration candidates, resulting in the final expiration tree structure. The approach was evaluated on a diverse set of 40 scan pairs and compared to the baseline method, which utilizes only the expiration CT scan.

RESULTS

The proposed method produced a significant (p < 0.05) increase in airway tree length by 13.35 cm, on average, which represents an 11.21% increase relative to the baseline result using only the expiration CT scan. A detailed analysis of all additionally identified airways resulted in a true and false positive rate of 94.8% and 5.2%, respectively. The true positive rate was found to be significantly higher than the false positive rate (p < 0.05).

CONCLUSIONS

The proposed method allowed increasing the number of found airways in expiration scans significantly. In addition, the algorithm establishes correspondence between inspiration and expiration airway trees, which can facilitate label transfer between airway trees and quantitative assessment of change in airways. The approach can be adapted to facilitate airway reconstruction in several longitudinal lung CT scans by means of mutual information transfer.

An approach for reducing the error rate in automated lung segmentation

Robust lung segmentation is challenging, especially when tens of thousands of lung CT scans need to be processed, as required by large multi-center studies. The goal of this work was to develop and assess a method for the fusion of segmentation results from two different methods to generate lung segmentations that have a lower failure rate than individual input segmentations. As basis for the fusion approach, lung segmentations generated with a region growing and model-based approach were utilized. The fusion result was generated by comparing input segmentations and selectively combining them using a trained classification system. The method was evaluated on a diverse set of 204 CT scans of normal and diseased lungs. The fusion approach resulted in a Dice coefficient of 0.9855±0.0106 and showed a statistically significant improvement compared to both input segmentation methods. In addition, the failure rate at different segmentation accuracy levels was assessed. For example, when requiring that lung segmentations must have a Dice coefficient of better than 0.97, the fusion approach had a failure rate of 6.13%. In contrast, the failure rate for region growing and model-based methods was 18.14% and 15.69%, respectively. Therefore, the proposed method improves the quality of the lung segmentations, which is important for subsequent quantitative analysis of lungs. Also, to enable a comparison with other methods, results on the LOLA11 challenge test set are reported.

Lung Segmentation in 4D CT Volumes Based on Robust Active Shape Model Matching

Dynamic and longitudinal lung CT imaging produce 4D lung image data sets, enabling applications like radiation treatment planning or assessment of response to treatment of lung diseases. In this paper, we present a 4D lung segmentation method that mutually utilizes all individual CT volumes to derive segmentations for each CT data set. Our approach is based on a 3D robust active shape model and extends it to fully utilize 4D lung image data sets. This yields an initial segmentation for the 4D volume, which is then refined by using a 4D optimal surface finding algorithm. The approach was evaluated on a diverse set of 152 CT scans of normal and diseased lungs, consisting of total lung capacity and functional residual capacity scan pairs. In addition, a comparison to a 3D segmentation method and a registration based 4D lung segmentation approach was performed. The proposed 4D method obtained an average Dice coefficient of 0.9773 ± 0.0254, which was statistically significantly better (p value ≪0.001) than the 3D method (0.9659 ± 0.0517). Compared to the registration based 4D method, our method obtained better or similar performance, but was 58.6% faster. Also, the method can be easily expanded to process 4D CT data sets consisting of several volumes.