support healpix decomposition for RA, Dec distribution

probably faster if this is implemented directly using healpix rather than constructing separate distributions for each pixel and then constucting a mixture model. However, having a separate object for each pixel and building a mixture model would allow us to more conveniently swap this with the IAU constellation implementation...