Multi-classifier-based identification of COVID-19 from chest computed tomography using generalizable and interpretable radiomics features.

PURPOSE: To investigate the efficacy of radiomics in diagnosing patients with coronavirus disease (COVID-19) and other types of viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19.
METHODS: Between 18 January 2020 and 20 May 2020, 110 SARS-CoV-2 positive and 108 SARS-CoV-2 negative patients were retrospectively recruited from three hospitals based on the inclusion criteria. Manual segmentation of pneumonia lesions on CT scans was performed by four radiologists. The latest version of Pyradiomics was used for feature extraction. Four classifiers (linear classifier, k-nearest neighbour, least absolute shrinkage and selection operator [LASSO], and random forest) were used to differentiate SARS-CoV-2 positive and SARS-CoV-2 negative patients. Comparison of the performance of the classifiers and radiologists was evaluated by ROC curve and Kappa score.
RESULTS: We manually segmented 16,053 CT slices, comprising 32,625 pneumonia lesions, from the CT scans of all patients. Using Pyradiomics, 120 radiomic features were extracted from each image. The key radiomic features screened by different classifiers varied and lead to significant differences in classification accuracy. The LASSO achieved the best performance (sensitivity: 72.2%, specificity: 75.1%, and AUC: 0.81) on the external validation dataset and attained excellent agreement (Kappa score: 0.89) with radiologists (average sensitivity: 75.6%, specificity: 78.2%, and AUC: 0.81). All classifiers indicated that "Original_Firstorder_RootMeanSquared" and "Original_Firstorder_Uniformity" were significant features for this task.
CONCLUSIONS: We identified radiomic features that were significantly associated with the classification of COVID-19 pneumonia using multiple classifiers. The quantifiable interpretation of the differences in features between the two groups extends our understanding of CT imaging characteristics of COVID-19 pneumonia.

Introduction

The number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is rapidly increasing. As of 1 December 2020, over 62.1 million cases and 1.4 million deaths were reported globally since the start of the pandemic [1]. The coronavirus strain SARS-CoV-2 causes the 2019 coronavirus disease (COVID-19), resulting in impaired lung function and decreased blood oxygen saturation in affected patients. The main symptoms of COVID-19 include fever, cough, fatigue, and body aches. In addition to the recommended viral nucleic acid real-time reverse transcriptase-polymerase chain reaction (RT-PCR) test, the value of radiological images for diagnosis of COVID-19 has been verified [2,3]. Based on radiologists' evaluations of the signs of pneumonia on chest computed tomography (CT), an average sensitivity of 80 % and specificity of 83 % of clinical diagnosis have been achieved [4].

The distinct characteristics of COVID-19 pneumonia on CT images include bilateral involvement, peripheral distribution, multifocality, mixed ground-glass opacity, consolidation, and vascular thickening [5,6]. However, as radiologists' evaluations of these specific characteristics may be affected by subjective experience, the accuracy of COVID-19 diagnosis varies [4], especially in cases of viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19. Previous studies have confirmed the value of a pre-designed artificial intelligence (AI) framework for task-oriented COVID-19 analysis based on radiological images [7,8]. However, this method is currently hampered by insufficient training data. Recent studies underscore the need for further verification of the sensitivity and stability of current AI-based solutions [9] and for more evidence to establish AI as a production-ready solution for COVID-19 diagnosis [10,11].

Based on the region of interest (ROI) for pneumonia lesions delineated by radiologists, radiomics may provide additional knowledge for survival prognosis and classification of illness severity for COVID-19 pneumonia [[12], [13], [14], [15]]. The latest analytical tool of radiomics, Pyradiomics, has paved the way for standardised radiomics analysis [16]. As an open-source Python package for extracting radiomic features from the ROI on medical images, the standardised radiomic feature extraction program provided by Pyradiomics can reduce programming bias in researchers, improve the reproducibility of radiomic features, and thus enhance the credibility of radiomics results. We hypothesised that radiomics features could be used to differentiate COVID-19 from other types of viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19. We evaluated the performance of radiomics features by using different feature selection and prediction classifiers and investigated standardised high-throughput CT image features extracted by Pyradiomics to clarify the value of radiomics features for the diagnosis of COVID-19.

