Python 年月日、儒略日、年纪日互转函数

import datetime

def time2mjd(dateT):
    t0=datetime.datetime(1858,11,17,0,0,0,0)#简化儒略日起始日
    mjd=(dateT-t0).days
    mjd_s=dateT.hour*3600.0+dateT.minute*60.0+dateT.second+dateT.microsecond/1000000.0
    return mjd+mjd_s/86400.0

#mjd转datetime类
def mjd2time(mjd):
    t0=datetime.datetime(1858,11,17,0,0,0,0)#简化儒略日起始日
    return t0+datetime.timedelta(days=mjd)

def date2doy(date):
    date_base = datetime.datetime(date.year,1,1)
    doy = (date - date_base).days + 1
    return doy