使用healpy从numpy数据集中获得功率谱图

【问题标题】：Power spectrum plot from numpy data set by using healpy使用healpy从numpy数据集中获得功率谱图
【发布时间】：2019-09-19 07:29:53
【问题描述】：

我刚开始在python中使用healpy。很快，我有一个位置数据让

x = np.array([7.47098722,  7.47805867,  7.47098779....12.57804461, 12.5809734])

y = np.array([58.32516929, 58.33223959, 58.33516828, ..., 63.41415155, 63.41707922])

其中 phi(x, y) 是物质在这些位置的潜力。现在，我想为此计算功率谱。但我不明白如何为healpy 设置我的数据以使用mollview 制作地图并使用anafast 计算功率谱。

【问题讨论】：

标签： python numpy healpy

【解决方案1】：

您只需使用以下方法之一将 thetas 和 phis 转换为像素索引 healpy 的像素函数，healpy.pixelfunc.ang2pix。您可以从https://healpy.readthedocs.io/en/latest/generated/healpy.pixelfunc.ang2pix.html了解更多信息

有了索引之后，你也有了这些索引的潜在值，换句话说，你有相应的 HEALPix 映射（使用mollview 来获得 mollview 此地图的投影）。最后，使用healpy的anafast、https://healpy.readthedocs.io/en/latest/healpy_spht.htm，你就有了功率谱。

有关healpy的更多信息，请参考- https://buildmedia.readthedocs.org/media/pdf/healpy/1.8.6/healpy.pdf

https://healpix.sourceforge.io/ 另外，请访问https://healpix.jpl.nasa.gov

【讨论】：