Soil erodibility (K) refers to the inherent ability of soil to withstand erosion. Accurate estimation and spatial prediction of K values are vital for assessing soil erosion and managing land resources. However, as most K-value estimation models are empirical, they suffer from significant extrapolation uncertainty, and traditional studies on spatial prediction focusing on individual empirical K values have neglected to explore the spatial pattern differences between various empirical models. This work proposed a universal framework for selecting an optimal soil-erodibility map using empirical models enhanced by machine learning. Specifically, three empirical models, namely, the erosion-productivity impact calculator model (K_EPIC), the Shirazi model (K_Shirazi), and the Torri model (K_Torri) were used to estimate K values. Random Forest (RF) and Gradient-Boosting Decision Tree (GBDT) algorithms were employed to develop prediction models, which led to the creation of three K-value maps. The spatial distribution of K values and associated environmental covariates were also investigated across varying empirical models. Results showed that RF achieved the highest accuracy, with R² of K_EPIC, K_Shirazi, and K_Torri increasing by 46%, 34%, and 22%, respectively, compared to GBDT. And distinctions among environmental variables that shape the spatial patterns of empirical models have been identified. The K_EPIC and K_Shirazi are influenced by soil porosity and soil moisture. The K_Torri is more sensitive to soil moisture conditions and terrain location. More importantly, our study has highlighted disparities in the spatial patterns across the three K-value maps. Considering the data distribution, spatial distribution, and measured K values, the K_Torri model outperformed others in estimating soil erodibility in the plateau lake watershed. This study proposed a framework that aimed to create optimal soil-erodibility maps and offered a scientific and accurate K-value estimation method for the assessment of soil erosion.