POSIXct毫秒解析异常求助：0.001毫秒丢失如何解决？

问题现象 #

使用as.POSIXct解析含特定毫秒/微秒的时间字符串时，会出现精度丢失：

• 解析2023-08-03T14:25:00.001Z时，毫秒位被舍入为0：

as.POSIXct("2023-08-03T14:25:00.001Z", format = "%Y-%m-%dT%H:%M:%OSZ", tz = "UTC")
# 输出："2023-08-03 14:25:00.000 UTC"

• 解析2023-08-03T14:25:00.001001Z并输出6位微秒时，最后一位被舍入为0：

library(magrittr)
as.POSIXct("2023-08-03T14:25:00.001001Z", format = "%Y-%m-%dT%H:%M:%OSZ", tz = "UTC") %>% 
  strftime(format = "%Y-%m-%dT%H:%M:%OS6Z")
# 输出："2023-08-03T14:25:00.001000Z"

使用环境：R 4.3.0 (x86_64-pc-linux-gnu)

原因分析 #

这并非POSIXct的Bug，而是二进制浮点数的精度限制导致的：
POSIXct以浮点数存储从1970-01-01 00:00:00 UTC以来的秒数，而十进制的0.001（毫秒）无法用二进制浮点数精确表示，存储时会产生微小的误差。当这个误差刚好处于舍入阈值时，就会被截断为0；其他毫秒值的误差未触及舍入阈值，因此表现正常。

规避方案 #

方案1：使用lubridate包的高精度时间解析函数 #

lubridate包的ymd_hms函数对小数秒的解析精度更稳定：

library(lubridate)
# 解析毫秒时间
ymd_hms("2023-08-03T14:25:00.001Z", tz = "UTC")
# 输出："2023-08-03 14:25:00.001 UTC"

# 解析微秒时间并输出6位精度
ymd_hms("2023-08-03T14:25:00.001001Z", tz = "UTC") %>% 
  strftime(format = "%Y-%m-%dT%H:%M:%OS6Z")
# 输出："2023-08-03T14:25:00.001001Z"

方案2：使用readr包的parse_datetime函数 #

readr包的解析函数采用了不同的精度处理逻辑，能避免浮点数舍入问题：

library(readr)
# 解析毫秒时间
parse_datetime("2023-08-03T14:25:00.001Z", format = "%Y-%m-%dT%H:%M:%OSZ", tz = "UTC")

# 解析微秒时间并输出6位精度
parse_datetime("2023-08-03T14:25:00.001001Z", format = "%Y-%m-%dT%H:%M:%OSZ", tz = "UTC") %>% 
  strftime(format = "%Y-%m-%dT%H:%M:%OS6Z")

方案3：手动拆分时间字符串并计算 #

如果不想依赖第三方包，可以手动拆分时间的整数秒和小数部分，再合并：

parse_precise_time <- function(time_str, tz = "UTC") {
  # 移除末尾的Z，拆分整数秒和小数部分
  cleaned_str <- sub("Z$", "", time_str)
  parts <- strsplit(cleaned_str, "\\.")[[1]]
  
  # 解析整数秒部分
  base_time <- as.POSIXct(parts[1], format = "%Y-%m-%dT%H:%M:%S", tz = tz)
  if (length(parts) == 1) {
    return(base_time)
  }
  
  # 计算小数秒并添加到基础时间
  decimal_sec <- as.numeric(paste0("0.", parts[2]))
  base_time + decimal_sec
}

# 测试毫秒时间
parse_precise_time("2023-08-03T14:25:00.001Z")
# 输出："2023-08-03 14:25:00.001 UTC"

# 测试微秒时间并输出6位精度
parse_precise_time("2023-08-03T14:25:00.001001Z") %>% 
  strftime(format = "%Y-%m-%dT%H:%M:%OS6Z")
# 输出："2023-08-03T14:25:00.001001Z"

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