python - 适合 WCS 来治愈

Question

假设我有一个包含 WCS header 的 FITS 文件，这样我就可以执行以下操作:

#import healpy as hp
#import astropy.io.fits as pyfits
#from astropy.wcs import WCS

listofhdus = pyfits.open(FITS)
wcs = WCS(listofhdus[0].header)

listofhdus[0].data 将包含一个 2D numpy 数组(NY 乘以 NX)，对应于银河坐标中完整天空的一小部分。

如果我想使用 Mollweide 投影将该天空图叠加到以下散点图中，将 2d numpy 转换为 helpix 格式(称为 WCS)的最佳方法是什么？

NPIX = hp.nside2npix(512)
m = 0*np.arange(NPIX)
hp.mollview(m, title="test")
hp.projscatter(longitude, latitude, lonlat=True, 
               coord='G',marker='+',color='black')
hp.graticule()

其中经度和纬度是两个二维数组(使用 numpy.meshgrid 生成)以及我感兴趣的一些天体物理源的银河坐标？

我想我可以尝试将我的healpix像素转换为坐标，以某种方式将它们与我的天空图中可用的像素相匹配，并从那里插入值，但一定有一些更简单，更优雅和精确的东西，对吧？ >

Answer 1

您应该查看 reproject ，它可以满足您与 HEALPix 之间的所有重投影需求。

特别是查看reproject_to_healpix .

如果您需要更好地控制插值例程(使用基于卷积的网格器，确保通量守恒和性能)，我建议您查看 cygrid ，其中有 some examples for HEALPix也是如此。