Konvertieren ein Datum-Vektor in julianischen Tag in R

Ich habe eine Spalte mit Datum im format:
16Jun10
und ich möchte zu extrahieren, die den julianischen Tag.
Ich habe verschiedenen Jahre.

Habe ich versucht die Funktionen julian und mdy.Datum, und es scheint nicht zu funktionieren.

InformationsquelleAutor user3166363 | 2014-01-28
  1. 22

    Versuchen Sie Folgendes zum konvertieren von Klasse character(d.h. text -) Klasse POSIXlt, und extrahieren Sie dann julianischen Tag (yday):

    tmp <- as.POSIXlt("16Jun10", format = "%d%b%y")
tmp$yday
# [1] 166

    Für mehr details über die Funktion Einstellungen:

    ?POSIXlt
?DateTimeClasses

    Andere Möglichkeit ist die Verwendung eines Date - Klasse, und verwenden Sie dann format zum extrahieren eines julianischen Tag (beachten Sie, dass diese Klasse definieren, die julianischen Tage zwischen 1:366, während POSIXlt 0:365):

    tmp <- as.Date("16Jun10", format = "%d%b%y")
format(tmp, "%j")
# [1] "167"
    • Diese top-Antwort spricht über Tag des Jahres/Ordnungszahl Datum und nicht über Julianische Datum. Julian Termine sind Tag zählt von 4713 v. Chr. (!!). Warum ist diese Antwort so viele Male von Ihnen positiv bewertet werden?.
    • Siehe Terminologie Unterabschnitt - "Der Begriff Julianisches Datum kann sich auch darauf beziehen, außerhalb der Astronomie, der Tag-des-Jahres-Anzahl (mehr richtig, die Ordnungszahl Datum) im Gregorianischen Kalender, vor allem in der computer-Programmierung...". Ich denke, die Programmierer haben bis-stimmten mehr als die Astronomen haben nach unten gestimmt.
    • Eine gute :). Hinweis: Andrei Antwort unten. Tatsächlich berechnet Astronomie-bezogenen Julian Termine. Beachten Sie auch, dass der tag julian bezieht sich auf die Astronomie bezogene Bedeutung. Und jetzt ich ergebe mich..
    InformationsquelleAutor Marc in the box
  2. 18

    Ähnlich:

    require(lubridate)
x = as.Date('2010-06-10')
yday(x)

[1] 161

    Beachten Sie auch, mit lubridate:

    > dmy('16Jun10')
[1] "2010-06-16 UTC"
    • Könnten Sie mir bitte helfen, mit immer ein Julianisches Datum, wie 112161 für 2012/06/09 (yyyy/mm/dd)?
    InformationsquelleAutor John
  3. 6

    Können Sie R insol - Paket hat eine JD(x, inverse=FALSE) Funktion, die konvertiert POSIXct zu Julian Day Number (JDN).

    insol Paket hat auch JDymd(year,month,day,hour=12,minute=0,sec=0) für benutzerdefinierte Termine.

    Zur Anzeige des ganzen Julianische Datum (JD), die Sie ggf. einstellen müssen options(digits=16).

    InformationsquelleAutor Andrei
  4. 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

    Gibt, das ist technisch falsch, denn Julian Termine nicht zurück zu 1 auf den ersten Tag eines jeden Januar:

    http://en.wikipedia.org/wiki/Julian_day

    Die Daten sind Ordinal Termine:

       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
    InformationsquelleAutor Mark Miller
  5. 1

    Hier sind meine R-Versionen der code ursprünglich geschrieben in APL und konvertiert J. Wir nennen diesen pseudo-Julian, denn es soll nur für Datumsangaben nach dem 15 Oktober 1582 was ist, wenn die Kalender-reform, die in einigen teilen der westlichen Welt, willkürlich verändert das Datum.

    #* 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
}


    Hier ist die inverse Funktion:

    #* 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)
}
    InformationsquelleAutor DevonMcC
Schreibe einen Kommentar