获得healpy map的积分

Question

我有一个 healpy 贴图，它是在 CMB 的特定透镜模拟过程中获得的偏转（透镜势的梯度）。如果可能，我想使用 healpy 获得透镜潜力图。我注意到有一个治疗功能alm2map_der1() 它将给我一个治疗地图及其给定地图施舍的一阶导数。我假设这个一阶导数是地图的梯度 - 如果我错了，请纠正我。本质上我想知道我是否可以使用healpy来做这个的倒退过程。我想移除渐变，我只想要透镜潜力。

到目前为止，我一直在尝试使用偏转和透镜势能谱之间的关系；偏转的 Cls = l(l+1) * 透镜势的 Cls，并将其重新排列为：透镜势的 Cls = 偏转的 Cls / l(l+1)，然后使用 synfast 将其转换回映射。我似乎没有得到正确的地图。

有没有更好的方法来做我想做的事情？也许甚至不使用healpy？

Answer 1

我对第一部分无能为力，但我知道转换为 Cls 会破坏方向信息。 'synfast' 提供了一个您输入的功率谱图，但方向随机。如果我在只有偶极子功率的 Cls 列表上运行“synfast”，我每次运行时都会得到一个随机偶极子方向：

import healpy
healpy.mollview(healpy.synfast([0,1],32,lmax=1))

如果你最终想要一张地图，我建议只使用施舍。