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

最新推荐文章于 2024-01-25 12:11:14 发布

有翅膀

最新推荐文章于 2024-01-25 12:11:14 发布

阅读量933

点赞数 2

分类专栏： R 文章标签： ggplot2

本文链接：https://blog.csdn.net/weixin_44921915/article/details/111149704

版权

本文介绍如何利用ggplot2包绘制病例分级统计地图。首先，讲解了导入地图数据时需要注意的事项，包括文件夹结构和解决中文乱码问题。接着，详细阐述了数据预处理步骤，包括去除重复数据、提取关键列、统计病例数量以及数据框的合并。最后，展示了绘图代码，使用scale_fill_gradient设置颜色梯度，coord_map进行投影变换，并添加指北针和比例尺。由于保密原因，文章未展示完整地图结果。

摘要由CSDN通过智能技术生成

会用到的包

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,