R地图系列(2)：maptools＋ggplot2包添加地图元素

接上一篇，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))

可以得到：

Paste_Image.png

气泡图形式的各个城市对应的数据就出来了，当然也可以在每个气泡标出城市名称。

文字图形式

执行以下代码：

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

Paste_Image.png

这样，我们以后就可以将任意维度的各个城市的数据（比如降雨量、人口、pm2.5等等）通过不同的方式（比如气泡图、问文字图、柱状图等ggplot2可以实现的图）投影到中国地图上了。