python-3.x - 使用 healpix 的 CMB 坐标旋转

Question

我有一张正在阅读的宇宙微波背景图healpy。我有兴趣提取北半球+10度纬度和南半球-10度纬度的像素条。hp.ang2pix()我可以使用函数轻松地从 +10 度北提取像素。但是在南半球，我发现很难定义角度，因为 theta 从 0 到 pi 不等。

我应该将球体的坐标系旋转 pi 弧度以便从南半球 10 度提取像素吗？

我正在使用以下程序从北半球提取像素条：

import numpy as np
import healpy as hp

fname = 'COM_CMB_IQU-070-fgsub-sevem-field-Pol_1024_R2.01_full.fits'

tmap = hp.read_map(fname)
nside = hp.get_nside(tmap)

x = hp.ang2pix(nside, np.deg2rad(10) , [0, 2*3.14])

print(x)

Answer 1

您可以使用 width 执行以下操作eps来定义恒定纬度环：

import healpy as hp
import numpy as np

nside = 128
npix = hp.nside2npix(nside)
x = np.arange(npix)

# All in degrees
glon, glat = hp.pix2ang(nside, np.arange(npix), lonlat=True)
eps = 0.5

# Set up the mask
mask = (glat < 10. + eps) & (glat > 10. - eps)
mask |= (glat > -10. - eps) & (glat < -10. + eps)

hp.mollview(x*mask)

这给了我以下信息：