Materials and methods

Between 18 January 2020 and 20 May 2020, patients with viral pneumonia were initially recruited from two hospitals in China and one hospital in the United States. All procedures performed in studies that involved human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration. Our institutional review board approved this retrospective study (IRB: no.2020-P-036 and no. 51059) and waived the need for written informed consent. There are no conflicts of interest to declare.

Patients

All patients included in this study were admitted during the COVID-19 pandemic. RT-PCR tests for COVID-19 were performed using respiratory secretions obtained by bronchoalveolar lavage, endotracheal aspirate, nasopharyngeal swab, or oropharyngeal swab. These COVID-19 RT-PCR tests were performed at least twice for each patient, to reduce the potential false-negative RT-PCR test.

For SARS-CoV-2 positive patients, the inclusion criteria were as follows: 1) patient over 18 years of age, 2) complete clinical data and CT scan available at admission, and 3) all RT-PCR test results positive. CT scans were reviewed by at least two radiologists to identify the outline of each pneumonia lesion. Patients whose pneumonia lesions could not be visually detected by radiologists were excluded from this study because the ROI of the pneumonia abnormalities could not be delineated for radiomics analysis. For SARS-CoV-2 negative patients, in addition to the points 1) and 2) of the above-mentioned inclusion criteria, 3) all RT-PCR test results negative, 4) patients with viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19, as reviewed by radiologists, and 5) relevant exposure history, was selected as an inclusion criterion for negative patients from China because these patients were at a high risk of COVID-19 infection during the outbreak. The related exposure histories were as follows: H1) history of travel to Wuhan within 14 days before admission, H2) history of contact with a confirmed COVID-19 patient and H3) history of contact with a dense crowd. All SARS-CoV-2 positive patients recruited from China were classified as having mild, common, severe, or critical illness according to the National Health Commission of the People's Republic of China, diagnosis and treatment of COVID-19 pneumonia (Tentative Standard 7).

CT scans

All patients underwent a chest CT scan at admission. The detailed CT scanning parameters of the three hospitals included in this study are listed in Supplementary Appendix A. We tabulated the number of days between symptom onset and date of the CT scan. The time from symptom onset to admission for CT examination was defined as early (0–2 days), intermediate (3–5 days), or late (6–12 days).

Manual segmentation for pneumonia lesions

Four radiologists from our local hospital, each with >5 years of radiological experience, completed the manual segmentation of all pneumonia lesions layer-by-layer using ITK-Snap software (v.3.6.0) [17,18]. Any potential signs of COVID-19 on chest CT were reviewed by radiologists, including bilateral involvement, peripheral distribution, multifocality, mixed ground-glass opacity, consolidation, and vascular thickening, as reported previously [5,6]. Examples of pneumonia lesion segmentation are provided in Supplementary Appendix B. All four radiologists underwent training to identify the characteristics of COVID-19 pneumonia on CT images, to distinguish the CT signs of COVID-19 pneumonia and other types of pneumonia, and to detect the boundary of the pneumonia lesions on CT images. All training sessions were completed online, and approximately eight lessons were provided by radiologists/respiratory experts with experience in the clinical diagnosis of COVID-19 pneumonia. All segmentations were summarised and reviewed by an expert with >10 years of radiological experience. Based on the expert evaluation, the four radiologists held an online consensus meeting to resolve disagreements as to lesion segmentation and to determine the boundaries of these lesions on the CT images.

