J'ai des données dont je veux faire une interpolation spatiale par krigeage avec le logiciel R.s'il vous au secours.
Ma variable régionale ici est la variable pop. Maintenant j'ai esseye de faire une grille qui prend les dimensions de ma carte
#grille
x.range<- as.integer(range(base@coords[,1]))
y.range<-as.integer(range(base@coords[,2]))
x=seq(from=x.range[1], to=x.range[2], by=1)
y=seq(from=y.range[1], to=y.range[2], by=1)
base.grid<-expand.grid(x=x, y=y)
plot(base.grid)
coordinates(base.grid) <- ~ x+y
gridded(base.grid)= TRUE
proj4string(base.grid) <-CRS("+init=epsg:28992")
Et quand j’essaie de faire le krigeage en tenant compte des limites de ma cartes,
g<-gstat(id="pop", formula=pop ~ 1, data=base,nmax = 1, set = list(idp = .5))
pop.kriged<- predict(g, model=pop.fit, newdata=base.grid)
pop.kriged
spplot(pop.kriged["pop.pred"],main = "ordinary kriging predictions")
c'est là que je n'arrive pas.svp aide moi c'est pour mon stage