我刚刚开始在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计算功率谱。
答案 0 :(得分：0)
您只需要使用以下方法之一将theta和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