R中同时绘制小时级与日度水文数据的日期转换及绘图报错问题

核心问题说明 #

你遇到的三类报错/异常本质是时间类型选择、解析函数使用不匹配导致的：

1. as.Date()仅存储年月日粒度的日期值，会强制丢弃时分秒信息，直接导致小时级数据被合并
2. ymd_hms()仅适用于解析文本格式的YYYY-MM-DD HH:MM:SS字符串，无法直接识别Excel存储的数值型日期序列号，因此会全部返回NA
3. 未转换为R标准时间类型时，ggplot2的日期变换函数无法识别输入值，触发date_trans works with objects of class Date only报错

分步解决流程 #

1. 正确转换为带时分秒的标准时间类型 #

注意：必须使用POSIXct类型存储时间列，该类型是R原生支持的、可精确到秒级的时间格式，既保留小时数据粒度，也可被绘图函数正常识别
根据你的数据导入场景，二选一即可：

• 方案A：导入时直接指定时间列类型（推荐，无后续转换误差）

library(readxl)
# 读取Excel时直接指定第一列为日期时间类型，其余为数值型
df <- read_excel("你的水文数据文件.xlsx", 
                 col_types = c("date", "numeric", "numeric", "numeric"))

• 方案B：对已经导入为数值型的Date列做转换

library(dplyr)
library(lubridate)
# 将Excel数值型日期转换为带时分秒的POSIXct格式
df <- df %>%
  mutate(Date = excel_numeric_to_date(Date, include_time = TRUE))

转换完成后可执行class(df$Date)校验，返回结果包含POSIXct即为转换成功。

2. 校验日度数据对齐逻辑 #

转换完成后可快速确认河水位、降水值是否正确对齐到每日12:00：

# 提取所有12点的记录校验非NA值
check_align <- df %>% filter(hour(Date) == 12)
head(check_align)

若所有非NA的River、Rain值均落在12:00时间点，说明对齐逻辑符合预期。

3. 绘图注意事项 #

由于x轴是秒级精度的POSIXct类型，绘图时必须搭配scale_x_datetime()做轴配置，不能使用仅支持日期粒度的scale_x_date()，基础绘图示例如下：

library(ggplot2)
ggplot(df, aes(x = Date)) +
  # 绘制小时级井水位折线
  geom_line(aes(y = `Well 1`, color = "井水位"), linewidth = 0.7) +
  # 绘制日度河水位点
  geom_point(aes(y = River, color = "河水位"), size = 2) +
  # 绘制日度降水柱形
  geom_col(aes(y = Rain*20, fill = "降水量"), alpha = 0.4, width = 3600*10) +
  # 配置时间轴，按天显示刻度
  scale_x_datetime(date_breaks = "1 day", date_labels = "%m-%d") +
  labs(x = "时间", y = "水位(m)", color = "", fill = "") +
  theme_bw()

避坑清单 #

• 不要用as.Date()处理需要保留小时粒度的时间列，该函数会自动截断时分秒部分
• 不要对数值型Excel日期直接调用ymd_hms()、dmy_hms()这类文本时间解析函数，必然返回NA
• 不要对POSIXct类型的时间轴使用scale_x_date()，会触发类型不匹配报错

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