如何提取xts对象中的每隔X个元素(生成指定步长时间序列)
对秒级xts时间序列降采样(保留每隔X秒的数据)
你有一个时间分辨率为秒级的xts对象,示例代码如下:
times <- seq.POSIXt(from=as.POSIXct("2023-04-01 00:00:00"), to=as.POSIXct("2023-04-04 00:01:00"), by = "sec") xdata <- runif(n=length(times)) series <- xts(xdata, order.by=times)
以下是几种实现保留每隔X秒数据的方法,可根据你的数据情况选择:
方法1:按元素间隔筛选(适用于严格等间隔秒级数据)
如果你的xts对象是严格每秒一条数据(无缺失),可以直接按元素索引间隔取值,比如保留每10秒的数据:
X <- 10 downsampled_series <- series[seq(1, nrow(series), by = X), ]
方法2:按时间戳条件筛选(适用于非严格等间隔数据)
如果数据存在缺失或时间间隔不固定,可通过时间戳的秒数条件筛选,精准保留秒数为X倍数的记录:
X <- 10 # 提取时间戳的秒数,筛选符合条件的行 downsampled_series <- series[seconds(index(series)) %% X == 0, ]
方法3:用xts内置降采样函数(支持周期聚合)
如果需要按时间周期处理(比如取每个X秒窗口的首个/末个值,或做均值、求和等聚合),可以使用endpoints配合period.apply:
X <- 10 # 按X秒分割时间序列,获取每个周期的端点索引 ep <- endpoints(series, on = "seconds", k = X) # 保留每个周期的最后一个点,若要首个点替换为first() downsampled_series <- period.apply(series, INDEX = ep, FUN = last)
内容的提问来源于stack exchange,提问作者smoff
相关产品推荐
相关产品推荐

