如何在R中按正负值填充曲线与阈值间的区域？

核心思路是补充曲线与y=0的交点数据，让填充区域在跨越零点时自动截断，无需大量插值即可解决geom_ribbon超出曲线范围的问题。

步骤1：准备可复现数据 #

先模拟一组类似的时间序列数据（替换为你的真实数据即可）：

library(ggplot2)
library(dplyr)

# 模拟北太平洋海温距平数据
set.seed(123)
date_seq <- seq(as.Date("2000-01-01"), as.Date("2020-12-01"), by = "month")
sst_anom <- cumsum(rnorm(length(date_seq), 0, 0.2))
df <- tibble(Date = date_seq, North_Pacific_SST_anom = sst_anom)

步骤2：计算曲线与y=0的交点 #

编写函数定位跨零点区间，通过线性插值计算交点坐标，再将交点整合到原始数据中：

# 找到曲线与y=0的交点并合并到数据
find_zero_crossings <- function(df, x_col, y_col) {
  x <- df[[x_col]]
  y <- df[[y_col]]
  
  # 定位所有跨越零点的相邻数据对
  cross_idx <- which(diff(sign(y)) != 0)
  
  # 计算每个跨零点区间的交点
  crossings <- lapply(cross_idx, function(i) {
    x1 <- x[i]
    x2 <- x[i+1]
    y1 <- y[i]
    y2 <- y[i+1]
    
    # 线性插值求交点的x值
    x_intersect <- x1 + (0 - y1) * (x2 - x1) / (y2 - y1)
    tibble(!!x_col := x_intersect, !!y_col := 0)
  }) %>% bind_rows()
  
  # 合并原始数据与交点，按x轴排序
  bind_rows(df, crossings) %>% arrange(!!sym(x_col))
}

# 处理数据
df_processed <- find_zero_crossings(df, "Date", "North_Pacific_SST_anom")

步骤3：分色填充绘图 #

使用geom_ribbon分别定义正负值区域的填充范围，即可得到无超出的分色填充效果：

ggplot(df_processed, aes(x = Date)) +
  # 填充正值区域（曲线到y=0之间）
  geom_ribbon(aes(ymin = 0, ymax = ifelse(North_Pacific_SST_anom > 0, North_Pacific_SST_anom, 0)),
              fill = "#FF6B6B", alpha = 0.5) +
  # 填充负值区域（曲线到y=0之间）
  geom_ribbon(aes(ymin = ifelse(North_Pacific_SST_anom < 0, North_Pacific_SST_anom, 0), ymax = 0),
              fill = "#4ECDC4", alpha = 0.5) +
  # 绘制原始曲线
  geom_line(aes(y = North_Pacific_SST_anom), color = "black", linewidth = 0.8) +
  labs(x = "日期", y = "北太平洋海温距平", title = "海温距平与y=0阈值间的分色填充") +
  theme_minimal()

替代方案：使用geom_area分组填充 #

也可以通过标记正负值分组，结合geom_area实现相同效果：

df_processed <- df_processed %>%
  mutate(anom_type = case_when(
    North_Pacific_SST_anom > 0 ~ "正距平",
    North_Pacific_SST_anom < 0 ~ "负距平",
    TRUE ~ "零"
  ))

ggplot(df_processed, aes(x = Date, y = North_Pacific_SST_anom)) +
  geom_area(aes(fill = anom_type), alpha = 0.5, position = "identity") +
  geom_line(color = "black", linewidth = 0.8) +
  scale_fill_manual(values = c("正距平" = "#FF6B6B", "负距平" = "#4ECDC4", "零" = "transparent")) +
  labs(x = "日期", y = "北太平洋海温距平", title = "海温距平与y=0阈值间的分色填充", fill = "距平类型") +
  theme_minimal()

这种方法仅补充必要的交点数据，比全局插值更高效，且能精准控制填充区域的边界，完美解决你遇到的geom_ribbon超出、geom_polygon路径错误的问题。

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