anafast 如何处理 healpy 中的掩蔽？

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