Координация вращения CMB с использованием healpix

Question

У меня есть космическая микроволновая фоновая карта, с которой я читаю healpy, Я заинтересован в извлечении полос пикселей на широте +10 градусов в северном полушарии и на широте -10 градусов в южном полушарии. Я мог легко извлечь пиксели с +10 градусов на север, используя hp.ang2pix() функция. Но в южном полушарии мне трудно определить угол, потому что тета варьируется от 0 до пи.

Должен ли я повернуть систему координат сферы на пи радианы, чтобы извлечь пиксель из 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)

Это дает мне следующее: