使用 Healpy 制作星图

Question

我正在使用 heappy 从 RA 和 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 创建一个包含每个像素的星系数量的数组，然后创建一个地图。