Purpose. To assess the role of multiple radiomic features of lymph nodes in the preoperative prediction of lymph node metastasis (LNM) in patients with esophageal squamous cell carcinoma (ESCC). Methods. Three hundred eight patients with pathologically confirmed ESCC were retrospectively enrolled (training cohort, n = 216; test cohort, n = 92). We extracted 207 handcrafted radiomic features and 1000 deep radiomic features of lymph nodes from their computed tomography (CT) images. The t-test and least absolute shrinkage and selection operator (LASSO) were used to reduce the dimensions and select key features. Handcrafted radiomics, deep radiomics, and clinical features were combined to construct models. Models I (handcrafted radiomic features), II (Model I plus deep radiomic features), and III (Model II plus clinical features) were built using three machine learning methods: support vector machine (SVM), adaptive boosting (AdaBoost), and random forest (RF). The best model was compared with the results of two radiologists, and its performance was evaluated in terms of sensitivity, specificity, accuracy, area under the curve (AUC), and receiver operating characteristic (ROC) curve analysis. Results. No significant differences were observed between cohorts. Ten handcrafted and 12 deep radiomic features were selected from the extracted features ( p < 0.05 ). Model III could discriminate between patients with and without LNM better than the diagnostic results of the two radiologists. Conclusion. The combination of handcrafted radiomic features, deep radiomic features, and clinical features could be used clinically to assess lymph node status in patients with ESCC.