matlab代码读电离层TEC文件

matlab代码读电离层TEC文件

function matrix_vtec = load_inx_file（ filename ）
% 读取一个INX文件
% 将每天的结果保存为一个mat文件
% 2013年4月5日22:13:47创建
% 2015年1月31日13:54:13修改

fid=fopen（ filename ）;
if fid==-1
    error（‘文件打开失败！‘）;
end

% 读INX文件的文件头
while feof（fid）==0
    tline=fgetl（fid）;
    if strfind（tline‘END OF HEADER‘）
        break;
    end
end

% 保存VTEC的矩阵，大小为71*73*12
matrix_vtec = nan（717312）;

epoch = 0;
lat_cnt = 0;

% 读INX文件记录部分
while ~feof（fid）
    tline=fgetl（fid）;
    if strfind（tline‘START OF RMS MAP‘）
        break;
    end
    if strfind（tline‘START OF TEC MAP‘）
        epoch = epoch + 1;
        if epoch > 12
            break;
        end
        lat_cnt = 0;
    end
    if strfind（tline‘LAT/LON1/LON2/DLON/H‘）
        vtec_lat = [];
        for i=1:5     % 一个时段的记录有5行
            tline=fgetl（fid）;
            vtec = str2num（ tline ）;  % 将一行的记录转换为向量
            vtec_lat = [vtec_lat vtec];
        end
        lat_cnt = lat_cnt + 1;
        matrix_vtec（lat_cnt:epoch） = vtec_lat;
    end 
end

fclose（fid）;

% INX文件的单位为0.1TECu，将单位变成TECu
matrix_vtec = matrix_vtec/10; % #ok

% len = length（filename）;
% save（[filename（len-11:len-5）‘.mat‘]‘matrix_vtec‘）;

 disp OK!! 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1314  2015-11-01 10:40  load_inx_file.m

----------- ---------  ---------- -----  ----

                 1314                    1