anafast如何处理健康中的遮罩？

Question

我正在尝试使用anafast in healpy从被遮罩的地图计算角功率谱。与F90 version（显式地将遮罩文件作为可选输入）相比，anafast的python版本如何处理遮罩的影响？

谢谢！

Answer 1

最简单的方法是将所有掩码值简单地设置为hp.UNSEEN或0（在内部UNSEEN期间，所有hp.anafast()值都将设置为0）。在下面的示例中，我使用numpy.where()。

import numpy as np
import healpy as hp

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

您还可以通过将输入映射变成masked array来传递遮罩，其效果与上述方法完全相同。

内部，为了计算功率谱，所有掩码值都设置为0，因此请注意没有零均值的映射。无论如何，您可能要在计算功率谱之前减去单极子。