Je souhaite faire un Krigeage simple limité à la France métropolitaine par voisinage unique mais je rencontre des difficultés.
Mes données représentent des entrepôts répartis sur toute la France avec des capacités de stockage.
Voici mon code :
- Code: Select all
#Création de la base de données :
data.db = db.create(data[, c("LB_LONG", "LB_LAT", "CAP_TOTALE")], flag.grid = F, ndim = 2, autoname = F)
- Code: Select all
#Création de la grille -
#je voudrais une grille qui maille la France avec des carrés de 10*10km - je ne sais pas trop comment déclarer ma grille :
grid.db=db.grid.init(data.db,nodes=20)
- Code: Select all
#Calcul du variogramme expérimental
v_exp=vario.calc(data.db, lag = 0.1, nlag = 20)
plot(v_exp)
- Code: Select all
#Ajustement du variogramme expérimental
v_fit=model.auto(v_exp)
- Code: Select all
#Krigeage par Voisinage unique
data.u.neigh = neigh.create(ndim = 2, type = 0)
grid.db = kriging(data.db, grid.db, v_fit, data.u.neigh, uc = "", mean = mean(data.db[, "CAP_TOTALE"], na.rm = T), radix = "KS")
Par ailleurs si je voulais faire un Krigeage prenant en compte non une grille, mais un périmètre autour de mes points, que devrais-je faire ?
J'ai vu dans la doc qu'on pouvait choisir le type de voisinage dans le neigh.create
0 : Unique Neighborhood
1 : Bench Neighborhood
2 : Moving Neighborhood
3 : Image Neighborhood
mais je ne comprends pas l'articulation avec la grille (si on fait un voisinage mouvant, on n'a plus besoin de grille je suppose ??)
Merci de votre aide