CMB 地图的 HEALPix FITS 文件如何转换为 ndarrays?坐标是多少?

How do HEALPix FITS files of CMB maps translate into ndarrays? What are the coordinates?

我正在使用 Healpy(在 Python 中开发的 HEALPix 工具)来读取和写入全天 CMB 地图。

我对天空地图上的像素坐标如何转换为 numpy ndarray 的条目感到困惑。 FITS 文件中的一个 HEALPix 像素如何转换为一个 ndarray 条目?

例如,假设我有一个生成的CMB图,Nside=64,lmax=64,使用默认的RING方案。像素总数由 Npix=12*Nside**2 给出。因此,对于我的示例,像素总数为 49152。

我的 FITS 文件格式为 1 列乘 48 行。 (第 1 列是所有温度值。)因此,我将此文件转换为 n-dim numpy ndarray。该数组的形状为 (48,1024)。

问题一:这个二维数组是不是所有地图像素在天空中的投影成"matrix"形式?我觉得是这样的。使用我的示例,48*1024 = 49152,像素总数。

问题二:1024标准从何而来?就 HEALPix 地图及其坐标而言,这意味着什么?一个 2 维 ndarray 条目(像素)如何与另一个相关(按角度、按 HEALPix 位置等)?

再比如map Nside=1024。当我将其转换为二维 ndarray 时,我得到了形状 (12288, 1024)。

编辑:为什么存在这种 1024 元素数组的约定?这如何对应于地图上的像素坐标?

有关 HEALPIX 像素化的最佳信息来源是 Healpix 入门:

http://healpix.jpl.nasa.gov/html/intro.htm

地图是一维向量,每个索引是天空中的一个像素。 您可以使用 hp.pix2ang.

将像素数转换为天空中的坐标

将其按 1024 行排列只是使 read/write 更快地写入磁盘的一个技巧。