r - fastPOSIXct 等效于将非 UTC 转换为 UTC

Question

嗨，我有一个长度为几百万的字符向量( rr )，它以 %Y-%m-%d %H:%M:%S 格式表示时间和日期戳。记录在澳大利亚/悉尼。

如何获得代表这一点的 POSIXct 对象(快速)。

我找到了 fastPOSIXct在 fasttime包，但为了准确，它要求原始字符串采用 GMT/UTC，(我的不是)，然后使用 tz 转换回正确的时区争论...

> head(rr)
[1] "2009-05-01 10:01:00" "2009-05-01 10:02:00" "2009-05-01 10:03:00" "2009-05-01 10:04:00"
[5] "2009-05-01 10:05:00" "2009-05-01 10:06:00"

> as.POSIXct(head(rr),tz="Australia/Sydney")
[1] "2009-05-01 10:01:00 EST" "2009-05-01 10:02:00 EST" "2009-05-01 10:03:00 EST"
[4] "2009-05-01 10:04:00 EST" "2009-05-01 10:05:00 EST" "2009-05-01 10:06:00 EST"

如果在完整的数据集上执行上述操作需要很长时间……因此任何速度改进都将不胜感激。谢谢。

Answer 1

受到 Dirk 对这个 qn 的回答的启发，我制作了这个包装器来处理一年中的一大堆日期:

    fastPOSIXct_generic <- function(x, mytz = "America/New_York")
{
    # Caution, read: ?DateTimeClasses
    stopifnot(is.character(x))
    times_UTC <- fastPOSIXct(x, tz='UTC')
    num_times <- as.numeric(times_UTC)
    t1 <- as.POSIXct(x[1], tz = mytz)
    t2 <- as.POSIXct(x[1], tz = "UTC")
    offset <- as.numeric(difftime(t1, t2, units = "secs"))
    daylightoffset <- as.POSIXlt(t1)$isdst
    # For this first 'time' in t1 and t2, remove possible impact of losing one hour by setting clocks one hour forward during summer months:
    offset <- offset + daylightoffset * 3600
    num_times <- num_times + offset
    new_num_times <- as.POSIXct(num_times, tz = mytz, origin = '1970-01-01')
    new_num_times2 <- new_num_times - as.POSIXlt(new_num_times)$isdst * 3600
    return(new_num_times2)
}

# Test Sydney time

mm <- as.POSIXct(c("2015-03-15 15:00:00", "2015-4-10 15:00:00", "2014-10-01 15:00:00", "2015-10-15 15:00:00"), tz = "Australia/Sydney")
# "2015-03-15 15:00:00 AEDT" "2015-04-10 15:00:00 AEST" "2014-10-01 15:00:00 AEST" "2015-10-15 15:00:00 AEDT"
aus_stamps <- as.character(mm)
aus_back <- fastPOSIXct_generic(x = aus_stamps, mytz = "Australia/Sydney")
#"2015-03-15 15:00:00 AEDT" "2015-04-10 15:00:00 AEST" "2014-10-01 15:00:00 AEST" "2015-10-15 15:00:00 AEDT"
identical(mm, aus_back)
# TRUE

我的用例几乎总是 UTC 到 America/New_York，到目前为止它似乎工作正常。我不知道它是否适用于其他时区；只是 dst 有时间前进一个小时的情况。