查看原文
其他

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

2017-06-18 杜雨 数据小魔方

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


  • 如何在一个版面上呈现不同行政级别的地理信息;

  • 如何给版面添加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() 



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

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() 



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

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导航窗口。


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


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

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



#拼贴合并:

vie<-viewport(width=0.30,height=0.25,x=0.25,y=0.25)

p1;print(p2,vp=vie)



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


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


如果你能获取带有九段线数据的地图素材,那么这是最佳的解决方案,倘若不行,我们就要对图表后期进行处理(大部分人想到的方法是使用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()



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


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




阅读原文进入本人博客,所有文章均会在博客同步更新,博客内可直接跳转GitHub主页,直达所有案例数据仓库。(也可进入交流群获取案例数据)


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存