FITS (Flexible Image Transport System)是天文学界常用的数据格式,由文件头和数据组成。
我需要分别处理来自EIT和AIA的fits文件,但因为Fits文件的标准并不统一,做实验的时候,对于EIT的fits文件在IDL中进行处理,AIA的fits文件在matlab中处理。很长一段时间都没有找到能同时处理这两种fits文件的方法(如果有,烦请不吝赐教),到现在也没有明白,为什么我用IDL中的fisread不能读取AIA的fits文件;matlab中的fitsread不能读取EIT的fits文件(在此也向大家请教)。
IDL中读取和显示Fits文件:
mreadfits,'D:\EIT_Dat\efz19970512.155012',index1,data1,header=header
data1=mrdfits('D:\test_data\AIA20130120_1428_0211.fits',1,header)
index2map,index1,data1,map1
plot_map,map1,/log
运行结果:
Matlab中读取和显示Fits文件:
data = fitsread('D:\Learning\Master\Project\AIAData\20100705\Test\AIA20100705_0252_0171.fits', 'image');
data = flipud(data);
fitsShow(data);
%%显示函数,将灰度值转换到[0,255]