今天合作者建议我把数据平滑一下重新计算。我用的那个程序的输入就是一个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
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)