Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 03 July 2025

Sec. Pulmonary Medicine

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1603472

CT radiomics combined with neural networks predict the malignant degree of pulmonary grinding glass nodules

Pengfei ChenPengfei ChenHuiyuan GongHuiyuan GongLei ZhangLei ZhangYang Geng
Yang Geng*
  • Department of Thoracic Surgery, The First Affiliated Hospital of Bengbu Medical University, Bengbu, China

Background: This study investigates the use of CT radiomics combined with convolutional neural networks (CNN) to predict the malignancy of lung ground glass nodules (GGN), which are challenging to diagnose due to their ambiguous boundaries. The goal is to improve diagnostic accuracy and support personalized treatment planning.

Methods: Retrospective data from 670 patients with pulmonary nodules (2019–2023) were analyzed. CT images were preprocessed using Gaussian filtering and manually segmented to define regions of interest (ROI). A CNN model was trained using MATLAB’s Deep Learning Toolbox, and its performance was compared to the Mayo and Brock models.

Results: Key predictors of malignancy included nodule diameter, volume, mean CT value, and consolidation-to-tumor ratio (CTR). The CNN-based model achieved an AUC of 0.887, with 82.4% sensitivity and 75.5% specificity, outperforming existing models (Mayo: AUC = 0.655; Brock: AUC = 0.574). Validation accuracy reached 85.07%.

Conclusion: In this single-center retrospective study, integrating CT radiomics with CNN depicted promising potential for GGN malignancy prediction, though external validation remains necessary. These findings warrant verification in multicenter prospective cohorts.

Introduction

Lung cancer is the second most common type of cancer all over the world, and according to the latest data published by GLOBALCAN, there will be about 2.48 million new cases of lung cancer and 18,187,000 deaths in 2022, which is the highest among all cancers (12.4 and 18.7%) (1). Lung nodules are the early manifestation of lung cancer, and according to their tissue composition, they can be classified as solid nodules and glassy nodules (GGN), and GGN are more challenging for physicians due to their blurred borders (2, 3). GGN is a non-specific manifestation of a cloudy, thin shadow of mildly increased density compared with that of normal lung tissue in chest computed tomography (CT) but does not obstruct the pulmonary vasculature or bronchiolar structures (4). The pathology of GGN may be benign, such as inflammatory infections and interstitial fibrosis, or malignant, such as lung adenocarcinoma (5, 6). Lung adenocarcinoma is the most common histological subtype of lung cancer (~50%) and has an overall poor prognosis due to its highly aggressive tumor, high metastatic rate and low therapeutic efficiency (7). Lung cancer usually has no obvious physical changes in the early stages, and about 70% of lung cancer cases are diagnosed in the late stages, thus missing the best time for treatment (8). Therefore, performing GGN screening is clinically important for early diagnosis of lung cancer.

CT has a critical part in the detection of disease and assessment of efficacy in the field of oncology (9). Studies reported that screening of people at high risk for lung cancer by spiral CT scanning reduced mortality by 20% (10). For GGN, CT scanning requires a wide span of follow-up time for lung cancer surveillance. As a result, it reduces patient compliance and delays diagnosis and treatment, thereby increasing treatment costs and reducing lung cancer survival. In addition, the high cost and potential radiotoxicity of CT make it limited in the prediction of GGN (11). Recently, with the rapid development of Artificial Intelligence (AI), people having used AI to detect CT images has become a major research direction in the field of GGN detection (12). Convolutional neural network (CNN) is the basic class of deep learning neural networks, which is a branch of machine learning technology based on regularized multilayer networks (13), and it is a subtype of AI systems, and currently CNN has been proved to be a promising tool for medical image interpretation and analysis (14). Based on the current domestic and international related research dynamics, CNN, as an AI model, shows certain academic value and application value in benign and malignant screening of lung nodules (15).

Here, we investigated the value of CNN and CT models of GGN in benign and malignant screening of lung nodules, and the benign and malignant prediction results of the CNN model can be used as an auxiliary tool for doctors’ judgement, which can provide more objective and reliable guidance for clinical decision-making, optimize patient management and formulate individualized treatment plans.

Materials and methods

Patient inclusion and allocation

This study intends to retrospectively include inpatients who were diagnosed with pulmonary nodules and underwent surgical treatment in the Department of Thoracic Surgery of the Chest Hospital of the First Affiliated Hospital of Bengbu Medical College from January 2019 to December 2023. Inclusion criteria: (1) lung nodules (diameter ≤3 cm) with a definite pathological diagnosis; (2) high-resolution CT images with good images (layer thickness ≤1.5 mm); (3) lung nodule images can be analyzed in AI-assisted diagnosis system; (4) images are clear enough to obtain specific CT image features; (5) patients with a confirmed diagnosis of GGN on CT examination. Exclusion criteria: (1) postoperative pathology suggesting peripheral metastases, not the primary tumor; (2) images suggesting that the size of the tumor diameter is more than 3 cm; (3) images with severe interstitial lung disease and motion artifacts and affecting the analysis and reading of the AI-assisted diagnostic system; (4) nodules with severe irregularities in shape, which prevented the correct segmentation of the tumor from the extratumoural tissues; and (5) unable to obtain the complete information from the electronic patient record system of our hospital.

