0

我想知道 healpy 中的 alm2fits 是什么。

假设我在 healpy 中有 3 个用于 T、E 和 B 的施舍,我如何导出它们以便通过 fit2alm 的 fortran 代码读取它们?

hp.fitsfunc.write_alm('alms.fits',[alm_T,alm_E,alm_B])

似乎不喜欢列表和

hp.fitsfunc.mwrfits(filename= 'alm.fits',data=[alm_T,alm_E,alm_B],colnames=['T','E','B'],keys=None)

似乎以错误的顺序保存它们,因此 fit2alm 引发 'Inconsistent l^2+l+m+1 -> l,m mapping' 错误

谢谢

一个。

4

0 回答 0