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

【问题标题】：Galactic and equatorial projected grids with matplotlib使用 matplotlib 的银河和赤道投影网格
【发布时间】：2018-03-26 15:26:27
【问题描述】：

我正在尝试使用两个坐标网格：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（下图），但如果可能的话，我宁愿避免使用它（我发现它使用起来相当复杂）。

【问题讨论】：

明显丑陋的最后解决方法是用点（标记='。'）绘制位于指定赤道子午线/平行线上的灰色“点云/点网格”，并在转换后绘制它们银河系坐标（？）
您链接到的包似乎基于 matplotlib。所以在他们的源代码中你会找到那个网格的定义。

标签： python matplotlib coordinates astronomy

【解决方案1】：

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

【讨论】：