我想知道healpy中alm2fits的等价物。
假设我在T,E和B中有3个施舍,如何导出它们以便可以通过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'
错误
由于
一个。