I'd like to plot the path of an object near L2 (like the Webb space telescope) in a rotating libration point coordinate system. That makes it is easy to distinguish halo orbits from more general Lissajous forms, etc. I expect to get Earth-centered coordinates from Horizons (see e.g. nealmcb/jwst_orbit: Explore the orbit of the James Webb Space Telescope in Python).

Here are details on the RLP coordinate system I'm looking for:: Technical Notes - Rotating Libration Points

Putting the origin at L2 as they suggesst seems good.

I was hoping to find an implementation in astropy, but it doesn't seem to be there. Has anyone implemented it for astropy or another open-source framework?

  • 1
    $\begingroup$ That's essentially what I did in the corotating frame plots here: space.stackexchange.com/a/57832/38535 $\endgroup$
    – PM 2Ring
    Commented Jan 28, 2022 at 11:21
  • $\begingroup$ I did something like this here but in a less-than-rigorous way, I rotated with the Earth/Moon barycenter's instantaneous position rather than at a steady rate. $\endgroup$
    – uhoh
    Commented Jan 29, 2022 at 23:38


You must log in to answer this question.

Browse other questions tagged .