Astropy support for selenocentric (Moon)reference frames and lunar surface observatories.
An extension to
astropy, providing selenocentric and topocentric reference frames
for the Moon and transformations of star positions from the ICRS system to these
frames. This is to describe the sky as observed from the surface of the Moon.
Non-relativistic transformations are calculated using the SPICE toolkit. Relativistic corrections (stellar aberration) will be added.
lunarsky may be installed with pip:
pip install lunarsky
or from the repository:
git clone https://github.com/aelanman/lunarsky python setup.py install
Definition of the MCMF and lunar topocentric frames, from Fig. 2 of [Ye et al.]
lunarsky provides the following classes:
MCMF– The "Moon-Centered-Moon-Fixed" frame, this is a cartesian reference frame that rotates with the moon. This is chosen to be the Mean Earth/Polar frame, with a Z axis defined by the mean rotation axis of the Moon and a prime meridian defined by the mean direction to the Earth's center (the X axis is through this meridian 90° from the Z axis, and Y is defined such that XYZ is a right-handed system). This is analogous to
- LunarTopo – A topocentric (East/North/Up) frame defined at a position on the Moon's surface. This is analogous to
MoonLocation– Analogous to the
astropy.coordinates.EarthLocationclass, this describes positions on the Moon in either selenocentric (x, y, z) or selenodetic (lat, lon, height) coordinates. The cartesian axes of the selenocentric system are those of the MCMF frame. In the selenodetic coordinates, "height" is defined relative to a sphere of radius 1737.1 km.
SkyCoord– A replacement for
astropy.coordinates.SkyCoord, with modifications that ensure compatibility with the
Time– A replacement for the
astropy.time.Timeclass, which can accept a
MoonLocationinstance as its location parameter, allowing for the calculation of local sidereal times on the Moon. The LST of an observer on the Moon is defined as the right ascension of the zenith.
This package makes use fo the
spiceypy wrapper  for the JPl SPICE Toolkit, produced by the NASA Navigation and Ancillary Information Facility (NAIF)  . The transformations are defined using data in kernel files
moon_assoc_me.tf. These may be found at the NAIF website, and were produced by Nat Bachman (NAIF/JPL) in March 2008. Further information may be found in the comments in these files in the
: Ye, Hanlin, et al. "Looking Vector Direction Analysis for the Moon-Based Earth Observation Optical Sensor." IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 11, no. 11, Nov. 2018, pp. 4488–99. IEEE Xplore, doi:10.1109/JSTARS.2018.2870247.
: Annex et al., (2020). SpiceyPy: a Pythonic Wrapper for the SPICE Toolkit. Journal of Open Source Software, 5(46), 2050, https://doi.org/10.21105/joss.02050
: Acton, C.H.; "Ancillary Data Services of NASA's Navigation and Ancillary Information Facility;" Planetary and Space Science, Vol. 44, No. 1, pp. 65-70, 1996.
: Charles Acton, Nathaniel Bachman, Boris Semenov, Edward Wright; A look toward the future in the handling of space science mission geometry; Planetary and Space Science (2017); DOI 10.1016/j.pss.2017.02.013; https://doi.org/10.1016/j.pss.2017.02.013
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.