Как сделать Обычный Кригинг с помощью gstat
Я пытаюсь написать код на 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 вместо обычного кригинга? Это может быть связано с типом координат, которые у меня есть? Если, например, у меня есть координаты, близкие друг к другу, или если область интереса слишком велика? Любая помощь будет действительно полезной.
Заранее спасибо димитрис