A CT-Based Automated Algorithm for Airway Segmentation Using Freeze-and-Grow Propagation and Deep Learning

Automated CT-based measurements of radial and longitudinal expansion of airways due to breathing-related lung volume change

BACKGROUND

Respiratory function is impaired in chronic obstructive pulmonary disease (COPD). Automation of multi-volume CT-based measurements of different components of breathing-related airway deformations will help understand multi-pathway impairments in respiratory mechanics in COPD.

PURPOSE

To develop and evaluate multi-volume chest CT-based automated measurements of breathing-related radial and longitudinal expansion of individual airways between inspiratory and expiratory lung volumes.

METHODS

We developed a method to compute breathing-related airway deformation metrics and applied it to total lung capacity (TLC) and functional residual capacity (FRC) chest CT scans. The computational pipeline involves: (1) segmentation of airways; (2) skeletonization of airways; (3) labeling of anatomical airway segments at TLC and FRC; and (4) computation of radial and longitudinal expansion metrics of individual airways across lung volumes. Radial expansion (∆CSA) of an airway is computed as the percent change of its cross-sectional area (CSA) between two lung volumes. Longitudinal expansion (∆L) of an airway is computed as the percent change in its airway path-length from the carina between lung volumes. These measures are summarized at different airway anatomic generations. Agreement of automated measures with their manually derived values was examined in terms of concordance correlation coefficient (CCC) of automated measures with those derived using manual outlining. Intra-class correlation coefficient (ICC) of automated measures from repeat CT scans (n = 37) was computed to assess repeatability. The method was also applied to a set of participants from the Genetic Epidemiology of COPD (COPDGene) Iowa cohort, distributed across COPD severity groups (n = 4 × 60).

RESULTS

The CCC values for the automated ∆CSA measure with manually derived values were 0.930 at the trachea, 0.898 at primary bronchi, and greater than 0.95 at pre-segmental and segmental airways; these CCC values were consistently greater than 0.95 for ∆L at all airway generations. ICC values for repeatability of ∆CSA were 0.974, 0.950, 0.943, and 0.901 at trachea, primary bronchi, pre-segmental, and segmental airways, respectively; these ICC values for ∆L were 0.973, 0.954, and 0.952 at primary bronchi, pre-segmental, and segmental airways, respectively. ∆CSA values were significantly reduced (p < 0.001) with increasing COPD severity at each of primary bronchi, pre-segmental, and segmental airways. Significantly lower ∆L values were observed for moderate (p = 0.042 at pre-segmental and p = 0.037 at segmental) and severe (p = 0.019 at pre-segmental and p < 0.001 at segmental) COPD groups as compared to the preserved lung function group. Body mass index (BMI) and smoking status were found to significantly associate with ∆CSA at segmental airways (r = 0.17 and -0.19, respectively; significance threshold = 0.13), while age and sex were significantly associated with ∆L (r = -0.21 and -0.17, respectively); COPD severity was significantly associated with both ∆CSA and ∆L (r = -0.35 and -0.22, respectively).

CONCLUSION

Our CT-based automated measures of breathing-related radial and longitudinal expansion of airways are repeatable and in agreement with manually derived values. Automation of different airway mechanical biomarkers and their observed significant associations with age, sex, BMI, smoking, and COPD severity establish an effective tool to investigate multi-pathway impairments of respiratory mechanics in COPD and other lung diseases.

Automated CT-based decoupling of the effects of airway narrowing and wall thinning on airway counts in chronic obstructive pulmonary disease

OBJECTIVE

We examine pathways of airway alteration due to wall thinning, narrowing, and obliteration in chronic obstructive pulmonary disease (COPD) using CT-derived airway metrics.

METHODS

Ex-smokers (N = 649; age mean ± std: 69 ± 6 years; 52% male) from the COPDGene Iowa cohort (September 2013-July 2017) were studied. Total airway count (TAC), peripheral TAC beyond 7th generation (TACp), and airway wall thickness (WT) were computed from chest CT scans using previously validated automated methods. Causal relationships among demographic, smoking, spirometry, COPD severity, airway counts, WT, and scanner variables were analysed using causal inference techniques including direct acyclic graphs to assess multi-pathway alterations of airways in COPD.

RESULTS

TAC, TACp, and WT were significantly lower (P < .0001) in mild, moderate, and severe COPD compared to the preserved lung function group. TAC (TACp) losses attributed to narrowing and obliteration of small airways were 4.59%, 13.29%, and 32.58% (4.64%, 17.82%, and 45.51%) in mild, moderate, and severe COPD, while the losses attributed to wall thinning were 8.24%, 17.01%, and 22.95% (12.79%, 25.66%, and 33.95%) in respective groups.

CONCLUSIONS

Different pathways of airway alteration in COPD are observed using CT-derived automated airway metrics. Wall thinning is a dominant contributor to both TAC and TACp loss in mild and moderate COPD while narrowing and obliteration of small airways is dominant in severe COPD.

ADVANCES IN KNOWLEDGE

This automated CT-based study shows that wall thinning dominates airway alteration in mild and moderate COPD while narrowing and obliteration of small airways leads the alteration process in severe COPD.

BCNet: Bronchus Classification via Structure Guided Representation Learning

