Add HEALPIX reprojection

Tool to convert HEALPIX images to Galactic or ICRS coordinates, based tightly on @cdeil's function in gammapy.  See https://github.com/healpy/healpy/issues/129.
Ah, I had the impression that this was used here fore the *output* header but I agree it wouldn't work on the input HEALPIX header.