如何让 NaN 在带有 astropy 的大米压缩 FITS 文件中工作？

问题描述

我在 astropy.fits 中遇到了 NaN 问题。当它们出现在文件中的任何位置时，我们会在整个图像中看到白色的 NaN 水平线，但前提是 HDU 被压缩。我们意识到这些线等于图块大小。

我们希望使用 zBLANK 属性来消除这些，但它似乎不起作用。例如：

from astropy.io import fits as pyfits
import numpy as np

# Create date to be output as data_rd.  Then...
data_rd[np.where(np.isnan(data_rd)) = -2147483648

pri_hdu = pyfits.PrimaryHDU()
sec_hdu = pyfits.CompImageHDU(data=data_rd,tile_size=[10,10],quantize_level=1000)
# Add other headers.  Then...
sec_hdu.header.append(("zBLANK",-2147483648,"The integer value to replace with NaNs"))

hdu = pyfits.HDUList([pri_hdu,sec_hdu])
hdu.writeto(product_filename,overwrite=True,checksum=True)

但是当我查看图像中的数据时，数据总是看起来像-2147483648，而不是NaN。

我做错了什么？我们也不能让它与整数缩放一起工作。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）