CT-based bronchial tree analysis is a key step for the diagnosis of lung and airway diseases. However, the topology of bronchial trees varies across individuals, which presents a challenge to the automatic bronchus classification. To solve this issue, we propose the Bronchus Classification Network (BCNet), a structure-guided framework that exploits the segment-level topological information using point clouds to learn the voxel-level features. BCNet has two branches, a Point-Voxel Graph Neural Network (PV-GNN) for segment classification, and a Convolutional Neural Network (CNN) for voxel labeling. The two branches are simultaneously trained to learn topology-aware features for their shared backbone while it is feasible to run only the CNN branch for the inference. Therefore, BCNet maintains the same inference efficiency as its CNN baseline. Experimental results show that BCNet significantly exceeds the state-of-the-art methods by over 8.0% both on F1-score for classifying bronchus. Furthermore, we contribute BronAtlas: an open-access benchmark of bronchus imaging analysis with high-quality voxel-wise annotations of both anatomical and abnormal bronchial segments. The benchmark is available at https://osf.io/pskr9/?viewonly=94fa3d87274b4095ac9a4b88cc9a1341.

Detail-sensitive 3D-UNet for pulmonary airway segmentation from CT images

The segmentation of airway from computed tomography (CT) images plays a vital role in pulmonary disease diagnosis, evaluation, surgical planning, and treatment. Nevertheless, it is still challenging for current methods to handle distal thin and low-contrast airways, leading to mis-segmentation issues. This paper proposes a detail-sensitive 3D-UNet (DS-3D-UNet) that incorporates two new modules into 3D-UNet to segment airways accurately from CT images. The feature recalibration module is designed to give more attention to the foreground airway features through a new attention mechanism. The detail extractor module aims to restore multi-scale detailed features by fusion of features at different levels. Extensive experiments were conducted on the ATM'22 challenge dataset composed of 300 CT scans with airway annotations to evaluate its performance. Quantitative comparisons prove that the proposed model achieves the best performance in terms of Dice similarity coefficient (92.6%) and Intersection over Union (86.3%), outperforming other state-of-the-art methods. Qualitative comparisons further exhibit the superior performance of our method in segmenting thin and confused distal bronchi. The proposed model could provide important references for the diagnosis and treatment of pulmonary diseases, holding promising prospects in the field of digital medicine. Codes are available at https://github.com/nighlevil/DS-3D-UNet/tree/master .

Ultra-low dose hip CT-based automated measurement of volumetric bone mineral density at proximal femoral subregions

BACKGROUND

Forty to fifty percent of women and 13%-22% of men experience an osteoporosis-related fragility fracture in their lifetimes. After the age of 50 years, the risk of hip fracture doubles in every 10 years. x-Ray based DXA is currently clinically used to diagnose osteoporosis and predict fracture risk. However, it provides only 2-D representation of bone and is associated with other technical limitations. Thus, alternative methods are needed.

PURPOSE

To develop and evaluate an ultra-low dose (ULD) hip CT-based automated method for assessment of volumetric bone mineral density (vBMD) at proximal femoral subregions.

METHODS

An automated method was developed to segment the proximal femur in ULD hip CT images and delineate femoral subregions. The computational pipeline consists of deep learning (DL)-based computation of femur likelihood map followed by shape model-based femur segmentation and finite element analysis-based warping of a reference subregion labeling onto individual femur shapes. Finally, vBMD is computed over each subregion in the target image using a calibration phantom scan. A total of 100 participants (50 females) were recruited from the Genetic Epidemiology of COPD (COPDGene) study, and ULD hip CT imaging, equivalent to 18 days of background radiation received by U.S. residents, was performed on each participant. Additional hip CT imaging using a clinical protocol was performed on 12 participants and repeat ULD hip CT was acquired on another five participants. ULD CT images from 80 participants were used to train the DL network; ULD CT images of the remaining 20 participants as well as clinical and repeat ULD CT images were used to evaluate the accuracy, generalizability, and reproducibility of segmentation of femoral subregions. Finally, clinical CT and repeat ULD CT images were used to evaluate accuracy and reproducibility of ULD CT-based automated measurements of femoral vBMD.

RESULTS

Dice scores of accuracy (n = 20), reproducibility (n = 5), and generalizability (n = 12) of ULD CT-based automated subregion segmentation were 0.990, 0.982, and 0.977, respectively, for the femoral head and 0.941, 0.970, and 0.960, respectively, for the femoral neck. ULD CT-based regional vBMD showed Pearson and concordance correlation coefficients of 0.994 and 0.977, respectively, and a root-mean-square coefficient of variation (RMSCV) (%) of 1.39% with the clinical CT-derived reference measure. After 3-digit approximation, each of Pearson and concordance correlation coefficients as well as intraclass correlation coefficient (ICC) between baseline and repeat scans were 0.996 with RMSCV of 0.72%. Results of ULD CT-based bone analysis on 100 participants (age (mean ± SD) 73.6 ± 6.6 years) show that males have significantly greater (p < 0.01) vBMD at the femoral head and trochanteric regions than females, while females have moderately greater vBMD (p = 0.05) at the medial half of the femoral neck than males.

CONCLUSION

Deep learning, combined with shape model and finite element analysis, offers an accurate, reproducible, and generalizable algorithm for automated segmentation of the proximal femur and anatomic femoral subregions using ULD hip CT images. ULD CT-based regional measures of femoral vBMD are accurate and reproducible and demonstrate regional differences between males and females.

AeroPath: An airway segmentation benchmark dataset with challenging pathology and baseline method

