上次说到用mrdfits读FITS文件，用/fscale关键字以后可以对数据进行重新标定，得到真实值。但是最近我画一幅图的时候发现有些地方总是画不出来，想了很久也没明白。

前几天写信给提供数据的人，这才知道，文件里有些地方是没有观测的，应该是空白的，而在FITS文件里是用BLANK关键字的值来替代这些空白的。mrdfits并没有处理这些空白，只是把它们当成数据读进来和普通数据一样进行处理，于是得出来的结果就莫名其妙了。

我想了一下，先把空白的地方的指标记下来
IDL> iblank=where(a eq blank,count)
随后待标定完了再另这些地方为零
IDL> if(count gt 0) then a(iblank)=0.0
如此之后就一切正常了。

转载本文请联系原作者获取授权，同时请注明本文来自钱磊科学网博客。
链接地址：https://m.sciencenet.cn/blog-117333-386619.html

上一篇：听徐遐生（Frank Shu）先生的报告
下一篇：小记IDL写FITS文件