MATLAB 等间隔(2h)读取RINEX气象文件(简洁版)

以darw站2014年数据为例

clear all;

clc;
day='缺失天数';         %定义变量名day,记录缺失的文件年积日
m=[];                   %定义空矩阵,存储数据
for n=1:366              %m文件天数1天到366天
    ns=num2str(n,'%03d');         % 数值转字符串
name=['darw',ns,'0.14m'];  % 字符串拼接得到文件名称
ffid=fopen(name,'rt');     % 打开文本文件,从第一行开始读取,数据存储在cell里
if ffid>0
i = 1;
while feof(ffid) == 0    
    tline{i,1} = fgetl(ffid);
    i = i+1;
end    
% 提取间隔2h的数据,从0时至22时,RINEX m文件采样间隔为30s,2h间隔为240行
for i=18:240:2658         
    a((i+222)/240,:)=str2num(tline{i,1});
end
% 提取压强p 和温度 t的值
for i=1:12
    b(i,1)=a(i,7);
    b(i,2)=a(i,8);
end
m=cat(1,m,b);
else 
    day=strcat(day,'-',ns);
end
 
end
xlswrite('14m.xlsx',m)
阅读更多
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页