Nothing Special   »   [go: up one dir, main page]

Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 13 December 2021
Sec. Head and Neck Cancer

Texture Analysis of Fat-Suppressed T2-Weighted Magnetic Resonance Imaging and Use of Machine Learning to Discriminate Nasal and Paranasal Sinus Small Round Malignant Cell Tumors

Chen ChenChen Chen1Yuhui QinYuhui Qin1Junying ChengJunying Cheng2Fabao Gao*Fabao Gao1*Xiaoyue ZhouXiaoyue Zhou3
  • 1Department of Radiology, West China Hospital, Sichuan University, Chengdu, China
  • 2Department of MRI, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • 3MR Collaboration, Siemens Healthineers Ltd., Shanghai, China

Objective: We used texture analysis and machine learning (ML) to classify small round cell malignant tumors (SRCMTs) and Non-SRCMTs of nasal and paranasal sinus on fat-suppressed T2 weighted imaging (Fs-T2WI).

Materials: Preoperative MRI scans of 164 patients from 1 January 2018 to 1 January 2021 diagnosed with SRCMTs and Non-SRCMTs were included in this study. A total of 271 features were extracted from each regions of interest. Datasets were randomly divided into two sets, including a training set (∼70%) and a test set (∼30%). The Pearson correlation coefficient (PCC) and principal component analysis (PCA) methods were performed to reduce dimensions, and the Analysis of Variance (ANOVA), Kruskal-Wallis (KW), and Recursive Feature Elimination (RFE) and Relief were performed for feature selections. Classifications were performed using 10 ML classifiers. Results were evaluated using a leave one out cross-validation analysis.

Results: We compared the AUC of all pipelines on the validation dataset with FeAture Explorer (FAE) software. The pipeline using a PCC dimension reduction, relief feature selection, and gaussian process (GP) classifier yielded the highest area under the curve (AUC) using 15 features. When the “one-standard error” rule was used, FAE also produced a simpler model with 13 features, including S(5,-5)SumAverg, S(3,0)InvDfMom, Skewness, WavEnHL_s-3, Horzl_GlevNonU, Horzl_RLNonUni, 135dr_GlevNonU, WavEnLL_s-3, Teta4, Teta2, S(5,5)DifVarnc, Perc.01%, and WavEnLH_s-2. The AUCs of the training/validation/test datasets were 1.000/0.965/0.979, and the accuracies, sensitivities, and specificities were 0.890, 0.880, and 0.920, respectively. The best algorithm was GP whose AUCs of the training/validation/test datasets by the two-dimensional reduction methods and four feature selection methods were greater than approximately 0.800. Especially, the AUCs of different datasets were greater than approximately 0.900 using the PCC, RFE/Relief, and GP algorithms.

Conclusions: We demonstrated the feasibility of combining artificial intelligence and the radiomics from Fs-T2WI to differentially diagnose SRCMTs and Non-SRCMTs. This non-invasive approach could be very promising in clinical oncology.

Introduction

Malignant tumors in the nasal and paranasal sinuses are rare, comprise less than 1% of all malignancies and about 3% of head and neck malignancies (1, 2), including small round cell malignant tumors (SRCMTs) and non-SRCMTs. SRCMTs form a specific group of malignancies in the nasal and paranasal sinuses based on neuroectodermal, soft tissue, and hematopoietic differentiation, such as seen in rhabdomyosarcoma (RMS), malignant melanoma (MM), olfactory neuroblastoma (ONB), neuroendocrine carcinoma (NEC), and lymphoma. In contrast, non-SRCMTs form another common group of malignant tumors in the nasal and paranasal sinuses based on epithelial differentiation, including squamous cell carcinomas (SCCs) and adenoid cystic carcinomas (ACCs) (3). The distinction between these two groups is crucial as tumors are variably managed with radiation, chemotherapy, conservative medical therapy, local surgery, exenterative surgery, and multimodal therapy, indicating that therapeutic decisions, surgical planning, and prognoses are very different for each tumor type (4).

Conventional magnetic resonance imaging (MRI) has limitations of its own when differentiating between SRCMTs and Non-SRCMTs. Under the circumstances, as texture analysis (TA) techniques, by using mathematically defined features, can analyze pixel distributions, intensities and dependencies, it can provide a wealth of information beyond what can be seen with the human eye and thus can be used to characterize SRCMTs and Non-SRCMTs, quantitatively (5). Other sequences, such as the apparent diffusion coefficient, have been used to discriminate benign and malignant nasal and paranasal sinus lesions or different histopathologic types of sinonasal malignancies (611). However, less attention has been given to the application of TA for fat-suppressed T2-weighted MR images (Fs-T2WI) collected as part of routine clinical practice.

As a branch of artificial intelligence, machine learning (ML) includes various algorithms that can enhance diagnosis, treatments and follow-up results in neuro-oncology medicine by analyzing huge complex datasets (12, 13). More importantly, not depending on user experience, ML is more objective than other conventional analyses and has good repeatability. To our knowledge, no studies using TA and ML for differentiating sinonasal SRCMTs from non-SRCMTs have been reported. To bridge this gap, this retrospective study was intended to evaluate the potential value of the ML-based Fs-T2WI texture analysis for distinguishing SRCMTs from non-SRCMTs. To achieve the optimal predictive ability and clinical utility, we compared two-dimensional reduction, four feature selection methods and ten ML algorithms.

Materials And Methods

Patients

We used the surgical pathology database from January 1, 2018, to January 1, 2021, in our hospital. Exclusion criteria were (1) patients who received treatments before MRI scans and (2) inadequate image quality. All methods were performed in accordance with the relevant guidelines and regulations, and the informed consent requirement was waived. This retrospective study was approved by the Institutional Ethics Review Committee of our hospital.

Image acquisition

Patients were examined with a 3T MR scanner (MAGNETOM Skyra, Siemens Healthcare, Erlangen, Germany) with a standard head coil. The MRI scan protocols included: axial fat-suppressed T2-weighted imaging (Fs-T2WI) (TR/TE= 5000/117 ms, matrix=256 x 256, field of view=24 x 24 cm2, slice thickness=5 mm, intersection gap =1mm).

Extraction of Textural Features

