使用Healpy制作星图

问题描述

我正在使用Healpy从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,180),cmap = 'binary')
hp.graticule()

解决方法

您可以使用numpy.bincount创建每个像素的星系数目的数组，然后创建一个映射。