def test_rho_roundtrip():
a = 10 * u.arcsec
eph = {'delta': 1 * u.au}
b = rho_as_angle(rho_as_length(a, eph), eph)
assert np.isclose(a.value, b.to(u.arcsec).value)
微信
公众号
扫码关注公众号