PDF document and several MATLAB functions and demonstration scripts that can be used to compute important selenocentric (moon-centered) coordinate information and transformations. Functions are provided for computing the orientation of the moon with respect to the Earth mean equator and equinox of J2000 (EME2000) system along with lunar libration angles extracted from a JPL binary ephemeris file.
Most of the information in this section was extracted from JPL D-32296, “Lunar Constants and Models Document” which is available at ssd.jpl.nasa.gov/dat/lunar_cmd_2005_jpl_d32296.pdf. Another useful reference is “Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009”, Celestial Mechanics and Dynamical Astronomy, 109: 101-135, 2011.
David Eagle (2020). Selenocentric Coordinates and Transformations (https://www.mathworks.com/matlabcentral/fileexchange/39216-selenocentric-coordinates-and-transformations), MATLAB Central File Exchange. Retrieved .