一个程序计算的结果会输出到一个.FIT文件中,以前用软件打开看过,就是一个数据表。但是用IDL尝试了很多此,都读不出来,总是提示没有图(image)。于是以前都是手动把数据表搞出来。
最近数据量越来越大,手动已经可以把人烦死了。于是迫不得已再开始尝试用IDL读这种文件。今天的尝试比较顺利,顿悟了一下,终于知道这个文件存的是二进制表,要用mrdfits读的时候要指明数据所在的位置(就是一个标号)。于是有如下脚本
PRO readfits
fitsname='two_cat.FIT'
a=mrdfits(fitsname,1)
print,a
END
这样可以把数据表里的数输出到屏幕,不过我没有搞明白怎么把一个二进制表转换成一个数组。所以就用sed和tr等处理了。
https://blog.sciencenet.cn/blog-117333-380421.html
上一篇:
CUPID FINDCLUMPS 拟合分子云核(三)下一篇:
运行简单的IDL脚本