# ajustement à un modèle de variogramme linéaire simple en utilisant les paramètres par défaut
> v.fit <- fit.variogram(variogram(g) ,vgm(model="Lin") )
# adaptation de l'objet gstat en fonction du modèle choisi
> g <- gstat(g, id="elev", model=v.fit )
# interpolation par krigeage ordinaire
> krig_OK <- predict(g, id='elev', newdata=grille)
# renvoi des résultats de l'interpolation dans Grass
> writeRAST6(krig_OK, 'elev.ok', zcol='elev.pred')
# carte de variance d'estimation minimale du krigeage ou variance de krigeage
> writeRAST6(krig_OK, 'elev.ok_var', zcol='elev.var')
> q()
GRASS 6.4.0RC5 (geol):~ >
