HealPix / Healpy:根据(MODIS)卫星数据创建 healpix 地图

HealPix / Healpy: Creat healpix map from (MODIS) satellite data

我打算将来自 MODIS 的卫星数据映射到一个球体上,我认为 Healpix 可能是正确的方法。但是,我不知道该怎么做:

输入地图是否需要为 FITS 格式,或者我可以将 HDF4 MODIS 数据文件读入数组(python),然后使用 healpy 将其映射到球体上吗?

我遇到了这个答案: 并希望在我的案例中也能有类似的东西。有什么建议吗?

是的,您可以使用标准 numpy 数组。

您需要将数据映射到正确的 HEALPix 像素中。 为此,您可以获得数据的余纬度 theta 和经度 phi 数组,然后使用 hp.ang2pix.

找到它们映射到的像素

然后您可以创建一维数组 healpy 地图并使用它。

查看示例 Jupyter Notebook https://gist.github.com/zonca/680c68c3d60697eb0cb669cf1b41c324