使用 ggplot2 绘制比较各省份及其区县的详细数据-阿里云开发者社区

使用 ggplot2 绘制比较各省份及其区县的详细数据

2023-05-25 119

版权

本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。

简介： 使用 ggplot2 绘制比较各省份及其区县的详细数据

简介

赵小编之前介绍过完整中国，及区域地图的画法。如：使用 ggplot2 绘制单个和多个省份地图；R 语言绘制十段线地图，给特定省份填色；今天小编将介绍如何比较各个省份及其区县的详细数据，本文参考《R 语言数据可视化之美》^[1]。

导入 R 包

安装并加载所需的 R 包：

library(rgdal)   
library(ggplot2)
library(dplyr)
library(grid)
library(RColorBrewer)

数据处理

这里需要用到包含县级信息的 BOUNT_poly.shp 文件，并利用 readOGR()函数来读取。

该文件可以在 GitHub^[2] 中可以下载。

dataProjected <- readOGR("China_adm_shp/bou4_4m/BOUNT_poly.shp") 
dataProjected@data$id <- rownames(dataProjected@data)

利用fortify()函数将对象转化为数据框形式，以便后续使用。

watershedPoints <- fortify(dataProjected)
df_China4 <- merge(watershedPoints, dataProjected@data, by = "id")

字段 ADCODE99 是行政区域代码，前两位表示所属的省。利用 substr() 函数将其取出，并存储在 province 变量中。

df_China4$province<-substr(as.character(df_China4$ADCODE99),1,2)

各县级的实际数据，由于没有真实数据，本文采用了随机数。读者可以使用自己的数据进行替换。数据形式如下：

##字段 NAME99 是各个县的中文名。
mydata<-data.frame(NAME99=unique(df_China4$NAME99),
                   value=round(runif(length(unique(df_China4$NAME99)),0,10)))
df_China<-full_join(df_China4,mydata)

由于 df_China$province 保存的是各省份的编号，为了更直观地展示信息，将各编号替换为省份名称。

# ProCode文件包含省份编号和名称
ProCode <- read.csv("PcCode.csv") 
ProCode$ID<-as.character(ProCode$ID)
for (i in 1:nrow(ProCode)){
  df_China$province<- gsub(pattern = ProCode$ID[i], x = df_China$province,replacement =ProCode$PcCode[i])
}
df_China$province<-factor(df_China$province,levels=as.character(ProCode$PcCode))

图形绘制

使用 ggplot 及facet_wrap()分面函数来完成图形绘制。

df_China <- na.omit(df_China)
ggplot(df_China, aes(x = long, y = lat,group=group,fill=value)) +
  geom_polygon(colour="black",size=0.1) +
  scale_fill_gradientn(colours=brewer.pal(9,"Spectral")) +
  facet_wrap(~province, ncol = 7,scales="free",strip.position = "top")+
  coord_cartesian()+
  theme_void()+
  theme(strip.background = element_rect(color="grey60",size=0.25),
        panel.background = element_rect(fill=NA,color="grey60",size=0.25))

由于南海地区维度跨度太大，不便于数据展示，可对其进行调整。添加下列语句：

df_China<-df_China[df_China$NAME99!='南沙群岛' &
                     df_China$NAME99!='西沙群岛' &
                     df_China$NAME99!='中沙群岛' &
                     df_China$NAME99!='中沙群岛的岛礁及其海域', ]
ggplot(df_China, aes(x = long, y = lat,group=group,fill=value)) +
  geom_polygon(colour="black",size=0.1) +
  scale_fill_gradientn(colours=brewer.pal(9,"Spectral")) +
  facet_wrap(~province, ncol = 7,scales="free",strip.position = "top")+
  coord_cartesian()+
  theme_void()+
  theme(strip.background = element_rect(color="grey60",size=0.25),
        panel.background = element_rect(fill=NA,color="grey60",size=0.25))

如果只想展示部分省份，可以添加下列语句：

df_China <- df_China[df_China$province == c('上海市','浙江省','江苏省','安徽省') ,]
ggplot(df_China, aes(x = long, y = lat,group=group,fill=value)) +
  geom_polygon(colour="black",size=0.1) +
  scale_fill_gradientn(colours=brewer.pal(9,"Spectral")) +
  facet_wrap(~province, ncol = 7,scales="free",strip.position = "top")+
  coord_cartesian()+
  theme_void()+
  theme(strip.background = element_rect(color="grey60",size=0.25),
        panel.background = element_rect(fill=NA,color="grey60",size=0.25))