数据下载
GIS 地图有很多种存储格式，其中 shapefile 格式（.shp）可以通过 R 的 maptools 包打开。[这里]可以(http://cos.name/wp-content/uploads/2009/07/chinaprovinceborderdata_tar_gz.zip)下载中国地图的GIS数据,完全解压后包含三个文件（bou2_4p.dbf、bou2_4p.shp和bou2_4p.shx.
library(maptools)
mydat = readShapePoly("../china-province-border-data/bou2_4p.shp")
plot(mydat)
给地图上色
plot(mydat,col=gray(924:0/924))
image.png
plot(mydat,col="green")
length(mydat)
## [1] 925
names(mydat)
## [1] "AREA" "PERIMETER" "BOU2_4M_" "BOU2_4M_ID" "ADCODE93"
## [6] "ADCODE99" "NAME"
mydat中有 925 条记录，每条记录中含有面积（AREA）、周长（PERIMETER）、各种编号、中文名（NAME）等字段。其中中文名（NAME）字段是以 GBK 编码的。
#统计每个省 / 直辖市的多边形数目。
#环境是UTF-8时，需要iconv函数转化一下才能正常显示。
table(iconv(mydat$NAME, from = "GBK"))
image.png
根据ADCODE99画提取某个省
Shanghai = mydat[mydat$ADCODE99 == 310000,]
plot(Shanghai)
image.png
参考：
https://cosx.org/2014/08/r-maps-for-china
https://cosx.org/2009/07/drawing-china-map-using-r
https://www.jianshu.com/p/40e14f6e0aa0