MaZda software (version 4.7, The Technical University of Lodz, Institute of Electronics, http://www.eletel.p.lodz.pl/mazda/) was used for the analyses. We applied the limitation of dynamics to μ± 3δ(μ: mean grey-level value, δ: standard deviation) (14) to achieve reliable results for the MRI texture classifications. Regions of interest (ROIs) on the Fs-T2WI images of the largest layer were selected. Two physicians delineated ROIs manually along the edge of the lesion and filled the lesion in with a red marker, excluding the various necrotic and cystic regions. In total, 271 features were extracted for each ROI. The number of radiomics features based on feature classes were as follows and shown in Table 1: (i) 9 histogram features based on the number of pixel counts in the image that possessed a certain grey-level value (15); (ii) 220 grey-level co-occurrence matrix (GLCM) features based on the extraction of statistical information about the distribution of pixel pairs (16); (iii) 20 grey-level run-length matrix (GLRLM) features based on searching the image for runs that have the same grey-level values in a pre-defined direction (17); (iv) a 5 auto-regressive model (ARM) based on the weights associated with four neighboring pixels and the variance of the minimized prediction error; (v) 12 wavelet transform (WAV) features on texture frequency components extracted from the energies computed within the channels (18); and (vi) 5 absolute gradient statistics (AGS) features based on smooth or steep variations, resulting in low or high gradient values (15). Multiple GLRLMs were computed along the 0°, 45°, 90°, 135°, and z-axis directions, and 1, 2, 3, and 4 pixels. Multiple GLCMs were computed along four different angles (horizontal, vertical, diagonal 45°, and diagonal 135°).

TABLE 1
www.frontiersin.org

Table 1 Texture analysis methods and the corresponding texture features.

Feature Selections

Computer-generated random datasets were used to assign 70% of datasets to the training set and 30% of the datasets to the independent test set. FeAture Explorer software (FAE, V 0.3.6) software on Python (3.7.6) (https://github.com/salan668/FAE) was used. Firstly, the synthetic minority oversampling technique (SMOTE) was used to balance the training dataset. This method works by taking each minority class sample and introducing synthetic examples along the line segments joining any/all of the k minority class nearest neighbors. The neighboring points were randomly chosen depending on the amount of over-sampling required. Secondly, we normalized the dataset by Z-score Normalization, which subtracts the mean value and divides the standard deviation for each feature. Lastly, we used a Pearson Correlation Coefficient (PCC) and principal component analysis (PCA) to reduce the dimensions. PCC is used for each pair of two features to reduce the row space dimensions of the feature matrix (19). If the PCC was larger than 0.99, one of them was randomly removed. PCA is an unsupervised feature reduction technique that explains the variance-covariance structure of a set of variables through linear combinations. Analysis of Variance (ANOVA) and Kruskal-Wallis (KW) and Recursive Feature Elimination (RFE) and Relief (20) were used for the feature selection. ANOVA was a common analytic method to explore the significant features corresponding to the labels. The KW is a non-parametric version of ANOVA, which hypothesizes that the population median of all groups is equal. The relief selects the sub-data set and finds the relative features according to label recursivity. The goal of the RFE is to select features based on a classifier by recursively considering a smaller set of features. The feature number range was set from 1 to 20.

Classification Performances

The classification performances were tested using 10 ML algorithms, including the support vector machine (SVM), linear discriminant analysis (LDA); auto-encoder (AE); random forests (RF); linear regression (LR); logistic regression using Lasso (LRLasso); ada-boost (AB); decision tree (DT); gaussian process (GP); and naive Bayes (NB)(Table 2).

TABLE 2
www.frontiersin.org

Table 2 The parameters of the algorithms.

Evaluations

The results were evaluated using leave-one-out cross-validation(LOOCV). Using LOOCV, learning sets were created by taking all samples but one, which was used as the validation set. The accuracy, sensitivity, and specificity were also calculated at a cutoff value that maximized the value of the Youden index. The area under the receiver operator characteristics curve (AUC) for the classification of results was calculated for each tested condition.

Results

Of the 171 consecutive patients with a pathologic diagnosis of SRCMTs or Non-SRCMTs over a 2-year period from January 2018 until January 2021, seven were excluded for poor MRI image quality, and 164 patients were finally selected for the study. There were 70 patients with SRCMTs and 94 patients with Non-SRCMTs; RMS (n=16), lymphoma(n=18), MM (n=14), NEC (n=14), ONB (n=8), SCC (n=66), and ACC (n=28). There were 94 males and 70 females in the entire cohort. The mean age of the patients was 55.22 years with a range of 13 to 87 years. After removing invalid cases automatically with FAE, 162 cases were included with 68 SRCMTs and 94 Non-SRCMTs. We assigned 70% of the datasets to the training set (114 patients with 48 SRCMTs and 66 Non-SRCMTs) and 30% of datasets to the independent test set (48 patients with 20 SRCMTs and 28 Non-SRCMTs).

The SMOTE technique was used to automatically create 18 synthetic SRCMTs samples in the training set by operating in the feature space. We compared the AUC of all pipelines on the validation dataset with FAE. The pipeline using PCC dimension reduction, Relief feature selection, and a GP classifier yielded the highest AUC using 15 features. When the “one-standard error” rule was used, FAE also produced a simpler model with 13 features (21), whose ROC curves are shown in Figure 1. The AUCs of the training/validation/test datasets achieved 1.000/0.965/0.979, and the accuracy, sensitivity, and specificity were 0.89, 0.88, and 0.92 (Supplementary Materials), respectively. Features selected by FAE were S(5,-5)SumAverg, S(3,0)InvDfMom, Skewness, WavEnHL_s-3, Horzl_GlevNonU, Horzl_RLNonUni, 135dr_GlevNonU, WavEnLL_s-3, Teta4, Teta2, S(5,5)DifVarnc, Perc.01%, and WavEnLH_s-2 (the weight=1.48, 1.29, 1.28, 1.25, 1.21, 1.19, 1.19, 1.13, 1.10, 1.02, 1.00, 0.99, and 0.98). The AUCs of the training/validation/test datasets by the two-dimensional reduction methods and four feature selection methods were greater than ~0.800 using GP algorithm (Figure 2). Especially, the AUCs of different datasets were more than about 0.900 using the PCC, RFE/Relief, and GP algorithm (Figures 3, 4).

FIGURE 1
www.frontiersin.org

Figure 1 Performance of models generated using Pearson correlation coefficient (PCC) analyses, and Relief and Gaussian process (GP) algorithms. (A) Receiver operating characteristic (ROC) curves of this model on different datasets, (B) FeAture Explorer software suggested a candidate 13-feature model according to the “one-standard error” rule, and (C) A contribution of features in the final model.

FIGURE 2
www.frontiersin.org

Figure 2 Areas under the curve (AUCs) of the different datasets using Pearson correlation coefficient (PCC) and principal component analysis (PCA) methods and Relief and Gaussian process (GP) algorithms. (A) Analysis of Variance (ANOVA), (B) Kruskal-Wallis (KW), (C) Relief, and (D) Recursive Feature Elimination (RFE).

FIGURE 3
www.frontiersin.org

Figure 3 Areas under the curve (AUCs) on different datasets using the Analysis of Variance (ANOVA), Kruskal-Wallis (KW), and Recursive Feature Elimination (RFE) and Relief using Gaussian process (GP). (A) Pearson correlation coefficient (PCC), (B) principal component analysis (PCA).

FIGURE 4
www.frontiersin.org

Figure 4 Areas under the curve (AUCs) of different datasets using 10 machine learning algorithms. (A) Pearson correlation coefficient (PCC), (B) principal component analysis (PCA). SVM, support vector machine; LDA, linear discriminant analysis; AE, auto-encoder; RF, random forests; LR, linear regression; LRLasso, logistic regression via Lasso; AB, ada-boost; DT, decision tree; GP, Gaussian process; NB, naive Bayes.

Discussion

This study investigated the potential value of the Fs-T2WI texture analysis of maximum tumor solid components for distinguishing SRCMTs from non-SRCMTs with ML. The key findings were as follows: (1) The pipeline using PCC dimension reduction, relief feature selection, and GP classifier yielded the highest AUC. (2) The best algorithm was GP whose AUCs of the training/validation/test datasets by the two-dimensional reduction methods and four feature selection methods were greater than approximately 0.800. Especially, the AUCs of different datasets were more than about 0.900 using the PCC, RFE/Relief, and GP algorithm. (3) TA with ML appears to be most helpful in tumor differentiation using standard Fs-T2WI routinely acquired with a high accuracy of 0.89.

Radiomics data contain first-, second-, and higher-order statistics (22). First-order statistics are described as the distribution of individual voxel values regardless of spatial relationships and are generally histogram-based features that refer to statistical parameters of pixel intensities within the ROI, such as mean, variance, skewness, kurtosis, and Perc.10%. Second-order statistics describe texture features. Specifically, they describe statistical interrelationships between voxels with similar or dissimilar contrast values, such as the co-occurrence matrix, which is calculated from the intensities of pixel pairing, with the spatial relationship of pixel pairing defined. Higher-order statistics impose filter grids on images to extract repetitive or nonrepetitive patterns, such as wavelets, which are data on texture frequency components extracted from the energies computed within the channels. In 2015, Fujima et al. (23) assessed the utility of the histogram analysis on tumor blood flow (TBF) obtained with pseudo-continuous arterial spin labeling to differentiate SCC and lymphoma in nasal or sinonasal cavities, achieving an accuracy of 0.87 for the mean TBF, the coefficient of variation, and kurtosis. In 2019, this group (24) also applied histograms and TAs for Fs-T2WI to differentiate SCC and lymphoma in the head and neck and found that the relative mean signal, contrast, and homogeneity could be useful. Another study (25) applied GLCM and WAV, using the ANOVA, rank-sum test, and SVM classification feature selections to diagnose osteosarcoma and achieved an accuracy of 0.82-0.96. Muramatsu et al. (26) applied GLCM and WAV using artificial neural network, SVM, and RF algorithms to classify and diagnose malignant and benign breast masses, achieving AUCs of 0.83-0.86. Finally, 440 radiomics features were applied, including histograms, GLCMs, GLRLMs, and WAVs, using 24 feature selection methods and three classification methods to predict lung cancer histologic subtypes (27)and they found that the Relief and NB algorithms achieved higher accuracies and the highest AUCs compared with other studied algorithms. Similarly, in our study, the Relief feature selection and GP classifier yielded the highest AUCs using the histogram, GLCM, GLRLM, and WAV. The Relief can avoid a heuristic search and select relevant features in linear time based on the given features and training instances regardless of the complexity of the target concept to be learned. The GP classification is a nonparametric method based on a Laplace approximation used for approximating the non-Gaussian posterior by a Gaussian method. It can easily handle a variety of problems, such as an insufficient capacity for the classical linear method, complex data types, and the curse of dimensions (28). Sovizi et al. reported that they achieved high sensitivity results using GP classification model (29). Consistent with their results, the AUCs of different datasets by GP algorithm were greater than approximately 0.900 with the RFE/Relief feature selection methods. RFE is recursively repeated in the pruned set, removing the least important features until the desired number of feature selections is eventually reached. Chatterjee et al. (30) have applied TAs to differentiate melanoma, dysplastic nevi, and basal cell carcinoma in dermoscopic images using SVM-RFE, achieving an accuracy of 0.952. Vamvakas et al. (31) applied SVM-RFE using 3D TA to classify low- and high-grade gliomas. With a LOOCV, they achieved an accuracy of 0.955 and an AUC of 0.955.

Some studies reported their diverse automated or semi-automated segmentation methods. For example, Muzzamil et al. proposed lung CT image segmentation using intensity thresholding (32). Husham et al. (33) and Malathi et al. (34), for another example, analysed between active contour and Otsu thresholding segmentation algorithms in segmenting brain tumor MRI and pleura diseases CT, respectively. In addition, Hussein et al. (35) proposed a new Viola–Jones model for the segmentation of ovarian and breast ultrasound images. Artificial neural networks and SVM have been tried for division of nasopharyngeal carcinoma, respectively (36, 37). The accuracy and consistency of the tumor delineation plays an important role in differential diagnosis. However, most of the tumors in the site of the nasal cavity and paranasal sinus adjacent to the air in Fs-T2 maps are without edema areas. Thus, we chose manual segmentation by experts not automated or semi-automatic segmentation to determine the boundary.

There were some limitations. First, as the SRCMTs studied were of various histologic types, subgroup analyses in more details should be performed in future studies after obtaining a larger sample size and a careful consideration of the study groups. Second, our model used manually delineated ROIs performed along the edge of the lesion. Segmenting precise tumor regions is the focus of future work. In our further studies, we will propose a multiparametric MRI investigation including ADC, T2-weighted MRI and dynamic contrast-enhanced MRI involving early and delayed phases to generate a robust model to differentially diagnose SRCMTs and Non-SRCMTs by segmenting precisely three-dimensional tumor regions in a larger sample.

Conclusions

We demonstrated the feasibility of combining artificial intelligence and radiomics using Fs-T2WI in the differential diagnosis of SRCMTs and Non-SRCMTs. This approach could be a very promising non-invasive method in clinical oncology.

Data Availability Statement

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

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Ethics Review Committee of West China Hospital. Written informed consent from the participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements. Written informed consent was not obtained from the individual(s), nor the minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.

Author Contributions

CC, YQ and JC initiated this study, participated in its design, performed study selection, data extraction, and data analysis. CC drafted the manuscript. FG supervised all aspects of the study. XZ revised the language. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (no. 81930046, 81771800, and 81829003).

Conflict of Interest

Author XZ was employed by Siemens Healthineers Ltd., Shanghai, China.

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

References

1. Weber AL, Stanton AC. Malignant Tumors of the Paranasal Sinuses: Radiologic, Clinical, and Histopathologic Evaluation of 200 Cases. Head Neck Surg (1984) 6(3):761–76. doi: 10.1002/hed.2890060310

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Silverberg E, Grant RN. Cancer Statistics, 1970. CA Cancer J Clin (1970) 20(1):11–23. doi: 10.3322/canjclin.20.1.10

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Simons SA, Bridge JA, Leon ME. Sinonasal Small Round Blue Cell Tumors: An Approach to Diagnosis. Semin Diagn Pathol (2016) 33(2):91–103. doi: 10.1053/j.semdp.2015.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Thompson LD. Small Round Blue Cell Tumors of the Sinonasal Tract: A Differential Diagnosis Approach. Mod Pathol (2017) 30(s1):S1–S26. doi: 10.1038/modpathol.2016.119

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Fetit AE, Novak J, Peet AC, Arvanitits TN. Three-Dimensional Textural Features of Conventional MRI Improve Diagnostic Classification of Childhood Brain Tumours. NMR BioMed (2015) 28(9):1174–84. doi: 10.1002/nbm.3353

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ozturk K, Gencturk M, Caicedo-Granados E, Li F, Cayci Z. Positron Emission Computed Tomography and Magnetic Resonance Imaging Features of Sinonasal Small Round Blue Cell Tumors. Neuroradiol J (2020) 33(1):48–56. doi: 10.1177/1971400919873895

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gencturk M, Ozturk K, Caicedo-Granados E, Li F, Cayci Z. Application of Diffusion-Weighted MR Imaging With ADC Measurement for Distinguishing Between the Histopathological Types of Sinonasal Neoplasms. Clin Imaging (2019) 55:76–82. doi: 10.1016/j.clinimag.2019.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang X, Zhang Z, Chen Q, Li J, Xian J. Effectiveness of 3 T PROPELLER DUO Diffusion-Weighted MRI in Differentiating Sinonasal Lymphomas and Carcinomas. Clin Radiol (2014) 69(11):1149–56. doi: 10.1016/j.crad.2014.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sasaki M, Sumi M, Eida S, Ichikawa Y, Sumi T, Yamada T, et al. Multiparametric MR Imaging of Sinonasal Diseases: Time-Signal Intensity Curve- and Apparent Diffusion Coefficient-Based Differentiation Between Benign and Malignant Lesions. AJNR Am J Neuroradiol (2011) 32(11):2154–9. doi: 10.3174/ajnr.A2675

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Razek AA, Sieza S, Maha B. Assessment of Nasal and Paranasal Sinus Masses by Diffusion-Weighted MR Imaging. J Neuroradiol (2009) 36(4):206–11. doi: 10.1016/j.neurad.2009.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

11. White ML, Zhang Y, Robinson RA. Evaluating Tumors and Tumorlike Lesions of the Nasal Cavity, the Paranasal Sinuses, and the Adjacent Skull Base With Diffusion-Weighted MRI. J Comput Assist Tomogr (2006) 30(3):490–5. doi: 10.1097/00004728-200605000-00023

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Villanueva-Meyer JE, Chang P, Lupo JM, Hess CP, Flanders AE, Kohli M, et al. Machine Learning in Neurooncology Imaging: From Study Request to Diagnosis and Treatment. AJR Am J Roentgenol (2019) 212(1):52–6. doi: 10.2214/AJR.18.20328

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Senders JT, Zaki MM, Karhade AV, Chang B, Gormley WB, Broekman ML, et al. An Introduction and Overview of Machine Learning in Neurosurgical Care. Acta Neurochir (Wien) (2018) 160(1):29–38. doi: 10.1007/s00701-017-3385-8

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Collewet G, Strzelecki M, Mariette F. Influence of MRI Acquisition Protocols and Image Intensity Normalization Methods On Texture Classification. Magn Reson Imaging (2004) 22(1):81–91. doi: 10.1016/j.mri.2003.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Castellano G, Bonilha L, Li LM, Cendes F. Texture Analysis of Medical Images. Clin Radiol (2004) 59(12):1061–9. doi: 10.1016/j.crad.2004.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Haralick R, Shanmugam K, Dinstein I. Textural Features for Image Classification. IEEE Trans Syst Man Cybern (1973) SMC-3:610–21. doi: 10.1109/TSMC.1973.4309314

CrossRef Full Text | Google Scholar

17. Galloway M. Texture Analysis Using Gray Level Run Lengths. Comput Graphics Image Process (1975) 4(2):172–9. doi: 10.1016/S0146-664X(75)80008-6

CrossRef Full Text | Google Scholar

18. Orphanidou-Vlachou E, Vlachos N, Davies NP, Arvanitis TN, Grundy RG, Peet AC, et al. Texture Analysis of T1 - and T2 -Weighted MR Images and Use of Probabilistic Neural Network to Discriminate Posterior Fossa Tumours in Children. NMR BioMed (2014) 27(6):632–9. doi: 10.1002/nbm.3099

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Benesty J, Chen J, Huang Y, et al. Pearson Correlation Coefficient. In: Cohen I, Huang Y, Chen J, et al, editors. Noise Reduction in Speech Processing. Berlin, Heidelberg: Springer Berlin Heidelberg (2009). p. 1–4.

Google Scholar

20. Kira K, Rendell LA. A Practical Approach to Feature Selection. In: Sleeman D, Edwards P, editors. Machine Learning Proceedings 1992. San Francisco (CA: Morgan Kaufmann (1992). p. 249–56.

Google Scholar

21. 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(8):e237587. doi: 10.1371/journal.pone.0237587

CrossRef Full Text | Google Scholar

22. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More Than Pictures, They Are Data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Fujima N, Kameda H, Tsukahara A, Yoshida D, Sakashita T, Homma A, et al. Diagnostic Value of Tumor Blood Flow and its Histogram Analysis Obtained With pCASL to Differentiate Sinonasal Malignant Lymphoma From Squamous Cell Carcinoma. Eur J Radiol (2015) 84(11):2187–93. doi: 10.1016/j.ejrad.2015.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Fujima N, Homma A, Harada T, Shimizu Y, Tha KK, Kano S, et al. The Utility of MRI Histogram and Texture Analysis for the Prediction of Histological Diagnosis in Head and Neck Malignancies. Cancer Imaging (2019) 19(1):5. doi: 10.1186/s40644-019-0193-9

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hu S, Xu C, Guan W, Tang Y, Liu Y. Texture Feature Extraction Based on Wavelet Transform and Gray-Level Co-Occurrence Matrices Applied to Osteosarcoma Diagnosis. BioMed Mater Eng (2014) 24(1):129–43. doi: 10.3233/BME-130793

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Muramatsu C, Hara T, Endo T, Fujita H. Breast Mass Classification on Mammograms Using Radial Local Ternary Patterns. Comput Biol Med (2016) 72:43–53. doi: 10.1016/j.compbiomed.2016.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wu W, Parmar C, Grossmann P, Quackenbush P, Lambin P, Bussink J, et al. Exploratory Study to Identify Radiomics Classifiers for Lung Cancer Histology. Front Oncol (2016) 6:71. doi: 10.3389/fonc.2016.00071

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Maniruzzaman M, Jahanur RM, Ahammed B, Abedin MM, Suri HS, Biswas M, et al. Statistical Characterization and Classification of Colon Microarray Gene Expression Data Using Multiple Machine Learning Paradigms. Comput Methods Programs BioMed (2019) 176:173–93. doi: 10.1016/j.cmpb.2019.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sovizi J, Mathieu KB, Thrower SL, Stefan W, Hazle JD, Fuentes D, et al. Gaussian Process Classification of Superparamagnetic Relaxometry Data: Phantom Study. Artif Intell Med (2017) 82:47–59. doi: 10.1016/j.artmed.2017.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chatterjee S, Dey D, Munshi S. Integration of Morphological Preprocessing and Fractal Based Feature Extraction With Recursive Feature Elimination for Skin Lesion Types Classification. Comput Methods Programs BioMed (2019) 178:201–18. doi: 10.1016/j.cmpb.2019.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Vamvakas A, Williams SC, Theodorou K, Kapsalaki E, Fountas K, Kappas C, et al. Imaging Biomarker Analysis of Advanced Multiparametric MRI for Glioma Grading. Phys Med (2019) 60:188–98. doi: 10.1016/j.ejmp.2019.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Javaid M, Javid M, Rehman MZ, Shah SI. A Novel Approach to CAD System for the Detection of Lung Nodules in CT Images. Comput Methods Programs BioMed (2016) 135:125–39. doi: 10.1016/j.cmpb.2016.07.031

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Husham S, Mustapha A, Mostafa S, Mohammed M, George ST. Comparative Analysis Between Active Contour and Otsu Thresholding Segmentation Algorithms in Segmenting Brain Tumor Magnetic Resonance Imaging. J Informat Technol Manag (2020). doi: 10.22059/jitm.2020.78889

CrossRef Full Text | Google Scholar

34. Malathi M, Sinthia P, Jalaldeen K. Active Contour Based Segmentation and Classification for Pleura Diseases Based on Otsu's Thresholding and Support Vector Machine (SVM). Asian Pac J Cancer Prev (2019) 20(1):167–73. doi: 10.31557/APJCP.2019.20.1.167

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hussein I, Burhanuddin MA, Mohammed M, Elhoseny M, Garcia-Zapirain B, Maashi MS, et al. Fully Automatic Segmentation of Gynaecological Abnormality Using a New Viola–Jones Model. Computers Materials \& Continua (2021) 66(3):3161–82. doi: 10.32604/cmc.2021.012691

CrossRef Full Text | Google Scholar

36. Abd Ghani MK, Mohammed MA, Arunkumar N, Mostafa SA, Ibrahim DA, Abdullah MK, et al. Decision-Level Fusion Scheme for Nasopharyngeal Carcinoma Identification Using Machine Learning Techniques. Neural Computing Appl (2020) 32(3):625–38. doi: 10.1007/s00521-018-3882-6

CrossRef Full Text | Google Scholar

37. Mohammed MA, Abd Ghani MK, Arunkumar N, Mostafa SA, Abdullah MK, Burhanuddin MA. Trainable Model for Segmenting and Identifying Nasopharyngeal Carcinoma. Comput Electrical Eng (2018) 71:372–87. doi: 10.1016/j.compeleceng.2018.07.044

CrossRef Full Text | Google Scholar

Keywords: Fs-T2WI, artificial intelligence, machine learning, texture analysis, radiomics, small round cell malignant tumors

Citation: Chen C, Qin Y, Cheng J, Gao F and Zhou X (2021) Texture Analysis of Fat-Suppressed T2-Weighted Magnetic Resonance Imaging and Use of Machine Learning to Discriminate Nasal and Paranasal Sinus Small Round Malignant Cell Tumors. Front. Oncol. 11:701289. doi: 10.3389/fonc.2021.701289

Received: 27 April 2021; Accepted: 18 November 2021;
Published: 13 December 2021.

Edited by:

Christophe Le Tourneau, Institut Curie, France

Reviewed by:

Mazin Mohammed, University of Anbar, Iraq
Yingwei Qiu, Third Affiliated Hospital of Guangzhou Medical University, China

Copyright © 2021 Chen, Qin, Cheng, Gao and Zhou. 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: Fabao Gao, gaofabao@wchscu.cn

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.