如何使用 PyEphem 计算正确的经度和行星星座答案

【问题标题】：How to calculate correct longitude and constellation of planets using PyEphem如何使用 PyEphem 计算正确的经度和行星星座
【发布时间】：2013-04-15 01:12:57
【问题描述】：

我尝试使用 PyEphem 计算太阳经纬度，但我不明白的是当我计算太阳经度和星座时

>>> sun = ephem.Sun()
>>> sun.compute('2011/05/04')
>>> print ephem.Ecliptic(sun).lon
43:02:58.8

所以 43:02:58.8 是 13.02 金牛座但是当我试图得到星座时

>>> print ephem.constellation(sun)
('Ari', 'Aries')

它说白羊座怎么了，我不明白请高人帮忙

【问题讨论】：

标签： python astronomy pyephem

【解决方案1】：

好问题！

有两个原因导致您得到的答案与您的预期不同。

首先，天文学家目前使用的星座边界是根据 1875 年黄道恰好所在的位置确定的。从那时起，黄道上的“零经度”点随着地轴的进动而不断移动，这样你得到的答案就会偏离几度。要真正了解 1875 坐标中的经度，您必须这样问：

import ephem
sun = ephem.Sun()
sun.compute('2011/05/04')
print "Today's coordinates:", ephem.Ecliptic(sun).lon
sun.compute('2011/05/04', epoch='1875')
print "In 1875 coordinates:", ephem.Ecliptic(sun).lon

运行此脚本时，我得到输出：

Today's coordinates: 43:02:58.8
In 1875 coordinates: 41:18:14.1

这样就可以从您的结果中消除一个错误来源。请注意，PyEphem constellation() 函数会在其星座边界表中查找天空位置之前自动将坐标转换为 1875，前提是您提供一个具有不同于 1875 的 epoch 的物体作为其参数。

您遇到的另一个问题是，官方星座边界标记了天空中实际的恒星群，这些恒星群形成了公羊（金牛座）、天秤座（天秤座）等的形象。碰巧的是，这些数字甚至没有接近每个横跨黄道的 30°（正如您可能已经预料到的，因为您从返回的角度减去 30°？）。如果你打印出沿着黄道的每个角度的星座，你会 (a) 发现星座的宽度变化很大，并且 (b) 你会发现有 13 个而不是 12 个星座穿过黄道，因为蛇夫座的一部分位于沿着 1875 年的黄道。

以下是一些可能有用的参考资料：

http://en.wikipedia.org/wiki/Constellation#IAU_constellations

http://en.wikipedia.org/wiki/Precession

【讨论】：

【解决方案2】：

在白羊座。这是我使用 Beta Ari 作为参考的计算：

>>> sun = e.Sun()
>>> sun.compute('2011/05/04')
>>> e.Ecliptic(sun).lon
43:02:58.8
>>> beta_ari = e.readdb('Sheratan,f|S|A5,01:54:38.5|98.74,28:48:28.9|-110.41,2.65,2000,0')
>>> beta_ari.compute('2011/05/04')
>>> e.Ecliptic(beta_ari).lon
36:53:55.2
>>> e.Ecliptic(sun).lon - e.Ecliptic(beta_ari).lon
0.10735523133236013
>>> e.degrees(e.Ecliptic(sun).lon - e.Ecliptic(beta_ari).lon)
6:09:03.6

黄道线上的太阳与 Beta Ari 相距 6:09 度，因此它仍在白羊座。我还使用该日期与其他程序仔细检查了这个结果，结果是一样的。

【讨论】：