In working on DM-33591 I noticed that the way that HealSparseInputMapTask converts bounding boxes to sky polygons will be inaccurate in regions of large distortion. This can be easily fixed by (a) adding a padding around the bounding box before converting to sky coordinates, and (b) confirming that all pixels are within the bounding box prior to this line: https://github.com/lsst/pipe_tasks/blob/38dd528b55dbde8a9bb3cbc6bd9163444e0d701d/python/lsst/pipe/tasks/healSparseMapping.py#L195