Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement