我是PyEphem的新手,我正在试图弄清楚它能做什么以及它是如何工作的.因为我不想把它用作黑盒子而且盲目地相信我得到的任何数字,我想重新创建一个解释为
here的例子.
该示例在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
RA 16:41:39.23 DEC 36:28:33.5
AZ 269:09:54.9 ALT 49:10:57.7
对于这个例子来说,AZ ALT的结果显然很棒,但远非相同.令我感到困惑的是RA和DEC在打印输出中略有修改,与我输入的相比.
如果有人能帮助我阐明为什么结果不同以及我能够或应该做些什么来复制结果,我将非常感激.谢谢.
编辑:纠正了下面的答案中指出的拼写错误.这个问题仍然有效.
编辑2:好的,我已经阅读(并且有点理解)为什么PyEphem从link开始调整正确的提升和赤纬.我不明白的是,是否有任何办法让PyEphem忽略调整相对论偏转,章动和光的像差与你可以忽略大气折射的方式相同吗?我假设方位角的差异是由于RA和DEC的调整,但确认很好.