我是新来的皮耶芬,我正试图弄清楚它能做什么以及它是如何工作的。因为我不想把它当作黑匣子,盲目地相信我得到的任何数字,所以我想重新创建一个解释为here的示例。在
这个例子计算1998年8月10日UT 23:10时给定观测者的目标方位角和高度。给出了以下参数:
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
运行该程序可以得到以下输出:
^{pr2}$
AZ+ALT的结果显然与示例大致相同,但相差甚远。我还感到困惑的是,与我输入的内容相比,RA和DEC在打印输出中略有修改。在
如果有人能帮我解释一下为什么结果不同,以及我可以或应该做些什么来复制结果,我将不胜感激。谢谢。在
编辑:更正了下面答案中指出的一个错误。这个问题仍然有效。在
好吧,我已经读过(并且有点理解)为什么赤经和赤纬是由PyEphem从这个link来调整的。我不明白的是,如果有什么方法可以让PyEphem忽略对相对论性偏转、章动和光像差的调整,就像你可以让它忽略大气折射一样?我假设方位角的差异是由于RA和DEC的调整,但最好能确认一下。在