我有一列日期格式: 16Jun10 我想提取朱利安日。 我有多年。
我已经尝试了julian和mdy.date这些函数,但它似乎没有用。
答案 0 :(得分:22)
尝试以下操作从类character
(即文本)转换为类POSIXlt
,然后提取朱利安日(yday
):
tmp <- as.POSIXlt("16Jun10", format = "%d%b%y")
tmp$yday
# [1] 166
有关功能设置的更多详细信息:
?POSIXlt
?DateTimeClasses
另一种选择是使用Date
类,然后使用format
提取朱利安日(请注意,此类定义1:366之间的julian日,而POSIXlt为0:365):
tmp <- as.Date("16Jun10", format = "%d%b%y")
format(tmp, "%j")
# [1] "167"
答案 1 :(得分:18)
类似地:
require(lubridate)
x = as.Date('2010-06-10')
yday(x)
[1] 161
另请注意,使用lubridate:
> dmy('16Jun10')
[1] "2010-06-16 UTC"
答案 2 :(得分:6)
您可以使用具有insol
功能的R JD(x, inverse=FALSE)
包将POSIXct转换为Julian Day Number(JDN)。
insol
包也有JDymd(year,month,day,hour=12,minute=0,sec=0)
。
要显示整个朱利安日期(JD),您可能需要设置options(digits=16)
。
答案 3 :(得分:3)
my.data = read.table(text = "
OBS MONTH1 DAY1 YEAR1
1 3 1 2012
2 3 31 2012
3 4 1 2012
4 4 30 2012
5 5 1 2012
6 5 31 2012
7 6 1 2012
8 6 30 2012
9 7 1 2012
10 7 31 2012
", header = TRUE, stringsAsFactors = FALSE)
my.data$MY.DATE1 <- do.call(paste, list(my.data$MONTH1, my.data$DAY1, my.data$YEAR1))
my.data$MY.DATE1 <- as.Date(my.data$MY.DATE1, format=c("%m %d %Y"))
my.data$my.julian.date <- as.numeric(format(my.data$MY.DATE1, "%j"))
my.data
退货,技术上不正确,因为Julian日期在每年一月的第一天不会返回1:
http://en.wikipedia.org/wiki/Julian_day
以下日期是序数日期:
OBS MONTH1 DAY1 YEAR1 MY.DATE1 my.julian.date
1 1 3 1 2012 2012-03-01 61
2 2 3 31 2012 2012-03-31 91
3 3 4 1 2012 2012-04-01 92
4 4 4 30 2012 2012-04-30 121
5 5 5 1 2012 2012-05-01 122
6 6 5 31 2012 2012-05-31 152
7 7 6 1 2012 2012-06-01 153
8 8 6 30 2012 2012-06-30 182
9 9 7 1 2012 2012-07-01 183
10 10 7 31 2012 2012-07-31 213
答案 4 :(得分:1)
这是我最初用APL编写并转换为J的R版本代码。我们称其为伪Julian,是因为它仅适用于1582年10月15日之后的日期,在西方世界的某些地区,这是日历改革的日期,任意更改日期。
#* toJulian: convert 3-element c(Y,M,D) timestamp into pseudo-Julian day number.
toJulian<- function(TS3)
{ mm<- TS3[2]
xx<- 0
if( mm<=2) {xx<- 1}
mm<- (12*xx)+mm
yy<- TS3[1]-xx
nc<- floor(0.01*yy)
jd<- floor(365.25*yy)+floor(30.6001*(1+mm))+TS3[3]+1720995+(2-(nc-floor(0.25*nc)))
return(jd)
#EG toJulian c(1959,5,24) -> 2436713
#EG toJulian c(1992,12,16) -> 2448973
}
这是反函数:
#* toGregorian: convert pseudo-Julian day number to timestamp in form c(Y,M,D)
# (>15 Oct 1582). Adapted from "Numerical Recipes in C" by Press,
# Teukolsky, et al.
toGregorian<- function(jdn)
{ igreg<- 2299161 # Gregorian calendar conversion day c(1582,10,15).
ja<- floor(jdn)
xx<- 0
if(igreg<=ja){xx<- 1}
jalpha<- floor((floor((xx*ja)-1867216)-0.25)/36524.25)
ja<- ((1-xx)*ja) + ((xx*ja)+1+jalpha-floor(0.25*jalpha))
jb<- ja+1524
jc<- floor(6680+((jb-2439870)-122.1)/365.25)
jd<- floor(365.25*jc)
je<- floor((jb-jd)/30.6001)
id<- floor((jb-jd)-floor(30.6001*je))
mm<- floor(je-1)
if(12<mm){mm<- mm-12}
iyyy<- floor(jc-4715)
if(mm>2){iyyy<- iyyy-1}
if(0>iyyy){iyyy<- iyyy-1}
gd<- c(iyyy, mm, id)
return(gd)
#EG toGregorian 2436713 -> c(1959,5,24)
#EG toGregorian 2448973 -> c(1992,12,16)
}