Myelodysplastic syndromes (MDSs) are a group of hematologic neoplasms accompanied by dysplasia of bone marrow (BM) hematopoietic cells with cytopenia. Recently, digitalized pathology and pathomics using computerized feature analysis have been actively researched for classifying and predicting prognosis in various tumors of hematopoietic tissues. This study analyzed the pathomic features of hematopoietic cells in BM aspiration smears of patients with MDS according to each hematopoietic cell lineage and dysplasia. We included 24 patients with an MDS and 21 with normal BM. The 12,360 hematopoietic cells utilized were to be classified into seven types: normal erythrocytes, normal granulocytes, normal megakaryocytes, dysplastic erythrocytes, dysplastic granulocytes, dysplastic megakaryocytes, and others. Four hundred seventy-six pathomic features quantifying cell intensity, shape, and texture were extracted from each segmented cell. After comparing the combination of feature selection and machine learning classifier methods using 5-fold cross-validation area under the receiver operating characteristic curve (AUROC), the quadratic discriminant analysis (QDA) with gradient boosting decision tree (AUROC = 0.63) and QDA with eXtreme gradient boosting (XGB) (AUROC = 0.64) showed a high AUROC combination. Through a feature selection process, 30 characteristics were further analyzed. Dysplastic erythrocytes and granulocytes showed lower median values on heatmap analysis compared to that of normal erythrocytes and granulocytes. The data suggest that pathomic features could be applied to cell differentiation in hematologic malignancies. It could be used as a new biomarker with an auxiliary role for more accurate diagnosis. Further studies including prediction survival and prognosis with larger cohort of patients are needed.