获得healpy地图的积分

Question

我有一个healpy贴图，它是偏转（透镜电位的梯度），这是在CMB的特定透镜模拟过程中获得的。如果可能的话，我想用healpy获得透镜电位图。我注意到有一个healpy函数alm2map_der1()，它会给我一张healpy map和它的第一个导数给出地图的施舍。我假设这个第一个衍生物是地图的渐变 - 如果我错了，请纠正我。基本上我想知道我是否可以使用healpy来做这个向后的过程。我想删除渐变，我只想要透镜效果。

到目前为止，我的尝试是利用偏转和透镜电位功率谱之间的关系;偏转的Cl = 1（l + 1）*透镜电位的Cls，并将其重新排列为：透镜电位的Cls =偏转的Cls / l（l + 1），然后使用synfast将其转换回图。我似乎没有得到正确的地图。

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

Answer 1

我无法帮助完成第一部分，但我知道转换为Cl会破坏方向信息。 ＆＃39; synfast＆＃39;提供一个地图，其中包含您输入的功率谱，但具有随机方向。如果我跑了'synfast＆＃39;在只有偶极子功率的Cl列表中，每次运行时都会得到一个随机偶极方向：

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

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