我已经使用Raster包上传各省的西班牙地图(2级),我想根据人均收入填充颜色。 这里有每个收入的文件和分配给每个省的数字。请注意,ID_2是包Raster分配的编号,变量PROV是西班牙政府指定的官方编号。
library(raster)
esp<-getData('GADM', country="ESP", level=2)
espPols <- unionSpatialPolygons(esp, esp$ID_2)
renta <- read.table("renta.csv",sep = ";", header=TRUE)
espMapRenta <- SpatialPolygonsDataFrame(espPols, renta)
plot(espMapRenta)
我遇到的第一个问题是包装中有一些省份重复,第二个问题是我不知道如何按收入水平以渐变颜色填充每个省份。
非常感谢你的帮助!! PS。人均收入数据的链接如下:https://www.dropbox.com/s/si6zpv7p2nap9zg/renta.csv?dl=0