In order to obtain the location of erosion gullies in Bayan County, understand the spatial relations between multi-features and erosion gullies, assessing the occurrence risk, and provide information and method reference to erosion gully management and precaution, high-resolution satellite imagery, digital elevation model, soil, and precipitation data were used to acquire erosion gully locations and topographical, soil, hydrological, meteorological features. Subsequently, the spatial relations between erosion gullies and multi-feature were analyzed, and erosion gully occurrence risk was assessed using the random forest method. As results, towns like Waxing, Dexiang, and Xinglong could be paid more attention on the erosion gully management. Elevation relief, slope, slope factor, slope length factor, water flow density, and distance to water flow tend to show more spatial relations with erosion gullies, and erosion gully would participate in the redistribution of soil nutrition. On the other hand, for natural or data spatial resolution reasons, slope aspect, terrain curvature, catchment quantity, rainfall erosivity factor, soil bulk density, and soil erodibility factor show little spatial relations with erosion gullies. Assessment indicates that the erosion gully risks are mainly in arable lands with slopes between 1.5° ~ 6°. Results indicates that topographical and hydrological features are closely related to erosion gullies, and diagonal ridge is an economical and effective measure.