Как Anafast заботится о маскировке в здоровых?

Question

Я пытаюсь вычислить угловой спектр мощности из маскированной карты, используя anafast в healpy .Как Python-версия anafast заботится о последствиях маскирования по сравнению, например, с F90-версией , которая явно принимает файл маски в качестве необязательного ввода?

Спасибо!

Answer 1

Проще всего установить все маскированные значения на hp.UNSEEN или на 0 (внутренне все значения UNSEEN в любом случае будут установлены на 0 в течение hp.anafast()).В приведенном ниже примере я использую numpy.where () .

import numpy as np
import healpy as hp

masked_map = np.where(mask, raw_map, hp.UNSEEN)
cl = hp.anafast(masked_map)

Вы также можете передать маску, превратив входную карту в маскированный массив ,который фактически действует так же, как подход, описанный выше.

Внутренне все маскированные значения установлены в 0 для вычисления спектра мощности, поэтому следует помнить о картах, которые не имеют нулевого среднего.В любом случае вы можете вычесть монополь перед вычислением спектра мощности.