In [76]:
import healpy as hp
In [77]:
import numpy as np
In [78]:
dip = np.array([1,0,0])
In [79]:
nside = 8
npix = hp.nside2npix(nside)
In [80]:
pix_vec_x, pix_vec_y, pix_vec_z = hp.pix2vec(nside, np.arange(npix))
In [81]:
m = np.array([np.dot(dip, vec) for vec in zip(pix_vec_x, pix_vec_y, pix_vec_z)])
In [82]:
%matplotlib inline
In [83]:
hp.mollview(m)
In [84]:
alm = hp.map2alm(m)
In [85]:
m2 = hp.alm2map(alm, nside)
Sigma is 0.000000 arcmin (0.000000 rad) -> fwhm is 0.000000 arcmin
In [87]:
hp.mollview(m-m2, title=" error")
In [ ]: