问题描述
我正在尝试计算0年之前巴比伦的古代现象。而且我能够获得与Pyephem的已知记录一致的东西,但我想使用skyfield,因为这似乎更现代图书馆。但是,当我尝试将工作代码从Pyephem转换为skyfield时,会得到截然不同的结果。我不确定日期本身是否有问题,或者计算是否错误。
可以在以下位置找到to程序的快速Jupyter Notebook:https://gist.github.com/willismonroe/ae49480cd4cb1c21c5a214a70eb6f3d6
解决方法
粗略的信封计算表明,在29天之内绕月飞行的月球,每天运动约360°÷29≈12½度。
因此,您的Skyfield脚本存在一个问题:仅当月亮相对于太阳的经度在6°至14°之间(仅14 − 6 = 8°宽)时,才打印日期。在某些月份中,它的12½°跳跃恰好落在了狭窄的8°范围内,而在几个月中,它完全跳了过来。这就是为什么Skyfield脚本只打印几个月而不能打印几个月的原因。
相比之下,PyEphem脚本以弧度表示的角度乘以30,将范围0…6.28变为范围0…188.4。粗略地说,您使用的度量单位是“双学位”,其中188个为整圆。月球每天仅跳跃约6¼“双度”,因此保证每月至少一次降落在6至14范围内,因为6¼的跳跃不足以完全跨越该范围。>
我的猜测是,您应该尝试选择一种您感兴趣的精确度角(例如6°)的方法,然后编写一些例程以准确找到每个月超过该阈值的日期和时间。这是一个希望可以帮助您入门的示例: