python - 如何从 Healpy 中的完整天空地图中获取方形地图（较小的补丁）

Question

我有一张完整的天空会聚图。我想从这张地图上剪下更小的方形斑块。

我想知道这怎么可能。

如果我阅读地图并绘制它，这就是它的外观：

map_kappa_lcdm = hp.read_map('lightcone.fits',nest=False);
import seaborn as sns
color_pal = sns.color_palette("colorblind", 6).as_hex()
fig, (ax1) = plt.subplots(figsize = (17,8),ncols=1)
################
fontsize = 16
plt.rcParams.update({'font.size':fontsize})
cax = hp.mollview(map_kappa_lcdm,title="convergence ($\kappa$) map, $\Lambda$CDM",xsize=4000, 
                  min=-0.02, max=0.025,hold=True,norm="%hist", format="%.01g",cmap='RdBu',cbar=True,nlocs=4)
hp.graticule()
plt.rcParams.update({'font.size':fontsize})

此后我不知道如何切割更小的方形补丁。

标签： pythonastronomyhealpy

Answer 1

您可以hp.gnomview用于绘制方形补丁。

如果您想拥有一片天空的 2D 数组，您可以进行设置return_projected_map=True，gnomview然后将其返回以进行进一步处理。