如何通过FIP编码在单张地图上绘制带渐变效果的指定美国郡县
解决方法
现有代码的核心问题是将选中状态处理为二分类因子,同时搭配离散手动色阶配置,无法实现数值对应的连续渐变效果,可按以下步骤修改:
步骤1:预处理自定义数据
首先将FIP编码补全为5位字符串,同时处理重复FIP的值(示例数据中FIP=12129出现了2次,可根据需求选择求和、取均值等聚合方式,这里以求和为例)
library(tidyverse) library(usmap) library(grDevices) # 原始数据 mydf <- data.frame(FIP=c(1017,1059,1061,1133,5103,12129,13123,12129), val=c(100,100,225,400,180,278,180,321)) # 清洗mydf mydf_clean <- mydf %>% # 补全为5位FIP编码 mutate(fips = stringr::str_pad(FIP, 5, pad = "0")) %>% # 重复FIP聚合,也可根据需求替换为mean等函数 group_by(fips) %>% summarise(val = sum(val), .groups = "drop")
步骤2:合并到郡县基础数据
将自定义数值和usmap自带的全量郡县数据关联,未匹配到的非目标郡县val默认设为NA,后续可统一填充为浅灰色
dt <- countypop %>% left_join(mydf_clean, by = "fips")
步骤3:绘制渐变填充地图
使用连续色阶函数配置低蓝高红的渐变效果
usmap::plot_usmap(data = dt, values = "val", color = "grey") + # 配置蓝到红的连续渐变,非目标郡县填充为浅灰 scale_fill_gradient( low = "blue", high = "red", na.value = "lightgrey", name = "数值" # 可自定义图例名称 ) + theme_minimal()
如果需要更丰富的渐变效果,也可以用scale_fill_gradient2设置中间过渡色,或scale_fill_gradientn传入自定义色板。
内容的提问来源于stack exchange,提问作者Mathica
相关产品推荐
相关产品推荐

