用Healpy制作星图

Question

我正在使用healpy从RAs和Decs列表中绘制出天空中星系的位置。到目前为止，我认为我已经能够正确地绘制星系图，但我希望改进最终的产品。有没有办法将每个健康的磁贴中出现的星系数量进行分类，而不是仅仅根据磁贴中是否有星表成员进行着色？

这里我展示了我目前正在制作的图像-

​

​

现在它只对告诉你银河系不在哪里有用。这是我正在使用的代码。

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创建每个像素的星系数组，然后创建该数组的地图。