ggplot geom_ribbon绘制水热平衡图填充区域异常问题求助

问题描述 #

你尝试用ggplot绘制水热平衡图，要求：

• 降水曲线（Precipitations/2）高于温度曲线（Temperature）时，用蓝色填充区域
• 温度曲线高于降水曲线时，用红色填充区域

但实际生成的图表存在问题：蓝色区域仅部分显示，红色区域完全不显示，部分区域出现空白。

问题分析 #

问题出在geom_ribbon的使用逻辑和分组设置上：

1. group=1的错误使用：设置group=1会让所有年份的月份数据被视为一个整体连接，而facet_wrap(~Year)是按年份分面，导致ribbon在分面内无法正确识别分组，填充区域断裂。
2. ifelse生成NA导致填充中断：当条件不满足时用NA赋值给ymin或ymax，geom_ribbon遇到NA会中断填充，无法形成连续的区域。

解决方案 #

修正思路 #

1. 移除全局group=1：利用facet_wrap(~Year)的自动分组，每个分面内的月份数据会被正确识别为独立序列。
2. 用条件赋值替代ifelse生成NA：当不满足填充条件时，让ymin和ymax相等（无高度差），避免NA中断填充；条件满足时设置正确的上下限。
3. 明确填充范围逻辑：分别为红色、蓝色填充定义专属的上下限规则，确保仅在目标区域生成有效填充。

修正后的代码 #

# 加载必要包
library(tidyverse)

# 模拟数据
data <- tibble(
  Year = rep(2021:2024, each = 12),
  Month = factor(rep(month.abb, 4), levels = month.abb, ordered = TRUE),
  Precipitations = runif(48, min = 30, max = 100), # 降水量（mm）
  Temperature = runif(48, min = 0, max = 25)       # 温度（°C）
) %>%
  # 预处理：计算降水半值，简化后续代码
  mutate(half_precip = Precipitations / 2)

# 绘制水热平衡图
data %>% 
  ggplot(aes(x = Month)) +
  
  # 红色填充：温度 > 降水半值时，填充降水半值到温度之间的区域
  geom_ribbon(aes(
    ymin = ifelse(Temperature > half_precip, half_precip, Temperature),
    ymax = ifelse(Temperature > half_precip, Temperature, Temperature)
  ), fill = "red", alpha = 0.3) +
  
  # 蓝色填充：降水半值 > 温度时，填充温度到降水半值之间的区域
  geom_ribbon(aes(
    ymin = ifelse(half_precip > Temperature, Temperature, half_precip),
    ymax = ifelse(half_precip > Temperature, half_precip, half_precip)
  ), fill = "blue", alpha = 0.3) +
  
  # 降水曲线
  geom_line(aes(y = half_precip, color = "Precipitations")) +
  
  # 温度曲线
  geom_line(aes(y = Temperature, color = "Temperature")) +
  
  # 坐标轴设置
  scale_y_continuous(
    expression("Monthly temperatures (°C)"),
    sec.axis = sec_axis(~ . * 2, name = "Monthly precipitation (mm)")
  ) +
  
  # 自定义颜色图例
  scale_colour_manual("", 
                      values = c("Temperature" = "red", "Precipitations" = "blue")) +
  
  # 图例样式调整
  guides(colour = guide_legend(override.aes = list(linewidth = 3))) +
  
  # 按年份分面
  facet_wrap(~Year) +
  
  # 主题设置
  theme_bw()

代码说明 #

• 数据预处理：新增half_precip列，避免重复计算Precipitations/2，提升代码可读性。
• geom_ribbon逻辑优化：当条件不满足时，让ymin与ymax相等，该区域不会产生填充；条件满足时设置正确的上下限，生成连续填充区域。
• 移除group=1：facet_wrap已按Year分组，每个分面内的Month是有序因子，geom_ribbon会自动按分面内的序列连接点，避免跨分面的错误连接。

如果你偏好更清晰的条件逻辑，可以用case_when替代ifelse：

# 红色填充的简化写法
geom_ribbon(aes(
  ymin = case_when(Temperature > half_precip ~ half_precip, TRUE ~ Temperature),
  ymax = case_when(Temperature > half_precip ~ Temperature, TRUE ~ Temperature)
), fill = "red", alpha = 0.3)

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