使用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

您可以使用宽度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)

这给了我以下内容：