我正计划将MODIS的卫星数据映射到一个球体上，我认为Healpix可能是正确的方法。然而，我不知道如何进行：
输入地图需要是FITS格式吗？或者我可以将HDF4 MODIS数据文件读取到数组中(使用python(，然后使用healpy将其映射到球体上吗？
我发现了这样的答案：https://stackoverflow.com/a/50495134/7157742并希望在我的情况下也能有类似的事情发生。有什么建议吗？
是的，您可以使用标准的numpy数组。
您需要将数据映射到正确的HEALPix像素中。为此，您可以获得数据的经度theta和经度phi的数组，然后使用hp.ang2pix找到它们映射到的像素。
然后您可以创建一个1D阵列healpy映射并使用它
查看Jupyter笔记本示例https://gist.github.com/zonca/680c68c3d60697eb0cb669cf1b41c324