To improve the prognosis of patients suffering from pulmonary diseases, such as lung cancer, early diagnosis and treatment are crucial. The analysis of CT images is invaluable for diagnosis, whereas high quality segmentation of the airway tree are required for intervention planning and live guidance during bronchoscopy. Recently, the Multi-domain Airway Tree Modeling (ATM'22) challenge released a large dataset, both enabling training of deep-learning based models and bringing substantial improvement of the state-of-the-art for the airway segmentation task. The ATM'22 dataset includes a large group of COVID'19 patients and a range of other lung diseases, however, relatively few patients with severe pathologies affecting the airway tree anatomy was found. In this study, we introduce a new public benchmark dataset (AeroPath), consisting of 27 CT images from patients with pathologies ranging from emphysema to large tumors, with corresponding trachea and bronchi annotations. Second, we present a multiscale fusion design for automatic airway segmentation. Models were trained on the ATM'22 dataset, tested on the AeroPath dataset, and further evaluated against competitive open-source methods. The same performance metrics as used in the ATM'22 challenge were used to benchmark the different considered approaches. Lastly, an open web application is developed, to easily test the proposed model on new data. The results demonstrated that our proposed architecture predicted topologically correct segmentations for all the patients included in the AeroPath dataset. The proposed method is robust and able to handle various anomalies, down to at least the fifth airway generation. In addition, the AeroPath dataset, featuring patients with challenging pathologies, will contribute to development of new state-of-the-art methods. The AeroPath dataset and the web application are made openly available.

Hunting imaging biomarkers in pulmonary fibrosis: Benchmarks of the AIIB23 challenge

Airway-related quantitative imaging biomarkers are crucial for examination, diagnosis, and prognosis in pulmonary diseases. However, the manual delineation of airway structures remains prohibitively time-consuming. While significant efforts have been made towards enhancing automatic airway modelling, current public-available datasets predominantly concentrate on lung diseases with moderate morphological variations. The intricate honeycombing patterns present in the lung tissues of fibrotic lung disease patients exacerbate the challenges, often leading to various prediction errors. To address this issue, the 'Airway-Informed Quantitative CT Imaging Biomarker for Fibrotic Lung Disease 2023' (AIIB23) competition was organized in conjunction with the official 2023 International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI). The airway structures were meticulously annotated by three experienced radiologists. Competitors were encouraged to develop automatic airway segmentation models with high robustness and generalization abilities, followed by exploring the most correlated QIB of mortality prediction. A training set of 120 high-resolution computerised tomography (HRCT) scans were publicly released with expert annotations and mortality status. The online validation set incorporated 52 HRCT scans from patients with fibrotic lung disease and the offline test set included 140 cases from fibrosis and COVID-19 patients. The results have shown that the capacity of extracting airway trees from patients with fibrotic lung disease could be enhanced by introducing voxel-wise weighted general union loss and continuity loss. In addition to the competitive image biomarkers for mortality prediction, a strong airway-derived biomarker (Hazard ratio>1.5, p < 0.0001) was revealed for survival prognostication compared with existing clinical measurements, clinician assessment and AI-based biomarkers.

Modeling Realistic Geometries in Human Intrathoracic Airways

Geometrical models of the airways offer a comprehensive perspective on the complex interplay between lung structure and function. Originating from mathematical frameworks, these models have evolved to include detailed lung imagery, a crucial enhancement that aids in the early detection of morphological changes in the airways, which are often the first indicators of diseases. The accurate representation of airway geometry is crucial in research areas such as biomechanical modeling, acoustics, and particle deposition prediction. This review chronicles the evolution of these models, from their inception in the 1960s based on ideal mathematical constructs, to the introduction of advanced imaging techniques like computerized tomography (CT) and, to a lesser degree, magnetic resonance imaging (MRI). The advent of these techniques, coupled with the surge in data processing capabilities, has revolutionized the anatomical modeling of the bronchial tree. The limitations and challenges in both mathematical and image-based modeling are discussed, along with their applications. The foundation of image-based modeling is discussed, and recent segmentation strategies from CT and MRI scans and their clinical implications are also examined. By providing a chronological review of these models, this work offers insights into the evolution and potential future of airway geometry modeling, setting the stage for advancements in diagnosing and treating lung diseases. This review offers a novel perspective by highlighting how advancements in imaging techniques and data processing capabilities have significantly enhanced the accuracy and applicability of airway geometry models in both clinical and research settings. These advancements provide unique opportunities for developing patient-specific models.

Artificial intelligence in COPD CT images: identification, staging, and quantitation

Chronic obstructive pulmonary disease (COPD) stands as a significant global health challenge, with its intricate pathophysiological manifestations often demanding advanced diagnostic strategies. The recent applications of artificial intelligence (AI) within the realm of medical imaging, especially in computed tomography, present a promising avenue for transformative changes in COPD diagnosis and management. This review delves deep into the capabilities and advancements of AI, particularly focusing on machine learning and deep learning, and their applications in COPD identification, staging, and imaging phenotypes. Emphasis is laid on the AI-powered insights into emphysema, airway dynamics, and vascular structures. The challenges linked with data intricacies and the integration of AI in the clinical landscape are discussed. Lastly, the review casts a forward-looking perspective, highlighting emerging innovations in AI for COPD imaging and the potential of interdisciplinary collaborations, hinting at a future where AI doesn't just support but pioneers breakthroughs in COPD care. Through this review, we aim to provide a comprehensive understanding of the current state and future potential of AI in shaping the landscape of COPD diagnosis and management.

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.

Fuzzy Attention Neural Network to Tackle Discontinuity in Airway Segmentation

Airway segmentation is crucial for the examination, diagnosis, and prognosis of lung diseases, while its manual delineation is unduly burdensome. To alleviate this time-consuming and potentially subjective manual procedure, researchers have proposed methods to automatically segment airways from computerized tomography (CT) images. However, some small-sized airway branches (e.g., bronchus and terminal bronchioles) significantly aggravate the difficulty of automatic segmentation by machine learning models. In particular, the variance of voxel values and the severe data imbalance in airway branches make the computational module prone to discontinuous and false-negative predictions, especially for cohorts with different lung diseases. The attention mechanism has shown the capacity to segment complex structures, while fuzzy logic can reduce the uncertainty in feature representations. Therefore, the integration of deep attention networks and fuzzy theory, given by the fuzzy attention layer, should be an escalated solution for better generalization and robustness. This article presents an efficient method for airway segmentation, comprising a novel fuzzy attention neural network (FANN) and a comprehensive loss function to enhance the spatial continuity of airway segmentation. The deep fuzzy set is formulated by a set of voxels in the feature map and a learnable Gaussian membership function. Different from the existing attention mechanism, the proposed channel-specific fuzzy attention addresses the issue of heterogeneous features in different channels. Furthermore, a novel evaluation metric is proposed to assess both the continuity and completeness of airway structures. The efficiency, generalization, and robustness of the proposed method have been proved by training on normal lung disease while testing on datasets of lung cancer, COVID-19, and pulmonary fibrosis.

Chest CT-based automated vertebral fracture assessment using artificial intelligence and morphologic features

BACKGROUND

Spinal degeneration and vertebral compression fractures are common among the elderly that adversely affect their mobility, quality of life, lung function, and mortality. Assessment of vertebral fractures in chronic obstructive pulmonary disease (COPD) is important due to the high prevalence of osteoporosis and associated vertebral fractures in COPD.

PURPOSE

We present new automated methods for (1) segmentation and labelling of individual vertebrae in chest computed tomography (CT) images using deep learning (DL), multi-parametric freeze-and-grow (FG) algorithm, and separation of apparently fused vertebrae using intensity autocorrelation and (2) vertebral deformity fracture detection using computed vertebral height features and parametric computational modelling of an established protocol outlined for trained human experts.

METHODS

A chest CT-based automated method was developed for quantitative deformity fracture assessment following the protocol by Genant et al. The computational method was accomplished in the following steps: (1) computation of a voxel-level vertebral body likelihood map from chest CT using a trained DL network; (2) delineation and labelling of individual vertebrae on the likelihood map using an iterative multi-parametric FG algorithm; (3) separation of apparently fused vertebrae in CT using intensity autocorrelation; (4) computation of vertebral heights using contour analysis on the central anterior-posterior (AP) plane of a vertebral body; (5) assessment of vertebral fracture status using ratio functions of vertebral heights and optimized thresholds. The method was applied to inspiratory or total lung capacity (TLC) chest scans from the multi-site Genetic Epidemiology of COPD (COPDGene) (ClinicalTrials.gov: NCT00608764) study, and the performance was examined (n = 3231). One hundred and twenty scans randomly selected from this dataset were partitioned into training (n = 80) and validation (n = 40) datasets for the DL-based vertebral body classifier. Also, generalizability of the method to low dose CT imaging (n = 236) was evaluated.

RESULTS

The vertebral segmentation module achieved a Dice score of .984 as compared to manual outlining results as reference (n = 100); the segmentation performance was consistent across images with the minimum and maximum of Dice scores among images being .980 and .989, respectively. The vertebral labelling module achieved 100% accuracy (n = 100). For low dose CT, the segmentation module produced image-level minimum and maximum Dice scores of .995 and .999, respectively, as compared to standard dose CT as the reference; vertebral labelling at low dose CT was fully consistent with standard dose CT (n = 236). The fracture assessment method achieved overall accuracy, sensitivity, and specificity of 98.3%, 94.8%, and 98.5%, respectively, for 40,050 vertebrae from 3231 COPDGene participants. For generalizability experiments, fracture assessment from low dose CT was consistent with the reference standard dose CT results across all participants.

CONCLUSIONS

Our CT-based automated method for vertebral fracture assessment is accurate, and it offers a feasible alternative to manual expert reading, especially for large population-based studies, where automation is important for high efficiency. Generalizability of the method to low dose CT imaging further extends the scope of application of the method, particularly since the usage of low dose CT imaging in large population-based studies has increased to reduce cumulative radiation exposure.

UNSUPERVISED AIRWAY TREE CLUSTERING WITH DEEP LEARNING: THE MULTI-ETHNIC STUDY OF ATHEROSCLEROSIS (MESA) LUNG STUDY

