Computing power spectrum of already masked map about namaster HOT 3 CLOSED

Hi @umilta
I think what you're doing should work. The main thing to remember is that, when one calls w.compute_coupling_matrix, the two fields one passes should have their true masks associated with them.
Note however, that if you try to use the systematics deprojection or purification functionality in NaMaster, it won't really work for f0, since these functions also need to know the mask.

To make things slightly less hackish, I'd do the following:

f_other = nmt.NmtField(apodized_mask, other_map)
w = nmt.NmtWorkspace()
# Divide by mask where it's not zero
map_unmasked = np.zeros_like(my_map)
map_unmasked[binary_mask>0.001] = my_map[binary_mask>0.001] / apodized_mask[binary_mask>0.001]
f0 = nmt.NmtField(apodized_mask, map_unmasked)
cl_decoupled = w.decouple_cell(nmt.compute_coupled_cell(f0, f0))

I agree that it'd be nice to be able to pass masked maps to NaMaster, so I'll open an issue to address that. Thanks!

from namaster.

@umilta this feature has now been implemented in the master branch.

from namaster.

Thanks!

from namaster.