Suite

Ajouter un point sur un spplot dans R

Ajouter un point sur un spplot dans R


Je dois ajouter un point sur un existantspplotdans R. J'en ai besoin uniquement à des fins visuelles. J'ai les coordonnées du point, mais je ne sais pas comment changer le code pour ajouter le point.

C'est le code duspplot:

spplot(pioggiaTP[1])

Et voici la carte obtenue en exécutant la commande :

Les coordonnées du point sont X=720731,8704 et Y=4861300.795 dans le système de référence WGS 84 / UTM zone 32N EPSG:32632, le même système de la parcelle.

Est-ce que quelqu'un peut m'aider?


spplotrenvoie un tracé en treillis de classe "treillis". Vous pouvez utiliser lecouchefonction de latreillisExtrapackage pour mettre à jour cet objet.

Exemple:

library("sp") library("latticeExtra") data(meuse) data(meuse.grid) coordonnées(meuse.grid) <- ~x+y proj4string(meuse.grid) <- CRS("+init=epsg : 28992") p <- spplot(meuse.grid, "dist") print(p)

p + layer(panel.points(x, y, col="green", pch=19), data=meuse)


Une autre option serait d'utiliser l'argument 'sp.layout' fourni avecspplot. Notez que contrairement à l'approche suggérée par @rcs, vous devez fournir un objet 'sp' (que vous avez peut-être déjà disponible puisque vous travaillez avecspplot). En utilisant l'ensemble de données 'meuse.grid', votre code pourrait ressembler à peu près à ce qui suit.

## test de point d'échantillon <- data.frame(x = 180000, y = 331500) coordonnées (test) <- ~ x + y proj4string(test) <- "+init=epsg:28992" ## afficher les données maillées, y compris l'échantillon point spplot(meuse.grid, "dist", scales = list(draw = TRUE), sp.layout = list("sp.points", test, pch = 16, cex = 2, col = "black"))


Voir la vidéo: Taso avaruudessa