HealPix / Healpy：从（MODIS）卫星数据创建 healpix 地图

【问题标题】：HealPix / Healpy: Creat healpix map from (MODIS) satellite dataHealPix / Healpy：从（MODIS）卫星数据创建 healpix 地图
【发布时间】：2020-03-02 12:37:39
【问题描述】：

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

输入映射是否需要采用 FITS 格式，或者我可以将 HDF4 MODIS 数据文件读入数组（使用 python），然后使用 healpy 将其映射到球体上吗？

我遇到了这个答案：https://stackoverflow.com/a/50495134/7157742 并希望在我的情况下也有类似的可能。有什么建议吗？

【问题讨论】：

标签： healpy

【解决方案1】：

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

您需要将数据映射到正确的 HEALPix 像素中。为此，您可以获取数据的 colatitude theta 和 longitude phi 数组，然后使用 hp.ang2pix 查找它们映射到的像素。

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

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

【讨论】：