Healpy (Python で開発された HEALPix ツール) を使用して、全天の CMB マップを読み書きしています。
スカイ マップ上のピクセルの座標が numpy ndarray のエントリにどのように変換されるかについて混乱しています。FITS ファイル内の 1 つの HEALPix ピクセルはどのように ndarray エントリに変換されますか?
たとえば、デフォルトの RING スキームを使用して、生成された CMB マップ、Nside=64、lmax=64 があるとします。ピクセルの総数は、Npix=12*Nside**2 で与えられます。したがって、私の例では、合計 49152 ピクセルです。
私が持っている FITS ファイルは、1 列 x 48 行の形式です。(1 列はすべて温度値です。) したがって、このファイルを n-dim numpy ndarray に変換します。この配列の形状は (48,1024) になります。
質問 1: この 2 次元配列は、空のすべてのマップ ピクセルを「行列」形式に投影したものですか? それはそうですね。私の例を使用すると、48*1024 = 49152、ピクセルの総数です。
質問 2: 1024 規格はどこから来たのですか? HEALPix マップとその座標に関して、これはどういう意味ですか? 1 つの 2-dim ndarray エントリ (ピクセル) は、別のエントリ (角度、HEALPix の位置など) とどのように関連していますか?
別の例として、Nside=1024 をマップします。これを 2 次元の ndarray に変換すると、(12288, 1024) という形になります。
編集: 1024 要素の配列の規則が存在するのはなぜですか? これは地図上のピクセル座標にどのように対応していますか?