3D computed tomography airway geometry for predicting bronchoscopic accessibility in peripheral pulmonary nodules: a prospective study
Original Article

3D computed tomography airway geometry for predicting bronchoscopic accessibility in peripheral pulmonary nodules: a prospective study

Byeong-Ho Jeong1#, Jonghoon Kim2,3#, Hwanho Cho4, Youjin Oh2,3, Ho Yun Lee2,3

1Division of Pulmonary and Critical Care Medicine, Department of Medicine, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea; 2Department of Health Sciences and Technology, Samsung Advanced Institute for Health Science & Technology (SAIHST), Sungkyunkwan University, Seoul, Republic of Korea; 3Department of Radiology and Center for Imaging Science, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea; 4Department of Electronics Engineering, Incheon National University, Incheon, Republic of Korea

Contributions: (I) Conception and design: BH Jeong, HY Lee; (II) Administrative support: HY Lee, BH Jeong; (III) Provision of study materials or patients: BH Jeong; (IV) Collection and assembly of data: BH Jeong, J Kim, Y Oh; (V) Data analysis and interpretation: J Kim, Y Oh, H Cho; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Ho Yun Lee, MD, PhD. Department of Health Sciences and Technology, Samsung Advanced Institute for Health Science & Technology (SAIHST), Sungkyunkwan University, Seoul, Republic of Korea; Department of Radiology and Center for Imaging Science, Samsung Medical Center, Sungkyunkwan University School of Medicine, 81 Irwon-Ro, Gangnam-Gu, Seoul 06351, Republic of Korea. Email: hoyunlee96@gmail.com.

Background: Bronchoscopic diagnosis of peripheral pulmonary nodules (PPNs) using radial probe endobronchial ultrasound (rEBUS) has been widely studied. However, the prospective utility of computed tomography (CT)-derived airway geometrical analysis for predicting rEBUS accessibility remains underexplored. This study aimed to determine whether quantitative CT airway geometric characteristics can improve preprocedural planning by predicting the bronchoscopic accessibility of PPNs.

Methods: The accessibility of 219 PPNs in 199 patients to rEBUS was prospectively evaluated as easily accessible or difficult/inaccessible. Preprocedural airway geometry was quantified from individual CT scans and used both to guide the procedure and develop predictive models using logistic-least absolute shrinkage and selection operator (LASSO) analyses. Model performance was assessed using the area under the curve (AUC).

Results: Of the 219 PPNs, 182 (83.1%) were easily accessible, and 37 (16.9%) were difficult or inaccessible. The mean age was 68.7 (±10.3) years, the average size of the PPNs was 24.5 (±16.1) mm, the mean distance from the pleural surface was 10.8 (±12.5) mm, and 83.1% showed a ‘within’ bronchus sign. Airway geometrical features—including more acute bifurcation angles, more sharply curved branches, and more elliptical and narrower lumen shapes—were significant risk factors for limited bronchoscopic access. A composite model incorporating these features with clinical variables achieved the highest performance (AUC =0.84), better than both the clinical model (AUC =0.68) and the airway geometrical model (AUC =0.79).

Conclusions: CT-based quantitative airway analysis enhances the prediction of PPN accessibility to rEBUS and could support more accurate procedural planning in clinical practice.

Keywords: Peripheral pulmonary nodule (PPN); radial probe endobronchial ultrasound (rEBUS); quantitative computed tomography (quantitative CT); airway geometry; accessibility


Submitted Jan 30, 2026. Accepted for publication Apr 08, 2026. Published online May 26, 2026.

doi: 10.21037/tlcr-2026-1-0133


Highlight box

Key findings

• Quantitative computed tomography-derived airway geometric features—specifically bifurcation angle, curvature, and luminal narrowing—are independent predictors of bronchoscopic accessibility for peripheral pulmonary nodules.

• A composite model integrating these 3D geometric metrics with clinical variables demonstrated improved predictive performance (area under the curve =0.84), offering incremental value beyond models based on conventional clinical or radiologic features alone.

What is known and what is new?

• Bronchoscopic accessibility is a prerequisite for successful radial probe endobronchial ultrasound-guided biopsy, yet its prediction has traditionally relied on subjective operator intuition or simple, non-quantitative bronchus sign assessments.

• This study establishes a data-driven framework that converts 3D airway morphology into quantifiable geometric biomarkers. By integrating metrics like bifurcation angles and curvature, this approach provides reproducible navigation mapping that replaces subjective qualitative assessments.

What is the implication, and what should change now?

• Preprocedural planning using quantitative airway analysis enables early identification of challenging pathways, allowing clinicians to optimize bronchoscopic strategies and select appropriate navigational tools.

• Integrating these objective geometric assessments into routine clinical practice can reduce operator-dependent bias and improve the overall efficiency and safety of peripheral lung nodule biopsies.


Introduction

Diagnosis of peripheral pulmonary nodules (PPNs) through bronchoscopy has been extensively investigated in recent literature (1). Bronchoscopic interventions, such as radial probe endobronchial ultrasound (rEBUS)-guided transbronchial lung biopsy (TBLB) based on computed tomography (CT)-guided approaches, have been extensively explored (2). Preprocedural chest CT imaging enables precise procedural planning by delineating key anatomical structures, such as lesion location and bronchial pathways, whereas rEBUS provides intraprocedural real-time localization of the lesion relative to the probe, enhancing the accuracy of imaging-guided bronchoscopic interventions when other forms of intraprocedural imaging are unavailable (3,4).

