如何從 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})

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

Answer 1

您可以使用hp.gnomview繪制方形補丁。

如果你想擁有一片天空的二維數組，你可以在gnomview設置return_projected_map=True ，這將返回它以供進一步處理。

Answer 2

這就是解決方案，感謝@Andrea Zonca 現在我們可以跳過 Ra 和 Dec（參數 = 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')