4

我是 PyEphem 的新手，我正试图弄清楚它可以做什么以及它是如何工作的。由于我不想将其用作黑匣子并盲目相信我得到的任何数字，我想重新创建一个示例，此处解释。

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

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

观察者的纬度为北纬 52 分 30 分，西经 1 分 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 略有修改，这也让我感到困惑。

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

编辑：更正了下面答案中指出的错字。这个问题仍然有效。

EDIT2：好的，我已经阅读（并且有点理解）为什么 PyEphem 从此链接调整了正确的提升和赤纬。我不明白的是，是否有任何方法可以让 PyEphem 忽略对相对论偏转、章动和光像差的调整，就像你可以让它忽略大气折射一样？我假设方位角的差异是由于 RA 和 DEC 的调整，但很高兴确认。

4

3 回答 3

2

观察者包括“temp=15.0C pressure=1010.0mBar”这一事实意味着计算将包括折射。你想关闭折射

如帮助中所述：

这些表观位置包括模拟观测者温度和压力的大气折射的调整；将观察者的压力设置为零以忽略折射。

于 2012-09-04T15:38:49.660 回答
2

PyEphem 底层的 C 库没有任何方法可以关闭偏转、畸变或章动——也许是因为大自然也不允许我们关闭这些效果，但我不确定！我会注意到，它不会对地球轨道卫星进行这些计算，但我想不出一个简单的方法让你把卫星放在一个精确的 RA 和你的位置之上，这样你就可以问 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()从第 263 行开始）被注释掉，那么您可能会得到更接近本文中产生的答案nut_eq()。ab_eq()您可以通过以下方式应用这些更改：

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

可能还有另一个障碍，如果岁差以某种方式开始发挥作用——在这种情况下，你可能需要将star对象的纪元设置为准确'1998/8/10 23:10:00'的——但请先尝试关闭三个光效调用，看看这是否能让你更接近！

于 2012-09-05T11:32:31.520 回答
1

您想输入：RA = 16 h 41.7 min

但你输入了：star._ra = '16:41.42.0'

而不是star._ra = '16:41:42.0'

于 2012-09-04T15:13:16.243 回答