将西班牙语 shapefile 与数据框合并

Question

我在将西班牙各省的 shapefile 与特定数据框合并时遇到问题。合并它们之后，形状以某种方式与虚假省份相关联（f.e。巴塞罗那 returns 特鲁埃尔的形状）。我搜索了例子，但我不明白我做错了什么。

所有必要的文件都存储在这个文件夹中： Dropbox folder

library(maptools)
library(dplyr)
library(data.table)
library(reshape2)

我将所有必要的文件输入 R:

gor=readShapeSpatial('prov_map.shp')
prov=read.csv("prov.csv",sep = ';')
prov=subset(prov,select=-X)

我更正了数据文件中的省份名称，使它们在实际名称之前没有数字，然后我将结果变量转换为一个因子：

provcorr=colsplit(prov$Province," ",c("Prov_num","Province"))
prov$prov_num=provcorr$Prov_num
prov$province_nonum=provcorr$Province 
prov$provfact=as.factor(prov$province_nonum)

由于西班牙口音，shapefile 中的省名有点奇怪，所以我更正了它们，使它们与数据框中的相匹配：

prov_nom=c("Melilla","Ceuta", "Zaragoza","Zamora", "Bizkaia",
       "Valladolid","Valencia/València","Toledo","Teruel",
       "Tarragona","Soria","Sevilla","Segovia","Cantabria",
       "Salamanca","Pontevedra","Palencia","Asturias","Ourense",
       "Navarra","Murcia","Málaga","Madrid","Lugo","Rioja, La",
       "Lleida", "León", "Jaén", "Huesca", "Huelva", "Gipuzkoa",
       "Guadalajara","Granada","Girona","Cuenca","Coruña, A",
       "Córdoba", "Ciudad Real", "Castellón/Castelló","Cádiz",
       "Cáceres", "Burgos","Barcelona", "Balears, Illes", "Badajoz",
       "Ávila", "Almería", "Alicante/Alcant","Albacete","Araba/Álava") 
gor@data$NAMEBIEN=as.factor(prov_nom)

至此，形状与省名完美对齐。例如，如果我输入：

plot(subset(gor,gor@data$NAMEBIEN=='Barcelona'))

巴塞罗那的形状出现了。但是在我尝试将它们合并为：

gor@data=merge(gor@data,prov,by.x='NAMEBIEN',by.y='provfact',all.x=T)

现在已经不是这样了。当我输入与之前相同的代码时：

plot(subset(gor,gor@data$NAMEBIEN=='Barcelona'))

特鲁埃尔的形状出现了

这个问题快把我逼疯了。有什么建议么？

Answer 1

问题是合并改变了数据框的顺序。在原始数据框中，特鲁埃尔与第 9 个多边形相关联。但合并后 data.frame 中的第 9 个条目是巴塞罗那，因此与特鲁埃尔的多边形相关联。添加参数 sort = FALSE 它应该可以工作。

gor@data=merge(gor@data,prov,by.x='NAMEBIEN',by.y='provfact',all.x=T, sort = FALSE)
plot(subset(gor,gor@data$NAMEBIEN=='Barcelona'))

将西班牙语 shapefile 与数据框合并

Merging spanish shapefile with dataframe

gis

merge

r

spatial

shapefile