géostatistique (krigeage)

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Armel Soubeiga
Messages : 12
Enregistré le : 30 Juin 2016, 09:17

géostatistique (krigeage)

Messagepar Armel Soubeiga » 04 Juil 2016, 17:27

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

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité