R语言小白，第一次画图，查了很多步骤文章，发现R语言做地图上的分析这篇，很详细，完全没有报错，自己follow做了一遍，很强大。
1.准备工作
SHP等文件可以到国家基础地理信息中心的网站（http://nfgis.nsdi.gov.cn）下载，但这个政府网站总是打不开！可到微信公众号 iamExcelPro 发送 shapefile 获得下载地址，解压到同一文件夹下，如E:/rstudy/chinaprovince目录，有3个文件bou2_4p.dbf、bou2_4p.shp、bou2_4p.shx ，都是需要的。
R中install并library了几个可能会用到的包：
library(mapdata)
library(maptools)
library(ggplot2)
library(plyr)
2.处理地图数据china_map=readShapePoly("E:/rstudy/chinaprovince/bou2_4p.shp") # 读取地图空间数据
结果报错：Error: '\c' is an unrecognized escape in character string starting ""E:\rstudy\c"
原来路径中应该是“/”不是“\”，不能直接复制路径粘贴！！！会报错。
重新读取地图文件：china_map<-readShapePoly("E:\\rstudy\\chinaprovince\\bou2_4p.shp")
警示信息：readShapePoly is deprecated; use rgdal::readOGR or sf::st_read
虽然但是右侧变量栏中还是读出来了
plot一下看看，还不赖，包括了重庆和南海的岛屿，截图是酱婶儿的
plot(china_map)
图1
据说投影方式不对？那就用ggplot绘制，并使用polyconic投影方式，显示一下。
ggplot(china_map,aes(x=long,y=lat,group=group)) +
geom_polygon(fill="white",colour="grey") +
coord_map("polyconic")
画出来是酱婶儿的
图2
emmm......确实，图1太扁了
现在地图是可用的了，但还需要加载和拼接行政信息，以便能与业务数据映射。
x <- china_map@data #读取行政信息
xs <- data.frame(x,id=seq(0:924)-1) #含岛屿共925个形状
china_map1 <- fortify(china_map) #转化为数据框
china_map_data <- join(china_map1, xs, type = "full") #合并两个数据框
提示信息：Joining by: id
2.处理业务数据
有的教程在命令行里输入数据，很坑爹，没有从Excel表格转存来得方便。
直接读取xlsx文件需要装的包比较麻烦，把要用的Excel数据转存为csv格式文件比较简单。！要注意csv业务数据中的第1列的字段名为NAME，各省名称也要与地图数据框里的省名一致，便于合并。各省名称可是用以下命令查看并记下。
unique(china_map@data$NAME)
结果显示：
[1] 黑龙江省
[2] 内蒙古自治区
[3] 新疆维吾尔自治区
[4] 吉林省
[5] 辽宁省
[6] 甘肃省
[7] 河北省
[8] 北京市
[9] 山西省
[10] 天津市
[11] 陕西省
[12] 宁夏回族自治区
[13] 青海省
[14] 山东省
[15] 西藏自治区
[16] 河南省
[17] 江苏省
[18] 安徽省
[19] 四川省
[20] 湖北省
[21] 重庆市
[22] 上海市
[23] 浙江省
[24] 湖南省
[25] 江西省
[26] 云南省
[27] 贵州省
[28] 福建省
[29] 广西壮族自治区
[30] 台湾省
[31] 广东省
[32] 香港特别行政区
[33] 海南省
然后读取业务指标数据，并与地图数据合并，这里的业务指标为te：
mydata<-read.csv("E:\\rstudy\\te.csv") #读取指标数据，csv格式
china_data <- join(china_map_data, mydata, type="full") #合并两个数据框
提示信息：Joining by: NAME
3.绘制地图
现在根据我的业务数据te画一下填色地图看看：
ggplot(china_data, aes(x = long, y = lat, group = group, fill = te)) +
geom_polygon(colour="grey") +scale_fill_gradient(low="white",high="steelblue") + coord_map("polyconic")
# 边界颜色选了"grey"，fill为要填色的业务数据指标名称，填充类型为两色渐变填充，低值选了"white"，高值选了"plum1"
#指定投影方式为polyconic，获得常见视角中国地图
这里贴一下百度文库中搜到的R语言颜色表，很多颜色可选。
截图出来是酱婶儿的：
图3
居然有点好看hhh，但图中的背景色、坐标轴、经纬线都是不需要的，图例也可以放到左下角，用theme命令清除重新画一下，这次te的颜色选ghostwhite和darkseagreen3：
ggplot(china_data, aes(x = long, y = lat, group = group,fill = te)) +
geom_polygon(colour="grey") +
scale_fill_gradient(low="ghostwhite",high="darkseagreen3") + #指定渐变填充色，可使用RGB
coord_map("polyconic") + #指定投影方式为polyconic，获得常见视角中国地图
theme( #清除不需要的元素，图例放在左下角
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank(),
legend.position = c(0.2,0.3)
)
再截图看一下叭：
图4
精致，完美，这个配色可以记下了！
4.添加省名标签
参考笔记中还涉及到了显示省名标签，暂时用不到，不过还是试一下吧。
用省会城市的经纬度数据标注省名，在最开始解压出来的文件中有chinaprovincecity.csv文件，先用它读取省会城市坐标：
province_city <- read.csv("E:\\rstudy\\chinaprovince\\chinaprovincecity.csv")
然后基本重复3的步骤画图，防止省份名称显示过大，size设置为了2
ggplot(china_data,aes(long,lat))+
geom_polygon(aes(group=group,fill=te),colour="grey")+
scale_fill_gradient(low="ghostwhite",high="darkseagreen") +
coord_map("polyconic") +
geom_text(aes(x = jd,y = wd,label = province), data =province_city,size=2)+
theme(
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank(),
legend.position = c(0.2,0.3)
)
图5
到此为止就画完了，效果还可以叭，以后只用复制粘贴代码就很方便啦~~
网友评论