【问题标题】:Healpy map2alm and alm2map errorsHealpy map2alm 和 alm2map 错误
【发布时间】:2016-08-15 07:50:06
【问题描述】:

我正在学习healpy。具体来说,我对 map2alm 和 alm2map 感兴趣:我发现 question 提供了将函数“应用于双方”(作为答案)的示例。如我们所见,在连续应用map2alm和alm2map后,与原始地图相比,我们得到了明显的错误。因此我想控制这些错误。所以我的问题是:

  1. 有什么方法可以减少错误,而不使用更高的 nside? (比如我想用某个nside)
  2. 另外,也许还有另一种方法可以得到地图的“光谱特征”,不会有这么大的错误?

【问题讨论】:

  • 能否提供具体的代码示例?
  • @AndreaZonca 我的意思是我在上面提到的 stackoverflow 问题中已经做了什么。例如:“m”的第一个元素是 0,但“new_map”的第一个元素是“0.15859344”。所以差异(误差)相当大,我的问题是关于如何减少它(例如到0.01)。

标签: python healpy


【解决方案1】:

不,没有别的办法。

从实空间到调和空间就像进行快速傅里叶变换。如果你没有足够的时间分辨率,你的 FFT 就不会很好,只是没有足够的信息来准确表示频域。

这里是一样的,在低 nside 下获得低误差的唯一方法是,如果你有一个很容易在谐波空间中表示的信号,例如纯偶极子,请参阅这个笔记本,即使在 nside 8、误差在1/1000:

http://nbviewer.jupyter.org/gist/zonca/2f87cb19ca8293539947b0d4aeca42b2

【讨论】:

    猜你喜欢
    • 2014-11-16
    • 2015-09-02
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2021-07-08
    • 1970-01-01
    • 2016-03-12
    • 2019-08-28
    相关资源
    最近更新 更多