首先，了解一下ESRI Shapefile(shp)格式文件：

用于描述几何体对象：点，折线与多边形。例如，Shapefile文件可以存储井、河流、湖泊等空间对象的几何位置。除了几何位置，shp文件也可以存储这些空间对象的属性，例如一条河流的名字，一个城市的温度等等。(百度百科)

#加载必要的包

library(mapproj)

library(ggplot2)

library(rgdal)

library(plyr)

#注意：由于重庆的行政编吗为500000，在R语言内会自动显示为5e5，后面利用行政编码进行数据联结会出错

options(scipen = 100000000) #取消科学记数法

#读取我国地图的边界信息

path="C:\\Users\\admin\\Desktop\\map\\dxt\\省.shp"

Sys.getlocale("LC_CTYPE") #找到数据编码的格式避免乱码

china_map

#注意要dxt文件夹下加载省.shp文件

View(china_map) #china_map长下面这样！！！

plot(china_map)

#直接使用经纬度作为平面坐标进行绘图(左图)是不合适的,需要更改投影模式(右图)

ggplot(china_map,aes(x=long,y=lat,group=group))+geom_polygon(fill="white",colour="grey")+coord_map("polyconic")左右

#添加行政区域信息

x

xs

china_map1

china_map_data" />

#读取展示用到的数据mydata和空气污染数据czh,并与地理数据联结

czh

czh_1" />

china_map_data$FDPAC

china_data

#利用FDPAC行政编码，联结地图数据和czh数据

china_data$group

#画图

p

#按GDP的值作为背景色

guides(fill=FALSE)+scale_fill_gradient(low="white",high="steelblue")+

theme(panel.grid.major =element_blank(), panel.grid.minor = element_blank(),

panel.background = element_blank(),axis.line = element_line(colour = "black"))+

#去除网格线

geom_errorbar(data=czh_1,aes(x=x-0.25,ymin=y,ymax=y+0.7*so2/max(so2,no,fc)*3),size=2,color="orange",width=0,alpha=0.8)+ #二氧化硫柱形图

geom_errorbar(data=czh_1,aes(x=x+0.35,ymin=y,ymax=y+0.7*no/max(so2,no,fc)*3),size=2,color="brown",width=0,alpha=0.8)+ #二氧化氮柱形图

geom_errorbar(data=czh_1,aes(x=x+1.05,ymin=y,ymax=y+0.7*fc/max(so2,no,fc)*3),size=2,color="black",width=0,alpha=0.8)+ #烟粉尘柱形图

geom_text(aes(x=x,y=y-0.6,label=paste0(czh_1$NAME)),data=czh_1,size=3)+ #添加省份文字标签

annotate("text", x=130, y=22, label="二氧化硫", color= "orange", size=4)+

annotate("text", x=130, y=20, label="氮氧化物", color= "brown", size=4) +

annotate("text", x=130, y=18, label="粉尘烟尘", color= "black", size=4)

当当当，最后的效果图，研究了好久，快夸我！！！还有一些细节下次修改后再更新~

代码：

要数据的私戳我