我有卫星经度,纬度和海拔高度。然后我有Ra& Decl和远处的物体。我需要知道在给定时间物体和月亮/太阳之间的角距离。
import ephem
sun = ephem.Sun()
scraft = ephem.Observer()
scraft.lon = lon
scraft.lat = lat
scraft.elevation = altit
scraft.date = time-15019.5 # my time is in MJD
sun.compute(scraft)
print float(sun.dec), float(sun.ra)
这是正确的做法吗? 有没有办法检查结果? (另一个带有Web界面的工具,我可以输入少数情况下的数字并查看是否存在差异)
答案 0 :(得分:2)
你可能想要创建一个固定的物体来表示远处物体的位置,然后要求从卫星观察它的位置与太阳的位置之间的分离:
f = ephem.FixedBody()
f._ra = '1:23:45.0'
f._dec = '6:78:90.0'
f.compute(scraft.date)
print ephem.separation(sun, f), 'degrees'
print float(ephem.separation(sun, f)), 'radians'