Today, soil erosion is considered as one of main subjects of watershed management at national and global levels. Assessment of decreasing risk of soil and it spatial distribution are one of key factors for successful evaluation of soil erosion. The objectives of research were calculation and modeling of rain erosivity factor in Ilam province using data of number of 28 pluviometry stations during of 23 years period. After that, rain erosivity factor was calculated, map of rain erosivity factor was provided using geostatistical algorithms. Spatial correlation between data of rain erosivity was drawn using Variogram. In order to, studying of erosivity index were used (IDW), (GPI), (RBF), (LPI) and kriging. The best algorithm was selected using statistical indices of Root Mean Square Error (RMSE) and MAE. Results showed that, GPI was the best algorithm with RMSE and MAE equal to 28.93 and 20.95, respectively.