Я пытаюсь написать код на R, который использует библиотеку gstat для создания интерполяции. Я уже прочитал руководство по gstat и, основываясь на некоторых примерах в Интернете, мне удалось написать этот код (это только часть):

 g <- gstat(id="tec", formula=TEC ~ 1, data=data)  ##I create an object
 v <- variogram(g) # plot the empirical variogram
 plot(v)
 mod<-vgm(sill=var(data$TEC),model="Sph",range=200,nugget=200) #create the variogram model

v.fit <- fit.variogram(v, model=mod,fit.method=1)  #fit the empirical variogram 
Theor_variogram=plot(variogram(g),v.fit,main="WLS Model") #plot the theoretical variogram
plot(Theor_variogram)
 ## Kriging interpolation
 p <- predict.gstat(g, model=v.fit, newdata=predGrid)

Моя проблема в том, что, когда я запускаю последнюю команду (предикат) вместо получения результата с обычной интерполяцией кригинга, я получаю команду с взвешенным обратным расстоянием (IDW). Я прочитал в руководстве gstat: «Когда вариограммы не указаны, взвешенная обратная интерполяция является действием по умолчанию. Когда вариограммы указаны, метод прогнозирования по умолчанию - обычный кригинг».

Но, как вы можете видеть в моем коде, я указываю как эмпирическую, так и теоретическую вариограмму. Вы знаете, почему я продолжаю получать IDW вместо обычного кригинга? Это может быть связано с типом координат, которые у меня есть? Если, например, у меня есть координаты, близкие друг к другу, или если область интереса слишком велика? Любая помощь будет действительно полезной.

Заранее спасибо димитрис