目标：用R绘制标准的中国地图，并且分省给图片上色（分别生成中文版和英文版）。

用到的分省数据如图1所示。

图1 本文用到的分省市自治区数据

在自然资源部的标准地图服务中，可以找到很规范的中国地图、世界地图。标准地图服务系统 (mnr.gov.cn)。利用该网站的自助地图服务，可以绘制精美的中国地图，如图2所示：

图2 利用自助上色服务生成的中国地图（我自己添加了指北针和比例尺）

用该网站绘制中国地图的好处是准确简单高效，缺点是无法绘制英文版的地图。我们可以利用R语言绘制更加符合个性化需要的中国地图。我们要用到的中国地图GIS数据如下（如果大家需要地图压缩包，我可以在评论区给出网盘地址）：

bou2_4p.shp：中国政区的面文件

bou2_4l.shp：中国政区的线文件，在这个线文件里包含了南海的九段线

九段线.shp：南海九段线的线文件。我自己数过，我们中国地图的下方有十段线，那为什么会多出来一段线呢？原因是南海断续线之外，又在台湾岛与琉球群岛之间增加了一段线，形成了南海九段线、台湾岛东侧一段线的“海疆十段线”基本格局。

#####本文需要的包

library(rgdal)

library(sf)

#####导入本文需要的中国地图、中国边界数据和南海九段线，并分别画图。rgdal包的readOGR函数或sf包的st_read函数都可以读shp文件。大家导入了三张基础地图后，可以打开数据文件，看一下数据的基本结构。

china_map<-rgdal::readOGR('F:\\地图\\bou2_4p.shp')#引号内为bou2_4p.shp文件路径

plot(china_map)

图3 利用bou2_4p文件绘制的基础中国地图数据（用于给省市数据填色）

bianjing<-rgdal::readOGR('F:\\地图\\中国省市县\\中国国界与省界-含九段线\\bou2_4l.shp')

plot(bianjing,color='grey30')

图4 利用bou2_4l文件绘制的中国地图边界数据

nine<-st_read('F:\\地图\\R中国地图模板\\SouthSea\\九段线.shp')

plot(nine)

图5 利用九段线文件单独展示南疆十段线

#绘制中国地图数据

unique(china_map@data$NAME)

图6 china_map中的省市名

#读取china_map中的行政信息，并引入id信息，方便后面与china_map1合并，-1为了与china_map1一致同为从0开始

x <-china_map@data

xs <-data.frame(x,id=seq(1:925)-1)

#将china_map转化为数据框，合并xs和china_map1 这两个数据框，joining by id

china_map1 <-fortify(china_map)

china_map_data <-plyr::join(china_map1, xs,by="id", type = "full")

#引入南海九段线。结合南海九段线的经纬度，在边界信息中只选择九段线相关的数据。合并china_map_data和china_map2这两个数据框，joining by id

china_map2<-fortify(bianjing)

china_map2<-subset(china_map2,id==1087|id==1336|id==1377|id==1475|id==1481|id==1784|id==1769|id==1764|id==1508|id==1488)

china_map_data <-plyr::join(china_map_data ,china_map2,by="id", type = "full")

#引入各省市自治区的数据，合并china_map_data和mydata这两个数据框，joining by NAME

mydata<- read_xlsx('F:\\地图\\地图尝试1.xlsx',sheet= 'Sheet1')

windowsFonts(A=windowsFont("Times New Roman"),B=windowsFont("Arial"))

china_data<-join(china_map_data,mydata,type='full')

#合并中国地图和九段线+东海中日领海分界线

china_data1<-merge(china_data,china_map2,by=c("order","piece","id","group"),all=TRUE)

china_data1$long<-ifelse(is.na(china_data1$long.x)==TRUE,china_data1$long.y,china_data1$long.x)#china_data和china_map2的经纬度

china_data1$lat<-ifelse(is.na(china_data1$lat.x)==TRUE,china_data1$lat.y,china_data1$lat.x)

midpos <- function(x) mean(range(x,na.rm=TRUE))

#找出各省市自治区地理位置的中心点，为了图片美观，我们调整了图中部分文字的位置。

centres <- ddply(china_data,.(eng),colwise(midpos,.(long,lat)))

centres<-mutate(centres,long=ifelse(eng=="Hainan",110.1,

                                    ifelse(eng=="Hebei",115.2,

                                          ifelse(eng=="Gansu",104,

                                                  ifelse(eng=="Inner Mongoria",107,long)))))

centres<-mutate(centres,lat=ifelse(eng=="Hainan",19.32,

                                  ifelse(eng=="Sichuan",30.8,

                                          ifelse(eng=="Shaanxi",34.8,

                                                ifelse(eng=="Gansu",36,

                                                        ifelse(eng=="Guangdong",23.2,

                                                              ifelse(eng=="Guangxi",24,

                                                                      ifelse(eng=="Guangxi",24,

                                                                            ifelse(eng=="Hebei",38.1,

                                                                                    ifelse(eng=="Tianjin",39.1,

                                                                                          ifelse(eng=="Inner Mongoria",41,lat)))))))))))

##绘制地图

ggplot(china_data1,aes(long,lat))+

geom_polygon(aes(group=group,fill=count),colour='grey41')+#以人口的值区分各省份的颜色，col即colour缩写，为边界颜色

scale_fill_gradient(low='white',high='aquamarine3')+

coord_map('polyconic') +

geom_text(aes(label=eng),family="A",data=centres )+#横坐标x对应经度，纵坐标y对应维度，经纬度确认的坐标为标签显示的位置，label即标签为图中显示的值

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.85,0.45)#调整图例位置

      )

属于我们中国的领土一点不能少！

图7 利用R语言绘制的中国地图(含南海诸岛和国界九段线