关于数据地图的几个遗留问题解决方案

今天跟大家分享三个数据地图的遗留问题,包括以前因为技术手段限制无法解决的问题,以及读者吐槽最多的问题。

  • 如何在一个版面上呈现不同行政级别的地理信息;
  • 如何给版面添加mini导航定位窗口;
  • 关于九段线问题。

问题一:

关于第一个问题的解决方案,其实很简单,就是通过两个不同行政级别的图层进行叠加来达到目的。

1
2
3
4
5
6
7
library(maps)
library("ggplot2")
library("ggthemes")
library("dplyr")
library("maptools")
library("grid")
library("ggimage")

1
2
3
4
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)

底层的区域地图:

1
2
3
4
ggplot(mymapdata)+
geom_polygon(aes(x=long,y=lat,group=group),fill="grey95",col="black")+
coord_map("polyconic") +
theme_map()

表层的中国省级行政区划地图:

1
2
china_map <- readShapePoly("D:/R/rstudy/CHN_adm/bou2_4p.shp")
china_map1 <- fortify(china_map)

1
2
3
4
ggplot(china_map1)+
geom_polygon(aes(x=long,y=lat,group=group),fill="white",col="black")+
coord_map("polyconic") +
theme_map()

将两个图层对象放在同一个图形对象中(注意先后顺序)。

1
2
3
4
5
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()

minimini导航窗口

第二个问题,关于minimini导航窗口。

其实这个问题多少跟第一个问题有些关联,也是旨在给我们要呈现的中心地图一个区域定位的作用。解决方案依然是拼图,不过这次是拼的图表(整个图表对象),而非图表中的图层。

位于中间的中国行政区划地图:

1
2
3
4
p1<-ggplot()+
geom_polygon(data=china_map1,aes(x=long,y=lat,group=group),fill="white",col="grey")+
coord_map("polyconic") +
theme_map();p1
1
2
3
4
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

拼贴合并:

1
2
vie<-viewport(width=0.30,height=0.25,x=0.25,y=0.25)
p1;print(p2,vp=vie)

第三个问题;也是被吐槽最多的九段线问题,这个问题,其实严格来说,我也没有办法搞定。

关于九段线问题:

一般的九段线标识方法有两种,一种是在实际位置框定九段线范围,另一种则是将九段线的实际范围单独移至版面的边缘位置。第一种方法需要精确的九段线经纬点数据,获取不易,第二种则九段线位置过于灵活,然后不好掌握数据。

如果你能获取带有九段线数据的地图素材,那么这是最佳的解决方案,倘若不行,我们就要对图表后期进行处理(大部分人想到的方法是使用PS等软件,当然 这也是一种解决方案,至少清晰度上可以保证)。

但是既然是在R里面写代码,突然跑去玩PS总感觉怪怪的,能用代码解决的问题就不算问题,没错拼图嘛,R依然搞的定。

首先你要有个九段线素材图片。

1
2
3
4
5
6
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()

我找的这个九段线素材可能跟中国地图风格搭配的不是很好,当然各位使用的时候可以自己去找相对准确、美观的素材使用。

最后提醒大家,关于九段线问题,一定要在报告或者文章末尾进行免责声明,涉及到领土争端问题,不仅仅要意识到,而且要尽可能的做好应对措施,避免不必要的麻烦。


关于数据源问题的说明:

点击博客左下侧的GitHub链接:

进入本人GitHub主页之后,打开名为datawarehouse的数据仓库;


联系方式:
wechat:ljty1991
Mail:578708965@qq.com
个人公众号:数据小魔方(datamofang)

qq交流群:[魔方学院]298236508

个人简介:
杜雨
财经专业研究僧;
伪数据可视化达人;
文科背景的编程小白;
喜欢研究商务图表与地理信息数据可视化,爱倒腾PowerBI、SAP DashBoard、Tableau、R ggplot2、Think-cell chart等诸如此类的数据可视化软件,创建并运营微信公众号“数据小魔方”。
Mail:578708965@qq.com


备注信息:
知识共享许可协议
本作品采用知识共享署名-非商业性使用 4.0 国际许可协议进行许可。

坚持原创技术分享,您的支持将鼓励我继续创作!