Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 18 November 2021
Sec. Cancer Imaging and Image-directed Interventions

Diagnostic Performance of 2D and 3D T2WI-Based Radiomics Features With Machine Learning Algorithms to Distinguish Solid Solitary Pulmonary Lesion

Qi Wan&#x;Qi Wan1†Jiaxuan Zhou&#x;Jiaxuan Zhou1†Xiaoying XiaXiaoying Xia1Jianfeng HuJianfeng Hu1Peng WangPeng Wang1Yu PengYu Peng1Tianjing ZhangTianjing Zhang2Jianqing SunJianqing Sun2Yang SongYang Song3Guang YangGuang Yang3Xinchun Li*Xinchun Li1*
  • 1Department of Radiology, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
  • 2Philips Healthcare, Guangzhou, China
  • 3Shanghai Key Laboratory of Magnetic Resonance, East China Normal University, Shanghai, China

Objective: To evaluate the performance of 2D and 3D radiomics features with different machine learning approaches to classify SPLs based on magnetic resonance(MR) T2 weighted imaging (T2WI).

Material and Methods: A total of 132 patients with pathologically confirmed SPLs were examined and randomly divided into training (n = 92) and test datasets (n = 40). A total of 1692 3D and 1231 2D radiomics features per patient were extracted. Both radiomics features and clinical data were evaluated. A total of 1260 classification models, comprising 3 normalization methods, 2 dimension reduction algorithms, 3 feature selection methods, and 10 classifiers with 7 different feature numbers (confined to 3–9), were compared. The ten-fold cross-validation on the training dataset was applied to choose the candidate final model. The area under the receiver operating characteristic curve (AUC), precision-recall plot, and Matthews Correlation Coefficient were used to evaluate the performance of machine learning approaches.

Results: The 3D features were significantly superior to 2D features, showing much more machine learning combinations with AUC greater than 0.7 in both validation and test groups (129 vs. 11). The feature selection method Analysis of Variance(ANOVA), Recursive Feature Elimination(RFE) and the classifier Logistic Regression(LR), Linear Discriminant Analysis(LDA), Support Vector Machine(SVM), Gaussian Process(GP) had relatively better performance. The best performance of 3D radiomics features in the test dataset (AUC = 0.824, AUC-PR = 0.927, MCC = 0.514) was higher than that of 2D features (AUC = 0.740, AUC-PR = 0.846, MCC = 0.404). The joint 3D and 2D features (AUC=0.813, AUC-PR = 0.926, MCC = 0.563) showed similar results as 3D features. Incorporating clinical features with 3D and 2D radiomics features slightly improved the AUC to 0.836 (AUC-PR = 0.918, MCC = 0.620) and 0.780 (AUC-PR = 0.900, MCC = 0.574), respectively.

Conclusions: After algorithm optimization, 2D feature-based radiomics models yield favorable results in differentiating malignant and benign SPLs, but 3D features are still preferred because of the availability of more machine learning algorithmic combinations with better performance. Feature selection methods ANOVA and RFE, and classifier LR, LDA, SVM and GP are more likely to demonstrate better diagnostic performance for 3D features in the current study.

Introduction

A solitary pulmonary lesion (SPL) is one of the most common findings on chest radiographs and computed tomography (CT). An increasing number of pulmonary nodules are detected by CT with the improvement in lung cancer screening. However, most of these positive detections are not cancerous (1). The high false-positive rate can lead to a waste of medical resources, additional radiation exposure, unnecessary patient anxiety, and so on. Recent advances in magnetic resonance imaging (MRI) techniques make it possible to use lung MRI in routine clinical practice. Published evidence showed that lung MRI could be a potentially effective screening tool because its performance was comparable with that of low-dose CT (2), even with a lower false-positive rate for nodule detection (3). A conventional sequence, such as T2 weighted imaging(T2WI), has the potential to detect pulmonary nodules no less than 6 mm in diameter (4), which is essential for screening. However, as a morphological sequence, it may have limited value in distinguishing malignant from benign SPLs.

Radiomics is an emerging field that extracts a large number of quantitative features from medical imaging and quantifies tumor heterogeneity related to cellularity, necrosis, and angiogenesis in the tumor microenvironment (5). Therefore, radiomics provides the possibility for early and accurate diagnosis of SPLs. Radiomics can increase the diagnostic accuracy of baseline CT (6). In addition, studies have shown the potential of radiomics based on CT and MRI in distinguishing pulmonary lesions (7, 8). However, MR radiomics research focusing on differentiating SPLs has not yet been reported. Also, the performance of 2D and 3D CT features in pulmonary lesions has been shown to be controversial in different studies (9, 10). However, the performance of 2D and 3D MR features as well as their corresponding optimal machine learning methods in distinguishing SPLs has not been discussed. These issues are crucial for further generalization of MR radiomics in clinical research and application in the lung.

The present study aimed to develop and validate a T2WI-based radiomics classifier to differentiate between malignant and benign SPLs. In addition, different machine learning methods were evaluated to achieve the best performance. Furthermore, 2D and 3D features and their combination with clinical features were compared.

Materials and Methods

Data Cohort

This retrospective study was approved by the local ethics committee of the hospital, which waived the need for patients’ informed consent. Preoperative MRI data of 231 patients with chest lesions from November 2015 to April 2018 were analyzed. The inclusion criteria were as follows: (a) lesions were measurable on previous CT scan or T2-weighted imaging; (b) no contraindication for MR examinations; and (c) patients received no therapies or anti-inflammatory therapies at least 2 weeks before the MRI scan, and lesions showed no shrinkage. The exclusion criteria were as follows: (a) operations or biopsies were not available (n = 27); (b) multiple lesions were reported (n = 42); and (c) mediastinal or pleural neoplasms were found (n = 30).

Finally, 132 patients (men and women; age range, 19–78 years; mean age, 54.9 years) were included in the study. The histopathological examination revealed 93 malignant lesions (46 adenocarcinoma, 23 squamous carcinoma, 13 small cell carcinoma, 6 lymphoepithelioid carcinoma, 3 mucoepidermoid carcinoma, 1 sarcomatoid carcinoma, and 1 pulmonary synoviosarcoma). Benign lesions included 7 organizing pneumonia, 6 tuberculosis, 6 granulomatous inflammation, 6 pulmonary cryptococcoses, 4 pulmonary hamartomas, 3 pulmonary aspergilloses, 2 sclerosing alveolar cell tumors, 2 bronchial cysts with infection, 2 focal inflammation, and 1 pulmonary adenofibroma. The patients were randomly divided into two independent groups in a ratio of 7:3. The training cohort included 92 patients (65/27 = positive/negative), whereas the independent test cohort included 40 patients (28/12 = positive/negative).

Image Data Acquisition

All patients were examined with a 3.0-T MRI (Achieva, Philips Healthcare, Best, The Netherlands) using a body phased-array coil. Turbo spin-echo T2-weighted (T2WI) imaging was performed using the following parameters: TR, 992 ms; TE 80 ms; field of view, 350 × 430 mm2; matrix, 640 × 640; thickness, 5 mm; gap, 0.5 mm.

Lesion Segmentation

