使用ggplot绘制栅格数据的两类问题：白边去除与NA值图例标记

嘿，我来帮你搞定这两个ggplot栅格绘图的问题！直接上代码和实用解释，包你解决痛点~

先基于你提供的数据模拟一份带NA值的示例数据（方便演示效果）：

library(ggplot2)

# 模拟你的栅格数据，特意添加几个NA值
df <- data.frame(
  X = rep(seq(75.875, 76.875, by = 0.25), 2),
  Y = rep(c(12.375, 12.125), each = 5),
  Z = c(2.6510951, 5.4310212, NA, -3.2694589, -1.6689374, 
        -9.1670256, NA, 1.1797042, 0.1593383, -0.56789)
)

问题1：彻底去除图像周围的白色边框 #

要消除白边，需要从两个关键维度入手：

• 让坐标轴严格匹配数据范围，不自动扩展留白
• 清除绘图区域的所有外边距

代码示例：

# 先搭建基础绘图框架
base_plot <- ggplot(df, aes(X, Y, fill = Z)) +
  geom_raster() +
  # 核心：关闭坐标轴自动扩展，让栅格填满整个绘图区域
  scale_x_continuous(expand = c(0, 0)) +
  scale_y_continuous(expand = c(0, 0)) +
  # 清除绘图区域的外边距
  theme(
    plot.margin = margin(0, 0, 0, 0),
    # 可选：如果不需要坐标轴，直接隐藏它们（按需选择）
    axis.text = element_blank(),
    axis.title = element_blank(),
    axis.ticks = element_blank(),
    # 可选：清除背景网格和默认面板背景
    panel.grid = element_blank(),
    panel.background = element_blank()
  )

base_plot

解释：expand = c(0,0)会强制x/y轴范围和你的数据完全重合，没有额外留白；plot.margin = margin(0,0,0,0)则去掉了整个绘图区域外的所有边距。如果需要保留坐标轴，只保留scale_*和plot.margin的设置即可。

问题2：在图例中把NA值标记为“缺失值” #

针对连续型的Z值（你的数据是连续型），我们可以通过色阶函数和图例自定义来实现：

base_plot +
  scale_fill_viridis_c(
    na.value = "gray50",  # 给NA值设置一个辨识度高的填充色（比如灰色）
    name = "Z值"  # 设置图例标题
  ) +
  guides(
    fill = guide_colorbar(
      # 保留自动生成的连续色阶标签，末尾加上NA对应的“缺失值”标签
      label = c(waiver(), "缺失值"),
      # 可选：调整色条高度，确保标签显示完整
      barheight = unit(10, "lines")
    )
  )

如果你的Z值是离散型（比如被分组），可以这样处理：

# 先把Z转成因子，保留NA作为因子水平
df$Z_factor <- factor(df$Z, exclude = NULL)

ggplot(df, aes(X, Y, fill = Z_factor)) +
  geom_raster() +
  scale_x_continuous(expand = c(0,0)) +
  scale_y_continuous(expand = c(0,0)) +
  scale_fill_discrete(
    na.value = "gray50",
    # 替换NA对应的因子标签为“缺失值”
    labels = c(levels(df$Z_factor)[-which(is.na(levels(df$Z_factor)))], "缺失值")
  ) +
  theme(plot.margin = margin(0,0,0,0))

最终整合代码 #

把两个问题的解决方案结合起来，完整代码如下：

library(ggplot2)

df <- data.frame(
  X = rep(seq(75.875, 76.875, by = 0.25), 2),
  Y = rep(c(12.375, 12.125), each = 5),
  Z = c(2.6510951, 5.4310212, NA, -3.2694589, -1.6689374, 
        -9.1670256, NA, 1.1797042, 0.1593383, -0.56789)
)

ggplot(df, aes(X, Y, fill = Z)) +
  geom_raster() +
  scale_x_continuous(expand = c(0, 0)) +
  scale_y_continuous(expand = c(0, 0)) +
  scale_fill_viridis_c(
    na.value = "gray50",
    name = "Z值"
  ) +
  guides(
    fill = guide_colorbar(
      label = c(waiver(), "缺失值"),
      barheight = unit(10, "lines")
    )
  ) +
  theme(
    plot.margin = margin(0, 0, 0, 0),
    axis.text = element_blank(),
    axis.title = element_blank(),
    axis.ticks = element_blank(),
    panel.grid = element_blank(),
    panel.background = element_blank()
  )

运行这段代码就能得到没有白边、图例明确标记缺失值的栅格图啦！

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