Add load_moon_relief to load Moon relief dataset in various resolutions and registrations

[seisman]

Description of proposed changes

Preview: https://pygmt-dev--2674.org.readthedocs.build/en/2674/api/generated/pygmt.datasets.load_moon_relief.html

import pygmt
from pygmt.datasets import load_moon_relief
grid = load_moon_relief(resolution="15m")

fig = pygmt.Figure()
fig.grdimage(grid, cmap="@moon_relief.cpt")
fig.show()

Address #2431.

Reminders

• Run make format and make check to make sure the code follows the style guide.
• Add tests for new features or tests that would have caught the bug that you're fixing.
• Add new public functions/methods/classes to doc/api/index.rst.
• Write detailed docstrings for all functions/methods.
• If wrapping a new module, open a 'Wrap new GMT module' issue and submit reasonably-sized PRs.
• If adding new functionality, add an example to docstrings or tutorials.
• Use underscores (not hyphens) in names of Python files and directories.

Slash Commands

You can write slash commands (/command) in the first line of a comment to perform
specific operations. Supported slash commands are:

• /format: automatically format and lint the code
• /test-gmt-dev: run full tests on the latest GMT development version

[seisman]

Please note that extra_attributes is empty here because I have no idea about the horizontal_datum attribute for the dataset.