Rainfall spatial analysis methods are very helpful since there are not enough rainfall gauge stations and watersheds are scattered in large extent. There are many different methods for estimating average precipitation such as; arithmetic method and Thiessen polygon. However, the arrangement and location of data and their correlations are not considered by classic methods. Thus, geostatistical techniques are applied instead. In the present article, 22 meteorological stations from within and around the basin with data collection period of 30 years were selected for the analysis. The geostatistic analysis methods including ordinary kriging, simple cokriging, ordinary cokriging, standardized ordinary kriging, moving average using inverse distance with powers of 1 to 5 were applied for spatial analysis of annual, monthly and 24 hourly maximum rainfall data in Hajighoshan watershed located in northeast of Iran. For this reason, rainfall data were fitted to different methods and compared using cross validation by removes rainfall values of each station, one at a time, and predicts the associated data value. The results of geostatistic analysis showed that ordinary kriging is the best method with MBE=34.26 for annual rainfall while moving average using inverse distance with power of 5 is the best method for monthly and 24 hourly maximum rainfall. According to the results obtained through analysis of variogram model, gaussian model are supposed as the best models for annual, monthly and 24 hourly maximum rainfall data.