Python PyEphem 计算方位角和高度

我是 PyEphem 的新手，我正在尝试弄清楚它可以做什么以及它是如何工作的。由于我不想将其用作黑匣子并盲目相信我得到的任何数字，因此我想重新创建一个解释的示例here http://www.stargazing.net/kepler/altaz.html.

该示例计算给定观察者在 1998 年 8 月 10 日 23:10 UT 的物体的方位角和高度。给出以下参数：

RA = 16 小时 41.7 分钟，DEC = 36 天 28 分钟

观察者的纬度为北纬 52 d 30 分，西经 1 d 55 分。

根据示例（我可以在 Excel 中重新创建），正确答案是 AZ = 269.14634 度，ALT = 49.169122 度。

我使用 pyephem 编写了以下代码来尝试达到相同的结果：

day = '1998/8/10 23:10:00'
longitude = ephem.degrees('-1.91667')
latitude = ephem.degrees('52.5')

star = ephem.FixedBody()
star._ra = '16:41:42.0'
star._dec = '36:28:00.0'

observer = ephem.Observer()
observer.date = day
observer.lon = longitude
observer.lat = latitude

star.compute(observer)

print 'Observer', observer
print 'RA', star.ra, 'DEC', star.dec
print 'AZ', star.az, 'ALT', star.alt

运行程序给我这个输出：

>>> 
Observer <ephem.Observer date='1998/8/10 23:10:00' epoch='2000/1/1 12:00:00' lon=-1:55:00.0 lat=52:30:00.0 elevation=0.0m horizon=0:00:00.0 temp=15.0C pressure=1010.0mBar>
RA 16:41:39.23 DEC 36:28:33.5
AZ 269:09:54.9 ALT 49:10:57.7

AZ + ALT 的结果显然与示例大致相符，但相差甚远。我还对打印输出中的 RA 和 DEC 与我输入的内容相比略有修改这一事实感到困惑。

如果有人可以帮助我阐明为什么结果不同以及我可以或应该做什么来复制结果，我将不胜感激。谢谢。

EDIT：更正了下面答案中指出的拼写错误。这个问题仍然有效。

EDIT2：好的，我已经阅读（并且有点理解）为什么 PyEphem 从此调整赤经和赤纬link http://rhodesmill.org/pyephem/radec.html。我不明白的是是否有任何方法可以让 PyEphem 忽略调整相对论偏转, nutation and 光像差就像你可以让它忽略一样大气折射？我假设方位角的差异是由于 RA 和 DEC 的调整造成的，但最好确认一下。

PyEphem 底层的 C 库没有任何方法来关闭偏转、像差或章动——也许是因为大自然也不让我们关闭这些效果，但我不确定！我要指出的是，它不会对地球轨道卫星进行这些计算，但我想不出一种简单的方法可以让您将卫星放置在您所在位置上方精确的赤经和赤经，以便您可以向 PyEphem 询问它的位置。

实际上，我这周将在 DjangoCon 上讨论 API，并思考 PyEphem 有朝一日如何使其内部工作更容易从 Python 访问，而不是将所有这些有趣的步骤锁定在 C 代码中。但是，在我准备好替代方案之前，完成您想要的任务的唯一方法是打开源文件circum.c并注释掉这些行：

/* allow for relativistic light bending near the sun */
deflect (mjed, lam, bet, lsn, rsn, 1e10, &ra, &dec);

/* TODO: correction for annual parallax would go here */

/* correct EOD equatoreal for nutation/aberation to form apparent 
 * geocentric
 */
nut_eq(mjed, &ra, &dec);
ab_eq(mjed, lsn, &ra, &dec);

如果那三个人打电话——deflect(), nut_eq(), and ab_eq()从第 263 行开始 - 被注释掉，那么您可能会得到更接近本文中给出的答案。您可以通过以下方式应用这些更改：

• 正在下载.tar.gz对于 PyEphem。
• 提取存档以生成文件。
• 进行我上面建议的编辑。
• Run python setup.py install安装您的自定义版本的软件。
• 试试看！

如果进动以某种方式发挥作用，可能还有另一个障碍 - 您可能需要设置您的纪元star完全反对'1998/8/10 23:10:00'在这种情况下 - 但请先尝试关闭三个光效调用，看看这是否能让您更接近！