如何高效为geom_sf多边形着色？ggplot分面地图性能优化

核心问题分析 #

渲染耗时剧增的本质是空间几何对象被重复复制：使用full_join将zones与包含多category的mydata合并后，每个区域的几何会被复制N次（N为category的数量），导致sf对象行数大幅膨胀，ggplot需要处理大量重复的几何数据，最终拖慢渲染速度。

优化方案 #

1. 避免重复复制空间数据，分面单独绘制 #

使用patchwork包将每个分面单独生成后拼接，每个分面仅加载对应category的数据，空间几何不会被重复复制：

library(patchwork)
library(tidyr)
library(dplyr)

# 按category分组嵌套数据
mydata_grouped <- mydata %>%
  group_by(category) %>%
  nest()

# 定义单个分面的绘图函数
plot_single_facet <- function(data, facet_title) {
  # 仅将当前分面的数据与空间数据关联
  zones_subset <- zones %>%
    left_join(data, by = "ZONE_NAME")
  
  ggplot() +
    # 绘制底图
    geom_polygon(data = basemap, aes(x = long, y = lat, group = group),
                 fill = "gray85", color = "gray80") +
    # 绘制当前分面的着色区域
    geom_sf(data = zones_subset, aes(fill = status)) +
    labs(title = facet_title) +
    theme_minimal()
}

# 生成所有分面图并拼接
all_facets <- pmap(mydata_grouped, plot_single_facet)
final_plot <- wrap_plots(all_facets)

2. 简化空间几何复杂度 #

如果zones shapefile的顶点过多，用sf::st_simplify减少顶点数，大幅降低渲染压力：

# 根据数据精度调整dTolerance（单位与空间数据的CRS一致，通常是米）
zones_simplified <- zones %>%
  sf::st_simplify(dTolerance = 100) # 示例值，可根据实际需求调整

3. 替换full_join为left_join #

full_join会保留mydata中不在zones里的ZONE_NAME，生成无几何信息的无效行，既无用又增加数据量，改用left_join仅保留有空间几何的区域：

# 原full_join替换为left_join
zones <- zones %>% left_join(mydata, by = "ZONE_NAME")

4. 关闭不必要的渲染选项 #

• 若不需要投影转换，添加coord_sf(datum = NA)跳过投影计算：

ggplot() +
  # ...其他图层...
  coord_sf(datum = NA)

• 若不需要区域边框，关闭边框绘制减少计算：

geom_sf(data = zones_subset, aes(fill = status), color = NA)

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