在R的ggplot中为无数据的北卡东北部县添加地图边界

你的问题核心是缺失的县不在当前过滤后的数据集里，或是这些县的空间几何未被正确加载，导致geom_sf无法绘制它们的轮廓与填充。以下是两种直接可行的解决方法：

方法1：确保数据集包含所有目标县（包括Rate为NA的） #

如果你的state_dat仅包含有有效Rate值的县，Rate为NA的县会被过滤掉，自然无法显示。你需要将所有东北部的县纳入数据集，哪怕它们的Rate是NA：

1. 准备一份包含北卡罗来纳州东北部所有县的完整空间数据集（比如从sf包自带的nc数据中筛选东北部区域，或使用官方地理数据）。
2. 将你的state_dat与这份完整数据集做左连接，确保每个县都保留空间几何，Rate值缺失的会自动填充为NA：

# 示例：用sf包的nc数据筛选东北部县（根据你的REGION规则调整）
full_ne_counties <- nc %>% filter(REGION == "1") # 替换成你的区域筛选逻辑
# 左连接保留所有县，匹配你的Rate数据
combined_dat <- full_ne_counties %>% 
  left_join(state_dat %>% select(COUNTY, Rate), by = "COUNTY") # 假设COUNTY是匹配字段

# 绘图
ggplot(combined_dat, aes(fill = Rate)) +
  geom_sf() +
  scale_fill_stepsn(colours = palmap, 
                    na.value = "#d0d1e6",
                    breaks = c(0,5, 10, 15, 20),
                    labels = c("0-5", "5-10", "10-15", "15-20", "20+"),
                    limits = c(0,25)) +
  theme_void()

这样Rate为NA的县会用你指定的#d0d1e6填充，所有县的轮廓都会正常显示。

方法2：单独添加边界图层 #

如果不想调整数据集，可以额外添加一个只绘制所有县边界的图层，确保缺失的县也有轮廓：

# 假设full_ne_counties是包含所有东北部县的空间数据
ggplot() +
  # 先绘制填充图层（你的原有数据）
  geom_sf(data = state_dat %>% filter(REGION=="1"), aes(fill = Rate)) +
  # 再绘制所有县的边界，填充设为NA只显示线条
  geom_sf(data = full_ne_counties, fill = NA, color = "gray30", size = 0.5) +
  scale_fill_stepsn(colours = palmap, 
                    na.value = "#d0d1e6",
                    breaks = c(0,5, 10, 15, 20),
                    labels = c("0-5", "5-10", "10-15", "15-20", "20+"),
                    limits = c(0,25)) +
  theme_void()

这个方法中，边界图层会覆盖所有县的轮廓，哪怕填充图层里没有这些县，也能保证边界连贯。

额外检查点 #

• 确认state_dat %>% filter(REGION=="1")是否真的包含了所有东北部的县，有些县可能REGION字段值不是"1"，需要核对地理数据的区域划分。
• 检查缺失县的空间几何是否有效：如果某个县的geometry字段是NA，geom_sf也无法绘制它，需要修复空间数据。

内容的提问来源于stack exchange，提问作者alocklear