转换站心坐标(方位角,仰角),以赤道坐标(RA,十二月)与PyEphem

这看起来相当的基本任务,但我不能梳理出来。

该PyEphem文档:

http://rhodesmill.org/pyephem/radec.html

介绍了如何围绕执行转换的另一种方式,从BodyObserver对象,以视站心位置,与仰角和方位角在.alt.az属性。

不过我应该怎么,相反,从仰角和方位角启动,并得到RA和十二月?

例如,下面是一组为此我想获得RA和DEC中赤道坐标系坐标:

az = 3.30084818 #rad el = 0.94610742 #rad lat = 34.64 #deg lon = -103.7 #deg alt = 35800.26 #m ut = 2455822.20000367 #julian date

谢谢!

--------------解决方案-------------

这里有两个微妙之处。 首先,你已经发生了用“高度”和“高度”的意思是什么这两个术语指的是在PyEphem库相反 - 所以你而不是调用点在天空中的“高度/方位”位置的“高空/方位“的位置;第二,看来PyEphem忘记提供一个简单的方法来朱利安日期转换为自己的格式。虽然是一个功能julian_date()将去另一个方向上,我们将不得不做一点而努力自己去搞清楚什么其他方向ephem的名字是。

考虑到这些规定,我认为这个脚本可能会回答你的问题:

import ephem

az = 3.30084818 #rad
el = 0.94610742 #rad
lat = 34.64 #deg
lon = -103.7 #deg
alt = 35800.26 #m
ut = 2455822.20000367 #julian date

# Which Julian Date does Ephem start its own count at?
J0 = ephem.julian_date(0)

observer = ephem.Observer()
observer.lon = str(lon) # str() forces deg -> rad conversion
observer.lat = str(lat) # deg -> rad
observer.elevation = alt
observer.date = ut - J0

print observer.date
print observer.radec_of(az, el)

它生产的答案看起来是正确的这一特定的观察? 下面是脚本打印对我来说:

2011/9/17 16:48:00
(9:16:24.95, -0:45:56.8)

让我知道,如果让身体感觉这个特殊的观察,或者如果一个号码是错在这里,仍然需要进行调整!

分类:蟒蛇 时间:2015-03-15 人气:4
分享到:

相关文章

Copyright (C) 55228885.com, All Rights Reserved.

55228885 版权所有 京ICP备15002868号

processed in 0.267 (s). 10 q(s)