У меня есть CSV-файл с именем seoul3112, содержащий концентрации PM10. пожалуйста, загрузите.. Я попытался построить образец вариограммы и подобрать модель в теме.
library(sp)
library(gstat)
library(rgdal)
seoul3112<-read.csv("seoul3112.csv",row.name=1)
seoul3112<-na.omit(seoul3112)
#назначить CRS и перепроецировать
coordinates(seoul3112)=~LON+LAT
proj4string(seoul3112) = "+proj=longlat +datum=WGS84"
seoul3112<-spTransform(seoul3112, CRS("+proj=utm +north +zone=52 +datum=WGS84"))
#построить полувариограмму
g<-gstat(id="PM10",formula=PM10~LON+LAT, data=seoul3112)
seoul3112.var<-variogram(g, cutoff=70000, width=6000)
seoul3112.var
plot(seoul3112.var, col="black", pch=16,cex=1.3,
xlab="Distance",ylab="Semivariance",
main="Omnidirectional Variogram for seoul 3112")
#Подходит для модели
model.3112<- fit.variogram(seoul3112.var,vgm(700,"Gau",40000,400),fit.method = 2)
plot(seoul3112.var,model=model.3112, col="black", pch=16,cex=1.3,
xlab="Distance",ylab="Semivariance",
main="Omnidirectional Variogram for seoul 3112")
После написания этого кода я получил вот такую полувариограмму.
Поскольку я очень новичок в геостатистике, я не понимаю, подходит ли моя вышеприведенная вариограмма для моего набора данных. Потому что в типичной вариограмме значение семидисперсии становится горизонтальным на пороге. Но эта вариограмма идет вверх! Должен ли я внести некоторые исправления в свой код?
Другое дело, что на самом деле моей конечной целью является проведение интерполяции методом кригинга на моем наборе данных (seoul3112). Я не понимаю, достаточно ли для проведения кригинга этой пробной вариограммы, или мне следует построить дирекционную вариограмму или что-то еще? Может ли кто-нибудь объяснить это подробно?