Suite

Comment ajouter des étiquettes, etc. avec les fichiers de formes EIA dans R ?


J'essaie de créer une carte des bassins de schiste américains et je suis très nouveau dans la création de cartes. J'ai téléchargé les fichiers de formes de l'US Energy Information Administration (EIA). J'ai également pu créer une carte de base avec les bases tracées (veuillez consulter le code ci-dessous). Je voudrais ajouter des étiquettes, etc. à cette carte. Je pense que ceux-ci se trouveraient dans les autres fichiers téléchargés en dehors du fichier de formes.

library(maps) map("state") shapLocation<- "C:/shapfile/shalegasbasin/" pcontorta <- readShapePoly(paste(shapLocation,"US_ShaleBasins_EIA_May2011.shp",sep="")) plot(pcontorta, add=TRUE , col=adjustcolor("noir", alpha.f = 0.6), border=FALSE)

Shapefile de : http://www.eia.gov/maps/layer_info-m.cfm

Description d'un shapefile : http://en.wikipedia.org/wiki/Shapefile


Jetez un œil à ?texte. Vous pouvez placer les étiquettes en utilisant les coordonnées du polygone qui renverront les centroïdes approximatifs du polygone.

bibliothèque(sp) Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2))) Sr2 = Polygon(cbind(c(5,4,2 ,5),c(2,3,2,2))) Sr3 = Polygone(cbind(c(4,4,5,10,4),c(5,3,2,5,5))) Srs1 = Polygones(list(Sr1), "1") Srs2 = Polygones(list(Sr2), "2") Srs3 = Polygones(list(Sr3), "3") SpP = SpatialPolygons(list(Srs1,Srs2,Srs3) , 1:3) SpP = SpatialPolygonsDataFrame(SpP, data.frame(ID=1:3)) # Tracer les polygones et placer le texte en fonction des centroïdes des polygones tracer (SpP) text(coordinates(SpP)[,1], coordinates(SpP )[,2], coller("p",SpP$ID,sep="-"))