我正试图在法国绘制不同面额的铅笔。使用以下代码,我可以实现我想要的地图,代表不同项目的不同颜色,以及每个部门/地区的每个项目的使用率。
ggplot() +
geom_polygon(data = plotDatafr, aes(x=long, y = lat, group = group, fill=item, alpha=score), colour = NA) +
scale_fill_manual(values = c("#009E73", "#F0E442", "#0072B2", "#D55E00"), name = "", na.value=NA) +
coord_map()
我现在正试图找到一种方法来平滑部门/地区之间的颜色过渡。地图here给出了我想要实现的目标。我的数据是here。我知道this post并且我尝试关注this tutorial,但是没有一个对我的问题有帮助。
修改
为了创建上面的地图,我根据每个部门的参与者总数(总共8272名参与者)计算每个部门给定项目(5个不同项目)的答案百分比。然后,我为每个部门选择了收到最大百分比的项目,并在地图上绘制这些项目。我在this folder中提供了调查的原始数据。
数据说明
“id_part”:参与者; “PLZ”:参与者城市的邮政编码 ; “Lat”/“Long”:参与者城市的长和纬度值 ; “NOM_DEPT”:参与者所在城市所在的部门/地区 ;其余列表示答案(5个选项)
Shapefile在文件夹中提供,我使用此代码连接数据(在用每个部门的最大百分比计算和提取项目之后)和shapefile:
library(rgdal)
mapa <- readOGR(dsn="France",layer="DEPARTEMENT")
mapa <- spTransform(mapa, CRS("+proj=longlat +datum=WGS84"))
mapa@data$id <- rownames(mapa@data)
mapa@data <- join(mapa@data, data, by="NOM_DEPT")
mapa.df <- fortify(mapa)
mapa.df <- join(mapa.df,mapa@data, by="id")
plotDatafr <- join(mapa.df, data)