更新: Skyfield 刚刚进行了重大修订,包括扩展文档和角度分离方法,请参见所接受的答案。
我正在使用Skyfield计算两个物体之间的视角分离。在该软件包中没有找到相应的方法,因此我通过计算两个视位置向量之间的点积“发明”了一种方法。
当前这种方法是最好的方法吗?在 Skyfield 的范围内基本正确吗?
我正在使用Skyfield计算两个物体之间的视角分离。在该软件包中没有找到相应的方法,因此我通过计算两个视位置向量之间的点积“发明”了一种方法。
当前这种方法是最好的方法吗?在 Skyfield 的范围内基本正确吗?
def separation(seconds, lat, lon):
lat, lon, seconds = float(lat), float(lon), float(seconds) # necessary it seems
place = earth.topos(lat, lon)
jd = JulianDate(utc=(2016, 3, 9, 0, 0, seconds))
mpos = place.at(jd).observe(moon).apparent().position.km
spos = place.at(jd).observe(sun).apparent().position.km
mlen = np.sqrt((mpos**2).sum())
slen = np.sqrt((spos**2).sum())
sepa = ((3600.*180./np.pi) *
np.arccos(np.dot(mpos, spos)/(mlen*slen)))
return sepa
from skyfield.api import load, now, JulianDate
import numpy as np
from scipy.optimize import minimize
data = load('de421.bsp')
sun = data['sun']
earth = data['earth']
moon = data['moon']
sep = separation(12000, 32.5, 215.1)
print "sun-moon aparent separation: ", sep, " arcsec"