使用 healpix 的 CMB 坐标旋转答案

【问题标题】：Coordinate rotation of CMB using healpix使用 healpix 的 CMB 坐标旋转
【发布时间】：2018-01-07 13:09:50
【问题描述】：

我正在使用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)

【问题讨论】：

请发布一个最小的工作示例。

标签： python-3.x healpy

【解决方案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)

这给了我以下信息：

【讨论】：