However, PPN identification with rEBUS depends heavily on operator knowledge of bronchial anatomy and procedural difficulty may vary according to the anatomical course of the target bronchus (5,6). Certain bronchial pathways may present greater mechanical challenges during scope advancement, underscoring the importance of additional preprocedural information for procedural planning. When a bronchoscope is advanced into an upper lobar bronchus, the direction of pushing force applied by the operator is opposite to the direction in which the bronchoscope needs to proceed (6). Consequently, accessibility to PPNs in the upper lobes can be lower than in other lobes. These anatomical constraints underscore the need for additional preprocedural information to support bronchoscopic navigation.

Lung airway morphometry using chest CT has emerged as a powerful tool for quantifying airway geometry and assessing structural abnormalities (7). Prior research has proposed CT-derived quantitative airway measures as predictors of lung function decline and disease progression (8-10). However, those findings have largely been limited to observational studies, and their utility in guiding real-world clinical procedures remains underexplored. In our previous retrospective study, we demonstrated the feasibility of CT-based airway analysis in preoperative planning (11). Although our results suggested potential clinical value, their application in prospective clinical settings had not been validated.

Therefore, in this study, we extend our previous work by implementing a prospective design to develop and clinically validate a predictive model for estimating the accessibility of PPNs during rEBUS-TBLB. We further assessed the utility of our predictive model as a preprocedural guide in routine bronchoscopic workflows. We present this article in accordance with the TRIPOD reporting checklist (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/rc).


Methods

Study participants

The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Institutional Review Board of Samsung Medical Center (No. 2018-03-021). Informed consent was taken from all the patients (or a statement that it was not required and why). Patient information was de-identified and anonymized prior to the analysis.

We collected data from 212 patients whose rEBUS-TBLB procedures were planned during a prospective trial at Samsung Medical Center, a single tertiary referral center, between December 2021 and February 2023. Since the procedures were performed using rEBUS alone—without additional guidance from fluoroscopy, cone-beam CT, guide sheaths, or navigation systems—patients were primarily referred for the procedure when a bronchus leading to the lesion was identifiable on pre-procedural CT to ensure procedural feasibility. During the initial clinical screening, and prior to any quantitative airway analysis, 10 patients were excluded for the following reasons: the operator determined that the lesion was clinically unsuitable for biopsy (n=6), the lesions were small pure ground-glass opacities (n=3), or the lesion improved on follow-up imaging before the procedure (n=1). Quantitative airway analysis was then attempted in the remaining 202 patients. However, three patients were further excluded because of processing failures due to technical or software issues. Consequently, a total of 219 PPNs from 199 patients were included in the final analysis. Because several patients presented with multiple PPNs, each lesion was prospectively assessed and included in the analysis to reflect the full clinical spectrum encountered in real-world practice. The overall study flow is illustrated in Figure 1.

Figure 1 Flow chart of patient inclusion. EBUS-TBLB, endobronchial ultrasound-guided transbronchial lung biopsy; GGN, ground-glass nodule.

CT acquisition and preoperative evaluation

All patients underwent chest CT scanning prior to rEBUS-TBLB as part of the prospective study protocol. CT scans were acquired with standard thin slices (≤1 mm) and high-resolution reconstruction; detailed scanning parameters are available in a supplementary appendix online (Appendix 1).

Several radiologic characteristics were visually assessed on chest CT: lesion size, lesion depth from the pleural surface, number of bronchial branches from the trachea to the lesion, and lesion solidity (solid, sub-solid, or cavitary). The lesion location was categorized into five anatomical regions: right upper, right middle, right lower, left upper, and left lower lobes. Furthermore, to account for anatomical and physiological factors affecting procedural difficulty, we categorized the segments into three bronchial clusters: (I) segments requiring acute angulation for access, specifically the apical and superior segments (RB1, LB1, RB6, LB6); (II) segments with high respiratory-induced mobility, including the basal segments (RB8–10, LB8–10); and (III) all other segments. This classification reflects established challenges in bronchoscopy, such as the acute angles of the upper lobes (6) and the significant respiratory displacement common in the lower lobes (12). All CT image assessments were performed by a board-certified thoracic radiologist (H.Y.L.) with 18 years of experience in thoracic imaging.

To assess preoperative airway information, we performed three-dimensional airway segmentation and lesion localization using chest CT data and commercial deep-learning software (AVIEW; Coreline Soft, Seoul, Korea). The segmented airway and lesion volumes were reviewed by an expert radiologist (H.Y.L.). Quantitative airway geometrical features were extracted using in-house software (AirGeo) that was developed and examined in a previous study (11). The extracted features included cross-sectional dimensions (e.g., minimum diameter, perimeter, circularity), spatial orientation (e.g., curvature, torsion), and branching characteristics (e.g., bifurcation angle). The detailed quantitative analysis process and a comprehensive list of airway geometry variables are provided in Appendix 1 and Table S1.

The extracted CT-based airway information was provided to the bronchoscopist prior to the procedure, enabling preoperative planning. It is important to note that the AirGeo-based analysis provided only raw quantitative data (e.g., subsegmental routes and bronchial angles) for preoperative planning. It did not offer any qualitative assessment of predicted difficulty. Thus, it served as an informative tool to guide the anatomical approach, rather than as a behavioral intervention. The proceduralist’s effort was standardized across all cases, and the study’s objective was to determine which of these objective preoperative parameters were significantly associated with the actual technical difficulty encountered during the procedure. Figure 2 shows the overall workflow of our study.

Figure 2 Overall workflow of the prospective study design. CT, computed tomography; EBUS-TBLB, endobronchial ultrasound-guided transbronchial lung biopsy; LA, luminal area; PPNs, peripheral pulmonary nodules.

Procedure for rEBUS-TBLB

All rEBUS-TBLB procedures were performed on patients under conscious sedation using midazolam and fentanyl, as previously reported (13). Briefly, a 4.1-mm bronchoscope (BF P260F; Olympus, Tokyo, Japan) was advanced to the sub-subsegmental bronchus nearest the lesion after reviewing the CT. Then, the rEBUS probe (1.4-mm, 20-MHz, UM S20-17S; Olympus) was inserted through the bronchoscope’s working channel. Upon confirming lesion visualization via ultrasound, a 1.8-mm biopsy forceps was inserted through the working channel for TBLB after the probe was removed. All procedures were performed by one operator (B.H.J.), a bronchoscopist with 11 years of clinical experience, without a guide sheath, fluoroscopy, or navigation system.

When the target was identified on ultrasound, the findings were classified as “within” (the radial probe was within the target lesion, surrounded by lesion material) or “adjacent” (the probe did not pierce the lesion but was nearby). When the nodule was not able to be visualized by the probe, it was defined as “invisible”, indicating failure of endobronchial navigation using thin bronchoscopy.

Independent of the relationship between the probe and the lesion, the difficulty in reaching the lesion was defined by the operator as “easily accessible”, “difficult but accessible”, or “inaccessible”. To provide a more objective assessment, “easily accessible” was defined as reaching the target lesion within one or two attempts after identifying the bronchial path on preprocedural CT. “Difficult” access (including “difficult but accessible” and “inaccessible”) was defined as cases requiring more than three attempts (>3 attempts) to successfully navigate the probe to the target.

We analyzed 219 procedures from 199 patients (some with multiple lesions) to validate the prediction model. Additionally, the procedure time was recorded from radial probe insertion to biopsy completion to assess procedural efficiency.

Statistical analysis

In this study, we focused on airway geometric characteristics as three levels of feature sets. First, six global level features were defined using five first-order statistics about the bifurcation angles of the total branches and the total length of each sectional branch. Next, 64 features were considered from the last two branches leading to the target lesion because they are likely most relevant for accessing the lesion. All cases were grouped by final accessibility (easily accessible vs. difficult or inaccessible).

For the prediction model, multicollinearity among the extracted features was assessed using the variance inflation factor (VIF), excluding features with VIF >4 (14). Univariate logistic regression identified variables associated with accessibility, and redundant variables (P<0.05) were removed. Feature selection and predictive scoring were performed using logistic-least absolute shrinkage and selection operator (LASSO) with 10-fold cross-validation to determine the optimal penalty term. The predictive score for each patient was calculated as a weighted sum of the selected features.

We developed three prediction models: clinical, airway geometrical, and composite (combining both feature sets). Model performance was evaluated using receiver operating characteristic (ROC) curves and corresponding area under the curve (AUC) values. DeLong’s test was applied to compare model classification performance (15). Calibration was assessed by plotting observed versus predicted probabilities, where an ideal nomogram would closely follow the 45-degree reference line (16).

Normality of continuous variables was assessed using the Shapiro-Wilk test. Because most continuous variables were not normally distributed, they are presented as median [interquartile range (IQR)] and were compared using the Mann-Whitney U test. Categorical variables were compared using the Chi-squared test. All statistical analyses were performed using both SPSS version 20.0.0.0 (IBM Corp., Armonk, NY, USA) and R 3.2.5 software (R Foundation for Statistical Computing, Vienna, Austria). A P value less than 0.05 was considered statistically significant.


Results

Characteristics of cohorts

Demographics for the enrolled patients are shown in Table 1. The patients [median age, 70.0 years (IQR, 62.0–77.0 years); 118 male and 81 female] had lesions with three types of solidity: solid (n=150, 68.5%), sub-solid (n=47, 21.5%), and cavitary (n=22, 10.0%). Compared with easily accessible cases, difficult or inaccessible cases had significantly smaller target lesions [19.2 mm (IQR, 15.6–22.8 mm) vs. 22.8 mm (IQR, 17.7–29.5 mm), P=0.008]. There was no significant difference between groups in the number of branch levels (P=0.56). Although lesion depth from the pleural surface did not reach statistical significance, it was numerically smaller in difficult or inaccessible cases [4.3 mm (IQR, 0–10.4 mm) vs. 6.9 mm (IQR, 0–16.9 mm), P=0.08]. The distribution of the PPNs in the five lobes did not differ statistically between the groups. However, in terms of bronchus clustering, more difficult or inaccessible cases than easily accessible cases tended to be found in the RB1/LB1/RB6/LB6 cluster (45.9% vs. 25.8%, P=0.059). The ‘invisible’ bronchus sign was significantly more frequent in the difficult or inaccessible group than the easily accessible group (16.2% vs. 3.3%; P=0.007). While the “within” was more frequently observed in the easily accessible group (155/182, 85.2%) than in the difficult or inaccessible group (27/37, 73.0%), its presence in a substantial proportion of the difficult or inaccessible group cases was notable.

