variogramme

# chargement de la librairie gstat
> library(gstat)
# lecture des points à traiter:
> mespoints = readVECT6('newcondZ@moi')
> coor <-coordinates(mespoints)
> coor
       coords.x   coords.y   coords.z 
  [1,]  169099.0  113900.0    239.00
  [2,]  169098.1  113895.8    214.38
  [3,]  169097.2  113891.5    189.76
  [...]
> z <- coor[,3]
> z 
 [1] 239.00 214.38 189.76 [...]
 
# l'interpolation nécessite une grille qui sera remplie par les valeurs interpolées
> grd <- gmeta2grd()
> grille <- SpatialGridDataFrame(grd, data=data.frame(k=rep(1,G$cols*G$rows)), proj4string=CRS(mespoints@proj4string@projargs))
 
# initialisation d'un objet gstat avec la valeur z à traiter 
> g <- gstat(id="elev", formula=z ~ 1, data=mespoints)
 
# vue des divers variogrammes (en fonctions de divers paramètres)
> plot(variogram(g,[...])