将儒略日转换为日期格式 yyyy-dd-mm HH:MM:SS

Converting julian days in to date format yyyy-dd-mm HH:MM:SS

借助以下函数(原创自R - convert POSIXct to fraction of julian day),我将不同的观察起点从日期格式转换为儒略日。

date <- as.POSIXct(c('2006-12-12 13:00:00', '2008-12-12 12:00:00', '2007-12-12 12:00:00'))

julian_conv <- function(x) { 
  if (is.na(x)) { 
    return(NA)
  }
  else {
    j <-julian(x, origin = as.POSIXlt(paste0(format(x, "%Y"),'-01-01')))
    temp <- unclass(j) 
    return(temp[1] + 1) 
  }
}

julian.days <- sapply(date, julian_conv)

结果:

print(julian.days)
[1] 346.5417 347.5417 346.5000

然后我取了那些起点的平均值。

mean <- mean(julian.days)
[1] 346.8611

现在我需要将平均起点转换回日期格式 (yyyy-dd-mm HH:MM:SS),一次用于平年,一次用于闰年。现在的问题是,这怎么可能?

您可以将儒略日添加到日期对象中:

2006 年是 365 天的一年:

as.Date("2006-01-01") + mean(julian.days)

2008 年是闰年:

as.Date("2008-01-01") + mean(julian.days)