使用Healpy制作星图

Question

我正在使用 healpy 从 RA 和 Dec 列表中绘制出天空中星系的位置。到目前为止，我认为我已经能够正确绘制星系，但我想改进成品。有什么方法可以对出现在每个 healpy 瓦片中的星系数量进行分类，而不是仅仅根据瓦片中是否存在目录成员来着色？

我在这里展示我正在制作的图像 —

目前，它仅在告诉您银河系不在何处时真正有用。这是我正在使用的代码。

phis = [np.deg2rad(ra) for ra in ra_list]
thetas = [np.pi / 2 - np.deg2rad(dec) for dec in dec_list]

pixel_indices = hp.ang2pix(NSIDE, thetas, phis)

m = np.zeros(hp.nside2npix(NSIDE))
m[pixel_indices] = np.ones(num_galaxies_to_plot)


hp.mollview(m, title = 'Sky Locations of GLADE Galaxies', cbar = False, rot=(180, 0, 180), cmap = 'binary')
hp.graticule()

Answer 1

您可以使用 numpy.bincount 创建每像素星系数的数组，然后创建它的地图。