使用matplotlib的银河系和赤道投影网格

Question

我试图使用两个坐标网格绘制投影在完整天空地图上的天文源的散点图：galactic和equatorial。生成单个网格图（使用其中一个预定义的投影）很简单：

import matplotlib.pyplot as plt
import numpy as np

# Random galactic coordinates
l_rad = np.random.uniform(-np.pi, np.pi, 100)
b_rad = np.random.uniform(-np.pi / 2., np.pi / 2., 100)

plt.figure()
plt.subplot(111, projection="aitoff")
plt.grid(True)
plt.plot(ra_rad, dec_rad, 'o', markersize=2, alpha=0.3)
plt.show()

我需要的是一个叠加了两个网格的地图：一个用于银河坐标（已经存在），另一个用于赤道。也就是说，第二个网格如下图所示（注意小时约定而不是右上角的十进制度数）应叠加到第一个网格上：

我认为Kapteyn包可以something similar（图片如下），但我希望尽可能避免使用它（我觉得使用起来相当复杂）。

Answer 1

我还没有使用过它，但是我注意到Astropy.WCSAxes包中有一个使用覆盖银河和赤道网格的地图绘制示例：https://docs.astropy.org/en/stable/visualization/wcsaxes/overlaying_coordinate_systems.html