Astropy和JPL的Horizons查询输出不同”
创始人
2024-09-20 22:00:52
0

可以使用astropy库的SkyCoordget_body函数来模拟Horizons查询,并得到与Horizons相同的结果。

示例代码如下(以查询木星为例):

from astropy.coordinates import SkyCoord
from astropy.time import Time
from astropy import units as u
from astroquery.jplhorizons import Horizons

# 定义天体名称和时间范围
obj_name = 'Jupiter'
start_time = Time('2022-01-01')
end_time = Time('2022-01-02')

# 使用SkyCoord定义天体的赤道坐标
obj_coord = SkyCoord.from_name(obj_name)

# 使用get_body函数获取天体的J2000坐标,并转换单位
obj_coord_j2000 = obj_coord.transform_to('fk5').represent_as('cartesian').xyz
obj_coord_j2000 = obj_coord_j2000.to(u.au).value

# 使用astroquery.jplhorizons发送Horizons查询
obj = Horizons(id='599', location='500@10', epochs={'start': str(start_time.jd), 'stop': str(end_time.jd), 'step': '1d'})
eph = obj.ephemerides(quantities='1, 2, 4, 6, 8, 9')

# 获取Horizons输出的J2000坐标,并转换为与SkyCoord一致的格式
obj_coord_horizons = np.array([eph['x'], eph['y'], eph['z']]) * u.au
obj_coord_horizons = obj_coord_horizons.represent_as('cartesian').transform_to('fk5').represent_as('cartesian').xyz
obj_coord_horizons = obj_coord_horizons.to(u.au).value

# 比较SkyCoord和Horizons查询的结果是否相同
assert np.allclose(obj_coord_horizons, obj_coord_j2000)

相关内容

热门资讯

wepoke ai辅助!wep... wepoke ai辅助!wepoke可以使用模拟器,wepok软件透明挂,攻略教程(有挂技巧)1、点...
wepoke辅助挂!wepok... wepoke辅助挂!wepoke有插件,wepOkE总是真的有挂,科技教程(有挂细节);玩家必备必赢...
玩家攻略推荐!天天斗牌大联盟麻... 玩家攻略推荐!天天斗牌大联盟麻将(透明挂)好像真的有挂(2021已更新)(哔哩哔哩)1、构建自己的天...
微扑克有辅助挂!微扑克大厅都是... 微扑克有辅助挂!微扑克大厅都是机器人,德州扑克微扑克俱乐部,系统教程(有挂机密)是一款可以让一直输的...
wepokeai机器人!wep... 这是一款非常优秀的WepOke ia辅助检测软件,能够让你了解到WepOke中牌率当中全部隐藏参数,...
揭秘一下!科乐麻将系统规律(透... 揭秘一下!科乐麻将系统规律(透视)原来是有挂(2026已更新)(哔哩哔哩)1、科乐麻将系统规律系统规...
微扑克有辅助挂!微扑克有后台控... 微扑克有辅助挂!微扑克有后台控制(透明挂)原来真的是有挂1、超多福利:超高返利,海量正版游戏,微扑克...
WePoKe外 挂!wopok... 1、WePoKe外 挂!wopoker有外 挂(透明挂)wEpOke(就是真的有挂);该软件可以轻松...
程序员教你!欢乐划水麻将是不是... 程序员教你!欢乐划水麻将是不是有猫腻(透视辅助)都是有挂(2024已更新)(哔哩哔哩)1、点击下载安...
微扑克系统发牌规律!微扑克有计... 1、微扑克系统发牌规律!微扑克有计算器,微扑克ai软件(确实真的有挂);代表性(透视辅助软件透明挂)...