To verify the credibility of the manual segmentation and feature extraction in this study, ten patients' CT scans were randomly selected and segmented by one radiologist from a different hospital to correct for the potential experience-related bias in feature extraction. Analysis of variance (ANOVA) and Mann-Whitney U tests were performed on the features extracted from the two rounds of segmenting the pneumonia lesion images to verify the inter-observer agreement.

Pyradiomics for feature extraction

The recommended standardised radiomics analysis workflow [19] was used in this study. Pyradiomics was used to extract high-throughput image features of the ROI of pneumonia lesions on CT images. The latest Pyradiomics version available online includes 120-dimensional features, including first-order statistics (19 features), 3D shape-based (16 features), 2D shape-based (10 features), Grey Level Co-occurrence Matrix (24 features), Grey Level Run Length Matrix (16 features), Grey Level Size Zone Matrix (16 features), Neighbouring Grey Tone Difference Matrix (5 features), and Grey Level Dependence Matrix (14 features). A detailed description has been presented previously [16,20].

Radiomics classifier

To verify the performance of radiomics features for the classification of COVID-19 pneumonia, the following classifiers were used in this study: linear classifier, k-nearest neighbour (KNN), least absolute shrinkage and selection operator (LASSO), and random forest (RF).

The pneumonia lesion images from China were randomly divided into training and test datasets with an 80:20 ratio based on the related radiomics studies [[21], [22], [23]]. After obtaining the classification model from the training dataset which executed on the server in China, the model was then provided to the dataset from the United States. The images from the United States were used as an external validation dataset to independently verify the method proposed in this study. The pneumonia lesion images manually delineated by the radiologists will be published for open-access after the review process.

Clinical use

To compare the diagnostic results of radiomics with radiologists, CT images in the two test datasets were additionally diagnosed by three radiologists from another hospital in China with 3, 5, and 10 years of radiological experience, respectively. The radiologists were blinded to the patients' RT-PCR results. We compared the sensitivity and specificity of the diagnosis obtained by radiologists and the radiomics methods. Also, a Kappa consistency analysis was used to evaluate the agreement of the classification result in the classifiers and radiologists on the external validation dataset.

Statistical analysis

The ROC curve, the area under the curve (AUC), sensitivity, and specificity were used to evaluate the diagnostic accuracy for COVID-19. All statistical analyses were performed using the R language (version 3.4.3, Vienna, Austria). The linear classifier, KNN, LASSO, and RF were implemented by the "lm", "kknn", "glmnet", and "randomforest" functions, respectively. Chi-square tests and ANOVA were used to evaluate the differences in demographics between the two groups. For Kappa analysis, the classification results of the two methods were considered to be excellent agreement when the Kappa score ≥0.85, and [0.6, 0.85) indicated good agreement, and [0.45, 0.6) indicated moderate agreement, and Kappa score <0.45 was considered poor agreement. P<0.05 was considered statistically significant.

Results

In total, 266 patients were initially recruited from the three hospitals. Based on the inclusion criteria, 218 patients were enrolled in this study, including 110 patients who were SARS-CoV-2 positive and 108 patients who were SARS-CoV-2 negative, but who had other types of viral pneumonia. There were 184 patients (112 males and 72 females) from China and 34 patients (18 males and 16 females) from the United States, as shown in Table 1 . No significant differences in demographic data were observed between male and female patients (P > 0.05). A flowchart of patient enrolment is shown in Fig. 1 . Regarding the time interval from initial symptom onset to admission CT scan, the ratio of individuals classified as early, intermediate, and late in the SARS-CoV-2 positive group was 21:30:49, and that in the SARS-CoV-2 negative group was 30:20:50, respectively. For all patients the second test was performed one day after the first test. In 26% patients the first RT-PCR test was performed immediately after (within the same day) the suspected pneumonia lesions were found on CT scan. In the other patients, the first RT-PCR test and CT scan were performed at the same time on admission. Statistics regarding the prevalence of the main characteristics of CT presentation in SARS-CoV-2 positive and SARS-CoV-2 negative patients in this study are shown in Table 2 .

