python-3.x - 这是如何在 Python 3 下使用 PyEphem 计算格林威治时角？

Question

一直在用python3和pyephem研究天文导航 主要是通过一些必须由以下人员完成的计算 无论如何都要交。
对于视力降低部分，我正在尝试使用 pyephem 输出数据 堪比航海年历，主要是格林威治时角和赤纬。 为了获得白羊座第一点的时角，我尝试添加一个 body 使用 ephem.readdb 在 0.0 RA 和 0.0 dec。但在进一步阅读文档后，我认为 这是工作。

#!/bin/python3

import math , ephem
#  zero longitude
gmt_long = '0:0:0'          #  deg, min, sec
myloc_date = ( '2012/12/02 22:00:00' )
# observer for greenwich  gst
utcz = ephem.Observer()
utcz.date = myloc_date
utcz.long = gmt_long
print ( utcz.date )
print (" gst hours", utcz.sidereal_time() )
print (" gst   deg", ephem.degrees( utcz.sidereal_time() ) )

将输出的最后一行与航海年历中的白羊座列进行比较 日期

航海年鉴在线版 www.tecepe.com.br

Answer 1

是的，这看起来确实是正确的方法；你的每一步都是有意义的，所以你的结果和年历之间的一致性是有意义的。对于任何其他想要进行比较的人，这里是您引用的年鉴中的特定页面:

http://www.tecepe.com.br/scripts/AlmanacPagesISAPI.dll/pages?date=12%2F02%2F2012