Pyephem和Skyfield在计算古代现象方面的区别

如何解决Pyephem和Skyfield在计算古代现象方面的区别

我正在尝试计算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°）的方法，然后编写一些例程以准确找到每个月超过该阈值的日期和时间。这是一个希望可以帮助您入门的示例：

https://rhodesmill.org/skyfield/searches.html