Table 1

Baseline characteristics

Variables Total cases (n=219) Easily accessible (n=182) Difficult or inaccessible (n=37) P
Age (years) 70.0 [62.0–77.0] 70.0 [62.0–77.0] 68.0 [62.0–75.0] 0.71
Male 118/199 (59.3) 98/166 (59.0) 20/33 (60.6) 0.58
Lesion size (mm) 22.4 [17.5–28.0] 22.8 [17.7–29.5] 19.2 [15.6–22.8] 0.008
Lesion depth from pleural surface (mm) 6.8 [0.0–15.8] 6.9 [0–16.9] 4.3 [0–10.4] 0.08
Solidity 0.91
   Solid 150 (68.5) 124 (68.1) 26 (70.3)
   Sub-solid 47 (21.5) 39 (21.4) 8 (21.6)
   Cavitary 22 (10.0) 19 (10.4) 3 (8.1)
Number of branch levels 7.0 [6.0–8.0] 7.0 [6.0–8.0] 7.0 [6.0–8.0] 0.56
Lesion location 0.12
   Right upper 47 (21.5) 37 (20.3) 10 (27.0)
   Right middle 16 (7.3) 16 (8.8) 0
   Right lower 47 (21.5) 35 (19.2) 12 (32.4)
   Left upper 69 (31.5) 59 (32.4) 10 (27.0)
   Left lower 40 (18.3) 35 (19.2) 5 (13.5)
Bronchus clustering 0.059
   RB1, LB1, RB6, LB6 64 (29.2) 47 (25.8) 17 (45.9)
   RB8–10, LB8–10 57 (26.0) 50 (27.5) 7 (18.9)
   Others 98 (44.7) 85 (46.7) 13 (35.1)
Bronchus sign on CT 0.007
   Within 182 (83.1) 155 (85.2) 27 (73.0)
   Adjacent 25 (11.4) 21 (11.5) 4 (10.8)
   Invisible 12 (5.5) 6 (3.3) 6 (16.2)

Data are reported as median [interquartile range] or number (%). P values for continuous variables were calculated using the Mann-Whitney U test, and P values for categorical variables were calculated using the Chi-squared test. We performed rEBUS-TBLB for 219 PPNs in 199 patients (19 patients had multiple PPNs). , age and sex were calculated for 199 patients. All remaining variables were analyzed for 219 PPNs. CT, computed tomography; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; LB8–10, anterior, lateral, and posterior segmental bronchi of the left lower lobar bronchus; PPN, peripheral pulmonary nodule; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus; RB8–10, anterior, lateral, and posterior segmental bronchi of the right lower lobar bronchus; rEBUS-TBLB, radial probe endobronchial ultrasound-guided transbronchial lung biopsy.

Procedure results are shown in Table S2. The procedure time was longer in the difficult or inaccessible cases than in the easily accessible cases (15.8±8.7 vs. 9.4±5.0 min, P<0.001). Radial probe positioning within the target lesion during the procedure was more common in the easily accessible cases than in the difficult or inaccessible cases (89.0% vs. 45.9%, P<0.001). Diagnostic results from pathology were obtained for 136/182 (74.7%) easily accessible cases and only 10/37 (27.0%) difficult or inaccessible cases (P<0.001).

Risk factors associated with difficult or inaccessible PPNs

We performed univariable logistic regression analyses to identify clinical variables (Table 2) associated with difficult or inaccessible PPNs. Larger lesion size [per 1 mm, odds ratio (OR) =0.931, P=0.007] and centrality (1 mm farther from pleura, OR =0.963, P=0.046) are likely to increase target lesion accessibility. For bronchus clustering, lower accessibility was observed if the targeted lesion was located at RB1, RB6, LB1, or LB6 (OR =2.288, P=0.044). When the bronchus sign on preprocedural CT was ‘invisible’ (OR =5.741, P=0.004) or the procedure time was longer (per 1 min, OR =1.156, P<0.001), bronchoscopic access to the target lesion was more likely to be difficult.

Table 2

Univariable analysis for clinical variables associated with difficult or inaccessible PPNs

Variables OR (95% CI) P
Age 0.995 (0.961–1.031) 0.79
Male (vs. female) 1.219 (0.596–2.494) 0.58
Lesion size 0.931 (0.883–0.981) 0.007
Lesion depth from pleural surface 0.963 (0.927–0.999) 0.046
Solidity
   Sub-solid (vs. solid) 0.978 (0.041–2.336) 0.96
   Cavitary (vs. solid) 0.753 (0.208–2.733) 0.66
Number of branch levels 1.040 (0.853–1.269) 0.69
Lesion location
   Left lower (vs. left upper) 0.843 (0.266–2.667) 0.77
   Right upper and middle (vs. left upper) 1.113 (0.430–2.884) 0.82
   Right lower (vs. left upper) 2.023 (0.792–5.166) 0.14
Bronchus clustering
   RB8–10, LB8–10 (vs. others) 0.905 (0.338–2.418) 0.84
   RB1, LB1, RB6, LB6 (vs. others) 2.288 (1.024–5.116) 0.044
