在 Python 中转换 Year/Month/Day to Day of Year

我使用的是 datetime模块,例如:

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

我想计算一年中考虑闰年的那一天。例如今天(2009年3月6日)是2009年的第65天。

我认为有两个选择:

  1. 创建一个 number_of_days_in_month = [31, 28, ...]数组,确定它是否是闰年,然后手动计算日期。

  2. 使用 datetime.timedelta进行猜测,然后进行二进制搜索,找出一年中正确的日期:

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

这两种方法都让人感觉很笨拙,而且我有一种直觉,那就是有一种更“ Python”的方法来计算一年中的某一天。有什么想法或建议吗?

209260 次浏览

你可以使用 strftime%j格式字符串:

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

但是如果您希望对这个数字进行比较或计算,则必须将其转换为 int(),因为 strftime()返回一个字符串。如果是这种情况,你最好使用 DzinX’s答案。

只要从日期中减去1月1日:

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

使用 datetime.timetuple()将您的 datetime对象转换为 time.struct_time对象,然后获取它的 tm_yday属性:

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

DZinX 的回答 是这个问题的绝佳答案。我发现了这个问题,在寻找反函数的时候使用了 DZinX 的答案: 将日期与朱利安年份转换为日期时间。

我发现这个有效:

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

或者从数字上来说:

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


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

或者使用在某些领域流行的基于分数1的 jdate:

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)

我不太懂这里的礼节,但我想一个指向反功能的指针可能对我这样的人有用。

我想在 Python 3.4,Linux x64上展示不同方法的性能:

      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
(...)

所以效率最高的是

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

如果您有理由避免使用 datetime模块,那么这些函数将起作用。

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

您可以简单地使用由“熊猫”提供的 dayofyear属性,该属性反过来为您提供特定年份的一天。 举例来说。

data["day_of_year"] = data.Datetime.apply(lambda x:x.dayofyear)