Table 1

Demographics of the patients enrolled in this study. Hospital 1 and 2 are in China, and Hospital 3 is in the United States of America.

	Hospital 1		Hospital 2		Hospital 3
	COVID-19	Non-COVID-19	COVID-19	Non-COVID-19	COVID-19	Non-COVID-19
Number	73	71	20	20	17	17
Age	44.1 (15.7)	39.8 (16.5)	39.1 (12.4)	35.5 (12.2)	56.6 (20.3)	56.9 (12.1)
Sex
Male	44	40	14	14	10	8
Female	29	31	6	6	7	9
Time Interval	6.1	4.9	7.4	4.2	10.5	11.3
History
H1	28	17	13	2	NA	NA
H2	11	10	6	5	NA	NA
H3	25	45	11	19	NA	NA

Note. Age is presented with mean and standard deviation. History denotes the related exposure history (H1-H3) defined for included patients in China in this study. NA: not applicable.

Fig. 1

Flowchart of patient enrolment in this study.

Table 2

Statistics of the prevalence of the main characteristics of CT presentation in SARS-CoV-2 positive and SARS-CoV-2 negative patients in this study.

	Main CT characteristics
	GGO	Crazy paving	Consolidation	Combined (2)	Combined (3)
Positive (%)	45.8	65.0	16.9	62.2	14.5
Negative (%)	32.6	6.7	88.8	7.0	2.2

Note. GGO: ground-glass opacity; Combined (2): two main characteristics that significantly appeared in a patient; Combined (3): three main characteristics that significantly appeared in a patient.

For manual segmentation, 12,001 images containing 24,216 pneumonia lesion ROIs were obtained from the Chinese datasets and were randomly divided into the training dataset (9601) and test dataset (2400). A total of 4052 images containing 8409 pneumonia lesion ROIs were manually delineated from the United States dataset. For a case with approximately 20 CT layers containing pneumonia lesions (single or multiple infection foci), the manual segmentation of the ROIs by one radiologist averaged approximately 10 min. Using Pyradiomics, 120 radiomic features were extracted from the original pneumonia lesion images. The statistical results of the ANOVA and Mann-Whitney U test revealed no significant differences between the features extracted from the two rounds of the segmented pneumonia images (P > 0.05).

For the training, test, and external validation datasets, the AUCs (sensitivity and specificity) of the linear classifier were 0.96 (92.7% and 87.8%), 0.96 (95.8% and 85.9%), and 0.79 (69.8% and 75.5%), respectively. The results of LASSO indicated that AUCs of 0.95 (93.1% and 84.1%) and 0.81 (72.2% and 75.1%) were obtained for the test and external validation datasets, respectively. The AUCs for RF were 0.99 (97.5% and 94.9%) and 0.60 (49.5% and 74.8%) for the test and external validation datasets, respectively. AUCs of 0.97 (94.2% and 95.5%) and 0.70 (58.8% and 74.1%) were obtained by the KNN for the test and external validation datasets, respectively. The ROC curves of the four classifiers are presented in Fig. 2 . The results of the radiologists' diagnosis indicated an average sensitivity of 75.6% and specificity of 78.2% for the images in the two test datasets, the detailed comparison is presented in Table 3 . For Kappa analysis, the Kappa score between the radiologists was 0.93, and results indicated that excellent agreement score was obtained between the LASSO and the radiologists (Kappa = 0.89). In contrast, the agreement between the radiologists and linear classifier (Kappa = 0.73), RF (Kappa = 0.40), and KNN (Kappa = 0.55) were lower. A significant difference was found when compared the classification results by the classifiers (P < 0.05).

Fig. 2

The receiver-operating characteristic curves of the four classifiers for the datasets in this study. Area under the curve with sensitivity and specificity are presented for each classifier. A: linear classifier, B: least absolute shrinkage and selection operator, C: k-nearest neighbour, and D: random forest.

