在软件使用方面,我很多时候都是参考别人的意见,不少事情没有亲自试验,于是脑子里就形成了一些过时的偏见,诸如,“Fortran的计算能力最强大,所以学天文的人要使用Fortran”、“IDL功能强大,学天文的人要会使用IDL”、“Matlab虽然好用,但是无法处理fits文件,所以天文中一般不用”。但是渐渐就会发现,这些说法是站不住脚的,多数只是习惯造成的,并不说明某种软件的强大或缺陷。
Matlab是我用过的计算工具(高级语言)中最人性化的,语法合理,各种帮助、示例都很齐全。但是一直以来我都很少使用Matlab,因为很长一段时间都没有明白如何使用Matlab读入FITS文件,而我也懒得去研究。
最近处理数据的时候发现某些操作用Matlab实在是方便,于是又开始考虑Matlab读fits文件的问题。多方查找,发现Matlab读fits文件其实非常简单
data=fitsread('文件路径/文件名')
就可以了,读入data的就是文件中所存储的数组。然后就可以按数组的方式对数据进行处理了。这相比IDL丝毫不逊色。而FITS文件头可以通过
[header dumy]=fitsload2('文件路径/文件名')
把文件头作为字符串数组读入header中,使用str2num就可以把字符串转换为数字为数据处理所用了。
fitsload2.m
https://m.sciencenet.cn/blog-117333-696838.html
上一篇:
python 笔记及示例(五)主分量分析与矩阵本征值下一篇:
Matlab处理分子云数据的可能性(二)积分强度图