tskymap: Calculates sky density maps
tskymap calculates a weighted density map
(or, to put it another way, a histogram) on the sky
from columns of an input table.
The sky is divided up into some discrete set of tiles according to
a specified tessellation scheme (currently HEALPix or HTM are supported),
and the required quantities are aggregated into bins corresponding to
The output table has a column giving the pixel index identifying each tile,
plus one or more columns each representing an aggregation of a quantity
from the input table.
By default the number of rows from the input table falling within
each tile is included as the first column in the output table.
But by specifying the
you can add more columns giving the sum, mean, median or other statistics
of input table columns or expressions as well.
The output table can then, for instance, be plotted using
healpix layer type
(though an alternative is to do that plot directly using a
In the case of HEALPix tiling, the result can also be output
in a FITS file suitable for use by external applications that
understand the semi-standard
Note in this case, for maximum compatibility,
fits-healpix output format should in general be used.