Healpy pix2ang如何读取像素索引？

Question

这是这个问题的延续：How do HEALPix FITS files of CMB maps translate into ndarrays? What are the coordinates?

CMB地图是FITS文件。地图是像素温度值的一维向量。在

我想使用Healpy的pix2ang函数来读取每个像素的位置，并知道哪个像素是哪个。在

http://healpy.readthedocs.org/en/latest/generated/healpy.pixelfunc.pix2ang.html

对于这个函数，“像素索引”的输入到底是什么？我必须先将每个温度值像素转换成球谐函数吗？在

我知道在球谐函数中，每个与创建的FITS文件扩展名相对应的aμm将包含一个索引为l^2+l+m+1的整数列。但是不使用l和m的像素数组呢？在

Answer 1

它与球谐函数无关。在

像素数是一维数组中的位置。所以第一个元素是像素0，第二个元素是像素1，依此类推。在

在nside 1的映射的最简单情况下（仅12像素）：

import healpy as hp
nside = 1
npix = hp.nside2npix(1)
print("Number of pixels:")
print(npix)
m = 100 + np.arange(npix)
print("Map:")
print(m)
print("Plot map")
%matplotlib
hp.mollview(m)
print("Pixel theta (colatitude) and phi (longitude) in radians:")
i_pix = np.arange(npix)
theta, phi = hp.pix2ang(nside, i_pix)
for i in range(npix):
    print("Pixel %d : [ %.2f, %.2f ]" % (i, theta[i], phi[i]))

打印：

^{pr2}$