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

Question

我正在尝试使用两个坐标网格绘制投影在全天 map 上的天文源的散点图: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()

但我需要的是一张叠加了两个网格的 map :一个用于银河坐标(已经存在)，另一个用于赤道。也就是说，下图中的第二个网格(注意赤经的小时约定而不是十进制度)应该叠加到第一个网格上:

我认为 Kapteyn包可以做something similar (下图)但我希望尽可能避免使用它(我发现它使用起来相当复杂)。

Answer 1

我没用过，但我注意到 Astropy.WCSAxes 包有一个绘制带有重叠银河和赤道网格的 map 的示例:https://docs.astropy.org/en/stable/visualization/wcsaxes/overlaying_coordinate_systems.html