今天跟大家分享三个数据地图的遗留问题,包括以前因为技术手段限制无法解决的问题,以及读者吐槽最多的问题。
如何在一个版面上呈现不同行政级别的地理信息;
如何给版面添加mini导航定位窗口;
关于九段线问题。
关于第一个问题的解决方案,其实很简单,就是通过两个不同行政级别的图层进行叠加来达到目的。
library(maps)
library("ggplot2")
library("ggthemes")
library("dplyr")
library("maptools")
library("grid")
library("ggimage")
world_map<-map_data("world")
data<-read.csv("D:/R/mapdata/Word_State.csv",stringsAsFactors = FALSE,check.names = FALSE)
data1<-filter(data,region=="Russia" | Address=="Asia")
mymapdata<-merge(world_map[,1:5],data1,type="right")%>%filter(long>=30 & lat>=10)%>%arrange(group,order)
#底层的区域地图:
ggplot(mymapdata)+
geom_polygon(aes(x=long,y=lat,group=group),fill="grey95",col="black")+
coord_map("polyconic") +
theme_map()
![image.png](https://ask.hellobi.com/uploads/article/20170619/fb065ed1c5fb1b09b2de60d074903421.png)
#表层的中国省级行政区划地图:
china_map <- readShapePoly("D:/R/rstudy/CHN_adm/bou2_4p.shp")
china_map1 <- fortify(china_map)
ggplot(china_map1)+
geom_polygon(aes(x=long,y=lat,group=group),fill="white",col="black")+
coord_map("polyconic") +
theme_map()
![image.png](https://ask.hellobi.com/uploads/article/20170619/092011fea5a21eb54657f7888bcf8caf.png)
#将两个图层对象放在同一个图形对象中(注意先后顺序)。
ggplot()+
geom_polygon(data=mymapdata,aes(x=long,y=lat,group=group),fill="grey95",col="grey")+
geom_polygon(data=china_map1,aes(x=long,y=lat,group=group),fill="white",col="black")+
coord_map("polyconic") +
theme_map()
![image.png](https://ask.hellobi.com/uploads/article/20170619/b392c1eb1692c3bd6bc15143ac65026b.png)
第二个问题,关于minimini导航窗口。
其实这个问题多少跟第一个问题有些关联,也是旨在给我们要呈现的中心地图一个区域定位的作用。解决方案依然是拼图,不过这次是拼的图表(整个图表对象),而非图表中的图层。
#位于中间的中国行政区划地图:
p1<-ggplot()+
geom_polygon(data=china_map1,aes(x=long,y=lat,group=group),fill="white",col="grey")+
coord_map("polyconic") +
theme_map();p1
p2<-ggplot()+
geom_polygon(data=world_map,aes(x=long,y=lat,group=group),fill="grey95",col="grey")+
annotate(geom="rect",xmin =30, xmax =155,ymin =10L, ymax = 55,colour="grey",fill="orange",alpha=.3)+
theme_map() %+replace% theme(plot.background=element_rect(colour="black"));p2
![image.png](https://ask.hellobi.com/uploads/article/20170619/8e9985d21cbd647be46378d2000173af.png)
#拼贴合并:
vie<-viewport(width=0.30,height=0.25,x=0.25,y=0.25)
p1;print(p2,vp=vie)
![image.png](https://ask.hellobi.com/uploads/article/20170619/d01e1f6672cf1c5198a733a18fad0649.png)
第三个问题;也是被吐槽最多的九段线问题,这个问题,其实严格来说,我也没有办法搞定。
一般的九段线标识方法有两种,一种是在实际位置框定九段线范围,另一种则是将九段线的实际范围单独移至版面的边缘位置。第一种方法需要精确的九段线经纬点数据,获取不易,第二种则九段线位置过于灵活,然后不好掌握数据。
如果你能获取带有九段线数据的地图素材,那么这是最佳的解决方案,倘若不行,我们就要对图表后期进行处理(大部分人想到的方法是使用PS等软件,当然 这也是一种解决方案,至少清晰度上可以保证)。
但是既然是在R里面写代码,突然跑去玩PS总感觉怪怪的,能用代码解决的问题就不算问题,没错拼图嘛,R依然搞的定。
首先你要有个九段线素材图片。
jiuduanxian<-"E:/数据可视化/R/R语言学习笔记/可视化/timg.jpg"
ggplot()+
geom_polygon(data=china_map1,aes(x=long,y=lat,group=group),fill="white",col="grey")+
geom_image(aes(image=jiuduanxian,x=130,y=20),size=.2)+
coord_map("polyconic") +
theme_map()
![image.png](https://ask.hellobi.com/uploads/article/20170619/7dc014a81217ab4597b3503f3db52f2d.png)
我找的这个九段线素材可能跟中国地图风格搭配的不是很好,当然各位使用的时候可以自己去找相对准确、美观的素材使用。
最后提醒大家,关于九段线问题,一定要在报告或者文章末尾进行免责声明,涉及到领土争端问题,不仅仅要意识到,而且要尽可能的做好应对措施,避免不必要的麻烦。
![](http://mmbiz.qpic.cn/mmbiz_png/x2S0icdguVd9OSChS1I0d2JFJVru08Vf4iaGdB5iaS7LJKJXAmOhEgPfz0zJYCIDicGxWtKT3wQeSApfe1zDkpmBDQ/640?wx_fmt=png&tp=webp&wxfrom=5&wx_lazy=1)