使用ggplot2的geom_contour_filled设置日期格式图例的技术问题

核心思路：geom_contour_filled基于数值型的z值生成等高线区间，我们需要将这些区间标签里的数值转换为%b %d格式的日期字符串，而非直接将z值转为日期对象（会导致离散水平重复等问题）。

方案一：自动提取等高线区间并修改标签 #

先让ggplot生成默认等高线，再解析自动生成的区间标签，替换为日期格式：

library(tidyverse)
library(stringr)

# 生成示例数据
predictions <- expand.grid(Var1 = 11437:11476, Var2 = 11475:11510)
predictions$Pred <- predictions[, 1] * (-1.60351)  + predictions[, 2] * (-0.9927131) + (41231.33)

# 绘制基础等高线图
p <- ggplot(predictions) +
  aes(x = as.Date(Var1), y = as.Date(Var2), z = Pred) +
  geom_contour_filled() +
  scale_x_date(date_labels = "%b %d") +
  scale_y_date(date_labels = "%b %d")

# 解析默认的区间标签，转换为日期格式
parse_level_label <- function(level_str) {
  # 提取区间内的数值
  num_vals <- str_extract_all(level_str, "\\d+\\.?\\d*")[[1]] %>% as.numeric()
  # 转为日期并格式化
  lower_date <- format(as.Date(num_vals[1]), "%b %d")
  upper_date <- format(as.Date(num_vals[2]), "%b %d")
  paste(lower_date, "–", upper_date)
}

# 获取原有填充层级并生成新标签
fill_levels <- levels(p$data$level)
new_labels <- sapply(fill_levels, parse_level_label)

# 替换图例标签，保留默认配色
p + 
  scale_fill_viridis_d(labels = new_labels, name = "预测日期") +
  guides(fill = guide_colorsteps(show.limits = TRUE))

方案二：手动指定等高线断点（更可控） #

如果需要自定义等高线的区间，可以先将日期转换为数值型断点，再生成对应标签：

library(tidyverse)
library(stringr)

predictions <- expand.grid(Var1 = 11437:11476, Var2 = 11475:11510)
predictions$Pred <- predictions[, 1] * (-1.60351)  + predictions[, 2] * (-0.9927131) + (41231.33)

# 转换Pred为日期，确定断点范围
pred_dates <- as.Date(predictions$Pred)
# 按周生成日期断点，可自定义间隔（如"10 days"）
date_breaks <- seq(floor_date(min(pred_dates), "week"), 
                   ceiling_date(max(pred_dates), "week"), 
                   by = "week")
# 转为数值型断点（适配geom_contour_filled的z值类型）
num_breaks <- as.numeric(date_breaks)

ggplot(predictions) +
  aes(x = as.Date(Var1), y = as.Date(Var2), z = Pred) +
  # 使用手动指定的断点生成等高线
  geom_contour_filled(breaks = num_breaks) +
  scale_x_date(date_labels = "%b %d") +
  scale_y_date(date_labels = "%b %d") +
  scale_fill_viridis_d(
    # 动态生成日期格式的标签
    labels = function(levels) {
      sapply(levels, function(lvl) {
        num_vals <- str_extract_all(lvl, "\\d+\\.?\\d*")[[1]] %>% as.numeric()
        paste(format(as.Date(num_vals[1]), "%b %d"), "–", format(as.Date(num_vals[2]), "%b %d"))
      })
    },
    name = "预测日期"
  ) +
  guides(fill = guide_colorsteps(show.limits = TRUE))

关键说明 #

• 不能直接将Pred转为日期对象传入z参数：geom_contour_filled需要数值型z值计算等高线，日期对象会被识别为离散值，导致重复水平报错。
• 使用scale_fill_viridis_d而非scale_fill_manual：保留默认的连续配色逻辑，避免手动指定颜色时的重复或配色混乱问题。
• stringr::str_extract_all用于解析等高线区间字符串（如[11413,11420)）中的数值，再转换为日期格式。

内容的提问来源于stack exchange，提问作者Grégwar