Правильный способ деградации карты healpix

Question

Я пытаюсь деградировать карту CMB в формате healpix от Nside = 4096 до Nside = 2048. Насколько мне известно, есть два способа достичь этого:

(a) конвертировать в alm с помощью map2alm, а затем преобразовать обратно в карту Nside = 2048, используя alm2map. (b) использовать ud_grade в healpix/healpy,

Я тестировал оба эти и получаю разные результаты. Для (a) я восстанавливаю правильный входной спектр, но вижу эффекты «звонка» вокруг пикселей, где градиент большой (т. Е. Вокруг точечных источников). Когда я использую (b), я не восстанавливаю входной спектр. Я попытался делить свой выходной спектр на pixwin (2048) ** 2, но это не устраняет проблему.

Каков правильный способ деградации реалистичной карты, сохраняющей правильный спектр и не вызывающий звуковых эффектов на карте?

Спасибо заранее,

Spectrum using method (b)

Ringing effect around point sources

Answer 1

Если основным требованием является сохранение спектра, происходит с помощью метода (а) является лучшим. Затем вам нужно замаскировать точечные источники, если ваше приложение имеет дело со спектрами, вы, вероятно, захотите это сделать. Резкие градиенты слабо представлены в области сферических гармоник.