如何使用geofacet与ggplot2为分面组添加矩形轮廓？

要实现仅用ggplot2给geofacet的指定分面添加轮廓或背景，核心思路是通过分面变量（code）关联矩形图层，让矩形自动匹配目标分面的位置，无需后期用magick修改图片。具体操作如下：

步骤1：准备高亮分面数据 #

先定义需要添加轮廓/背景的分面编码，再创建包含这些编码和矩形范围的数据框：

library(ggplot2)
library(geofacet)

# 指定要高亮的分面编码
highlight_codes <- c("NSW", "VIC")

# 创建高亮用数据框，设置矩形覆盖整个分面区域
highlight_df <- data.frame(
  code = highlight_codes,
  xmin = -Inf,
  xmax = Inf,
  ymin = -Inf,
  ymax = Inf
)

步骤2：修改原绘图代码，添加矩形图层 #

在原ggplot代码中插入geom_rect图层，通过inherit.aes = FALSE避免继承主图的映射关系，确保矩形仅出现在目标分面：

ggplot(aus_pop, aes(age_group, pop / 1e6, fill = age_group)) +
  # 可选：添加分面背景（浅黄半透明）
  geom_rect(data = highlight_df, 
            aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax),
            fill = "lightyellow", alpha = 0.3, inherit.aes = FALSE) +
  # 绘制主柱形图
  geom_col() +
  # 添加分面轮廓（红色粗边框，无填充）
  geom_rect(data = highlight_df, 
            aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax),
            color = "red", size = 1.5, fill = NA, inherit.aes = FALSE) +
  facet_geo(~ code, grid = "aus_grid1") +
  coord_flip() +
  labs(
    title = "Australian Population Breakdown",
    caption = "Data Source: ABS Labour Force Survey, 12 month average",
    y = "Population [Millions]") +
  theme_bw()

关键说明 #

• 用-Inf和Inf作为矩形范围，确保不管分面内的数据如何调整，矩形始终覆盖整个分面区域，位置自动对齐。
• 图层顺序很重要：背景矩形要放在柱形图之前，轮廓矩形要放在柱形图之后，避免被遮挡。
• 可根据需求调整轮廓的color、size，或背景的fill、alpha参数优化样式。

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