Abstract: Kriging interpolation is a common method applied to spatial data. Correct fitting of the variogram mode is the key to kriging interpolation. At present, experience is the guide when fitting variograms in analysis tools and development kits. An intuitional method is used to determine the parameters of a variogram and the parameters are adjusted reiteratively. These artificial operational steps obstruct the automation of kriging interpolation but also are not based on scientific theories relevant to setting parameters and option modes As a result, the predictive results are not exact. Aiming to address this problem and with the goal of automating the process, in this paper we calculate the sampling variogram and an automatic fitting algorithm for variograms is studied. Automatic kriging interpolation is realized and automatic cross validation is conducted based on the automatic fitting function in R language. The results of experiments show that a automatic method for kriging interpolation to spatial data is feasible.