python - PyEphem:RA 和分离问题

Question

我正在尝试使用 PyEphem 库进行天文学，并尝试重新创建/理解其中的一些基本计算。

有一个称为分离的函数，给定两个行星和日期/时间，它可以计算这两个行星之间相对于它们在平面上的“x”投影的分离(角度)。

因此，对于 2018 年 1 月 1 日的水星和火星，我们有:

import ephem
import math 

mercury = ephem.Mercury('2018/1/1')
mars = ephem.Mars('2018/1/1')

s1 = ephem.separation(mercury, mars)
print(math.degrees(s1))

哪个返回

33.792384499568264

但是如果我想在没有分离函数的情况下计算这个，那么计算就会像 1 的“赤经”减去另一个的“赤经”一样简单:

math.degrees(mercury.ra) - math.degrees(mars.ra)

哪个返回

35.114532008671574

为什么角度不同？根据 PyEphem 的说法，由于我没有包含观察者的纬度和经度，因此所有计算都应该以地为中心。

是否有人熟悉 PyEphem 或另一个具有内置星历的库背后进行的计算，可以产生一致的分离结果？

Answer 1

好问题!您需要查阅有关球面三角学的教科书。正如您无法仅根据两个地点的经度来预测它们之间的英里数一样，因为它们可能都位于赤道上且相距较远，或者都位于北极附近且彼此相距仅几英尺，因此您也无法预测它们之间的距离。两个物体之间的距离仅来自赤经(=“天空经度”)，但必须添加(a)两个物体的赤纬，以及(b)使用特殊公式来计算球体上两个位置之间的弧度沿着一个大圆路线。