应用ggplot2绘制病例分级统计地图

会用到的包

library(magrittr) #管道函数
library(stringr) #字符处理
library(ggplot2) #绘图
library(rgdal) #读取shp文件
library(plyr) #合并数据
library(legendMap) #加比例尺
library(ggspatial) #加指北针
library(tidyfst) # count_dt 
library(dplyr) # left_join 
library(patchwork) #拼图
library(RColorBrewer) # 挑颜色的

##我不记得哪些是GitHub的包了，反正这个是
devtools::install_github("3wen/legendMap")

四处找了很多文章，但别人做的跟自己的需求总是不太一样，东拼西凑加上自己码代码搞出了这篇文章。

导入地图数据

setwd("E:/数据分析")
shp_city <- readOGR("C:/Users/地图/quxian",encoding = "UTF-8")
shp_city <- spTransform(shp_city, CRS("+proj=longlat +datum=WGS84"))
dtframe <- fortify(shp_city) # 转化为数据框
data1 <- shp_city@data #提取行政区划信息
data2 <- data.frame(id=row.names(data1),data1) #添加序号
map_data <- join(dtframe,data2, type