The authors describe a methodology of automated making soil erosion class maps using ordinal logistic regression, where relief characteristics are predictors. For mapping, a probabilistic approach is used, in which it is proposed to evaluate the possibility of belonging to several arbitrary erosion classes for each point in space