gse_plot.py

import matplotlib.pyplot as plt
from astropy import units as u
from astropy.time import Time

from poliastro.bodies import Earth
from poliastro.frames import GeocentricSolarEcliptic
from poliastro.plotting.static import StaticOrbitPlotter
from poliastro.twobody import Orbit

EPOCH = Time("2018-02-19")

ss_2018cb = Orbit.from_sbdb("2018cb").propagate(EPOCH)
ss_earth = Orbit.from_body_ephem(Earth, epoch=EPOCH)

plotter = StaticOrbitPlotter()
plotter.set_frame(GeocentricSolarEcliptic(obstime=EPOCH))
plotter.plot(ss_2018cb, label="2018cb", color="green")
plotter.plot(ss_earth, label="Earth", color="navy")
plt.show()