Table 3

Comparison of sensitivity (%) and specificity (%) of classification among the four machine learning classifiers and radiologists on different datasets.

	Linear		LASSO		RF		KNN		Radiologists
	Sen (%)	Spe (%)	Sen (%)	Spe (%)	Sen (%)	Spe (%)	Sen (%)	Spe (%)	Sen (%)	Spe (%)
Training	92.7	87.8	91.3	84.4	100	100	NA	NA	NA	NA
Test	95.8	85.9	93.1	84.1	97.5	94.9	94.2	95.5	75.6	78.2
External test	69.8	75.5	72.2	75.1	49.5	74.8	58.8	74.1

NA: not applicable. Sen: sensitivity, Spe: specificity.

Different radiomics features were selected as the key features for classification when using different machine learning classifiers. For the linear classifier, a total of 26 features were identified as statistically significant for the classification task in this study (P < 0.0001). For the LASSO, 17 features were identified as the key features for the classification when the minimum objective function was obtained. For the RF, the top-20 features of importance for the classification were acquired. Of the four classifiers, LASSO had the best classification performance (highest AUC) for external validation dataset. A detailed list of the significant features selected by the three classifiers (ranking of significant features was not available in KNN) and the classifier constructed by the LASSO are shown in Supplementary Appendix C. Also, the source code of model training and the structure of the model is shown in Supplementary Appendix D.

The features of "original_firstorder_RootMeanSquared" (RMS) and "original_firstorder_Uniformity" were significant features validated by all classifiers. Specifically, the values of the RMS feature in COVID-19 lesion images were significantly lower than those in non−COVID-19 lesion images in all datasets (P < 0.05, t-test) in this study (Table 4 ). The values of the "original_firstorder_Uniformity" feature in the COVID-19 lesion images were significantly lower than those in the non−COVID-19 lesion images in the three datasets (P < 0.05, t-test). Evaluation of the differences in this feature reflected in the CT images revealed that the structure of COVID-19 lesion images was more chaotic than that of non−COVID-19 lesion images with high feature values. Specifically, the internal structure of COVID-19 lesions and boundary structure were more likely to demonstrate a sharp contrast, and the internal structure tended to be more heterogeneous than that of non−COVID-19 lesions. Further, the uniformity of COVID-19 pneumonia lesion images with low feature values tended to be worse (Fig. 3 ). Additionally, although the discrimination ability of the "diagnostics_Image-original_Mean" feature was significant in both the training and test datasets, the performance was decreased in the external validation dataset.

Table 4

The mean and standard deviation of the three significant features selected by the classifiers based on the three datasets.

	Mean		RMS		Uniformity
	COVID-19	Non-COVID-19	COVID-19	Non-COVID-19	COVID-19	Non-COVID-19
Training	−707.98 ± 138.98	−994.61 ± 170.22	584.00 ± 100.45	662.61 ± 127.87	0.03 ± 0.01	0.06 ± 0.03
Test	−706.98 ± 135.03	−991.52 ± 170.70	586.28 ± 99.27	662.99 ± 127.02	0.03 ± 0.01	0.06 ± 0.03
External	−513.18 ± 105.32	−561.83 ± 77.57	552.90 ± 90.07	638.77 ± 125.23	0.03 ± 0.01	0.06 ± 0.02

Note: Mean: diagnostics_Image-original_Mean. RMS: original_firstorder_RootMeanSquared. Uniformity: original_firstorder_Uniformity.

Fig. 3

