儒略日 儒略历

儒略历(Julian Calendar)、儒略日((julian Day))与儒略年这三个词字面上我并不陌生,但对于具体内涵却不知就里。

   最近无意检索发现(因研究生态水文模型IHA,里面采用了儒略日计算),原来此儒略非彼儒略,一为威风八面的君主,另一为子报父恩的学者。

   儒略历是格里历的前身,由罗马共和国独裁官儒略·恺撒(Julius Caesar)采纳埃及亚历山大希腊数学家天文学家索西琴尼计算的历法,在公元前45年1月1日起执行,取代旧罗马历法的一种历法。一年设12个月,大小月交替,四年一闰,平年365日,闰年于二月底增加一闰日,年平均长度为365.25日。(来自维基百科)

   儒略历一直用到大约 1582 年,这时各国开始使用罗马历法。在儒略历里面, 一年是近似 365 1/4 天 = 365.25 天。 这样大约每 128 年就有一天的误差。由于累积误差随着时间越来越大,促使教皇格里高利十三世(Gregory XIII)按照与弥撒议会 ( Council of Trent)一致的精神改革了历法,变为格里历,即沿用至今的西历。在罗马历法里,一年是近似 365 + 97 / 400 天 = 365.2425 天。因此对应于罗马历法,大约要 3300 年,才会积累一天的误差。(本段根据维基百科与其他来源修改)

儒略日(Julian day)是指由公元前4713年1月1日,协调世界时中午12时开始所经过的天数,多为天文学家采用,用以作为天文学的单一历法,把不同历法的年表统一起来。 

   儒略日是一种不用年月的长期纪日法,简写为JD。是由法国学者Joseph Justus Scliger(1540–1609)在1583年所创,这名称据说是为了纪念他的父亲——意大利学者Julius Caesar Scaliger(1484–1558)。 

   儒略日的起点订在公元前4713年(天文学上记为 -4712 年)1月1日格林威治时间平午(世界时12:00),即JD 0 指定为 4713 B.C. 1月1日12:00 UT到4713 B.C. 1月2日12:00 UT的24小时。每一天赋予了一个唯一的数字,顺数而下,如:1996年1月1日12:00:00的儒略日是2450084。这个日期是考虑了太阳、月亮的运行周期,以及当时收税的间隔而订出来的。Joseph Scliger定义儒略周期为7980年,是因28、19、15的最小公倍数为28×19×15=7980。其中: 

   28年为一太阳周期(solar cycle),经过一太阳周期,则星期的日序与月的日序会重复。 

   19年为一太阴周期,或称默冬章(Metonic cycle),因235朔望月=19回归年,经过一太阴周期则阴历月年的日序重复。 

  15年为一小纪(indiction cycle),此为罗马皇帝君士坦丁(Constantine)所颁,每15年评定财产价值以供课税,成为古罗马用的一个纪元单位, 

   故以7980年为一儒略周期,而所选的起点公元前4713年,则是这三个循环周期同时开始的最近年份。 

   以儒略日计日是为方便计算年代相隔久远或不同历法的两事件所间隔的日数。

儒略日(Julian Day)是一种表示日期和时间的系统,它将所有日期和时间统一为一个连续的计数。儒略日的起点是公元前4713年1月1日中午(以儒略历计)。儒略日转换是将公历日期转换为儒略日计数的过程,反之亦然。在Python中,可以使用第三方库如`ephem`或者内置模块如`datetime`来实现儒略日的转换。 使用Python内置的`datetime`模块进行转换时,你需要先计算从儒略日起点到目标日期的天数差,然后加上儒略日的起点值(2440587.5,因为起点是中午时分)。以下是一个简单的例子: ```python import datetime # 设定儒略日起点 JD_START = 2440587.5 # 给定一个公历日期(datetime对象) date = datetime.datetime(2023, 1, 1) # 计算自儒略日起点以来的天数(考虑到中午12点) jd = JD_START + date.toordinal() + (date.hour + date.minute / 60 + date.second / 3600) / 24 print("儒略日为:", jd) ``` 要将儒略日转换回公历日期,可以按照儒略日起点逆向计算: ```python from datetime import datetime, timedelta # 设定儒略日起点 JD_START = 2440587.5 # 给定儒略日数值 jd = 2459602.5 # 示例儒略日 # 转换回日期时间 # 注意:由于jd包含了时间部分,需要先减去起点的小数部分以获取整天数 jd -= JD_START ordinal_date = int(jd) # 转换为整数部分(整天数) day_fraction = jd - ordinal_date # 小数部分(一天中的时间) # 计算公历日期 date = datetime.fromordinal(ordinal_date) + timedelta(days=day_fraction) print("公历日期为:", date) ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值