A total of 670 patients were included in this study, including 244 benign and 426 malignant cases. Basic information (e.g., age, gender, smoking history, etc.), mean CT value, lobulation sign, burr sign, vascular cluster sign, pleural depression sign, and vacuole sign were collected in all included cases according to the same inclusion and exclusion criteria. This research got approved by the Ethics Committee of the hospital. Each patient and their families knew their consent and voluntarily participated in research.

CNN-based modelling of lung nodules

The MATLAB (R2023a) Deep Learning Toolbox was employed for initial model development due to native DICOM image handling capabilities matching our PACS system, simplified deployment for clinical validation within our institution’s MATLAB-enabled infrastructure and rapid prototyping advantages for multidisciplinary collaboration. The architecture follows standard Keras-style layer definitions ensuring straightforward conversion to TensorFlow/PyTorch.

CT image acquisition and pre-processing

The establishment of the CNN lung nodule model is shown in Figure 1. First, we retrieved the preoperative CT images of the study subject from our hospital using an image archiving and communication system (PACS), and then input them into the MATLAB (MathWorks, Inc., USA) system for preprocessing (Figure 2). To enhance the quality of the images, a Gaussian filter (Gaussian filter) was used for preprocessing. Gaussian filter is a traditional linear filter that has been widely used in image denoising. It works by giving different weights to pixels based on the distance between the pixel and the center of the filter to achieve image smoothing.

Figure 1

Figure 1. CNN based data processing process. The full pipeline from CT image acquisition, preprocessing, ROI segmentation, radiomics feature extraction, to CNN model training and prediction is illustrated.

Figure 2

Figure 2. CT image preprocessing (Gaussian filtering). Gaussian filtering was applied to enhance image quality, reduce noise, and improve the clarity of CT images before segmentation and radiomics analysis.

Image preprocessing is necessary to improve the integrity of lung images as low-quality images can affect subsequent analyses and the effectiveness of the system. By using Gaussian filter, we can reduce the noise in the image and improve the clarity and detail of the image. This helps to better observe and analyze lung images and makes subsequent processing steps more accurate and reliable. It provides a good foundation for subsequent research and analysis.

Radiomic feature extraction

CT radiomic features were extracted using PyRadiomics (v3.0) following Image Biomarker Standardization Initiative (IBSI) guidelines. Features included:

1. First-order statistics: Mean, median, skewness, kurtosis of intensity values.

2. Shape-based: Volume, sphericity, surface area (3D).

3. Texture features:

Gray-level co-occurrence matrix (GLCM): Contrast, correlation, entropy. Gray-level run-length matrix (GLRLM): Run emphasis, non-uniformity.

4. Higher-order: Wavelet-filtered features (LoG sigma = 1.0, 3.0). All features were Z-score normalized and extracted from isotropic-resampled (1 × 1 × 1 mm3) images with fixed bin width (25 HU).

Segmentation of preprocessed CT images

Ground-glass nodules (GGNs) exhibit low contrast and heterogeneous density gradients, particularly in subsolid lesions, which current automated algorithms (e.g., threshold-based or U-Net) frequently missegment. Manual delineation by experienced radiologists (≥5 years in thoracic imaging) was prioritized to preserve morphological details critical for radiomics analysis, account for adjacent anatomical interference (vessels, pleura). Inter-reader consistency received rigorous evaluation. To perform region of interest (ROI) segmentation (Figure 3), we used a manual approach for each preprocessed in the MATLAB (MathWorks, Inc., USA) system for nodules in CT images. In this process, we needed to perform segmentation at each level involving the nodule and manually exclude regions not related to lung tissue, such as bronchi, large vessels, bones and mediastinum. All segmentation was done on thin-layer lung window images.

Figure 3

Figure 3. Segmentation of the region of interest (lung nodule). Manual delineation was performed for each nodule to exclude unrelated structures (e.g., bronchi, vessels, bones), ensuring accuracy in ROI extraction for radiomics.

To assess the agreement between the two investigators, a certain number of patients were randomly selected from the entire sample for comparison. The two investigators segmented the lesions independently and assessed them using the intragroup correlation coefficient (region of interest, ROI). During segmentation, the nodule site, maximum diameter and nodule type were recorded, and a judgement was made as to whether the nodule had a burr sign, whether calcification was present, and whether it was combined with emphysema. If disputes arose, we would discuss and resolve them. Through this segmentation method and evaluation process, we can obtain detailed information about the nodule for subsequent analysis and study.

CNN-based classification model training and validation