Bronchus sign on CT before the procedure
   Adjacent (vs. within) 1.093 (0.348–3.435) 0.87
   Invisible (vs. within) 5.741 (1.724–19.120) 0.004
Procedure time 1.156 (1.090–1.226) <0.001

CI, confidence interval; CT, computed tomography; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; LB8–10, anterior, lateral, and posterior segmental bronchi of the left lower lobar bronchus; OR, odds ratio; PPNs, peripheral pulmonary nodules; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus; RB8–10, anterior, lateral, and posterior segmental bronchi of the right lower lobar bronchus.

Several airway geometrical features were significantly associated with difficult or inaccessible PPNs (Table S3). Among bifurcation angle metrics at the total branch level, higher sum (OR =1.505, P=0.02), maximum (OR =1.595, P=0.01), and average (OR =1.636, P=0.006) values were linked to reduced accessibility. At the branch preceding the final segment, a higher minimum of the min-max diameter ratio indicated greater accessibility (OR =0.641, P=0.02). At the final branch level, increased values of the average max inscribed sphere R (OR =0.623, P=0.01), maximum hydraulic luminal diameter (OR =0.349, P=0.040), and average hydraulic luminal diameter (OR =0.496, P=0.008) were associated with better accessibility. Conversely, greater maximum and average curvature were linked to inaccessibility (OR =1.512, P=0.01; OR =1.474, P=0.03, respectively). Representative cases illustrating these associations are shown in Figure 3.

Figure 3 Representative cases of peripheral pulmonary nodules demonstrating airway pathway visualization. Left: cavitary lesion in LB3 of a 66-year-old male patient, located close to the pleural surface. Right: solid lesion in LB4 of a 75-year-old male patient, positioned farther from the pleural surface within the lung parenchyma. Both cases showed a “within” bronchus sign on preprocedural computed tomography. Airway geometry maps display the reconstructed pathway from the trachea to each lesion, with segmental coloring indicating sequential airway branches and an enlarged view highlighting the distal approach route. LB3, anterior segmental bronchus of the left upper lobe; LB4, superior lingular bronchus of the left upper lobe.

Performance of the prediction models

To identify patients with PPNs that were difficult to access with rEBUS, we applied the logistic-LASSO method to develop a prediction model and assessed the performance metrics using 10-fold cross-validation (Table 3). In the final composite model, smaller lesion size; closeness to the pleural surface; location in RB1, RB6, LB1, or LB6; invisible bronchus sign; larger bifurcation angle-in; smaller min-max diameter ratio; larger curvature; and smaller hydraulic luminal diameter all predicted difficult PPN access for rEBUS.

Table 3

Multivariable logistic-LASSO analysis of the three predictive models associated with difficult or inaccessible PPNs

Proposed model Variables OR (95% CI)
Clinical model
   Clinical variables Lesion size 0.280 (0.231–0.340)
Lesion depth from pleural surface 0.467 (0.439–0.498)
Bronchus clustering (RB1, RB6, LB1, LB6) 1.372 (1.294–1.454)
Bronchus sign (invisible) 1.347 (1.308–1.387)
Airway geometrical model
   Total branches Sum of bifurcation angle-in 1.195 (1.118–1.277)
Max of bifurcation angle-in 1.621 (1.490–1.763)
Average of bifurcation angle-in 1.073 (1.022–1.127)
   Previous branch to final branch Min of min–max diameter ratio 0.634 (0.591–0.680)
   Final branch Average of max inscribed sphere R 0.789 (0.732–0.850)
Max curvature 1.585 (1.474–1.704)
Max hydraulic luminal diameter 0.528 (0.486–0.575)
Composite model
   Clinical variables Lesion size 0.363 (0.299–0.440)
Lesion depth from pleural surface 0.471 (0.415–0.533)
Bronchus clustering (RB1, RB6, LB1, LB6) 1.278 (1.183–1.381)
Bronchus sign (invisible) 1.355 (1.288–1.425)
   Total branches Sum of bifurcation angle-in 1.214 (1.131–1.302)
Max of bifurcation angle-in 1.436 (1.313–1.570)
   Previous branch to final branch Min of min–max diameter ratio 0.563 (0.518–0.611)
   Final branch Max curvature 1.377 (1.294–1.464)
Average curvature 1.105 (1.044–1.169)
Max hydraulic luminal diameter 0.481 (0.417–0.554)

CI, confidence interval; LASSO, least absolute shrinkage and selection operator; LB1, apical segmental bronchus of the left upper lobar bronchus; LB6, superior segmental bronchus of the left lower lobar bronchus; OR, odds ratio; PPNs, peripheral pulmonary nodules; RB1, apical segmental bronchus of the right upper lobar bronchus; RB6, superior segmental bronchus of the right lower lobar bronchus.

The composite model, which combines clinical variables with airway geometrical variables, performed better than either model alone [AUC (95% CI); composite model =0.84 (0.76–0.92), clinical model =0.68 (0.60–0.81), airway geometrical model =0.79 (0.74–0.84)] (Table 4, Figure 4). The Delong test showed significant model differences (P=0.05). Calibration curves showed good agreement between predicted and actual accessibility, indicating well-developed models (Figure S1). The composite model showed greater consistency between observed and predicted probabilities.

Table 4

Performance evaluation and statistics for the three proposed models

