Based on ASTER GDEM V2 (30 meter) and the corresponding Google Earth high-resolution images data, small-catchments artificial interaction and random points sample method were used to estimate gully density of five selected catchments in the Losses Plateau and analysis the relationship between the difference for gully density from two methods and the areas of catchments. Besides the gully density, terrain factors such as slope, flow accumulation as well as the distribution probability of gully were calculated by using spatial analysis method. The statistical analysis method was used to analyze the correlation between gully density and terrain factors and obtain gully characteristic distribution prediction model. The results of the study indicate the following: (1) the average difference of gully density extracted by random point method and artificial interactive method is 7.86%, and the difference appears to increase first and then decrease with the watershed area. When the area of catchment is 76 km2, the difference is closed to 0; when the area is 131 km2, the difference is closed to maximum value, which is 0.182. The time-saving and labor-saving random point method can replace the tedious artificial interaction method. (2) The gully slope on Loess Plateau is mainly distributed at 10°∼30° (74%), and the flow accumulation of gully distribution is mainly distributed at 0∼0.07 km2 (95%); 10°∼30° and 0∼0.07 km2 can be used as critical reference values for, respectively. (3) The gully distribution probability showed a when slope gradient was lower than 30 degree, while a negative correlation of gully density and the flow accumulation were found in this study.