Konwertuj rok / miesiąc / dzień na dzień roku w Pythonie

Korzystam z datetimemodułu tj .:

>>> import datetime
>>> today = datetime.datetime.now()
>>> print(today)
2009-03-06 13:24:58.857946

i chciałbym obliczyć dzień w roku uwzględniający lata przestępne. np. dzisiaj (6 marca 2009) jest 65 dzień 2009.

Widzę dwie opcje:

1. Utwórz number_of_days_in_month = [31, 28, ...]tablicę, zdecyduj, czy jest to rok przestępny i ręcznie zsumuj dni.

2. Użyj, datetime.timedeltaaby zgadnąć, a następnie wyszukać binarny właściwy dzień roku:

   >>> import datetime
   >>> YEAR = 2009
   >>> DAY_OF_YEAR = 62
   >>> d = datetime.date(YEAR, 1, 1) + datetime.timedelta(DAY_OF_YEAR - 1)
   

Oba wydają się dość niezgrabne i mam przeczucie, że istnieje bardziej „Pythonowy” sposób obliczania dnia w roku. Jakieś pomysły / sugestie?

Użyj, datetime.timetuple()aby przekonwertować datetimeobiekt na time.struct_timeobiekt, a następnie pobierz jego tm_ydaywłaściwość:

from datetime import datetime
day_of_year = datetime.now().timetuple().tm_yday  # returns 1 for January 1st

Możesz użyć strftimez %jciągiem formatu :

>>> import datetime
>>> today = datetime.datetime.now()
>>> today.strftime('%j')
'065'

ale jeśli chcesz robić porównania lub obliczenia z tą liczbą, musisz ją przekonwertować na, int()ponieważ strftime()zwraca ciąg. W takim przypadku lepiej jest skorzystać z odpowiedzi DzinX .

Odpowiedź DZinX jest świetną odpowiedzią na to pytanie. Znalazłem to pytanie i użyłem odpowiedzi DZinX, szukając funkcji odwrotnej: zamień daty z juliańskim dniem roku na czasy dat.

Znalazłem to do pracy:

import datetime
datetime.datetime.strptime('1936-077T13:14:15','%Y-%jT%H:%M:%S')

>>>> datetime.datetime(1936, 3, 17, 13, 14, 15)

datetime.datetime.strptime('1936-077T13:14:15','%Y-%jT%H:%M:%S').timetuple().tm_yday

>>>> 77

Lub liczbowo:

import datetime
year,julian = [1936,77]
datetime.datetime(year, 1, 1)+datetime.timedelta(days=julian -1)

>>>> datetime.datetime(1936, 3, 17, 0, 0)

Lub z ułamkowymi datami Jdates opartymi na 1, popularnymi w niektórych domenach:

jdate_frac = (datetime.datetime(1936, 3, 17, 13, 14, 15)-datetime.datetime(1936, 1, 1)).total_seconds()/86400+1
display(jdate_frac)

>>>> 77.5515625

year,julian = [1936,jdate_frac]
display(datetime.datetime(year, 1, 1)+datetime.timedelta(days=julian -1))

>>>> datetime.datetime(1936, 3, 17, 13, 14, 15)

Nie jestem pewien co do etykiety, ale pomyślałem, że wskaźnik do odwrotnej funkcjonalności może być przydatny dla innych takich jak ja.

Chcę zaprezentować wydajność różnych podejść, w Pythonie 3.4, Linux x64. Fragment z profilera linii:

      Line #      Hits         Time  Per Hit   % Time  Line Contents
      ==============================================================
         (...)
         823      1508        11334      7.5     41.6          yday = int(period_end.strftime('%j'))
         824      1508         2492      1.7      9.1          yday = period_end.toordinal() - date(period_end.year, 1, 1).toordinal() + 1
         825      1508         1852      1.2      6.8          yday = (period_end - date(period_end.year, 1, 1)).days + 1
         826      1508         5078      3.4     18.6          yday = period_end.timetuple().tm_yday
         (...)

Więc najbardziej wydajne jest

yday = (period_end - date(period_end.year, 1, 1)).days

Jeśli masz powody, aby unikać korzystania z datetimemodułu, te funkcje będą działać.

def is_leap_year(year):
    """ if year is a leap year return True
        else return False """
    if year % 100 == 0:
        return year % 400 == 0
    return year % 4 == 0

def doy(Y,M,D):
    """ given year, month, day return day of year
        Astronomical Algorithms, Jean Meeus, 2d ed, 1998, chap 7 """
    if is_leap_year(Y):
        K = 1
    else:
        K = 2
    N = int((275 * M) / 9.0) - K * int((M + 9) / 12.0) + D - 30
    return N

def ymd(Y,N):
    """ given year = Y and day of year = N, return year, month, day
        Astronomical Algorithms, Jean Meeus, 2d ed, 1998, chap 7 """    
    if is_leap_year(Y):
        K = 1
    else:
        K = 2
    M = int((9 * (K + N)) / 275.0 + 0.98)
    if N < 32:
        M = 1
    D = N - int((275 * M) / 9.0) + K * int((M + 9) / 12.0) + 30
    return Y, M, D

Po prostu odejmij 1 stycznia od daty:

import datetime
today = datetime.datetime.now()
day_of_year = (today - datetime.datetime(today.year, 1, 1)).days + 1