在R中如何将宽格式降雨txt数据转换为可用于分析的时间序列数据

前置依赖：如果未安装用到的第三方包，先运行以下命令安装：

install.packages(c("tidyr", "dplyr", "lubridate", "xts", "zoo"))

步骤1：读取原始txt文件 #

根据你的文件实际分隔符调整读取参数，示例适配空格+逗号混合分隔的原始格式：

# header设为TRUE识别第一行表头，strip.white自动清除列值前后多余空格
rain_wide <- read.table("你的降雨数据文件路径.txt", 
                        header = TRUE, 
                        sep = "", 
                        strip.white = TRUE, 
                        stringsAsFactors = FALSE)
# 处理列名：清除日序列后的逗号，统一为1-366的数值格式
colnames(rain_wide) <- c("Year", gsub(",", "", colnames(rain_wide)[-1]))

如果你的文件是纯逗号分隔，将sep = ""替换为sep = ","即可。

步骤2：宽表转长表，生成逐日降水对应关系 #

把按年横向排列的降水值，转为每行对应单天降水的长表结构：

library(tidyr)
library(dplyr)

rain_long <- rain_wide %>%
  pivot_longer(cols = -Year, 
               names_to = "day_of_year", 
               values_to = "rainfall") %>%
  # 把日序列转为数值型，方便后续生成日期
  mutate(day_of_year = as.integer(day_of_year))

步骤3：生成标准日期列，清洗无效数据 #

自动识别闰年，删除非闰年366天的无效数据：

library(lubridate)

rain_ts <- rain_long %>%
  # 用年份和日序自动生成标准日期，闰年非闰年自动适配
  mutate(date = make_date(year = Year, day = day_of_year)) %>%
  # 过滤掉非闰年366天生成的NA无效日期
  filter(!is.na(date)) %>%
  # 按日期升序排列，保证时间序列连续性
  arrange(date)

步骤4：校验时间序列连续性 #

检查是否存在缺测日期，可根据需求补全：

# 生成观测时间范围内的完整逐日日期序列
full_date_seq <- seq.Date(from = min(rain_ts$date), 
                          to = max(rain_ts$date), 
                          by = "1 day")
# 查找缺测日期
missing_dates <- full_date_seq[!full_date_seq %in% rain_ts$date]
length(missing_dates) # 输出为0说明时间序列无缺测

如果存在缺测，可根据分析需求选择插值方法补全，示例为线性插值：

library(zoo)

rain_ts <- rain_ts %>%
  # 先补全所有日期行
  complete(date = full_date_seq) %>%
  # 线性插值填充缺测降水值
  mutate(rainfall = na.approx(rainfall))

步骤5：转换为R常用时间序列对象，适配后续分析 #

根据你的分析需求选择对应格式：

• 转xts对象（适配大多数高频逐日时间序列分析场景）：

library(xts)
rain_xts <- xts(rain_ts$rainfall, order.by = rain_ts$date)

• 转ts对象（适配ARIMA等传统时序模型分析）：

# frequency设为365对应逐日序列，start参数指定序列起始的年份和当年日序
rain_ts_obj <- ts(rain_ts$rainfall, 
                  start = c(year(min(rain_ts$date)), yday(min(rain_ts$date))), 
                  frequency = 365)

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