Interpolation approaches underlie most attempts to create surface maps based on point samples or observations. To the best of my knowledge, Kriging associated with cross-validation is at least one of the most powerful methods. However it is too slow and we have to set a lot of initial values in order to get the best results. I would be thankful if anyone could introduce me a better and faster method (if it exists).