Description

Loosely related to #9652 and #9685, is there a specific reason why io.fits does not support lzma compression? From the way .xz files are handled in data.py it looks like it could be implemented quite analogously in the FITS file handlers, but of course it would involve somewhat more significant modifications.