Я думаю, что вы можете сделать что-то вроде следующего:
from astropy import coordinates as coord
from astropy import units as u
from astropy.time import Time
now = Time('2018-03-14 23:48:00')
xyz = [-6340.40130292, 3070.61774516, 684.52263588]
vxvyvz = [-10.90, 56.4, -74.6]
# put velocities into cartesian differential
cartdiff = coord.CartesianDifferential(*vxvyvz, unit='km/s')
cartrep = coord.CartesianRepresentation(*xyz, unit=u.km, differentials=cartdiff)
gcrs = coord.ITRS(cartrep, obstime=now)
itrs = gcrs.transform_to(coord.GCRS(obstime=now))
# print position
print(itrs.cartesian.xyz)
# print velocity
print(itrs.cartesian.differentials)
Но я не совсем уверен, что он делает то, что вам нужно, или нет. Альтернативно, в astropy v. 3.0.1 ITRS класс, похоже, может принимать значения скорости, поэтому вы можете использовать
now = Time('2018-03-14 23:48:00')
pos = [-6340.40130292, 3070.61774516, 684.52263588]*u.km
vel = [-10.90, 56.4, -74.6]*u.km/u.s
gcrs = coord.ITRS(x=pos[0], y=pos[1], z=pos[2], v_x=vel[0], v_y=vel[1], v_z=vel[2], representation_type='cartesian', differential_type='cartesian', obstime=now)
itrs = gcrs.transform_to(coord.GCRS(obstime=now))
# print position
print(itrs.cartesian.xyz)
# print velocity
print(itrs.cartesian.differentials)
Обе версии дают один и тот же ответ, но вторая немного аккуратнее.
person
Matt Pitkin
schedule
16.03.2018