R语言绘制中国各省填色图，SHP文件

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(