R ggplot2 fusionner avec shapefile et csv des données pour remplir les polygones

[NB: Cette question a été posée plus d'un mois, si l'OP a sans doute trouvé un autre moyen de résoudre leur problème. Je suis tombé sur elle tout en travaillant sur cette question relative à la. Cette réponse est incluse dans l'espoir qu'il sera bénéfique pour quelqu'un d'autre.]

Cela semble être ce que l'OP demande...

... et a été produit avec le code suivant:

require("rgdal")
require("maptools")
require("ggplot2")
require("plyr")
# read temperature data
setwd("<location if your data file>")
temp.data        <- read.csv(file = "levels.dat", header=TRUE, sep=" ", na.string="NA", dec=".", strip.white=TRUE)
temp.data$CODINE <- str_pad(temp.data$CODINE, width = 5, side = 'left', pad = '0')
# read municipality polygons
setwd("<location of your shapefile")
esp     <- readOGR(dsn=".", layer="poligonos_municipio_etrs89")
muni    <- subset(esp, esp$PROVINCIA == "46" | esp$PROVINCIA == "12" | esp$PROVINCIA == "3")
# fortify and merge: muni.df is used in ggplot
muni@data$id <- rownames(muni@data)
muni.df <- fortify(muni)
muni.df <- join(muni.df, muni@data, by="id")
muni.df <- merge(muni.df, temp.data, by.x="CODIGOINE", by.y="CODINE", all.x=T, a..ly=F)
# create the map layers
ggp <- ggplot(data=muni.df, aes(x=long, y=lat, group=group)) 
ggp <- ggp + geom_polygon(aes(fill=LEVEL))         # draw polygons
ggp <- ggp + geom_path(color="grey", linestyle=2)  # draw boundaries
ggp <- ggp + coord_equal() 
ggp <- ggp + scale_fill_gradient(low = "#ffffcc", high = "#ff4444", 
space = "Lab", na.value = "grey50",
guide = "colourbar")
ggp <- ggp + labs(title="Temperature Levels: Comunitat Valenciana")
# render the map
print(ggp)

Explication:

Les fichiers de formes importées dans R avec readOGR(...) sont de type SpacialDataFrame et ont deux sections principales: une ploygon section qui contient les coordonnées de tous les points sur chaque polygone, et un données section qui contient des informations sur chaque polygone (donc, une ligne par polygone). Ceux-ci peuvent être référencés, par exemple, à l'aide de muni@polygons et muni@data. La fonction d'utilité fortify(...) convertit le polygone de la section d'une trame de données organisées pour comploter avec ggplot. De sorte que le flux de production de base est:

[1] Import temperature data file (temp.data)
[2] Import polygon shapefile of municipalities (muni)
[3] Convert muni polygons to a data frame for plotting (muni.df <- fortify(...))
[4] Join columns from muni@data to muni.df
[5] Join columns from temp.data to muni.df
[6] Make the plot

Le joint doit être fait sur les champs communs, et c'est là que la plupart des problèmes viennent. Chaque polygone dans le shapefile original a un attribut ID unique. L'exécution de fortify(...) sur le shapefile crée une colonne, id, qui est basée sur ce. Mais il n'y a pas de colonne ID dans la section de données. Au lieu de cela, le polygone Identifiants sont stockés en tant que noms de ligne. Alors d'abord, nous devons ajouter une id colonne de muni@data comme suit:

muni@data$id <- rownames(muni@data)

Maintenant, nous avons une id champ dans muni@data et un correspondant id champ dans muni.df, afin que nous puissions faire la jointure:

muni.df <- join(muni.df, muni@data, by="id")

Pour créer la carte, nous aurons besoin de définir les couleurs de remplissage en fonction de la température. Pour ce faire nous avons besoin de joindre la LEVEL colonne de temp.data à muni.df. Dans temp.data il y a un champ CODINE qui identifie la municipalité. Il y a aussi, maintenant, un champ correspondant CODIGOINE dans muni.df. Mais il y a un problème: CODIGOINE est char(5), avec les zéros non significatifs, alors que CODINE est entier qui signifie que les zéros non significatifs sont manquantes (importé à partir d'Excel, peut-être?). Donc, juste de joindre ces deux champs ne produit pas de matchs. Il faut d'abord convertir CODINE en char(5) avec les zéros non significatifs:

temp.data$CODINE <- str_pad(temp.data$CODINE, width = 5, side = 'left', pad = '0')

Maintenant, nous pouvons rejoindre temp.dat à muni.df basée sur les champs correspondants.

muni.df <- merge(muni.df, temp.data, by.x="CODIGOINE", by.y="CODINE", all.x=T, a..ly=F)

Nous utilisons merge(...) au lieu de join(...) parce que les champs de jointure ont des noms différents et join(...) exige qu'ils aient le même nom. (Notez, cependant, que join(...) est plus rapide et devrait être utilisé si possible). Donc, finalement, nous avons un bloc de données qui contient toutes les informations pour tracer les polygones et de la température LEVEL qui peut être utilisé pour définir la couleur de remplissage pour chaque polygone.

Quelques notes sur les OP le code original:

1. OP première carte (le vert en haut) identifie les "30 zones distinctes pour notre région...". Je n'ai trouvé aucun fichier de formes (shapefile l'identification de ces zones. La municipalité fichier identifie 543 municipalités, et j'ai pu voir aucun moyen de groupe dans ces zones 30. En outre, le niveau de température fichier a 542 lignes, une pour chaque municipalité (plus ou moins).

2. OP est de l'importation de fichiers pour la municipalité de tracer les limites. Vous n'en avez pas besoin car geom_polygon(...) allez dessiner (et de remplir) les polygones et les geom_path(...) permettra de tracer les limites.

OriginalL'auteur jlhoward