Question

我现在开始使用healpix，并且试图从普朗克目录中绘制源。我有一个FITS文件，其中包含诸如银河系坐标之类的信息，并且我想用mollweide投影生成地图。

首先，我使用函数“ healpy.pixelfunc.ang2pix”将坐标从像素转换为像素。我仍然不了解的参数之一：NSIDE，所以我只输入了2的随机数。其次，我使用了“ healpy.visufunc.mollview”，其中的map参数是从上一个像素获得的像素矢量功能。但是我收到此错误消息“像素数错误”。我不知道我在做什么错..

    import healpy as hp
    from astropy.io import fits
    import numpy as np
    from astropy.table import Table
    import copy as copy

    catalogue=fits.open("/home/luana/Documents/Mestrado/Planck_data/COM_PCCS_070_R2.04.fits")
    COORD=Table(catalogue[1].data)
    nside=np.power(2,4)
    GLON=copy.copy(COORD['GLON'])
    GLAT=copy.copy(COORD['GLAT'])
    hp_pixel_numbers=hp.pixelfunc.ang2pix(nside,GLON,GLAT,nest=True,lonlat=True)     
    hp.visufunc.mollview(map=hp_pixel_numbers, fig=None, rot=None, nest=True,cbar=True, cmap=None, notext=False, norm=None, hold=False, margins=None, sub=None, nlocs=2, return_projected_map=False)

使用Healpix的问题

0 个答案: