Tutorial

Let's make an Enmap, the primary structure in this package.

using Pixell, Plots
shape, wcs = fullsky_geometry(300.0 * Pixell.arcminute)  # set up the map geometry
m = Enmap(randn(shape), wcs)  # generate a random map with the shape and WCS
plot(m)
Example block output

Let's compute a spherical harmonic transform with Libsharp, and then compute the power spectrum $C_{\ell}$.

cl = alm2cl(map2alm(m))
plot(cl, ylabel=raw"$C_{\ell}$", xlabel=raw"Multipole moment, $\ell$")
Example block output

Reading and writing maps

m = read_map("map.fits")

Noteworthy differences from Python

The Julia language has a number of differences from Python, and the Pixell.jl has a number of important differences from the Python package pixell.