HEALPix conventions

This page contains many equations that may not show up correctly in HTML. We recommend that the Postscript document be used instead.
A bandlimited function f on the sphere can be expanded in spherical harmonics, Ylm, as \begin{displaymath}f({ \gamma})=\sum_{l=0}^{l_{max}}\sum_{m}a_{lm}Y_{l
m}(\gamma),
\end{displaymath}
where $ \gamma$ denotes a unit vector pointing at polar angle $ \theta$ $ \in$ [0,$ \pi$] and azimuth $ \phi$ $ \in$ [0, 2$ \pi$). Here we have assumed that there is insignificant signal power in modes with l > lmax and introduce the notation that all sums over m run from - lmax to lmax but all quantities with index lm vanish for m > l. Our conventions for Ylm are defined in subsection A.4 below.

Pixelising f ($ \gamma$) corresponds to sampling it at Npix locations $ \gamma_{p}^{}$, p $ \in$ [0, Npix - 1]. The sample function values fp can then be used to estimate alm. A straightforward estimator is \begin{displaymath}\hat{a}_{lm}=\sum_{p=0}^{N_{pix}-1} Y^\ast_{lm}(\gamma_p) f(\gamma_p)
\end{displaymath}
where the superscript star denotes complex conjugation. This zeroth order estimator, as well as higher order estimators, are implemented in the Fortran90 facility anafast, included in the package. Publications discussing these estimators and the generalised quadrature problem on the sphere (Górski; Wandelt et al. ) are in preparation.



Subsections
Eric Hivon 2003-02-07