接上一篇,R地图系列(1):maptools包绘制中国地图,如果没有读过第一篇的读者建议回去浏览一遍。
由上一篇文章得知,我们的数据已经保存在map中了。那么,如何将将包含经纬度的数据投影到地图上去呢?
数据准备
画图之前需要准备的数据包括(随便举个例子):
城市名,经度,纬度,人口数(自行车数等等),随便编一些数据,大家可以看看数据格式。
name latitude longitude num
杭州 30.27415 120.15515 100
上海 31.23037 121.4737 70
成都 30.5702 104.06476 60
OK,按照以上格式准备好后我们就可以绘图了。
我这边用的的数据涉及隐私,不便公布,只是最后用来展现一下效果。
将数据存在city 中,待会绘图会用到。
city=read_excel("citydata.xlsx",col_names = TRUE)
气泡图形式
执行以下代码:
mymap = ggplot(data = fortify(map)) +
geom_polygon(aes(x = long, y = lat, group=group), colour = "gray",fill = "white") +
labs(title="中国地图",x="",y="") +
geom_point(data=city,aes(x=longitude,y=latitude,size=num),colour="red",pch=1) +
geom_text(aes(x=longitude,y=latitude,label=name), data=city,color="black",size=1.5) +
theme(text=element_text(family="STKaiti",size=14))
可以得到:
气泡图形式的各个城市对应的数据就出来了,当然也可以在每个气泡标出城市名称。
文字图形式
执行以下代码:
mymap1 = ggplot(data = fortify(map)) +
geom_polygon(aes(x = long, y = lat, group=group), colour = "lightblue",fill = "white") +
labs(title="中国地图",x="",y="") +
geom_text(aes(x=longitude,y=latitude,label=name,size=num,colour=factor(res_num)), data=city) +
scale_color_manual("name",values = mycolors)
文字图形式的各个城市对应的数据就出来了,数据中我用了双重映射,其中文字的大小影射到个城市某个数据的大小,颜色的深浅也映射到数据的大小。
备注:颜色部分建议先读RColorBrewer与ggplot2
这样,我们以后就可以将任意维度的各个城市的数据(比如降雨量、人口、pm2.5等等)通过不同的方式(比如气泡图、问文字图、柱状图等ggplot2可以实现的图)投影到中国地图上了。