Toujours dans l'optique de me former à R, j'essaie de faire un krigeage.
J'ai récurépé des données de prélévement de sol (environ 1000). Avec seulement 3 informations: lattitude, Longitude et une quantité de polluant.
Et j'utilise les packages gstat et sp.
Pour pouvoir faire mon krigeage, je dois faire une grille. et c'est mon problème.
Je n'arrive pas à créer la grille en reprénant les informations de mes données.
Voila ce que j'ai écris:
Code : Tout sélectionner
seqx <- seq(min(spatial$long), max(spatial$long), by=0.05)
seqy <- seq(min(spatial$lat), max(spatial$lat), by=0.05)
Jusque là tout va bien.
puis:
Code : Tout sélectionner
Mseqx <- rep(seqx, length(seqy))
Mseqy <- rep(seqy, length(seqx))
Mais ici j'ai un message d'erreur me disant que la taille du vecteur est trop importante.
Je ne peux donc pas continuer.
Existe-il une solution? comment réduire la taille sans impacter ?
Voila ce que je voulais écrire ensuite
Code : Tout sélectionner
Mseqy <- sort(Mseqy, decreasing=F)
grid <- data.frame(x=Mseqx, y=Mseqy)
coordinates(grid)=~x+y
gridded(grid) <- TRUE
Dites moi si vous avez une idée.
Merci d'avance.
(j'ai hate d'être le jour, où je pourrais enfin aider quelqu'un sur ce forum).