Models Performance and statistics
AUC Sensitivity Specificity Delong test
Clinical Airway geometrical
Clinical 0.68 (0.60–0.81) 0.15 (0.02–0.28) 0.99 (0.98–1.00) Reference
Airway geometrical 0.79 (0.74–0.84) 0.36 (0.22–0.50) 0.99 (0.97–1.00) 0.43 Reference
Composite 0.84 (0.76–0.92) 0.51 (0.28–0.74) 0.98 (0.96–1.00) 0.05 0.28

AUC, sensitivity, and specificity are represented as mean (95% confidence interval). Model performance was evaluated using 10-fold cross validation. The Delong test was performed to compare the model performance (e.g., AUC). AUC, area under the curve.

Figure 4 ROC curves of the three predictive models. AUC values are presented as AUC (95% confidence interval, 95% CI). AUC, area under the curve; ROC, receiver operating characteristic.

Discussion

In this study, we prospectively evaluated the accessibility of 219 PPNs (median size, 22.4 mm; median distance from the pleura, 6.8 mm) to rEBUS. To create a model for predicting accessibility, we analyzed both clinical and airway geometrical factors. The composite model showed superior predictive performance (AUC =0.84) compared with the clinical model (0.68). These results suggest that objective quantitative assessment using CT-based airway geometrical data may be more useful than clinical information alone.

Several meta-analyses have reported pooled diagnostic yields of approximately 70% for PPNs measuring 15–40 mm when using rEBUS (17-19). In our study, the median lesion size was 22.4 mm and the diagnostic yield was 66.7%, suggesting comparable lesion characteristics. Unlike previous studies focusing on diagnostic yield, we analyzed the accessibility of PPNs to rEBUS. Because diagnostic yield varies substantially depending on malignancy status and biopsy technique (20), we focused on anatomical accessibility using CT-based airway information. Clinical factor analysis showed that larger lesion size, farther distance from the pleura, and presence of the bronchus sign were associated with easier access—findings consistent with previous meta-analyses (18,21,22). Moreover, we assessed accessibility by bronchus clustering rather than lobe-wise locations and found that apical segments of the upper lobes (RB1, LB1) and posterior segments of the lower lobes (RB6, LB6) were significantly less accessible, likely due to their greater curvature (5).

In the composite model, six airway geometrical features were selected as significant risk factors for difficult or inaccessible cases (Table S4). These included larger bifurcation angles, lower min-max diameter ratio, greater curvature, and smaller hydraulic luminal diameter. Clinically, these features indicate that more acutely angled, irregularly shaped, and narrowed distal airways present barriers to successful rEBUS navigation—findings that are consistent with our previous pilot study (11). The present results not only validate these risk factors in a prospective cohort but also extend prior work by integrating both clinical and quantitative airway data within a unified predictive model.

Unlike earlier studies, our analysis directly compared the predictive value of clinical variables, airway geometrical features, and their combination. The composite model achieved a higher AUC (0.84) than either the clinical (0.68) or airway geometrical (0.79) models alone, with notably improved sensitivity. These findings support the clinical utility of combining quantitative CT-based airway measurements with conventional factors to better predict bronchoscopic accessibility.

Notably, 73.0% of the lesions classified as difficult or inaccessible exhibited a ‘within’ bronchus sign on preprocedural CT, highlighting the limitations of conventional visual assessment. The bronchus sign can often fail to capture underlying anatomical complexities such as sharp bifurcation angles, severe curvature, or airway narrowing. In contrast, our quantitative CT-based airway analysis could provide a more objective and comprehensive assessment, enabling early identification of difficult cases—even when CT appearances seemed favorable. This approach has the potential to enhance preprocedural planning and procedural success, supporting its integration into real-world bronchoscopic workflows.

The proposed CT-based airway geometric analysis is not intended to replace existing bronchoscopic navigation systems, but rather to complement them by providing objective preprocedural information. Unlike many navigation technologies that are applied primarily during bronchoscopy, our approach uses quantitative CT-derived airway measurements to estimate procedural difficulty before the intervention, thereby supporting earlier and more structured decision-making. Because this framework relies on reproducible geometric parameters rather than visual bronchial assessment alone, it may reduce subjectivity in procedural planning. Furthermore, by identifying lesions likely to be difficult to access in advance, this method may help clinicians selectively allocate adjunctive navigation or robotic systems to cases in which they are most likely to provide incremental benefit, thereby improving procedural efficiency and potential cost-effectiveness.

The procedures in this study were performed under conscious sedation without the use of a guide sheath, fluoroscopy, or navigation systems. While these adjuncts have been shown to enhance diagnostic yield in several studies (5,23-25), their routine use across all cases may be limited by medical costs and resource availability (26). Recent meta-analyses of high-cost robotic-assisted bronchoscopy (RAB) reported diagnostic yields approaching 80% (27,28). However, as noted in previous literature, RAB studies often involve a higher proportion of challenging targets, including smaller nodules, bronchus sign-negative lesions, and subsolid opacities, which are inherently more difficult to sample than the lesions typically included in standalone rEBUS studies. Therefore, rather than a direct comparison of diagnostic yields, we propose that our CT-driven airway geometrical analysis (AirGeo) could serve as a tool for procedural triage. By preoperatively identifying lesions predicted to be “inaccessible” via standalone rEBUS, clinicians could strategically refer these high-difficulty cases for advanced platforms like RAB or navigation-guided bronchoscopy under general anesthesia. Such a risk-stratified approach could optimize both clinical outcomes and cost-effectiveness by matching the complexity of the lesion with the appropriate level of technological intervention.

