HealPix / Healpy: создание карты healpix со спутниковых данных (MODIS)

Question

Я планирую отобразить спутниковые данные из MODIS на сферу, и я подумал, что Healpix может быть правильным способом сделать это. Тем не менее, я не знаю, как go об этом:

Должна ли входная карта быть в формате FITS или я могу прочитать файл данных HDF4 MODIS в массив (с python) и затем использовать healpy, чтобы отобразить его на сфере?

Я наткнулся на этот ответ: { ссылка } и надеялся, что что-то подобное будет возможно в моем случае. Есть предложения?

Answer 1

Да, вы можете использовать стандартный массив numpy.

Вам необходимо отобразить ваши данные в правильный пиксель HEALPix. Для этого вы можете получить массивы значений theta и долготы phi ваших данных, а затем определить, на какой пиксель они отображаются, используя hp.ang2pix.

Затем вы можете создать одномерный массив healpy карта и работа с ней.

См. пример Jupyter Notebook https://gist.github.com/zonca/680c68c3d60697eb0cb669cf1b41c324