R package：ggplot2画中国地图

1.安装包

install.packages("rgdal")
install.packages("ggplot2")

2.读取地图

2.1读取省界地图

m1 <-rgdal::readOGR("../bou2_4p.shp")

2.2读取市界地图

m2 <-rgdal::readOGR("市界_region.shp")

2.3读取县界地图

m3 <-rgdal::readOGR("县界_region.shp")

length(m1)
names(m1)
table(iconv(m1$NAME, from = "GBK"))
iconv(m1$NAME, from = "GBK")

面积（AREA），周长（PERIMETER），各种编号、中文名（NAME）等字段。
其中中文名（NAME）字段是以GBK编码的。利用iconv 格式转换函数来转换各省名称

3. 画图

library(ggplot2)
ggplot(data = m1, aes(x = long, y = lat, group = group)) +
  geom_polygon(fill="white",col = "grey") +
  coord_map("polyconic")#指定投影方向

省界.png

ggplot(data = m2[c(101:111),], aes(x = long, y = lat, group = group)) +
  geom_polygon(fill="white",col = "grey") +
  coord_map("polyconic")

山西.png

ggplot(data = m3[c(310:322),], aes(x = long, y = lat, group = group)) +
  geom_polygon(fill="white",col = "grey") +
  coord_map("polyconic")

吕梁.png

另一种画法

m = sf::st_read("../bou2_4p.shp")
m = sf::st_read("市界_region.shp")
m = sf::st_read("县界_region.shp")
plot(m)