Our study has several limitations. First, although the study was conducted prospectively, it was performed at a single center, and all bronchoscopic procedures were carried out by one experienced bronchoscopist. This design may limit the generalizability of the findings to other institutions, procedural settings, and operators with different levels of experience. However, the use of a single operator also reduced procedural heterogeneity and enabled a more consistent assessment of the relationship between CT-derived airway geometrical features and bronchoscopic accessibility. Future multicenter studies including bronchoscopists with varying levels of expertise are warranted to confirm the external validity and reproducibility of the proposed approach. Another limitation is that the operational definition of bronchoscopic accessibility—although based on the number of procedural attempts, and therefore more structured than a purely subjective assessment—may still be affected by operator-dependent factors such as procedural technique, persistence, and intraprocedural judgment. Accordingly, this definition should be interpreted as a pragmatic procedural surrogate rather than a fully anatomy-specific measure of lesion accessibility. In this study, the use of a single experienced bronchoscopist helped reduce inter-operator variability; however, residual operator dependence cannot be excluded. Future studies are needed to validate this concept using more standardized and less operator-sensitive criteria. Second, our proposed model was not externally validated due to the limited dataset. However, to minimize overfitting and improve model robustness, we conducted internal validation using cross-validation during logistic-LASSO regression modeling. Third, we evaluated only accessibility and not the diagnostic yield of rEBUS for PPNs. As noted in recent literature, diagnostic yield is heavily influenced by confounding variables such as the type of biopsy instruments used (e.g., forceps vs. cryobiopsy), the prevalence of malignancy, and the pathological nature of the lesion (19,20). By focusing strictly on accessibility—defined as the ability to navigate the probe to the target—we aimed to isolate the impact of airway geometric factors from these variables and provide a purer assessment of navigational predictability. Despite those limitations, our study identified the potential of CT-based quantitative assessment for enhancing bronchoscopic procedures and improving PPN accessibility to rEBUS.


Conclusions

In conclusion, this study prospectively evaluated the bronchoscopic accessibility of PPNs by constructing predictive models that incorporate CT-derived airway geometrical features. By integrating these imaging-derived variables with clinical information, we demonstrated their clinical significance in preprocedural planning. Our findings validated the utility of a quantitative airway analysis as an objective tool for assessing procedural difficulty and support its potential role in real-world bronchoscopic workflows. The prospective study design further strengthens the evidence that CT-based airway geometry holds practical relevance and promising potential for improving bronchoscopic strategy and decision-making in clinical practice.


Acknowledgments

None.


Footnote

Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/rc

Data Sharing Statement: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/dss

Peer Review File: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/prf