To assess model stability, we performed 5-fold cross-validation on the entire cohort (n = 670). The dataset was randomly partitioned into 5 equal subsets, with each subset used once as the validation set while the remaining 4 subsets formed the training set. The final performance metrics were averaged across all folds. After the segmentation process, the segmented region of interest (ROI) is analyzed by MATLAB’s Deep Learning Toolbox to classify the image as benign or malignant based on the postoperative pathology results. The basic principle of CNN is that it extracts the local features from the input and passes them to the lower layers to obtain more complex features. Figure 4 shows the main organization of CNN. The basic structure of CNN can be divided into five parts, i.e., Input Layer, Convolutional Layer, Convergence Layer, Fully Connected Layer and Output Layer.

Figure 4

Figure 4. CNN main organizational structure. The CNN consisted of input, convolutional, pooling, fully connected, and output layers, enabling multilevel feature extraction and classification of nodules as benign or malignant.

With CNN, feature extraction is performed, including first order statistics features (first order statistics), shape features (shape), texture features and higher order statistics features, using minimum redundancy maximum relevance (max-relevance and min-redundancy, mRMR) and minimum absolute shrinkage and selection operator (The basic principle of the LASSO algorithm is to shrink the regression coefficients of all variables by constructing the sum of the absolute values of all regression coefficients with the addition of a first-order paradigm penalty term, i.e., the sum of the absolute values of all regression coefficients). The basic principle of the LASSO algorithm is to shrink the regression coefficients of all variables by constructing the sum of the absolute values of all regression coefficients by adding a first-order paradigm penalty term, i.e., the sum of the absolute values of all regression coefficients and then shrinking the coefficients.

Mayo and Brock model analysis

The Mayo and Brock models are currently recognized nationally and internationally as predictive models for the probability of malignancy of lung nodules, with the probability of malignancy P = ex/(1 + ex), where x takes on different values.

1. Mayo model (16) x = −6.827 2 + (0.0391 × age) + (0.7917 × history of smoking) + (1.338 8 × history of malignancy) + (0.1274 × diameter of nodule) + (1.0407 × hairbrush sign) + (0.7838 × position of upper lobe).

2. Brock model (17) x = −6.614 4 + (0.6467 x sex) + (−5.5537 x diameter) + (0.9309 x burr sign) + (0.6009 x upper lobe position).

Using the pathological diagnosis results as the gold standard, the receiver operating characteristic (ROC) curves of subjects for these 2 models were plotted separately.

Statistical analysis

1. SPSS 25.0 and R4.2. 2 were applied to perform statistical analysis; the mean ± standard deviation was used to express the measurement data that satisfied normal distribution, otherwise, the median and interquartile spacing were used to express the measurement data; frequency, rate and percentage were used for the count data; and the chi-square test or Fisher’s exact probability method and t-test were used for the comparison between groups.

2. The infiltrative nature of lung nodules was used as the dependent variable, the general information of the patients and the imaging characteristics of the nodules were used as the independent variables, and statistically significant (p < 0.05) independent variables were distinguished by univariate analysis. The above independent variables were included in the multifactorial Logistic regression analysis to construct a clinical prediction model for predicting the degree of infiltration of lung nodules with CT image features extracted by the AI-assisted diagnosis system.

3. The malignancy probability collected by the CNN system was used to draw the receiver operating characteristic (ROC) curve, the area under the ROC curve (AUC) was used to evaluate the diagnostic ability assessment of the risk prediction model, and the cut-off value was determined to indicate the critical value for clinical application.

Data availability

The radiomics processing code and synthetic test data are available at GitHub. Due to institutional restrictions and patient privacy protections under ethics approval, the trained model weights and original clinical data cannot be shared. However, all feature extraction procedures are fully described in Methods, all key extracted features are reported in Tables 13, and model architecture details were provided in Figure 4. A Docker image containing all dependencies is provided for replication. Due to IRB restrictions, clinical data cannot be shared but synthetic examples demonstrating the processing pipeline are included.

Table 1
www.frontiersin.org

Table 1. Univariate analysis of infiltrative GGN(mean ± SD) or n (%).

Table 2
www.frontiersin.org

Table 2. Multifactorial logstic regression analysis of infiltrative GGN.

Table 3
www.frontiersin.org

Table 3. Threshold results and test efficacy for independent risk factors.

Results

Patient characteristics

A total of 670 patients were included in this study, including 469 patients in the training set and 201 patients in the validation set. The baseline characteristics of included patients showed no significant differences from excluded cases (n = 535) in age or gender distribution (p > 0.05). Primary exclusion reasons were unanalyzable image quality (53.8%), size criteria (29.3%), and pathological confirmation unavailable (16.9%). The mean age was (58.43 ± 10.12) years, 46.72% were male and 53.28% were female. The imaging characteristics of nodules such as nodule location, diameter, percentage of solid component, and pulmonary nodule characteristics were not significantly different in both groups (p < 0.05, Table 4). Therefore, the selection of patients for the training and validation sets met the requirements for subsequent analyses (Figure 5).

Table 4
www.frontiersin.org

Table 4. Analysis of general clinical data(mean ± SD) or n (%).

Figure 5

Figure 5. Research sample screening flowchart.

Inter-reader agreement analysis

