Healpy ошибки map2alm и alm2map

Я учусь на здоровье. В частности, меня интересуют map2alm и alm2map: я нашел вопрос, где пример применения функций "в обе стороны" (в ответ) предусмотрено. Как мы видим, мы получаем значительные ошибки по сравнению с исходной картой после последовательного применения map2alm и alm2map. Поэтому я хочу контролировать эти ошибки. Итак, мой вопрос:

  1. Есть ли способ уменьшить ошибку, не используя более высокий nside? (Например, если я хочу использовать определенный nside)
  2. Кроме того, может быть, есть другой способ получить "спектральную характеристику" карты, в котором не было бы таких больших ошибок?

python healpy
person Oiale    schedule 21.04.2016    source источник
comment
Можете ли вы привести конкретный пример с кодом?   -  person Andrea Zonca    schedule 22.04.2016
comment
@AndreaZonca Я имею в виду именно то, что было сделано в вопросе stackoverflow, о котором я упоминал выше. Например: первый элемент m равен 0, а первый элемент new_map равен 0,15859344. Так что разница (ошибка) довольно большая, и мой вопрос был о том, как ее уменьшить (например, до 0,01).   -  person Oiale    schedule 22.04.2016

Ответы (1)


arrow_upward
1
arrow_downward

Нет, другого пути нет.

Переход от реального пространства к гармоническому подобен быстрому преобразованию Фурье. Если у вас недостаточное временное разрешение, ваше БПФ будет не очень хорошим, просто недостаточно информации для точного представления в частотной области.

Здесь то же самое, единственный способ получить низкую ошибку на низком уровне nside — это если у вас есть сигнал, который легко представить в гармоническом пространстве, например, чистый диполь, см. этот блокнот, даже только при nside 8 ошибка составляет 1/1000:

http://nbviewer.jupyter.org/gist/zonca/2f87cb19ca8293539947b0d4aeca42b2

person Andrea Zonca    schedule 24.04.2016