Mass segmentation was performed to select the entire tumor using open-source software (ITK-SNAP v. 3.6.0, http://www.itksnap.org). Regions of interest (ROIs) of lesions were segmented manually by the consensus of two radiologists with 3 and 8 years of experience (Figure 1). The ROIs included the whole tumor and excluded visible air regions.

FIGURE 1
www.frontiersin.org

Figure 1 Segmentation of a nodule in the upper right lobe on T2-weighted images.

Extraction of Features

The radiomics features were extracted by the Philips Radiomics tool (Philips Medical Systems, Shanghai, China) based on pyRadiomics (11). The hyper-parameters were set to default parameters of the PyRadiomic. The details were described on the website: https://pyradiomics.readthedocs.io/en/latest/features.html. For each ROI, a total of 1,692 3D and 1231 2D radiomic features, including direct features, indirect features, Wavelet transform features, and Laplacian of Gaussian filtered features, were extracted as described in a previous study (12). The 2D features were generated using the slice with the maximum area in 3D ROI. The basic clinical data, sex and age, were included as clinical features. The flow chart for the data processing is displayed in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2 Flow chart for the data processing.

Radiomics Feature Selection and Classifier Building

Building one machine learning model usually consisted of the following steps (1): normalizing each feature to avoid the effect of the scale (2); reducing the dimension of the feature space to remove the information of no use (3); selecting features from the remained features according to the label; and (4) training a classifier to map the selected features onto the diagnosis. In each step, different methods could be selected to make a machine learning pipeline for the final diagnosis.

Three normalization methods [Min-max Normalization (norm-unit), Z-Score Normalization (norm0center), Mean Normalization (norm0centerunit)], two dimension reduction algorithms [principal component analysis (PCA) and Pearson correlation coefficient (PCC)], three feature selection methods [analysis of variance (ANOVA), relief, recursive feature elimination (RFE)], and 10 different classifiers [support vector machine (SVM), auto-encoder (AE), linear discriminant analysis (LDA), random forest (RF), logistic regression (LR), LR-least absolute shrinkage and selection operator (Lasso), Adaboost (AB), decision tree (DT), Gaussian process (GP), naive Bayes (NB)] with 7 different feature numbers (confined to 3–9) were used to build 1260 models for the diagnosis. These methods were chosen due to their popularity in the literature. More details of the combination of the pipelines are illustrated in Figure 3. The number of features was constrained to less than 10% of the training sample size to build a robust machine learning model (13). The ten-fold cross-validation on the training dataset was applied to choose the candidate machine-learning pipeline. The pipeline was selected as the optimal final model when the area under the receiver operating characteristic (ROC) curve (AUC) difference between the training and validation sets on cross-validation was less than 0.1, and the AUC in the test set was the highest. Grid search is used for hyperparameter tuning. The search grid of each classifier was documented in Supplementary Table 1.

FIGURE 3
www.frontiersin.org

Figure 3 Combination of the pipelines for radiomics analysis.

Statistical Analysis

The clinical characteristics between the training and testing sets were compared using the Student t test for continuous variables and the chi-squared (χ2) test for categorical variables. If the counting variable had a theoretical number <10, it was obtained by Fisher’s exact probability test. The aforementioned processes were executed in R software version 3.0.2 (The R Project for Statistical Computing, Vienna, Austria; http://www.r-project.org). A P value of <0.05 indicated statistical significance.

The performance of models was evaluated using the ROC curve, precision-recall(PR) plot, and Matthews Correlation Coefficient. The AUC-ROC and AUC-PR were calculated for quantification. The accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were also calculated at a cutoff value that maximized the value of the Youden index. Also, the estimation was boosted 1000 times to give the 95% confidence intervals (CIs) of AUC-ROC. The aforementioned processes were implemented with FeAture Explorer (14) (FAE, v0.2.5, https://github.com/salan668/FAE) on Python (3.6.8, https://www.python.org/), which is an open-source software based on scikit-learn (v0.22) (15).

Results

Clinical Characteristics

In general, a statistically significant difference was found in age (57.54 ± 10.24 vs. 49.23 ± 14.90, P < 0.001) and sex (male/female ratio 69:24 vs. 20:19, P = 0.01) between the malignant and benign groups. The clinical features of training and test cohorts are summarized in Table 1. No significant differences in lesion diameter were found between the training and test sets. Malignant tumors were more common in upper lobes, men, and elderly people in the training set (P = 0.039, 0.04, and 0.001, respectively) and showed a similar trend in the test set in an insignificant manner (P = 0.667, 0.121, and 0.13, respectively).

TABLE 1
www.frontiersin.org

Table 1 Clinical features of training and test cohorts.

Comparison of Different Machine Learning Classification Models

Figure 4 shows the AUC heat map of 2D and 3D features with different machine learning methods. The number of models with AUC greater than 0.7 in both the validation and test groups was used as an evaluation index. A total of 140 models based on 3D features (n = 129) and 2D features (n = 11) showed AUC greater than 0.7 in both groups. For dimension reduction algorithms, PCA (n = 111) showed higher performance than PCC (n = 29). In terms of normalization, min-max, Z-score, and mean normalization had similar performance in 3D features, while models using only the Z-score (n = 11) showed AUC >0.7 for 2D features. For feature selection, ANOVA (n = 80) performed the best followed by RFE (n = 60), while relief had poor performance (n = 0) in the dataset. As for classifiers, SVM, LDA, LR, GP, and NB performed better for 3D features, while RF, AB, and GP performed better for 2D features (Table 2).

FIGURE 4
www.frontiersin.org

Figure 4 AUC heat map in each dataset showed the performance of 2D and 3D features combined with different machine learning methods in distinguishing solitary pulmonary lesions. It can be clearly seen that the 3D feature group has much more machine learning combinations with higher AUC than 2D feature group in the test dataset. AB, Adaboost; AE, auto-encoder; ANOVA, analysis of variance; DT, decision tree; FN, feature numbers; GP, Gaussian process; LASSO, least absolute shrinkage and selection operator; LDA, linear discriminant analysis; LR, logistic regression; NB, naive Bayes; PCA, principal component analysis; PCC, Pearson correlation coefficient; RF, random forest; RFE, recursive feature elimination; SVM, support vector machine; Unitnorm, Min-max Normalization; Unit with Zerocenter, Mean Normalization; Zscorenorm, Z-score normalization.

TABLE 2
www.frontiersin.org

Table 2 The number of models with AUC greater than 0.7 in both validation and test groups.

Model Performance

The AUC-ROC and AUC-PR of different models are shown in Figures 5, 6. For 2D features, the model based on six features with Z-score + PCA + RFE + GP could achieve stable (the difference in the AUC of the training and validation sets was less than 0.1) and the highest AUC in the test dataset (training: 0.858; validation: 0.810; test: 0.740) with a sensitivity of 0.607 and specificity of 0.833. The model obtained an AUC-PR of 0.846 and MCC of 0.404. After adding clinical features, the AUC-ROC, AUC-PR, and MCC in the test dataset were improved to 0.780, 0.900, and 0.574 (sensitivity: 0.893, specificity: 0.667) with a machine learning pipeline of Z-score + PCA + ANOVA + SVM.

FIGURE 5
www.frontiersin.org

Figure 5 Receiver operating characteristic curves for 2D features, 2D + clinical features, 2D + 3D features, 3D features, and 3D + clinical features in distinguishing malignant from benign solitary pulmonary lesions.

FIGURE 6
www.frontiersin.org

Figure 6 The precision-recall plots of optimal models based on different features.

For 3D features, the model based on nine features with Min-max + PCA + ANOVA + GP could achieve stable and the highest AUC (training: 0.858; validation: 0.797; test: 0.824) with a sensitivity of 0.643 and specificity of 0.917. The model obtained an AUC-PR of 0.927 and MCC of 0.514. Incorporating clinical features with 3D radiomic features slightly improved the AUC-ROC of the testing dataset (AUC = 0.836, sensitivity: 0.821, specificity: 0.833), with an AUC-PR of 0.918 and MCC of 0.620.

For the joint 2D and 3D features, the model based on nine features with Z-score + PCC + ANOVA + SVM could achieve stable and the highest AUC (training: 0.770; validation: 0.737; test: 0.813) with a sensitivity of 0.607 and specificity of 1.000. The model obtained an AUC-PR of 0.926 and MCC of 0.563. The clinical statistics of the independent test dataset are summarized in Table 3. The hyperparameter settings of each classifier used in final models were provided in Supplementary Table 2.

TABLE 3
www.frontiersin.org

Table 3 Clinical statistics in the independent test dataset.

Discussion

The identification of optimal machine learning methods is essential for stable and clinical application (16). The present study provided a comprehensive and detailed assessment of machine learning approaches and explored the diagnostic value of multiple models including 2D, 3D radiomics models, and combinations of clinical and radiomics models based on MR T2WI in noninvasively differentiate SPLs. Our results demonstrated that the T2WI-based radiomics model showed potential in differentiating malignancy from benign SPLs. The 3D radiomics features were better than the 2D features in differentiating SPLs. The optimal machine learning methods were not consistent in different scenarios or with different features included.

Both 2D and 3D features have been employed in previous radiomics researches. A previous CT study suggested that 2D features were superior to 3D features in predicting the prognosis of non-small cell lung cancer (10). However, in some other studies, 3D features demonstrated better predictive performance (17, 18). In the present study, the number of extracted 2D radiomics features was less than that of 3D features because 2D features were extracted based on a single slice, thus losing the spatial information within lesions. Therefore, some features reflecting the spatial distribution of voxels became unavailable. The study found that the radiomics signature derived from 3D features outperformed the signature from 2D features, indicating that the 3D volumetric ROI contained more comprehensive information than 2D ROI and therefore had a better diagnostic performance. Although joint 2D and 3D features showed a higher AUC in a previous study (9), they failed to show superiority over 3D features in the present study. Instead, the performance of joint features was similar to that of 3D features. This finding suggested that the joint features contained information on both 2D and 3D features, while 2D features could not provide new information to 3D features in the present cohort. Accordingly, the classification performance of the joint model failed to show improvement.

In this study, the dimensionality reduction method of PCA was better than PCC, probably because each feature was linearly independent after PCA. Therefore, information could be expressed with fewer features, and thus the performance was stable. However, PCA made variables less interpretable. Therefore, PCA was not suitable for cases where feature “interpretability” was emphasized. In addition, In this study, the three normalization methods (min-max, Z-score, and mean normalization) showed little difference in 3D features, indicating that all these three can be used as effective normalization methods.

Our results showed that the feature selection methods ANOVA and RFE, and the classifier LR, LDA, SVM and GP yield relatively better diagnostic performance for 3D features compared with other methods. This explains the optimal machine learning approach reported in some previous studies (14, 19). Song et al. (14) reported that ANOVA feature selection and an LDA classifier yielded the highest AUC in classifying the clinical-significant prostate cancer (CS PCa) and non-CS Pca. Wang et al. (19) reported that RFE combined with SVM performed the best in distinguishing benign and malignant pulmonary lesions. These also suggested that the optimal machine learning methods were not consistent in various scenarios. Besides, we found that the optimal machine learning strategy was nonunique and different among 2D, 3D, and joint features, indicating that the optimal method might vary depending on the features included.

The lesion diameter did not differ between benign and malignant groups. This is due to the inclusion of inflammatory lesions in the benign group, which could be patchy and thus have a large diameter. Besides, the results showed that elderly patients and men were predisposed to lung malignancies, which was consistent with the findings of a previous study (20). Integrating these clinical data into the radiomics model further increased the accuracy of the model, indicating that clinical and radiomics features contained complementary information needed for differential diagnosis. However, the improvement after adding clinical data was not significant in the present cohort. This might be because the difference in sex and age between malignant and benign groups was not significant in the test dataset. The radiomics models (especially with 3D features) still performed well under such circumstances, which suggested that the radiomics models had the potential to differentiate SPLs in patients with an atypical clinical history (e.g., lung cancer in young patients).

The present study found that some models had pretty high AUC in the training group but relatively lower AUC in the validation or test group (e.g., the classifier RF and AB achieved AUC of 1.0 in the training group). This could be attributed to overfitting and should be avoided. This finding confirmed that it was necessary to have independent datasets to test the real performance of the developed radiomics model. Based on the experience, it was preferred that the difference in AUC between the training and validation or test groups was less than 0.1, indicating that the model was more stable. In the present study, the final models had good discrimination efficiency in the training cohort and showed similar performance in the validation and test groups, suggesting that the models developed in this study were robust.

Previous studies conducted on CT demonstrated that the radiomics model was helpful in distinguishing pulmonary lesions (8, 2123). Chen et al. (23) found that the accuracy of the radiomics signature in benign or malignant classification was 84% with a sensitivity of 92.85% and a specificity of 72.73%. Choi et al. (22) developed a radiomics model with an accuracy of 84.6%, which was 12.4% higher than that of lung-RADS. The performance of the proposed model in the present study was relatively satisfactory and promising, which was similar to that reported in CT radiomics studies. MRI showed great potential in lung application with the advantage of radiation-free and multi-parametric imaging. Recently, pulmonary nodule characterization using MR is recommended for clinical use (24). Nonetheless, the application of MR radiomics in assessing lung diseases is in the initial stage. Therefore, the potential clinical significance of our research is that, for one thing, it provides a non-invasive method that may help to increase the accuracy of MR routine sequence in the differentiation of SPL, for another, it lays a theoretical basis for more clinical applications of MR radiomics in the lung in the future.

This study had several limitations. First, the retrospective study design was subject to potential selection bias. Second, PCC threshold 0.9 has been used in previous studies to screen radiomics features (25, 26). In this study, we used a threshold value lower than 0.9 to more rigorously filter out redundant features. However, a seemingly arbitrarily chosen value of 0.86 was used since in the software we used, PCC can only be set as 0.86 by default. Third, the study included only T2WI to reduce the impact of parametric changes because its scanning parameter remained consistent in all patients while other sequences did not. However, other sequences such as T1W contrast enhancement, diffusion-weighted imaging, and ultrashort TE MRI (27) could provide valuable information and should be included in future studies. Fourth, lesion segmentation was not automatic in the present study and thus could be susceptible to potential human error. Finally, the sample size was relatively small. Therefore, a study with larger sample size and external validation from another institute is needed.

In conclusion, the present study developed and validated a radiomics model based on T2WI that might serve as a promising tool for noninvasive discrimination of SPLs. The 3D features were better than 2D features in differentiating SPLs and performed well in populations with different clinical characteristics. Therefore, 3D segmentations are recommended for further MR radiomics researches. Combining radiomics features with clinical data could further improve model performance. Nonetheless, the optimal machine learning method might not be consistent in different scenarios or with different features involved.

Data Availability Statement

The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by the First Affiliated Hospital of Guangzhou Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

QW and XL conceived the project. QW and JZ completed the article writing. XX, JH, PW, and YP completed data collection and statistical analysis. QW, TZ, and JS completed the Radiomics analysis and chart making. YS and GY provided technique support. All authors contributed to the article and approved the submitted version.

Funding

This study is supported by National Natural Science Foundation of China  (81601457) the Foundation of Guangzhou Municipal Science and Technology Bureau (202102010253), and Open Project Fund of the Sixth Affiliated Hospital of Guangzhou Medical University (2020–11–370).

Conflict of Interest

Authors TZ and JS were employed by Philips Healthcare.

The remaining 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.

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/fonc.2021.683587/full#supplementary-material

References

1. National Lung Screening Trial Research T, Aberle DR, Adams AM, Berg CD, Black WC, Clapp JD, 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

2. Meier-Schroers M, Homsi R, Gieseke J, Schild HH, Thomas D. Lung Cancer Screening With MRI: Evaluation of MRI for Lung Cancer Screening by Comparison of LDCT- and MRI-Derived Lung-RADS Categories in the First Two Screening Rounds. Eur Radiol (2019) 29:898–905. doi: 10.1007/s00330-018-5607-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Sommer G, Tremper J, Koenigkam-Santos M, Delorme S, Becker N, Biederer J, et al. Lung Nodule Detection in a High-Risk Population: Comparison of Magnetic Resonance Imaging and Low-Dose Computed Tomography. Eur J Radiol (2014) 83:600–5. doi: 10.1016/j.ejrad.2013.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Meier-Schroers M, Homsi R, Schild HH, Thomas D. Lung Cancer Screening With MRI: Characterization of Nodules With Different Non-Enhanced MRI Sequences. Acta Radiol (Stockholm Sweden 1987) (2019) 60:168–76. doi: 10.1177/0284185118778870

CrossRef Full Text | Google Scholar

5. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding Tumour Phenotype by Noninvasive Imaging Using a Quantitative Radiomics Approach. Nat Commun (2014) 5:4006. doi: 10.1038/ncomms5006

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Mao L, Chen H, Liang M, Li K, Gao J, Qin P, et al. Quantitative Radiomic Model for Predicting Malignancy of Small Solid Pulmonary Nodules Detected by Low-Dose CT Screening. Quant Imaging Med Surg (2019) 9:263–72. doi: 10.21037/qims.2019.02.02

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Wu H, Sun T, Wang J, Li X, Wang W, Huo D, et al. Combination of Radiological and Gray Level Co-Occurrence Matrix Textural Features Used to Distinguish Solitary Pulmonary Nodules by Computed Tomography. J Digit Imaging (2013) 26:797–802. doi: 10.1007/s10278-012-9547-6

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Jun W, Xia L, Di D, Jiangdian S, Min X, Yali Z, et al. Prediction of Malignant and Benign of Lung Tumor Using a Quantitative Radiomic Method. Conf Proc IEEE Eng Med Biol Soc (2016) 2016:1272–5. doi: 10.1109/EMBC.2016.7590938

CrossRef Full Text | Google Scholar

9. Yang L, Yang J, Zhou X, Huang L, Zhao W, Wang T, et al. Development of a Radiomics Nomogram Based on the 2D and 3D CT Features to Predict the Survival of non-Small Cell Lung Cancer Patients. Eur Radiol (2019) 29:2196–206. doi: 10.1007/s00330-018-5770-y

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Shen C, Liu Z, Guan M, Song J, Lian Y, Wang S, et al. 2D and 3D CT Radiomics Features Prognostic Performance Comparison in Non-Small Cell Lung Cancer. Transl Oncol (2017) 10:886–94. doi: 10.1016/j.tranon.2017.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

11. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res (2017) 77:e104–7. doi: 10.1158/0008-5472.CAN-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Liu X, Yang Q, Zhang C, Sun J, He K, Xie Y, et al. Multiregional-Based Magnetic Resonance Imaging Radiomics Combined With Clinical Data Improves Efficacy in Predicting Lymph Node Metastasis of Rectal Cancer. Front Oncol (2020) 10:585767. doi: 10.3389/fonc.2020.585767

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Park JE, Park SY, Kim HJ, Kim HS. Reproducibility and Generalizability in Radiomics Modeling: Possible Strategies in Radiologic and Statistical Perspectives. Korean J Radiol (2019) 20:1124–37. doi: 10.3348/kjr.2018.0070

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Song Y, Zhang J, Zhang YD, Hou Y, Yan X, Wang Y, et al. FeAture Explorer (FAE): A Tool for Developing and Comparing Radiomics Models. PloS One (2020) 15:e0237587. doi: 10.1371/journal.pone.0237587

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-Learn: Machine Learning in Python. J Mach Learn Res (2011) 12:2825–30.

Google Scholar

16. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: The Bridge Between Medical Imaging and Personalized Medicine. Nature Reviews. Clin Oncol (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141

CrossRef Full Text | Google Scholar

17. Xu L, Yang P, Yen EA, Wan Y, Jiang Y, Cao Z, et al. A Multi-Organ Cancer Study of the Classification Performance Using 2D and 3D Image Features in Radiomics Analysis. Phys Med Biol (2019) 64:215009. doi: 10.1088/1361-6560/ab489f

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ortiz-Ramon R, Larroza A, Ruiz-Espana S, Arana E, Moratal D. Classifying Brain Metastases by Their Primary Site of Origin Using a Radiomics Approach Based on Texture Analysis: A Feasibility Study. Eur Radiol (2018) 28:4514–23. doi: 10.1007/s00330-018-5463-6

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Wang X, Wan Q, Chen H, Li Y, Li X. Classification of Pulmonary Lesion Based on Multiparametric MRI: Utility of Radiomics and Comparison of Machine Learning Methods. Eur Radiol (2020) 30:4595–605. doi: 10.1007/s00330-020-06768-y

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Muller DC, Johansson M, Brennan P. Lung Cancer Risk Prediction Model Incorporating Lung Function: Development and Validation in the UK Biobank Prospective Cohort Study. J Clin Oncol Off J Am Soc Clin Oncol (2017) 35:861–9. doi: 10.1200/JCO.2016.69.2467

CrossRef Full Text | Google Scholar

21. Xiang Y, Sun Y, Liu Y, Han B, Chen Q, Ye X, et al. Development and Validation of a Predictive Model for the Diagnosis of Solid Solitary Pulmonary Nodules Using Data Mining Methods. J Thorac Dis (2019) 11:950–8. doi: 10.21037/jtd.2019.01.90

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Choi W, Oh JH, Riyahi S, Liu CJ, Jiang F, Chen W, et al. Radiomics Analysis of Pulmonary Nodules in Low-Dose CT for Early Detection of Lung Cancer. Med Phys (2018) 45:1537–49. doi: 10.1002/mp.12820

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Chen CH, Chang CK, Tu CY, Liao WC, Wu BR, Chou KT, et al. Radiomic Features Analysis in Computed Tomography Images of Lung Nodule Classification. PloS One (2018) 13:e0192002. doi: 10.1371/journal.pone.0192002

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hatabu H, Ohno Y, Gefter WB, Parraga G, Madore B, Lee KS, et al. Expanding Applications of Pulmonary MRI in the Clinical Evaluation of Lung Disorders: Fleischner Society Position Paper. Radiology (2020) 297:286–301. doi: 10.1148/radiol.2020201138

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhuo Z, Qu L, Zhang P, Duan Y, Cheng D, Xu X, et al. Prediction of H3K27M-Mutant Brainstem Glioma by Amide Proton Transfer–Weighted Imaging and its Derived Radiomics. Eur J Nucl Med Mol Imaging (2021) 48:4426–36. doi: 10.1007/s00259-021-05455-4.

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yu T-T, Lam S-K, To L-H, Tse K-Y, Cheng N-Y, Fan Y-N, et al. Pretreatment Prediction of Adaptive Radiation Therapy Eligibility Using MRI-Based Radiomics for Advanced Nasopharyngeal Carcinoma Patients. Front Oncol (2019) 9:1050. doi: 10.3389/fonc.2019.01050

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ohno Y, Koyama H, Yoshikawa T, Kishida Y, Seki S, Takenaka D, et al. Standard-, Reduced-, and No-Dose Thin-Section Radiologic Examinations: Comparison of Capability for Nodule Detection and Nodule Type Assessment in Patients Suspected of Having Pulmonary Nodules. Radiology (2017) 284:562–73. doi: 10.1148/radiol.2017161037

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: algorithms, area under the curve, lung neoplasms, machine learning, magnetic resonance imaging

Citation: Wan Q, Zhou J, Xia X, Hu J, Wang P, Peng Y, Zhang T, Sun J, Song Y, Yang G and Li X (2021) Diagnostic Performance of 2D and 3D T2WI-Based Radiomics Features With Machine Learning Algorithms to Distinguish Solid Solitary Pulmonary Lesion. Front. Oncol. 11:683587. doi: 10.3389/fonc.2021.683587

Received: 21 March 2021; Accepted: 29 October 2021;
Published: 18 November 2021.

Edited by:

Pavithra Viswanath, University of California, San Francisco, United States

Reviewed by:

Vishwa S. Parekh, Johns Hopkins University, United States
Zhongxiang Ding, Zhejiang University, China

Copyright © 2021 Wan, Zhou, Xia, Hu, Wang, Peng, Zhang, Sun, Song, Yang and Li. 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: Xinchun Li, xinchunli@163.com

These authors have contributed equally to this work

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.