微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

Astropy 和 JPL 的 Horizo​​ns 查询有不同的输出

如何解决Astropy 和 JPL 的 Horizo​​ns 查询有不同的输出

我正在尝试将地球在 Astropy 给出的日心系统中的位置与 jpl 的 Horizo​​ns 的位置进行比较。我发现了显着的差异,但我不知道它们可能是由于什么......

我想这样做是因为我的目的是将地理坐标转换为日心坐标。

    from astropy.coordinates import SkyCoord
    from astropy.time import Time
    from astropy import units as u
    from astroquery.jplhorizons import Horizons
    
    time = '2015-03-30T21:33:52.000'
    JD = Time(time).jd
    
    observing_time = Time(time,format='isot',scale='utc')
    
    #ASTROPY
    c = SkyCoord(x=0,y=0,z=0,unit='au',representation_type='cartesian',frame='gcrs',obstime=time)
    
    cc = c.transform_to("heliocentriceclipticiau76")
    print(cc.cartesian.x,cc.cartesian.y,cc.cartesian.z)
    
    cc = c.transform_to("heliocentricmeanecliptic")
    print(cc.cartesian.x,cc.cartesian.z)
    
    cc = c.transform_to("heliocentrictrueecliptic")
    print(cc.cartesian.x,cc.cartesian.z)
    
    #jpl'S HORIZONS
    obj = Horizons(id="Geocenter",location="@sun",epochs=JD,id_type='id').vectors()
    print(obj['x'],obj['y'],obj['z'])

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。