healpy.sphtfunc.map2alm¶
- healpy.sphtfunc.map2alm(maps, lmax=None, mmax=None, iter=3, pol=True, use_weights=False, datapath=None, gal_cut=0, use_pixel_weights=False, verbose=True)¶
Computes the alm of a Healpix map. The input maps must all be in ring ordering.
For recommendations about how to set lmax, iter, and weights, see the Anafast documentation
Pixel values are weighted before applying the transform:
when you don’t specify any weights, the uniform weight value 4*pi/n_pix is used
with ring weights enabled (use_weights=True), pixels in every ring are weighted with a uniform value similar to the one above, ring weights are included in healpy
with pixel weights (use_pixel_weights=True), every pixel gets an individual weight
Pixel weights provide the most accurate transform, so you should always use them if possible. However they are not included in healpy and will be automatically downloaded and cached in ~/.astropy the first time you compute a trasform at a specific nside.
If datapath is specified, healpy will first check that local folder before downloading the weights. The easiest way to setup the folder is to clone the healpy-data repository:
git clone –depth 1 https://github.com/healpy/healpy-data cd healpy-data bash download_weights_8192.sh
and set datapath to the root of the repository.
- Parameters:
- mapsarray-like, shape (Npix,) or (n, Npix)
The input map or a list of n input maps. Must be in ring ordering.
- lmaxint, scalar, optional
Maximum l of the power spectrum. Default: 3*nside-1
- mmaxint, scalar, optional
Maximum m of the alm. Default: lmax
- iterint, scalar, optional
Number of iteration (default: 3)
- polbool, optional
If True, assumes input maps are TQU. Output will be TEB alm’s. (input must be 1 or 3 maps) If False, apply spin 0 harmonic transform to each map. (input can be any number of maps) If there is only one input map, it has no effect. Default: True.
- use_weights: bool, scalar, optional
If True, use the ring weighting. Default: False.
- datapathNone or str, optional
If given, the directory where to find the pixel weights. See in the docstring above details on how to set it up.
- gal_cutfloat [degrees]
pixels at latitude in [-gal_cut;+gal_cut] are not taken into account
- use_pixel_weights: bool, optional
If True, use pixel by pixel weighting, healpy will automatically download the weights, if needed
- verbosebool, optional
Deprecated, has not effect.
- Returns:
- almsarray or tuple of array
alm or a tuple of 3 alm (almT, almE, almB) if polarized input.
Notes
The pixels which have the special UNSEEN value are replaced by zeros before spherical harmonic transform. They are converted back to UNSEEN value, so that the input maps are not modified. Each map have its own, independent mask.