The expression of the differences in COVID-19 (a1-a3) and non−COVID-19 (b1-b3) images by the feature of "original_firstorder_Uniformity". Figure a(1) represents the lung of a 63-year-old male with fever, chest tightness, and anorexia for 9 days. CT manifested as bilateral involvement and multifocality. Figure a(2) represents the lungs of a 88-year-old male with anorexia for 6 days and fever for 2 days. CT manifestation indicates peripheral distribution, diffuseness, and mixed ground-glass opacity. Figure a(3) represents the lungs of a 53-year-old male with weakness and muscle aches for more than 10 days. CT manifestation included bilateral involvement, consolidation, and vascular thickening. Figure b(1) represents the lungs of a 29-year-old male with cough for 6 days. CT manifested as pneumonia lesions similar to the visual characteristics of Figure a(1) in the right lung. Figure b(2) represents the lungs of a 30-year-old male with cough and chest pain for 2 days. CT manifestation included bilateral involvement with multiple consolidation lesions. Figure b(3) represents the lungs of a 37-year-old female with a cough for 3 days. CT manifested as mixed ground-glass opacity, consolidation, and vascular thickening in the left lung.

Discussion

We verified the value of radiomics for differentiating COVID-19 from other types of viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19 based on the standardised radiomics workflow. Using the features extracted by Pyradiomics, our study identified the radiomics features that were significantly associated with the classification of pneumonia patients with and without COVID-19 using multiple classifiers. We also clarified the differences in CT images reflected by different key features.

Image recognition technologies have been demonstrated to be effective methods for the clinical diagnosis of COVID-19 [24,25]. Deep learning-based methods have developed rapidly due to the unnecessary manual segmentation of pneumonia lesions on CT images [26,27]. However, owing to the black-box nature of AI [28,29], neither the operator nor the clinician can intuitively understand the working process of AI for distinguishing patients with COVID-19 pneumonia from those with non−COVID-19 pneumonia. Although high diagnostic accuracy was reported by AI, we require a deeper understanding in clinical practice to perceive the differences in imaging characteristics between COVID-19 and non−COVID-19 pneumonia. Our study demonstrated that the radiomics method can provide feasible solutions. No significant statistical differences in radiomics features obtained by manual segmentation between radiologists ensured the accuracy of the analysis of pneumonia lesions in this study. Based on the features extracted by Pyradiomics (and although the results of screening for distinctive features were different in the four classifiers), there were three measurements on the evaluation of the intensity and homogeneity of pneumonia lesions that were significantly associated with the classification task and were demonstrated cross multiple classifiers. As a measure of image homogeneity, the feature of "original_firstorder_Uniformity" demonstrated that higher feature values were observed in non−COVID-19 pneumonia lesion images. In addition, based on the measurement of the CT grey level magnitude, RMS indicated that the range of grey level of COVID-19 pneumonia was higher. These findings suggested that COVID-19 pneumonia lesions were more likely to be heterogenous, and the alveoli fluid of COVID-19 pneumonia contributed to the more-chaotic structure of infections foci on images. Furthermore, the significantly worse uniformity and the sharper contrast of peripheral and internal structure of COVID-19 pneumonia lesion on CT images indicated that the proposed radiomics features could be used to assist radiologists in improving the performance of visually differentiating CT manifestations of COVID-19 from other types of viral pneumonia. On the other hand, Pyradiomics-based analysis of the pneumonia lesions could further quantify the difference between COVID-19 and other types of viral pneumonia more intuitively, which provide a reference for more accurate clinical diagnosis.

Results of our study indicated that the performance of different classifiers varied substantially in the discrimination of pneumonia lesion features. Other classifiers, such as linear classifier, KNN, and RF achieved good performance in the training and test datasets. However, their performance in the external validation dataset was poorer (P < 0.05) than that of the LASSO. Therefore, studies that use a single machine learning classifier for radiomics feature selection and signature construction may have potential defects. The use of multiple classifiers to evaluate radiomic features in future radiomics studies may play an important role to ensure the credibility of radiomics results. Although the main purpose of this study was not to propose a method beyond the radiologists' diagnosis, we found that the diagnostic performance and Kappa score of the LASSO were closest to that of the radiologists' based on the diagnosis of CT slices. This finding supports the superiority of the LASSO in radiomics workflow, which has been reported in previous studies such as human oncological diseases and immunotherapy-induced pneumonitis [18,[30], [31], [32], [33], [34]].

