An accurate analysis of spatial rainfall distribution is of great importance for managing watershed water resources, in addition to giving support to meteorological studies and agricultural planning. This work compares the performance of two interpolation methods: Inverse distance weighted (IDW) and Kriging, in the analysis of annual rainfall spatial distribution. We use annual rainfall data for the state of Rio Grande do Sul (Brazil) from 1961 to 2017. To determine which proportion of the sample results in more accurate rainfall distribution maps, we use a certain amount of points close to the estimated point. We use mean squared error (MSE), coefficient of determination (R2), root mean squared error (RMSE) and modified Willmott's concordance index (md). We conduct random fields simulations study, and the performance of the geostatistics and classic methods for the exposed case was evaluated in terms of precision and accuracy obtained by Monte Carlo simulation to support the results. The results indicate that the co-ordinary Kriging interpolator showed better goodness of fit, assuming altitude as a covariate. We concluded that the geostatistical method of Kriging using nine closer points (50% of nearest neighbors) was the one that better represented annual rainfall spatial distribution in the state of Rio Grande do Sul.