R语言地图复现遇内存不足报错求助

这种内存耗尽问题通常是因为分面渲染时，每个子图都会重复加载完整的高分辨率空间几何数据，加上长格式的carbon观测数据，导致内存占用急剧上升。试试下面这些针对性的解决方案：

• 简化空间几何精度：用rmapshaper包的ms_simplify()降低几何顶点数量，大幅减少内存占用。调整keep参数平衡地图精度和内存需求：

  library(rmapshaper)
  fips_us_simplified <- ms_simplify(fips_us, keep = 0.1) # 保留10%顶点，可按需调整
  us_outline_simplified <- ms_simplify(us_outline, keep = 0.05)
  

  之后用简化后的数据集重新运行绘图代码。

• 清理冗余数据：检查fips_us是否存在重复的几何或观测记录，用dplyr去重：

  library(dplyr)
  fips_us_clean <- fips_us %>% distinct(geometry, month, carbon, .keep_all = TRUE)
  

• 扩大R内存限制（仅Windows）：如果你用的是64位R，可以手动调高内存分配上限：

  memory.limit(size = 64000) # 分配64GB内存，根据自己机器的配置调整数值
  

  Linux和macOS系统下R会自动利用可用内存，无需手动设置。

• 手动拼接子图替代分面：用patchwork包逐个生成单月份地图再拼接，避免ggplot分面时的内存累积：

  library(patchwork)
  # 生成每个月份的地图列表
  plot_list <- lapply(unique(fips_us$month), function(m) {
    ggplot(fips_us[fips_us$month == m, ]) +
      geom_sf(aes(fill = carbon), color = NA) +
      geom_sf(data = us_outline, fill = NA, size = 0.25) +
      scale_fill_scico(palette = "bilbao") +
      labs(title = paste("Month:", m))
  })
  # 按3列拼接所有地图
  wrap_plots(plot_list, ncol = 3)
  

• 换用更高效的绘图包：tmap对大空间数据的内存处理更友好，分面渲染效率更高：

  library(tmap)
  tm_shape(fips_us) +
    tm_fill("carbon", palette = "bilbao") +
    tm_shape(us_outline) +
    tm_borders(size = 0.25) +
    tm_facets(by = "month")
  

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