【问题标题】:How exactly does healpy's anafast compute the angular power spectrum?healpy 的 anafast 究竟是如何计算角功率谱的?
【发布时间】:2019-12-13 16:07:00
【问题描述】:

我试图理解为什么随机丢弃我的一半数据会使 Cls 降低。我在 healpix 文档或谷歌上找不到它。

Angular power spectrum

编辑:我以前从未使用过这样的东西。我拥有的数据覆盖了整个球体的 1/4,红移在 0.69 和 1.3 之间,质量光晕在 log10M = 13 和 log10M = 15 之间。我有角位置、红移、质量和粒子数。 我用零填充了一张地图,并在地图上的位置添加了粒子的数量。这是正确的方法吗,anafast 是否考虑到这一点并为部分地图计算正确的 Cls?

【问题讨论】:

    标签: healpy


    【解决方案1】:

    HEALPix 中的球面谐波变换始终在全天空,如果地图被遮盖,这些像素将设置为 0。建议在运行 anafast 之前从地图中移除单极子以减少边界效应。

    我将此文本添加到 anafast 的文档字符串中。

    【讨论】:

    • 您好,感谢您指出这一点。但我仍然不明白为什么会这样。 anafast 是否只是总结了我的地图中两个点的每个组合的角功率谱,没有任何归一化?如果我想比较它们,我是否需要在不同的红移 bin 中使用相同数量的数据点?
    • 最简单的方法是制作一个联合蒙版并将它们与相同的天空覆盖范围进行比较。否则,您至少可以校正天空分数。您必须将光谱除以天空分数或其平方,我不记得了。
    • 不确定我的问题与掩码有多大关系:当我使用 anafast 从 Planck 2018 数据图中提取 Cl 时,(wget -O COM_CMB_IQU-commander_2048_R3.00_full "pla.esac.esa.int/pla-sl/data-action?MAP.MAP_OID=13470") 我得到了一个 Cl 大约是 Planck 2018 发布的 CL 的六倍。我应该注意哪些可能导致这种不匹配的陷阱?
    猜你喜欢
    • 1970-01-01
    • 2019-07-13
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2016-08-12
    • 2015-05-24
    相关资源
    最近更新 更多