Удаление режимов низкого уровня с карты CMB с помощью HealPix ⇐ Python

[Anonymous]

Существует 2 подхода к удалению режимов низкого уровня с карты CMB.

• Имея доступ к милостыне в качестве милостыни;

Код: Выделить всё

    lmax = hp.Alm.getlmax(len(alm))
ls, ms = hp.Alm.getlm(lmax)
alm[ls < 10] = 0  #Removing first 10 ls
isw_map_2 = hp.alm2map(alm, nside)

Код: Выделить всё

def remove_low_ell(map,ell,nside):
alm = hp.map2alm(map,lmax=3*nside-1,mmax=3*nside-1)
k = []
for i in range(len(alm)):
if i < ell+1:
k.append(0.0)
else:
k.append(1.0)
alm2=hp.almxfl(alm, k)
map = hp.alm2map(alm2,nside)
return map

По сути они одинаковы, но результирующая карта будет немного отличаться.
Что вы порекомендуете?
Я попробовал оба метода выше

Подробнее здесь: https://stackoverflow.com/questions/790 ... th-healpix