Advertisement
Guest User

Untitled

a guest
Jul 21st, 2019
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.60 KB | None | 0 0
  1. import matplotlib.pyplot as plt
  2. from astropy import units as u
  3. from astropy.time import Time
  4.  
  5. from poliastro.bodies import Earth
  6. from poliastro.frames import GeocentricSolarEcliptic
  7. from poliastro.plotting.static import StaticOrbitPlotter
  8. from poliastro.twobody import Orbit
  9.  
  10. EPOCH = Time("2018-02-19")
  11.  
  12. ss_2018cb = Orbit.from_sbdb("2018cb").propagate(EPOCH)
  13. ss_earth = Orbit.from_body_ephem(Earth, epoch=EPOCH)
  14.  
  15. plotter = StaticOrbitPlotter()
  16. plotter.set_frame(GeocentricSolarEcliptic(obstime=EPOCH))
  17. plotter.plot(ss_2018cb, label="2018cb", color="green")
  18. plotter.plot(ss_earth, label="Earth", color="navy")
  19. plt.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement