带有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