The intraclass correlation coefficient (ICC) for ROI segmentation between two radiologists was 0.82 (95% CI: 0.76–0.87), indicating excellent agreement. For radiomic feature extraction, ICC values ranged from 0.75 to 0.91 across different feature categories (Table 5).

Table 5
www.frontiersin.org

Table 5. Verification of the accuracy of modelling.

Univariate analysis of GGN infiltrability

In order to verify whether the infiltrative nature of GGN was correlated with general clinical data and CT imaging features, a one-way ANOVA was performed on the indicators included in this study. Statistical results revealed that patient age, nodule maximum diameter, volume, mean CT value, consolidation tumor ratio (CTR), burr sign and pleural depression sign were statistically significant between non-infiltrating and infiltrating nodules (p < 0.01).

Feature selection

Of 1,310 initial features, 42 received selection via minimum redundancy maximum relevance (mRMR) and LASSO regression (λ = 0.01). Dominant predictors included: texture: GLCM_Entropy (β = 0.32, p < 0.01); shape: Sphericity (β = −0.21, p < 0.05).

Multifactor regression analysis

The seven indicators screened in the above univariate analysis were included in the multifactorial logistic regression analysis, in which the maximum diameter of the nodule, volume, mean CT value, and CTR were independent predictors (Table 2).

Analysis of the correlation between predictive factors and pulmonary nodal infiltrates

The differential expression of independent factors in infiltrating and non-infiltrating nodules was further analyzed and their correlations were statistically determined. As shown in Figure 6, the maximum diameter of the nodule, volume, mean CT value, and CTR were all significantly highly expressed in infiltrating nodules (p < 0.05). Correlation analysis revealed that the predictive factors were positively correlated with the infiltrative nature of the nodules, with correlation coefficients of 0.642 for the maximum diameter, 0.53 for the volume, 0.59 for the mean CT value, and 0.503 for the proportion of solid components.

Figure 6

Figure 6. Correlation of maximum nodule diameter, volume, mean CT value, CTR and infiltrative nature. Four independent predictive features showed significantly higher values in infiltrative GGNs and were positively correlated with nodule invasiveness.

Clinical prediction modelling

Combining the above independent prognostic factors, the ROC curve was drawn to calculate the prognostic value of GGN. The areas under the curve, 95% confidence intervals, thresholds, sensitivities and specificities of the four independent factors are shown below, and the thresholds of these four indexes were used as the thresholds to construct the composite scores of the prediction model of GGN infiltration. The scoring rule is as follows: prediction model = maximum diameter × 0.642/(0.642 + 0.53 + 0.59 + 0.503) + volume × 0.53/(0.642 + 0.53 + 0.59 + 0.503) + mean CT value × 0.59/(0.642 + 0.53 + 0.59 + 0.503) + CTR × 0.503/(0.642 + 0.53 + 0.59 + 0.503) (maximal diameter ≥12.63 mm was recorded as 1, otherwise 0; volume ≥663.1 mm3 was recorded as 1, otherwise 0; mean CT value ≥-445Hu was recorded as 1, otherwise 0; CTR ≥ 12.60% was recorded as 1, otherwise 0), and infiltrative nodules were considered when the composite score was ≥0.55, and when less than 0.55, infiltrative nodules were considered. When the composite score was ≥0.55, it was considered non-infiltrating nodules. The ROC curve was plotted according to the prediction model, in which the scoring threshold was 0.55. CNN model demonstrated discriminative ability for malignant GGNs (sensitivity of 82.4%, specificity of 75.5%), with AUC of 0.887 and 95% Cl of (0.871–0.898) (Figure 7).

Figure 7

Figure 7. ROC curve analysis. (a) ROC curve analysis of nodule maximum diameter, volume, mean CT value, consolidation tumor ratio versus infiltration. (b) ROC curve analysis for modelling. All four features showed good discrimination (AUC > 0.78), and the composite prediction model achieved higher predictive accuracy (AUC = 0.887).

Correctness validation of predictive models

The accuracy of the modelling was assessed using the validation set, in which the accuracy of the validation set data for maximum diameter volume, mean CT value and CTR were judged correctly at 82.9, 84.07, 79.1, and 77.11%, respectively. The accuracy of the predictive model was 85.07%, which is more effective compared to the test of independent factors.

In the 5-fold cross-validation, the model achieved a mean AUC of 0.872 (95% CI: 0.858–0.886) with sensitivity and specificity of 80.1% ± 3.2 and 74.6% ± 2.8%, respectively, indicating consistent performance across subsets.

Evaluation of model building and AI performance in predicting the benign and malignant nature of lung nodules

The Mayo and Brock model is a nationally and internationally recognized model for predicting the probability of malignancy of lung nodules, and its performance in predicting the benign and malignant nature of nodules was verified by comparing the prediction models constructed in this study. The AUCs of the Mayo and Brock models were 0.655 and 0.574, respectively, whereas that of our prediction model was 0.826, demonstrating a stronger goodness-of-fit (Figure 8).

Figure 8

Figure 8. ROC curve analysis of Mayo models, Brock models and Combinatorial model. Compared to Mayo and Brock models, the combinatorial CNN-based model demonstrated superior predictive performance in malignancy assessment of GGNs.

Discussion

In China, the incidence and mortality rates of lung cancer are at a high level, and it is the first among malignant tumors (18). Recently, adenocarcinoma has gradually replaced squamous carcinoma as the pathology with the highest incidence rate among lung cancers. GGN, as a possible form of early-stage lung adenocarcinoma, has also been gradually emphasized with the early screening of lung cancer (19). It has been found that patients with invasive lung adenocarcinoma also have a significantly lower 5-year survival rate than those with non-invasive adenocarcinoma (20). Therefore, correctly identifying the benign or malignant nature of GGN is important for improving the prognosis of lung cancer patients.

The screening of lung nodules requires the support of high-resolution chest CT, which can clearly present the diameter, CT value, burr sign, lobular sign, pleural depression sign and other imaging features of lung nodules and is an important means of diagnosis for clinicians (21). However, the diagnosis of nodules is a difficult task when the number of nodules is unknown, and time is short. AI is expected to be integrated with medicine at this time to reduce the diagnostic burden of clinicians. Ciompi et al. developed an AI algorithm that can differentiate between different types of nodules with certain reliability (22). Chen et al. based on big data, combined AI software with CT impact characterization can improve its diagnostic and screening ability (23). Ardila et al. (24) used 42,290 CT data from the NLST dataset and constructed a composite convolutional neural network with a three-dimensional Inception network as its core. It identified and extracted both global and local features of chest CT images, as well as analyzed chest CT images at different time points, and verified that the AUC reached 0.944, achieving the highest existing accuracy of artificial intelligence in predicting the risk of lung cancer. A large number of studies have revealed the application value of CNN model in benign and malignant screening of lung nodules. Based on the existing studies, the CNN model constructed in this study has higher efficiency and accuracy in feature extraction of nodule regions and can classify lung nodules as benign or malignant through accurate classification and judgement, and give the corresponding prediction probability, which strengthens the validity and feasibility of CNN in the clinical practice. Feasibility of CNN in clinical practice.

In this study, we firstly investigated the relevant factors affecting the infiltration of GGN, and the results revealed that there was a significant difference in age, maximum diameter of the nodule, volume, mean CT value, CTR, burr sign and pleural depression sign between patients with invasive and non-invasive GGN (p < 0.05) The age of patients with invasive GGN (60.04 ± 8.87) was higher than that of the non-invasive group (57.65 ± 11.2), which may be due to the fact that increasing age increases the risk of immune escape of tumor cells, thus increasing the incidence of lung cancer. Zhang et al. revealed that age was one of the independent risk factors for invasive and non-invasive GGN using CT (25), which is in line with the results of our study. Fleischner’s guideline states that nodule size is one of the most important factors in differentiating the benign and malignant nature of GGN (26), and the probability of malignancy increases with every 2 mm increase in the maximum diameter of a subsolid nodule (27). Among the 469 patients included in the training set of this study, the maximum diameter of infiltrating GGN (16.23 ± 4.65) was significantly larger than that of the non-infiltrating group (11.02 ± 5.11), similar to Jiang et al. who reported that the diameter of the nodule correlated with the degree of infiltration, which is a predictor of infiltrating adenocarcinoma (28). Further correlation analysis between the maximum diameter and the infiltrative nature of GGN showed that they were positively correlated to some extent and the correlation coefficient was 0.642. In addition, some studies revealed that the volume of the nodule was also related to the infiltrative nature of GGN (29), and our study similarly revealed a positive correlation between the volume and the infiltrative nature of GGN as well, and it could be used as one of the independent predictive factors of GGN. The density of the material is linearly related to the CT value, which has been shown by several studies to correlate with the degree of tumor infiltration (30), where tumor cells continuously proliferate along the alveolar wall, resulting in a local increase in alveolar density. The present study revealed that CT value is an independent predictor of infiltrative GGN, which is in line with the findings of Junji Ichinose et al. that maximum CT value is an independent predictor of tissue invasiveness (31).

Several studies have pointed out that the area of the solid component reflects the degree of infiltration of tumor cells and that CTR and its CT value correlate with tumor invasiveness (32, 33). In addition, burr sign and pleural depression are common features of malignant GGN22-23. In this study, burr sign and pleural depression and age were found to be risk factors for the degree of infiltration of GGN by univariate analysis, but the inclusion of multivariate analysis of roar was found not to have an effect on the degree of infiltration of GGN.

The novelty of this study is the combination of four variables: maximum diameter, volume, mean CT value, and CTR to infer the infiltrative nature of the nodule, determine its threshold value and develop a prediction model to predict the probability of malignant nodules. The accuracy of the determined thresholds was explored using a validation set, and the prediction model showed a certain degree of improvement in sensitivity, specificity, and correctness compared with the individual variables. The prediction model constructed in this study also has certain advantages over the Mayo and Brock models, which are recognized nationally and internationally for predicting the malignancy probability of lung nodules. While the CNN model demonstrated superior discriminative ability in imaging-based diagnosis, Mayo and Brock models retain clinical utility for rapid risk stratification without advanced imaging analysis. The hybrid approach’s performance gain suggests potential synergies between radiomic and clinical features-a direction for future multimodal AI development.

Compared to existing radiomics studies, our work specifically addresses two key literature gaps: integration of deep learning with handcrafted radiomic features for GGN analysis, and systematic comparison against clinical risk models. However, like multiple prior studies, our current validation remains institution-specific. While we employed rigorous cross-validation to mitigate overfitting, we fully acknowledge the necessity of external testing—a limitation shared by approximately 78% of similar AI-radiomics studies according to recent reviews. Our group will initiate a multicenter trial to formally assess generalizability across different scanner types and populations. While MATLAB served our research needs effectively, we recognized industry preference for TensorFlow/PyTorch in production environments. We have open-sourced a PyTorch version achieving identical performance and containerized the model using Docker for cloud deployment. The single-center retrospective design may introduce selection bias, and external validation is required to confirm generalizability. The term “effective” in this context specifically refers to internal validation results, not established clinical utility. Future multicenter prospective studies with larger cohorts will be essential to translate this model into clinical practice. Moreover, the method used to validate the predictive performance of the model in this study is limited to logistic regression, and more methods such as random forests and elasticity network regression should be added for comparison to develop a model with the best predictive performance. Additionally, while we openly share the implementation code and processing pipeline, institutional policies prevent distribution of the trained model weights and original clinical data. Researchers may replicate the approach using the provided architecture details and synthetic data or contact the corresponding author for collaborative verification. While manual segmentation ensures precision in this exploratory study, its labor-intensive nature limits clinical scalability. Recent advances in deep learning show promise for GGN segmentation but require large, annotated datasets-a focus of our ongoing work. A semi-automated pipeline with manual correction may balance efficiency and accuracy in future implementations. Furthermore, while our evaluation focused on AUC, sensitivity and specificity-metrics most clinically relevant for malignancy screening-we recognize that additional measures like precision, F1-score and confusion matrix analysis could provide deeper model characterization. These metrics will be prioritized in our upcoming multicenter validation study where class imbalance management becomes more critical.

The benign and malignant diagnosis of GGN, as well as the aggressiveness of the nodule is the focus and difficulty in clinical work, and in the imaging diagnosis and follow-up of lung adenocarcinoma, the CT imaging features have a high degree of accuracy in determining the degree of tumor infiltration. In this study, the correlation between clinical features and GGN infiltration was analyzed by CNN, and a prediction model of GGN infiltration probability with high specificity and sensitivity was constructed, so as to provide a strong imaging basis for the clinicians to grasp the patient’s condition quickly and to make the choice of appropriate surgical plan.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

The studies involving humans were approved by ethics committee of Bengbu Medical University (Approval number 2025-123. Approval date: 2025.4.10). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

PC: Writing – original draft, Writing – review & editing. HG: Writing – original draft, Writing – review & editing. LZ: Writing – original draft, Writing – review & editing. YG: Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Bengbu Medical College Key Project of Natural Science Research (2023byzd048) and the Key Project of Natural Science Research in Anhui Provincial Institutions of Higher Learning (2023AH051999).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1603472/full#supplementary-material

References

1. Bray, F, Laversanne, M, Sung, H, Ferlay, J, Siegel, RL, Soerjomataram, I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834

PubMed Abstract | Crossref Full Text | Google Scholar

2. Mazzone, PJ, and Lam, L. Evaluating the patient with a pulmonary nodule: a review. JAMA. (2022) 327:264–73. doi: 10.1001/jama.2021.24287

PubMed Abstract | Crossref Full Text | Google Scholar

3. Ma, H, Guo, H, Zhao, M, Qi, S, Li, H, Tian, Y, et al. Automatic pulmonary ground-glass opacity nodules detection and classification based on 3D neural network. Med Phys. (2022) 49:2555–69. doi: 10.1002/mp.15501

PubMed Abstract | Crossref Full Text | Google Scholar

4. Zhang, Y, Fu, F, and Chen, H. Management of Ground-Glass Opacities in the lung Cancer Spectrum. Ann Thorac Surg. (2020) 110:1796–804. doi: 10.1016/j.athoracsur.2020.04.094

PubMed Abstract | Crossref Full Text | Google Scholar

5. Gao, JW, Rizzo, S, Ma, LH, Qiu, XY, Warth, A, Seki, N, et al. Pulmonary ground-glass opacity: computed tomography features, histopathology and molecular pathology. Transl Lung Cancer Res. (2017) 6:68–75. doi: 10.21037/tlcr.2017.01.02

PubMed Abstract | Crossref Full Text | Google Scholar

6. Huang, KX, and Gibney, BC. Pure ground-glass opacities (GGO) lung adenocarcinoma: surgical resection is curative. J Thorac Dis. (2024) 16:3518–21. doi: 10.21037/jtd-23-1983

PubMed Abstract | Crossref Full Text | Google Scholar

