search

simonsobs / pspy

Python power spectrum code
Other
7 stars 5 forks source link

Add function to compute and remove monopole & dipole from map #8

Closed xgarrido closed 3 years ago

xgarrido commented 3 years ago

Removing monopole and dipole works for both HEALPIX and CAR pixellisation and all the job is done within the specialized subtract_mono_dipole function. We also create a so_map method in order to remove mono/dipole within a so_map.

Masking the maps can be done in several ways :

  1. Both maps and masks have the number of components (I, Q, U). Then each map component is masked with the corresponding mask component
  2. Maps have more component than mask. For instance, the maps are I, Q, U and only component is provided for the mask. Then the unique mask component is used for every maps.

Any other combination will raise an error.

xgarrido commented 3 years ago

Latest dev. now use a tuple of mask like it is done in the mode coupling computation. In order to remove mono/dipole for I, Q, U map, one should provide a mask as a tuple like (mask_T, mask_pol).