Check
out the new, improved version of HEALPix
(Hierarchical Equal Area isoLatitude Pixelisation of
the sphere). This routine pixelates the maps that will
be made by TopHat. It facilitate fast and accurate statistical
and astrophysical analysis of full-sky data sets.