How to get NaNs to work in rice compressed FITS files with astropy?

Question

I am having trouble with NaNs in astropy.fits. When they are present anywhere in the file, we get white horizontal lines of NaNs across the entire image, but only if the HD