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

【问题标题】：How to obtain square maps (smaller patches) from a full sky map in Healpy如何从 Healpy 中的完整天空地图中获取方形地图（较小的补丁）
【发布时间】：2021-11-15 16:35:08
【问题描述】：

我有一张完整的天空会聚图。我想从这张地图中切出更小的方形补丁。

我想知道这怎么可能。

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

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})

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

【问题讨论】：

标签： python astronomy healpy

【解决方案1】：

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

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

【讨论】：

【解决方案2】：

感谢@Andrea Zonca，这就是解决方案现在我们可以跳过 Ra 和 Dec (argument = rot[ra, dec]) 来覆盖整个天空

map_cut=hp.gnomview(map_kappa_lcdm,rot=[280.5,32.9],reso=0.8,min=-3e-4,max=3e-4,cmap='RdYlBu_r')

【讨论】：