如何用ggplot绘制含交替NA值的反向时间偏移多时间序列？

针对你的Harz 4雪温和气象数据，我会一步步带你完成数据整理到绘图的全流程，解决NA值和时间偏移的问题：

1. 先做好数据准备（加载包+规整数据） #

首先确保你加载了必要的包，然后把两个数据集的时间列统一命名，方便后续对齐：

library(tidyverse)
library(readxl)

# 读取雪温数据，把第一列（时间）重命名为time，去掉空列
snow_temp <- read_excel("20180219_Harz-4_Tsnow.xlsx.xlsx", skip = 5) %>%
  rename(time = 1) %>%
  select(-matches("^X__")) # 清理可能存在的空列

# 读取气象数据，把时间列重命名为time
weat_stat <- read_excel("20180219_Harz-4_T_RH.xlsx.xlsx", skip = 5) %>%
  rename(time = `Datum Zeit`) %>%
  select(-matches("^X__"))

2. 处理时间偏移与NA值 #

情况1：时间序列是反向的（比如从新到旧排列） #

如果其中一个数据集的时间是倒序的，先把它转成正序：

# 按时间正序排列，确保时间轴连贯
snow_temp <- snow_temp %>% arrange(time)
weat_stat <- weat_stat %>% arrange(time)

情况2：两个数据集有固定时间偏移（比如雪温数据滞后气象数据） #

如果需要手动调整时间偏移（比如让雪温时间提前2小时对齐气象数据）：

# 给雪温时间加上2小时（滞后就用hours(-2)）
snow_temp <- snow_temp %>% mutate(time = time + hours(2))

合并数据并处理交替NA值 #

用full_join合并两个数据集，这样会保留所有时间点，缺失值自动显示为NA；再转成长格式（ggplot最喜欢的格式）：

combined_data <- full_join(snow_temp, weat_stat, by = "time") %>%
  pivot_longer(cols = -time, names_to = "参数", values_to = "数值")
# 可选：如果不想显示NA对应的断点，用drop_na(数值)去掉NA行；保留NA的话，geom_line会自动在NA处断开

3. 用ggplot绘制反向时间序列 #

现在可以画图了，要实现反向时间（x轴从新到旧），只需要反转x轴的时间范围；同时处理NA值：

ggplot(combined_data, aes(x = time, y = 数值, color = 参数)) +
  geom_line(na.rm = TRUE) + # na.rm=TRUE会跳过NA，避免线条中断报错
  # 反转x轴时间，实现"反向时间偏移"的效果
  scale_x_datetime(limits = rev(range(combined_data$time, na.rm = TRUE))) +
  labs(title = "Harz 4 雪温与气象参数时间序列",
       x = "时间",
       y = "数值",
       color = "参数类型") +
  theme_minimal()

进阶优化：突出数据点与NA位置 #

如果想更清晰看到每个数据点的位置（尤其是交替NA的地方），可以加上散点：

ggplot(combined_data, aes(x = time, y = 数值, color = 参数)) +
  geom_line(na.rm = TRUE, linewidth = 0.8) +
  geom_point(size = 0.5, na.rm = TRUE, alpha = 0.6) + # 半透明散点，避免遮挡线条
  scale_x_datetime(limits = rev(range(combined_data$time, na.rm = TRUE))) +
  labs(title = "Harz 4 雪温与气象参数时间序列",
       x = "时间",
       y = "数值",
       color = "参数类型") +
  theme_minimal()

这样处理后，你就能得到一个包含所有时间序列、自动处理交替NA值、且时间轴反向的可视化图啦！

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