我正在尝试将healpix格式的CMB映射从Nside=4096降级到Nside=2048。据我所知,有两种方法可以做到这一点:
(a)使用map2alm转换为alm,然后使用alm2map转换回Nside=2048映射。 (b) 在healpix/healpy中使用ud_等级
我测试了这两个,得到了不同的结果。对于(a),我恢复了正确的输入光谱,但是在梯度较大的像素周围(即点源周围),我看到了“振铃”效应。当我使用(b)时,我不能恢复输入光谱。我试过用pixwin(2048)**2除以我的输出频谱,但这并不能解决问题。在
什么是正确的方法来降级一个真实的地图,保持正确的频谱,而不诱导振铃效应在地图?在
提前谢谢你
如果你的主要要求是保存光谱,那么使用方法(a)是最好的。然后你需要屏蔽点源,如果你的应用程序处理的是光谱,你可能还是想这么做。在球谐函数域中,锐利梯度的表达并不理想。在
相关问题 更多 >
编程相关推荐