High-resolution full lung CT scans now enable the detailed segmentation of airway trees up to the 6th branching generation. The airway binary masks display very complex tree structures that may encode biological information relevant to disease risk and yet remain challenging to exploit via traditional methods such as meshing or skeletonization. Recent clinical studies suggest that some variations in shape patterns and caliber of the human airway tree are highly associated with adverse health outcomes, including all-cause mortality and incident COPD. However, quantitative characterization of variations observed on CT segmented airway tree remain incomplete, as does our understanding of the clinical and developmental implications of such. In this work, we present an unsupervised deep-learning pipeline for feature extraction and clustering of human airway trees, learned directly from projections of 3D airway segmentations. We identify four reproducible and clinically distinct airway sub-types in the MESA Lung CT cohort.

Semiautomated Segmentation and Analysis of Airway Lumen in Pediatric Patients Using Ultra Short Echo Time MRI

RATIONALE AND OBJECTIVES

Ultra short echo time (UTE) magnetic resonance imaging (MRI) pulse sequences have shown promise for airway assessment, but the feasibility and repeatability in the pediatric lung are unknown. The purpose of this work was to develop a semiautomated UTE MRI airway segmentation pipeline from the trachea-to-tertiary airways in pediatric participants and assess repeatability and lumen diameter correlations to lung function.

MATERIALS AND METHODS

A total of 29 participants (n = 7 healthy, n = 11 cystic fibrosis, n = 6 asthma, and n = 5 ex-preterm), aged 7-18 years, were imaged using a 3D stack-of-spirals UTE examination at 3 T. Two independent observers performed airway segmentations using a pipeline developed in-house; observer 1 repeated segmentations 1 month later. Segmentations were extracted using region-growing with leak detection, then manually edited if required. The airway trees were skeletonized, pruned, and labeled. Airway lumen diameter measurements were extracted using ray casting. Intra- and interobserver variability was assessed using the Sørensen-Dice coefficient (DSC) and intra-class correlation coefficient (ICC). Correlations between lumen diameter and pulmonary function were assessed using Spearman's correlation coefficient.

RESULTS

For airway segmentations and lumen diameter, intra- and interobserver DSCs were 0.88 and 0.80, while ICCs were 0.95 and 0.89, respectively. The variability increased from the trachea-to-tertiary airways for intra- (DSC: 0.91-0.64; ICC: 0.91-0.49) and interobserver (DSC: 0.84-0.51; ICC: 0.89-0.21) measurements. Lumen diameter was significantly correlated with forced expiratory volume in 1 second and forced vital capacity (P < .05).

CONCLUSION

UTE MRI airway segmentation from the trachea-to-tertiary airways in pediatric participants across a range of diseases is feasible. The UTE MRI-derived lumen measurements were repeatable and correlated with lung function.

Multi-site, Multi-domain Airway Tree Modeling

