ggplot循环绘制不同RCP场景空间地图及ggsave保存异常排查

循环生成的地图着色一致，核心原因是ggplot的延迟求值特性：循环里的rcp变量不会立刻被解析，直到所有循环结束后绘图才会读取变量值，这时rcp已经是最后一个循环的取值，导致所有图都用了最后一组RCP的ecodisrup数据。

解决方法1：修复for循环 #

在循环内部用local()包裹绘图代码，强制ggplot即时解析当前循环的变量值：

# 假设你的空间数据框是df，包含rcp、ecodisrup和geometry列
rcp_list <- unique(df$rcp)

for (rcp in rcp_list) {
  local({
    current_rcp <- rcp  # 把当前循环值存入局部变量，避免延迟求值问题
    p <- ggplot(df[df$rcp == current_rcp, ]) +
      geom_sf(aes(fill = ecodisrup)) +
      scale_fill_viridis_c(limits = c(0, 1)) +
      labs(title = paste("RCP场景:", current_rcp)) +
      theme_minimal()
    
    ggsave(paste0("rcp_", current_rcp, "_map.png"), p, width = 8, height = 6)
  })
}

解决方法2：用purrr包代替for循环（适配后续物种拓展） #

如果后续要拓展到RCP+物种的双重循环，用purrr的遍历函数更简洁，还能彻底避免求值问题：

library(purrr)
library(ggplot2)
library(sf)

# 假设数据新增species列，先按rcp和物种分组拆分
plot_groups <- df %>%
  group_by(rcp, species) %>%
  group_split()

# 遍历每个分组绘图
walk(plot_groups, function(group) {
  current_rcp <- unique(group$rcp)
  current_species <- unique(group$species)
  
  p <- ggplot(group) +
    geom_sf(aes(fill = ecodisrup)) +
    scale_fill_viridis_c(limits = c(0, 1)) +
    labs(title = paste("RCP:", current_rcp, "| 物种:", current_species)) +
    theme_minimal()
  
  ggsave(paste0("rcp_", current_rcp, "_species_", current_species, "_map.png"), 
         p, width = 8, height = 6)
})

关键注意点 #

• 无论用哪种方法，都要确保每次绘图用的是当前循环分组的子集数据，而非整个数据集，避免ggplot自动匹配最后一组的变量值
• 如果你的shapefile是分开存储的，要提前把shapefile路径和RCP（或物种）对应存入列表，遍历的时候同步读取对应空间文件

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