Feature request: FITS header shape to region

[keflavich]

It should be possible, though not entirely trivial, to turn any celestial FITS header into a RectangleSkyRegion. That would be very useful for showing instrument/image footprints on the sky. (see also #457)