healpy - anafast 如何处理 healpy 中的掩蔽？

Question

我正在尝试使用healpy 中的 anafast从蒙版地图中计算角功率谱。与F90 版本相比，anafast 的 python 版本如何处理掩码的效果，后者明确将掩码文件作为可选输入？

谢谢！

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 以计算功率谱，因此请注意不具有零均值的地图。在任何情况下，您都可能希望在计算功率谱之前减去单极子。