Open international challenges are becoming the de facto standard for assessing computer vision and image analysis algorithms. In recent years, new methods have extended the reach of pulmonary airway segmentation that is closer to the limit of image resolution. Since EXACT'09 pulmonary airway segmentation, limited effort has been directed to the quantitative comparison of newly emerged algorithms driven by the maturity of deep learning based approaches and extensive clinical efforts for resolving finer details of distal airways for early intervention of pulmonary diseases. Thus far, public annotated datasets are extremely limited, hindering the development of data-driven methods and detailed performance evaluation of new algorithms. To provide a benchmark for the medical imaging community, we organized the Multi-site, Multi-domain Airway Tree Modeling (ATM'22), which was held as an official challenge event during the MICCAI 2022 conference. ATM'22 provides large-scale CT scans with detailed pulmonary airway annotation, including 500 CT scans (300 for training, 50 for validation, and 150 for testing). The dataset was collected from different sites and it further included a portion of noisy COVID-19 CTs with ground-glass opacity and consolidation. Twenty-three teams participated in the entire phase of the challenge and the algorithms for the top ten teams are reviewed in this paper. Both quantitative and qualitative results revealed that deep learning models embedded with the topological continuity enhancement achieved superior performance in general. ATM'22 challenge holds as an open-call design, the training data and the gold standard evaluation are available upon successful registration via its homepage (https://atm22.grand-challenge.org/).

A Survey on Artificial Intelligence in Pulmonary Imaging

Over the last decade, deep learning (DL) has contributed a paradigm shift in computer vision and image recognition creating widespread opportunities of using artificial intelligence in research as well as industrial applications. DL has been extensively studied in medical imaging applications, including those related to pulmonary diseases. Chronic obstructive pulmonary disease, asthma, lung cancer, pneumonia, and, more recently, COVID-19 are common lung diseases affecting nearly 7.4% of world population. Pulmonary imaging has been widely investigated toward improving our understanding of disease etiologies and early diagnosis and assessment of disease progression and clinical outcomes. DL has been broadly applied to solve various pulmonary image processing challenges including classification, recognition, registration, and segmentation. This paper presents a survey of pulmonary diseases, roles of imaging in translational and clinical pulmonary research, and applications of different DL architectures and methods in pulmonary imaging with emphasis on DL-based segmentation of major pulmonary anatomies such as lung volumes, lung lobes, pulmonary vessels, and airways as well as thoracic musculoskeletal anatomies related to pulmonary diseases.

Reproducibility of a combined artificial intelligence and optimal-surface graph-cut method to automate bronchial parameter extraction

OBJECTIVES

Computed tomography (CT)-based bronchial parameters correlate with disease status. Segmentation and measurement of the bronchial lumen and walls usually require significant manpower. We evaluate the reproducibility of a deep learning and optimal-surface graph-cut method to automatically segment the airway lumen and wall, and calculate bronchial parameters.

METHODS

A deep-learning airway segmentation model was newly trained on 24 Imaging in Lifelines (ImaLife) low-dose chest CT scans. This model was combined with an optimal-surface graph-cut for airway wall segmentation. These tools were used to calculate bronchial parameters in CT scans of 188 ImaLife participants with two scans an average of 3 months apart. Bronchial parameters were compared for reproducibility assessment, assuming no change between scans.

RESULTS

Of 376 CT scans, 374 (99%) were successfully measured. Segmented airway trees contained a mean of 10 generations and 250 branches. The coefficient of determination (R2) for the luminal area (LA) ranged from 0.93 at the trachea to 0.68 at the 6th generation, decreasing to 0.51 at the 8th generation. Corresponding values for Wall Area Percentage (WAP) were 0.86, 0.67, and 0.42, respectively. Bland-Altman analysis of LA and WAP per generation demonstrated mean differences close to 0; limits of agreement (LoA) were narrow for WAP and Pi10 (± 3.7% of mean) and wider for LA (± 16.4-22.8% for 2-6th generations). From the 7th generation onwards, there was a sharp decrease in reproducibility and a widening LoA.

CONCLUSION

The outlined approach for automatic bronchial parameter measurement on low-dose chest CT scans is a reliable way to assess the airway tree down to the 6th generation.

STATEMENT ON CLINICAL RELEVANCE

This reliable and fully automatic pipeline for bronchial parameter measurement on low-dose CT scans has potential applications in screening for early disease and clinical tasks such as virtual bronchoscopy or surgical planning, while also enabling the exploration of bronchial parameters in large datasets.

KEY POINTS

• Deep learning combined with optimal-surface graph-cut provides accurate airway lumen and wall segmentations on low-dose CT scans. • Analysis of repeat scans showed that the automated tools had moderate-to-good reproducibility of bronchial measurements down to the 6th generation airway. • Automated measurement of bronchial parameters enables the assessment of large datasets with less man-hours.

Two-stage contextual transformer-based convolutional neural network for airway extraction from CT images

Accurate airway segmentation from computed tomography (CT) images is critical for planning navigation bronchoscopy and realizing a quantitative assessment of airway-related chronic obstructive pulmonary disease (COPD). Existing methods face difficulty in airway segmentation, particularly for the small branches of the airway. These difficulties arise due to the constraints of limited labeling and failure to meet clinical use requirements in COPD. We propose a two-stage framework with a novel 3D contextual transformer for segmenting the overall airway and small airway branches using CT images. The method consists of two training stages sharing the same modified 3D U-Net network. The novel 3D contextual transformer block is integrated into both the encoder and decoder path of the network to effectively capture contextual and long-range information. In the first training stage, the proposed network segments the overall airway with the overall airway mask. To improve the performance of the segmentation result, we generate the intrapulmonary airway branch label, and train the network to focus on producing small airway branches in the second training stage. Extensive experiments were performed on in-house and multiple public datasets. Quantitative and qualitative analyses demonstrate that our proposed method extracts significantly more branches and longer lengths of the airway tree while accomplishing state-of-the-art airway segmentation performance. The code is available at https://github.com/zhaozsq/airway_segmentation.

3D-2D image registration in the presence of soft-tissue deformation in image-guided transbronchial interventions

Purpose. Target localization in pulmonary interventions (e.g. transbronchial biopsy of a lung nodule) is challenged by deformable motion and may benefit from fluoroscopic overlay of the target to provide accurate guidance. We present and evaluate a 3D-2D image registration method for fluoroscopic overlay in the presence of tissue deformation using a multi-resolution/multi-scale (MRMS) framework with an objective function that drives registration primarily by soft-tissue image gradients.Methods. The MRMS method registers 3D cone-beam CT to 2D fluoroscopy without gating of respiratory phase by coarse-to-fine resampling and global-to-local rescaling about target regions-of-interest. A variation of the gradient orientation (GO) similarity metric (denotedGO') was developed to downweight bone gradients and drive registration via soft-tissue gradients. Performance was evaluated in terms of projection distance error at isocenter (PDE_iso). Phantom studies determined nominal algorithm parameters and capture range. Preclinical studies used a freshly deceased, ventilated porcine specimen to evaluate performance in the presence of real tissue deformation and a broad range of 3D-2D image mismatch.Results. Nominal algorithm parameters were identified that provided robust performance over a broad range of motion (0-20 mm), including an adaptive parameter selection technique to accommodate unknown mismatch in respiratory phase. TheGO'metric yielded median PDE_iso= 1.2 mm, compared to 6.2 mm for conventionalGO.Preclinical studies with real lung deformation demonstrated median PDE_iso= 1.3 mm with MRMS +GO'registration, compared to 2.2 mm with a conventional transform. Runtime was 26 s and can be reduced to 2.5 s given a prior registration within ∼5 mm as initialization.Conclusions. MRMS registration via soft-tissue gradients achieved accurate fluoroscopic overlay in the presence of deformable lung motion. By driving registration via soft-tissue image gradients, the method avoided false local minima presented by bones and was robust to a wide range of motion magnitude.

Airway Detection in COPD at Low-Dose CT Using Deep Learning and Multiparametric Freeze and Grow

PURPOSE

To present and validate a fully automated airway detection method at low-dose CT in patients with chronic obstructive pulmonary disease (COPD).

MATERIALS AND METHODS

In this retrospective study, deep learning (DL) and freeze-and-grow (FG) methods were optimized and applied to automatically detect airways at low-dose CT. Four data sets were used: two data sets consisting of matching standard- and low-dose CT scans from the Genetic Epidemiology of COPD (COPDGene) phase II (2014-2017) cohort (n = 2 × 236; mean age ± SD, 70 years ± 9; 123 women); one data set consisting of low-dose CT scans from the COPDGene phase III (2018-2020) cohort (n = 335; mean age ± SD, 73 years ± 8; 173 women); and one data set consisting of low-dose, anonymized CT scans from the 2003 Dutch-Belgian Randomized Lung Cancer Screening trial (n = 55) acquired by using different CT scanners. Performance measures for different methods were computed and compared by using the Wilcoxon signed rank test.

RESULTS

At low-dose CT, 56 294 of 62 480 (90.1%) airways of the reference total airway count (TAC) and 32 109 of 37 864 (84.8%) airways of the peripheral TAC (TAC_p), detected at standard-dose CT, were detected. Significant losses (P < .001) of 14 526 of 76 453 (19.0%) airways and 884 of 6908 (12.8%) airways in the TAC and 12 256 of 43 462 (28.2%) airways and 699 of 3882 (18.0%) airways in the TAC_p were observed, respectively, for the multiprotocol and multiscanner data without retraining. When using the automated low-dose CT method, TAC values of 347, 342, 323, and 266 and TAC_p values of 205, 202, 289, and 141 were observed for those who have never smoked and participants at Global Initiative for Chronic Obstructive Lung Disease stages 0, 1, and 2, respectively, which were superior to the respective values previously reported for matching groups when using a semiautomated method at standard-dose CT.

CONCLUSION

A low-cost, automated CT-based airway detection method was suitable for investigation of airway phenotypes at low-dose CT.Keywords: Airway, Airway Count, Airway Detection, Chronic Obstructive Pulmonary Disease, CT, Deep Learning, Generalizability, Low-Dose CT, Segmentation, Thorax, LungClinical trial registration no. NCT00608764 Supplemental material is available for this article. © RSNA, 2022.

Human treelike tubular structure segmentation: A comprehensive review and future perspectives

Various structures in human physiology follow a treelike morphology, which often expresses complexity at very fine scales. Examples of such structures are intrathoracic airways, retinal blood vessels, and hepatic blood vessels. Large collections of 2D and 3D images have been made available by medical imaging modalities such as magnetic resonance imaging (MRI), computed tomography (CT), Optical coherence tomography (OCT) and ultrasound in which the spatial arrangement can be observed. Segmentation of these structures in medical imaging is of great importance since the analysis of the structure provides insights into disease diagnosis, treatment planning, and prognosis. Manually labelling extensive data by radiologists is often time-consuming and error-prone. As a result, automated or semi-automated computational models have become a popular research field of medical imaging in the past two decades, and many have been developed to date. In this survey, we aim to provide a comprehensive review of currently publicly available datasets, segmentation algorithms, and evaluation metrics. In addition, current challenges and future research directions are discussed.

A deep learning-based post-processing method for automated pulmonary lobe and airway trees segmentation using chest CT images in PET/CT

Background

The proposed algorithm could support accurate localization of lung disease. To develop and validate an automated deep learning model combined with a post-processing algorithm to segment six pulmonary anatomical regions in chest computed tomography (CT) images acquired during positron emission tomography/computed tomography (PET/CT) scans. The pulmonary regions have five pulmonary lobes and airway trees.

Methods

Patients who underwent both PET/CT imaging with an extra chest CT scan were retrospectively enrolled. The pulmonary segmentation of six regions in CT was performed via a convolutional neural network (CNN) of DenseVNet architecture with some post-processing algorithms. Three evaluation metrics were used to assess the performance of this method, which combined deep learning and the post-processing method. The agreement between the combined model and ground truth segmentations in the test set was analyzed.

Results

A total of 640 cases were enrolled. The combined model, which involved deep learning and post-processing methods, had a higher performance than the single deep learning model. In the test set, the all-lobes overall Dice coefficient, Hausdorff distance, and Jaccard coefficient were 0.972, 12.025 mm, and 0.948, respectively. The airway-tree Dice coefficient, Hausdorff distance, and Jaccard coefficient were 0.849, 32.076 mm, and 0.815, respectively. A good agreement was observed between our segmentation in every plot.

Conclusions

The proposed model combining two methods can automatically segment five pulmonary lobes and airway trees on chest CT imaging in PET/CT. The performance of the combined model was higher than the single deep learning model in each region in the test set.

Origins of and lessons from quantitative functional X-ray computed tomography of the lung

Functional CT of the lung has emerged from quantitative CT (qCT). Structural details extracted at multiple lung volumes offer indices of function. Additionally, single volumetric images, if acquired at standardized lung volumes and body posture, can be used to model function by employing such engineering techniques as computational fluid dynamics. With the emergence of multispectral CT imaging including dual energy from energy integrating CT scanners and multienergy binning using the newly released photon counting CT technology, function is tagged via use of contrast agents. Lung disease phenotypes have previously been lumped together by the limitations of spirometry and plethysmography. QCT and its functional embodiment have been imbedded into studies seeking to characterize chronic obstructive pulmonary disease, severe asthma, interstitial lung disease and more. Reductions in radiation dose by an order of magnitude or more have been achieved. At the same time, we have seen significant increases in spatial and density resolution along with methodologic validations of extracted metrics. Together, these have allowed attention to turn towards more mild forms of disease and younger populations. In early applications, clinical CT offered anatomic details of the lung. Functional CT offers regional measures of lung mechanics, the assessment of functional small airways disease, as well as regional ventilation-perfusion matching (V/Q) and more. This paper will focus on the use of quantitative/functional CT for the non-invasive exploration of dynamic three-dimensional functioning of the breathing lung and beating heart within the unique negative pressure intrathoracic environment of the closed chest.

Semi-Supervised Segmentation of Radiation-Induced Pulmonary Fibrosis From Lung CT Scans With Multi-Scale Guided Dense Attention

Computed Tomography (CT) plays an important role in monitoring radiation-induced Pulmonary Fibrosis (PF), where accurate segmentation of the PF lesions is highly desired for diagnosis and treatment follow-up. However, the task is challenged by ambiguous boundary, irregular shape, various position and size of the lesions, as well as the difficulty in acquiring a large set of annotated volumetric images for training. To overcome these problems, we propose a novel convolutional neural network called PF-Net and incorporate it into a semi-supervised learning framework based on Iterative Confidence-based Refinement And Weighting of pseudo Labels (I-CRAWL). Our PF-Net combines 2D and 3D convolutions to deal with CT volumes with large inter-slice spacing, and uses multi-scale guided dense attention to segment complex PF lesions. For semi-supervised learning, our I-CRAWL employs pixel-level uncertainty-based confidence-aware refinement to improve the accuracy of pseudo labels of unannotated images, and uses image-level uncertainty for confidence-based image weighting to suppress low-quality pseudo labels in an iterative training process. Extensive experiments with CT scans of Rhesus Macaques with radiation-induced PF showed that: 1) PF-Net achieved higher segmentation accuracy than existing 2D, 3D and 2.5D neural networks, and 2) I-CRAWL outperformed state-of-the-art semi-supervised learning methods for the PF lesion segmentation task. Our method has a potential to improve the diagnosis of PF and clinical assessment of side effects of radiotherapy for lung cancers.

Computed Tomography Image Texture under Feature Extraction Algorithm in the Diagnosis of Effect of Specific Nursing Intervention on Mycoplasma Pneumonia in Children

To evaluate the effect of specific nursing intervention in children with mycoplasma pneumonia (MP), a feature extraction algorithm based on gray level co-occurrence matrix (GLCM) was proposed and combined with computed tomography (CT) image texture features. Then, 98 children with MP were rolled into the observation group with 49 cases (specific nursing) and the control group with 49 cases (routine nursing). CT images based on feature extraction algorithm of optimized GLCM were used to examine the children before and after nursing intervention, and the recovery of the two groups of children was discussed. The results showed that the proportion of lung texture increase, rope shadow, ground glass shadow, atelectasis, and pleural effusion in the observation group (24.11%, 3.86%, 8.53%, 15.03%, and 3.74%) was significantly lower than that in the control group (28.53%, 10.23%, 13.34%, 21.15%, and 8.13%) after nursing (P < 0.05). There were no significant differences in the proportion of small patchy shadows, large patchy consolidation shadows, and bronchiectasis between the observation group and the control group (P > 0.05). In the course of nursing intervention, in the observation group, the disappearance time of cough, normal temperature, disappearance time of lung rales, and absorption time of lung shadow (2.15 ± 0.86 days, 4.81 ± 1.14 days, 3.64 ± 0.55 days, and 5.96 ± 0.62 days) were significantly shorter than those in the control group (2.87 ± 0.95 days, 3.95 ± 1.06 days, 4.51 ± 1.02 days, and 8.14 ± 1.35 days) (P < 0.05). After nursing intervention, the proportion of satisfaction and total satisfaction in the experimental group (67.08% and 28.66%) was significantly higher than that in the control group (40.21% and 47.39%), while the proportion of dissatisfaction (4.26%) was significantly lower than that in the control group (12.4%) (P < 0.05). To sum up, specific nursing intervention was more beneficial to improve the progress of characterization recovery and the overall recovery effect of children with MP relative to conventional nursing. CT image based on feature extraction algorithm of optimized GLCM was of good adoption value in the diagnosis and treatment of MP in children.

Automatic airway segmentation from computed tomography using robust and efficient 3-D convolutional neural networks

This paper presents a fully automatic and end-to-end optimised airway segmentation method for thoracic computed tomography, based on the U-Net architecture. We use a simple and low-memory 3D U-Net as backbone, which allows the method to process large 3D image patches, often comprising full lungs, in a single pass through the network. This makes the method simple, robust and efficient. We validated the proposed method on three datasets with very different characteristics and various airway abnormalities: (1) a dataset of pediatric patients including subjects with cystic fibrosis, (2) a subset of the Danish Lung Cancer Screening Trial, including subjects with chronic obstructive pulmonary disease, and (3) the EXACT'09 public dataset. We compared our method with other state-of-the-art airway segmentation methods, including relevant learning-based methods in the literature evaluated on the EXACT'09 data. We show that our method can extract highly complete airway trees with few false positive errors, on scans from both healthy and diseased subjects, and also that the method generalizes well across different datasets. On the EXACT'09 test set, our method achieved the second highest sensitivity score among all methods that reported good specificity.