Our results support the application of radiomics to assist with diagnosis of COVID-19 pneumonia. The recent development of radiomics has provided a new research paradigm in clinical studies [35,36], and there have already been radiomics studies published for COVID-19 survival prognosis and illness severity identification [[12], [13], [14], [15]]. However, current radiomics studies on the differentiation of COVID-19 from other types of viral pneumonia with clinical symptoms and CT signs similar to those of COVID-19, and the evaluation of radiomic feature among different classifiers on COVID-19, are scarce. Our study differed from previous studies [[37], [38], [39]] in several ways. First, we used a variety of machine learning classifiers to prove the effectiveness of radiomics for COVID-19 lesion feature analysis, and the performance of the current mainstream classifiers for COVID-19 classification was determined. Next, the radiomic features that were significant for COVID-19 classification were identified by multi-classifiers. This further strengthened the conclusion that radiomics can be used to assist in the diagnosis of COVID-19 pneumonia in previous studies, which provides a reference for more accurate COVID-19 diagnosis. Finally, we found that different classifiers produced significantly different classification results for the same radiomic features. This finding is of significance for the future radiomics studies to validate multiple machine learning classifiers to improve the credibility of the radiomics results.

In addition to the current specific CT manifestations of COVID-19 observed by radiologists, there is a clinical need for interpretable features that reflect the differential expression of COVID-19 and non−COVID-19 pneumonia on CT images. Radiomic features can be used as bridges to identify differences in pneumonia lesions on CT images in a manner that can be quantitatively interpreted and demonstrated the performance of different classifiers to this task. Based on the proposed significant radiomics features, the experiment could be reproduced using the open access Pyradiomics tool alongside other COVID-19 datasets. Furthermore, use of multiple machine learning classifiers in future radiomics studies will help identify more interpretable and reproducible radiomic features for clinical diagnosis of COVID-19 pneumonia.

This study has several limitations. First, only patients from China and the United States were included. Future studies should employ datasets from other countries to provide a more detailed classification of pneumonia subtypes, such as eosinophilic pneumonia, and other interstitial pneumonia to improve the robustness of the radiomics results. Next, four current mainstream classifiers were used in this study to differentiate COVID-19 and non−COVID-19 pneumonia. Still, the results indicated that the key radiomics features among the classifiers varied, and the results from the external validation dataset indicated that all four classifiers produced overfitting. More approaches should be tested to determine the optimal classifier. Also, even though each patient underwent at least two RT-PCR tests, the possibility of double false-negative RT-PCR tests cannot be eliminated. Finally, this study did not limit CT scanning parameters. Although robustness was improved, the accuracy of the external validation dataset was reduced. Future studies should design an optimised image normalisation method to mitigate the decrease in accuracy caused by different CT scanning parameters.

In conclusion, this study provides new evidence supporting the use of radiomics for the diagnosis of COVID-19 pneumonia. Our study identified the CT imaging features with significant discriminatory potential for patients with COVID-19 pneumonia and non−COVID-19 viral pneumonia using multiple classifiers. We also clarified the differences in CT images reflected by these key features. Our results expand knowledge regarding CT characteristics of COVID-19 pneumonia, which will improve the clinical diagnosis of COVID-19 pneumonia in the current radiological workflow.

Funding

This study received funding from (2018M630310), and (201908210051), and Irish Clinical Academic Training (ICAT) Programme, supported by the Wellcome Trust and the Health Research Board (Grant Number 203930/B/16/Z), the Health Service Executive National Doctors Training and Planning and the Health and Social Care, Research and Development Division, Northern Ireland and the Faculty of Radiologists, Royal College of Surgeons in Ireland. This work was completed as part of a Fulbright-HRB Health Impact award 2020. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Declaration of Competing Interest

None.