Funding: This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (No. NRF-2022R1A2C1003999) and by the Future Medicine 20*30 Project of Samsung Medical Center (No. SMO1250061) and partly supported by an Institute of Information & Communications Technology Planning & Evaluation grant funded by the Korean government (No. RS-2021-II212068, Artificial Intelligence Innovation Hub).

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2026-1-0133/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Institutional Review Board of Samsung Medical Center (No. 2018-03-021). Informed consent was taken from all the patients (or a statement that it was not required and why).

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Shen YC, Chen CH, Tu CY. Advances in Diagnostic Bronchoscopy. Diagnostics (Basel) 2021;11:1984. [Crossref] [PubMed]
  2. Steinfort DP, Khor YH, Manser RL, et al. Radial probe endobronchial ultrasound for the diagnosis of peripheral lung cancer: systematic review and meta-analysis. Eur Respir J 2011;37:902-10. [Crossref] [PubMed]
  3. Chen A, Chenna P, Loiselle A, et al. Radial probe endobronchial ultrasound for peripheral pulmonary lesions. A 5-year institutional experience. Ann Am Thorac Soc 2014;11:578-82.
  4. Pritchett MA, Bhadra K, Calcutt M, et al. Virtual or reality: divergence between preprocedural computed tomography scans and lung anatomy during guided bronchoscopy. J Thorac Dis 2020;12:1595-611. Erratum in: J Thorac Dis 2020;12:4593-5.
  5. Kurimoto N, Miyazawa T, Okimasa S, et al. Endobronchial ultrasonography using a guide sheath increases the ability to diagnose peripheral pulmonary lesions endoscopically. Chest 2004;126:959-65. [Crossref] [PubMed]
  6. Oki M, Saka H, Imabayashi T, et al. Guide sheath versus non-guide sheath method for endobronchial ultrasound-guided biopsy of peripheral pulmonary lesions: a multicentre randomised trial. Eur Respir J 2022;59:2101678. [Crossref] [PubMed]
  7. Kirby M, Smith BM, Quantitative CT. Scan Imaging of the Airways for Diagnosis and Management of Lung Disease. Chest 2023;164:1150-8. [Crossref] [PubMed]
  8. Diaz AA, Rahaghi FN, Ross JC, et al. Understanding the contribution of native tracheobronchial structure to lung function: CT assessment of airway morphology in never smokers. Respir Res 2015;16:23. [Crossref] [PubMed]
  9. Nakano Y, Muro S, Sakai H, et al. Computed tomographic measurements of airway dimensions and emphysema in smokers. Correlation with lung function. Am J Respir Crit Care Med 2000;162:1102-8.
  10. Washko GR, Diaz AA, Kim V, et al. Computed tomographic measures of airway morphology in smokers and never-smoking normals. J Appl Physiol 1985;2014:668-73. [Crossref] [PubMed]
  11. Cho HH, Choe J, Kim J, et al. 3D airway geometry analysis of factors in airway navigation failure for lung nodules. Cancer Imaging 2024;24:84. [Crossref] [PubMed]
  12. Salahuddin M, Sarkiss M, Sagar AS, et al. Ventilatory Strategy to Prevent Atelectasis During Bronchoscopy Under General Anesthesia: A Multicenter Randomized Controlled Trial (Ventilatory Strategy to Prevent Atelectasis -VESPA- Trial). Chest 2022;162:1393-401. [Crossref] [PubMed]
  13. Moon SM, Choe J, Jeong BH, et al. Diagnostic Performance of Radial Probe Endobronchial Ultrasound without a Guide-Sheath and the Feasibility of Molecular Analysis. Tuberc Respir Dis (Seoul) 2019;82:319-27. [Crossref] [PubMed]
  14. Miles J. Tolerance and variance inflation factor. In: Everitt BS, Howell DC. editors. Encyclopedia of Statistics in Behavioral Science. John Wiley & Sons, Ltd.; 2005. Available online: https://onlinelibrary.wiley.com/doi/epdf/10.1002/9781118445112.stat06593
  15. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 1988;44:837-45.
  16. Lagendijk J, Steyerberg EW, Daalderop LA, et al. Validation of a prognostic model for adverse perinatal health outcomes. Sci Rep 2020;10:11243. [Crossref] [PubMed]
  17. McGuire AL, Myers R, Grant K, et al. The Diagnostic Accuracy and Sensitivity for Malignancy of Radial-Endobronchial Ultrasound and Electromagnetic Navigation Bronchoscopy for Sampling of Peripheral Pulmonary Lesions: Systematic Review and Meta-analysis. J Bronchology Interv Pulmonol 2020;27:106-21. [Crossref] [PubMed]
  18. Lee J, Song JU. Diagnostic yield of radial probe endobronchial ultrasonography-guided transbronchial biopsy without fluoroscopy in peripheral pulmonary lesions: A systematic review and meta-analysis. Thorac Cancer 2023;14:195-205. [Crossref] [PubMed]
  19. Sryma PB, Mittal S, Madan NK, et al. Efficacy of Radial Endobronchial Ultrasound (R-EBUS) guided transbronchial cryobiopsy for peripheral pulmonary lesions (PPL...s): A systematic review and meta-analysis. Pulmonology 2023;29:50-64. [Crossref] [PubMed]
  20. Ali MS, Trick W, Mba BI, et al. Radial endobronchial ultrasound for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis. Respirology 2017;22:443-53. [Crossref] [PubMed]
  21. Tian S, Li X, Liu J, et al. Radial endobronchial ultrasound - guided bronchoscopy for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis of prospective trials. Heliyon 2024;10:e29446. [Crossref] [PubMed]
  22. Kim SH, Chung HS, Kim J, et al. Development of the Korean Association for Lung Cancer Clinical Practice Guidelines: Recommendations on Radial Probe Endobronchial Ultrasound for Diagnosing Lung Cancer - An Updated Meta-Analysis. Cancer Res Treat 2024;56:464-83. [Crossref] [PubMed]
  23. Jiang S, Xie F, Mao X, et al. The value of navigation bronchoscopy in the diagnosis of peripheral pulmonary lesions: A meta-analysis. Thorac Cancer 2020;11:1191-201. [Crossref] [PubMed]
  24. Tamiya M, Okamoto N, Sasada S, et al. Diagnostic yield of combined bronchoscopy and endobronchial ultrasonography, under LungPoint guidance for small peripheral pulmonary lesions. Respirology 2013;18:834-9. [Crossref] [PubMed]
  25. Eberhardt R, Anantham D, Ernst A, et al. Multimodality bronchoscopic diagnosis of peripheral lung lesions: a randomized controlled trial. Am J Respir Crit Care Med 2007;176:36-41. [Crossref] [PubMed]
  26. Vakil E, Jackson N, Sainz-Zuñega PV, et al. Optimizing Diagnostic and Staging Pathways for Suspected Lung Cancer: A Decision Analysis. Chest 2021;160:2304-23. [Crossref] [PubMed]
  27. Zhang C, Xie F, Li R, et al. Robotic-assisted bronchoscopy for the diagnosis of peripheral pulmonary lesions: A systematic review and meta-analysis. Thorac Cancer 2024;15:505-12. [Crossref] [PubMed]
  28. Ali MS, Ghori UK, Wayne MT, et al. Diagnostic Performance and Safety Profile of Robotic-assisted Bronchoscopy: A Systematic Review and Meta-Analysis. Ann Am Thorac Soc 2023;20:1801-12. [Crossref] [PubMed]
Cite this article as: Jeong BH, Kim J, Cho H, Oh Y, Lee HY. 3D computed tomography airway geometry for predicting bronchoscopic accessibility in peripheral pulmonary nodules: a prospective study. Transl Lung Cancer Res 2026;15(5):132. doi: 10.21037/tlcr-2026-1-0133

Download Citation