Predicting glioma grade plays a pivotal role in treatment and prognosis. However, several current methods for grading depend on the characteristics of the whole tumor. Predicting grade by analyzing tumor subregions has not been thoroughly investigated, which aims to improve the prediction performance. To predict glioma grade via analysis of tumor heterogeneity with features extracted from tumor subregions, it is mainly divided into four magnetic resonance imaging (MRI) sequences, including T2-weighted (T2), fluid-attenuated inversion recovery (FLAIR), pre-gadolinium T1-weighted (T1), and post-gadolinium T1-weighted methods. This study included the data of 97 patients with glioblastomas and 42 patients with low-grade gliomas before surgery. Three subregions, including enhanced tumor (ET), non-enhanced tumor, and peritumoral edema, were obtained based on segmentation labels generated by the GLISTRBoost algorithm. One hundred radiomic features were extracted from each subregion. Feature selection was performed using the cross-validated recursive feature elimination with a support vector machine (SVM) algorithm. SVM classifiers with grid search were established to predict glioma grade based on unparametric and multiparametric MRI. The area under the receiver operating characteristic curve (AUC) was used to evaluate the performance of the classifiers, and the performance of the subregions was compared with the results of the whole tumor. In uniparametric analysis, the features from the ET subregion yielded a higher AUC value of 0.8697, 0.8474, and 0.8474 than those of the whole tumor of FLAIR, T1, and T2. In multiparametric analysis, the ET subregion achieved the best performance (AUC = 0.8755), which was higher than the uniparametric results. Radiomic features from the tumor subregion can potentially be used as clinical markers to improve the predictive accuracy of glioma grades.