7. Jia, Y, Ji, Q, Zhang, L, She, Y, Su, M, and Shi, Z. Prognosis of early-stage lung adenocarcinoma in young patients. Clin Exp Pharmacol Physiol. (2023) 50:826–32. doi: 10.1111/1440-1681.13806

PubMed Abstract | Crossref Full Text | Google Scholar

8. Iqbal, MA, Arora, S, Prakasam, G, Calin, GA, and Syed, MA. MicroRNA in lung cancer: role, mechanisms, pathways and therapeutic relevance. Mol Asp Med. (2019) 70:3–20. doi: 10.1016/j.mam.2018.07.003

PubMed Abstract | Crossref Full Text | Google Scholar

9. DiCorpo, D, Tiwari, A, Tang, R, Griffin, M, Aftreth, O, Bautista, P, et al. The role of Micro-CT in imaging breast cancer specimens. Breast Cancer Res Treat. (2020) 180:343–57. doi: 10.1007/s10549-020-05547-z

PubMed Abstract | Crossref Full Text | Google Scholar

10. Aberle, DR, Adams, AM, Berg, CD, Black, WC, Clapp, JD, Fagerstrom, RM, et al. Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. (2011) 365:395–409. doi: 10.1056/NEJMoa1102873

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ren, S, Zhang, S, Jiang, T, He, Y, Ma, Z, Cai, H, et al. Early detection of lung cancer by using an autoantibody panel in Chinese population. Onco Targets Ther. (2018) 7:e1384108. doi: 10.1080/2162402X.2017.1384108

PubMed Abstract | Crossref Full Text | Google Scholar

12. Jungmann, F, Müller, L, Hahn, F, Weustenfeld, M, Dapper, AK, Mähringer-Kunz, A, et al. Commercial AI solutions in detecting COVID-19 pneumonia in chest CT: not yet ready for clinical implementation? Eur Radiol. (2022) 32:3152–60. doi: 10.1007/s00330-021-08409-4

PubMed Abstract | Crossref Full Text | Google Scholar

13. Lim, JJ, Kim, DW, Hong, WH, Kim, M, Lee, DH, Kim, SY, et al. Application of convolutional neural network (CNN) to recognize ship structures. Sensors (Basel). (2022) 22:3824. doi: 10.3390/s22103824

PubMed Abstract | Crossref Full Text | Google Scholar

14. Poonkuzhali, P, and Prabha, HK. Deep convolutional neural network based hyperspectral brain tissue classification. J Xray Sci Technol. (2023) 31:777–96. doi: 10.3233/XST-230045

PubMed Abstract | Crossref Full Text | Google Scholar

15. Yuan, H, Fan, Z, Wu, Y, and Cheng, J. An efficient multi-path 3D convolutional neural network for false-positive reduction of pulmonary nodule detection. Int J Comput Assist Radiol Surg. (2021) 16:2269–77. doi: 10.1007/s11548-021-02478-y

PubMed Abstract | Crossref Full Text | Google Scholar

16. Papalampidou, A, Papoutsi, E, and Katsaounou, PA. Pulmonary nodule malignancy probability: a diagnostic accuracy meta-analysis of the Mayo model. Clin Radiol. (2022) 77:443–50. doi: 10.1016/j.crad.2022.01.055

PubMed Abstract | Crossref Full Text | Google Scholar

17. Owens, C, Hindocha, S, Lee, R, Millard, T, and Sharma, B. The lung cancers: staging and response, CT, (18)F-FDG PET/CT, MRI, DWI: review and new perspectives. Br J Radiol. (2023) 96:20220339. doi: 10.1259/bjr.20220339

PubMed Abstract | Crossref Full Text | Google Scholar

18. Li, Y, Wu, X, Yang, P, Jiang, G, and Luo, Y. Machine learning for lung Cancer diagnosis, treatment, and prognosis. Genom Proteom Bioinform. (2022) 20:850–66. doi: 10.1016/j.gpb.2022.11.003

PubMed Abstract | Crossref Full Text | Google Scholar

19. Liang, X, Zhang, C, and Ye, X. Overdiagnosis and overtreatment of ground-glass nodule-like lung cancer. Asia Pac J Clin Oncol. (2024) 21:14042. doi: 10.1111/ajco.14042

PubMed Abstract | Crossref Full Text | Google Scholar

20. Zhang, G, Dong, R, Kong, D, Liu, B, Zha, Y, and Luo, M. The effect of GLUT1 on the survival rate and immune cell infiltration of lung adenocarcinoma and squamous cell carcinoma: a meta and bioinformatics analysis. Anti Cancer Agents Med Chem. (2022) 22:223–38. doi: 10.2174/1871520621666210708115406

PubMed Abstract | Crossref Full Text | Google Scholar

21. Groheux, D, Quere, G, Blanc, E, Lemarignier, C, Vercellino, L, de Margerie-Mellon, C, et al. FDG PET-CT for solitary pulmonary nodule and lung cancer: literature review. Diagn Interv Imaging. (2016) 97:1003–17. doi: 10.1016/j.diii.2016.06.020

