如上图所示,用
fv打开fits文件看到,有四个扩展(Extension)
读取第二个扩展,也就是RMS,然后写到一个新的fits文件中去IDL命令如下:
首先第一步,读取原文件:/root/data/*.fits,读取的是原文件的第二个扩展:EXTEN_NO=2,第二个扩展的数据放在了数组im中,第二个扩展的头文件放在了字符串数组:Header中。
im = readfits("/root/data/*.fits",Header, EXTEN_NO=2)
第二部将数据im和头文件Header写到文件/root/data/RMS.FITS中
WRITEFITS,"/root/data/RMS.FITS", im, Header
这样在目录/root/data/就可以看到新生成的fits文件
第三步转化为ndf格式,这种格式文适用于starlink。
convert
fits2ndf /root/data/RMS.FITS /root/data/RMS
这样在目录/root/data就可以看到RMS.sdf
ndf 转化为fits
ndf2fits /root/data/RMS.ndf /root/data/RMS
https://m.sciencenet.cn/blog-783377-641296.html
上一篇:
starlink中的cupid的使用方法(一)下一篇:
将fits数据中的NaN替换为随机数