Skip to Content

variogramme

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


Theme by Dr. Radut.