从开始做这个课题到现在就没少用IDL读FITS文件。这个方面用mrdfits比较容易，基本就是一行搞定数据，几行搞定文件头，用了不知多少次。其实在读FITS的时候就在想，把写FITS也搞明白吧，不过惰性太大，一直都回避这个问题。

今天合作者建议我把数据平滑一下重新计算。我用的那个程序的输入就是一个FITS文件，这就意味着我需要重新写一个平滑后的FITS文件，于是今天不得不去看看怎么写FITS文件了。原来知道和mrdfits对应的有mwrfits，专门写FITS文件的。我有一个数组a和文件头head，于是按照说明里写
IDL> mwrfits,a,'out.fits',head
这样倒是可以生成一个FITS文件，查看了也正常，可是我用来处理的那个程序就是不认。不得已，参考了一下别人的程序，用writefits
IDL> writefits, 'out.fits', a, head
这样生成的FITS文件就能被识别了。原因为何，有待研究。

在文件头某些信息改变的情况下还需要改一下文件头里的参数，可以用fxaddpar，例如改变参数'NAXIS'的值
IDL> fxaddpar,'NAXIS',2

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

上一篇：小记IDL读FITS文件（外一则）
下一篇：运行COSMOMC