Sentinel-1 orbits#
[1]:
import numpy as np
import matplotlib.pyplot as plt
from pymgp.satorb import eci2ecef
from pymgp.tleorb import tleread, tle2vec
from pymgp.orbplot import plh2eci, obj2sat, xyz2zassp
from pymgp.orbplot import pltposvel, pltsatview, pltsattrack, pltgroundtrack, skyplot, plt3dorbit
[3]:
tle = tleread('resource-20240802.tle', verbose=0)
[4]:
xsat, vsat, t, satid = tle2vec(tle, ['2024-08-02 0:00', 24*60 ,1], 'SENTINEL-1A')
#xsat, vsat, t, satid = tle2vec(tle, ['2024-08-02 0:00', 24*60 ,1], 'JASON-3')
xsat.shape
[4]:
(1441, 3)
[5]:
xsate, vsate = eci2ecef(t, xsat, vsat)
xsate.shape
[5]:
(1441, 3)
pltposvel#
[6]:
pltposvel(t,xsat,vsat, satid=satid, rstyle={'linewidth':1, 'color':'k'})
[6]:
array([<Axes: title={'center': 'SENTINEL-1A position and velocity in ECI (2024-08-02)'}, ylabel='Position [km]'>,
<Axes: xlabel='Date 2024-08-02', ylabel='Velocity [km/s]'>],
dtype=object)
