Я делаю нормальный кригинг с пакетом gstat. У меня есть координаты x и y и переменная для интерполяции. Координаты x и y представляют собой координаты долготы и широты в проекции 4326. Но когда я создаю экспериментальную вариограмму с этими данными, ось расстояния вариограммы не согласуется с моими данными x, y. Вот пример:
library(sp)
library(gstat)
# coordinates x and y, between 0 and 1.
x = runif(1000)
y = runif(1000)
# variable to interpolate
variable = runif(1000)
# coordinates are put in a data frame.
table = data.frame(x,y)
names(table) = c('coord_x','coord_y')
# dataframe is transformed to a SpatialPoints object
sp::coordinates(table) <- ~ coord_x + coord_y
# I need to change the projection (lattiude longitude coordinates)
sp::proj4string(table) = "+init=epsg:4326"
# creation of the experimental variogram
v1 <- gstat::variogram(variable~1, data=table)
plot(v1)
Как видите, вариограмма показывает расстояние, которое не находится в интервале [0,1]. 
Когда я пробую тот же сценарий без линии проекции, у меня нормальные расстояния:
Но мне нужен прогноз для дальнейшего анализа ...