PubMed Abstract | Crossref Full Text | Google Scholar

22. Ciompi, F, Chung, K, van Riel, SJ, Setio, AAA, Gerke, PK, Jacobs, C, et al. Towards automatic pulmonary nodule management in lung cancer screening with deep learning. Sci Rep. (2017) 7:46479. doi: 10.1038/srep46479

PubMed Abstract | Crossref Full Text | Google Scholar

23. Chen, L, Gu, D, Chen, Y, Shao, Y, Cao, X, Liu, G, et al. An artificial-intelligence lung imaging analysis system (ALIAS) for population-based nodule computing in CT scans. Comput Med Imaging Graph. (2021) 89:101899. doi: 10.1016/j.compmedimag.2021.101899

PubMed Abstract | Crossref Full Text | Google Scholar

24. Ardila, D, Kiraly, AP, Bharadwaj, S, Choi, B, Reicher, JJ, Peng, L, et al. End-to-end lung cancer screening with three-dimensional deep learning on low-dose chest computed tomography. Nat Med. (2019) 25:954–61. doi: 10.1038/s41591-019-0447-x

PubMed Abstract | Crossref Full Text | Google Scholar

25. Zhang, WZ, Zhang, YY, Yao, XL, Li, PL, Chen, XY, He, LY, et al. Computed tomography radiomics study of invasion and instability of lung adenocarcinoma manifesting as ground glass nodule. J Thorac Dis. (2024) 16:3828–43. doi: 10.21037/jtd-24-27

PubMed Abstract | Crossref Full Text | Google Scholar

26. Wang, Q. Management strategies of pulmonary ground Galss nodule. Zhongguo Fei Ai Za Zhi. (2018) 21:160–2. doi: 10.3779/j.issn.1009-3419.2018.03.06

PubMed Abstract | Crossref Full Text | Google Scholar

27. Cao, H, Liu, H, Song, E, Ma, G, Xu, X, Jin, R, et al. A two-stage convolutional neural networks for lung nodule detection. IEEE J Biomed Health Inform. (2020) 24:2006–15. doi: 10.1109/JBHI.2019.2963720

PubMed Abstract | Crossref Full Text | Google Scholar

28. Jiang, C, Zhao, M, and Zhang, W. Application of 18 F-fluorodeoxyglucose PET/computed tomography in the diagnosis of infiltrative subsolid nodules in lung adenocarcinoma. Nucl Med Commun. (2024) 45:1082–91. doi: 10.1097/MNM.0000000000001908

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kou, J, Gu, X, and Kang, L. Correlation analysis of computed tomography features and pathological types of multifocal ground-glass nodular lung adenocarcinoma. Comput Math Methods Med. (2022) 2022:1–10. doi: 10.1155/2022/7267036

PubMed Abstract | Crossref Full Text | Google Scholar

30. He, W, Guo, G, Du, X, Guo, S, and Zhuang, X. CT imaging indications correlate with the degree of lung adenocarcinoma infiltration. Front Oncol. (2023) 13:1108758. doi: 10.3389/fonc.2023.1108758

PubMed Abstract | Crossref Full Text | Google Scholar

31. Ichinose, J, Kawaguchi, Y, Nakao, M, Matsuura, Y, Okumura, S, Ninomiya, H, et al. Utility of maximum CT value in predicting the invasiveness of pure ground-glass nodules. Clin Lung Cancer. (2020) 21:281–7. doi: 10.1016/j.cllc.2020.01.015

PubMed Abstract | Crossref Full Text | Google Scholar

32. Ono, Y, Tagawa, T, Kinoshita, F, Haratake, N, Takada, K, Kohno, M, et al. Relationship between consolidation tumor ratio and tumor-infiltrating lymphocytes in small-sized lung adenocarcinoma. Thorac Cancer. (2022) 13:2134–41. doi: 10.1111/1759-7714.14524

PubMed Abstract | Crossref Full Text | Google Scholar

33. Yoon, DW, Kim, CH, Hwang, S, Choi, YL, Cho, JH, Kim, HK, et al. Reappraising the clinical usability of consolidation-to-tumor ratio on CT in clinical stage IA lung cancer. Insights Imaging. (2022) 13:103. doi: 10.1186/s13244-022-01235-2

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: convolutional neural network, lung ground glass nodules, CT radiomics, infiltration, predictive modelling

Citation: Chen P, Gong H, Zhang L and Geng Y (2025) CT radiomics combined with neural networks predict the malignant degree of pulmonary grinding glass nodules. Front. Med. 12:1603472. doi: 10.3389/fmed.2025.1603472

Received: 31 March 2025; Accepted: 22 May 2025;
Published: 03 July 2025.

Edited by:

Sandeep Kumar Mishra, Yale University, United States

Reviewed by:

Palash Ghosal, Sikkim Manipal University, India
Junxiang Huang, Boston College, United States

Copyright © 2025 Chen, Gong, Zhang and Geng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yang Geng, R2VuZ3lhbmd4d2tAaG90bWFpbC5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.