matlab读MICAPS第一类数据格式

daveryjl

function [stdata]=ReadfirstMicaps(filename,fminute,fsecond)
%stdata 为结构数据
%stdata.datetime 时间年月日时分秒
%stdata.num  站点数
%stdata.lon 站点经度
%stdata.lat 站点纬度
%stdata.code 站点代号
%stdata.altitude
%stdata.stlevel站点级别
%stdata.totolCoverCloud;
%stdata.windd
%stdata.windv
%stdata.sealevelpressure
%stdata.varypressure3h
%stdata.passweather1
%stdata.passweather2
%stdata.rainfallsixh
%stdata.lowcloudshape
%stdata.lowcovercloud
%stdata.lowhigh
%stdata.dewpoint
%stdata.visibility
%stdata.nowweather
%stdata.temp  温度
%stdata.midcloudshape
%stdata.highcloudshape
%stdata.indicate1
%stdata.indicate2
%stdata.varytemp24h
%stdata.varypressure24h
[inf,errmsg]=fopen(filename,'r');
if inf==-1
    disp(errmsg);
    stdata=NaN;
    return;
end
aheadscript=fscanf(inf,'diamond 1 %s',1);
stdata.datetime=fscanf(inf,'%d %d %d %d',[1,4]);
stdata.datetime=[stdata.datetime,fminute,fsecond];
stdata.num=fscanf(inf,'%d',1);
element=fscanf(inf,'%ld %f %f %f %d %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f',[26,stdata.num]);
fclose(inf);
stdata.code=element(1,:);
stdata.lon=element(2,:);
stdata.lat=element(3,:);
stdata.altitude=element(4,:);
stdata.stlevel=element(5,:);
stdata.totolCoverCloud=element(6,:);
stdata.windd=element(7,:);
stdata.windv=element(8,:);
stdata.sealevelpressure=element(9,:);
stdata.varypressure3h=element(10,:);
stdata.passweather1=element(11,:);
stdata.passweather2=element(12,:);
stdata.rainfallsixh=element(13,:);
stdata.lowcloudshape=element(14,:);
stdata.lowcovercloud=element(15,:);
stdata.lowhigh=element(16,:);
stdata.dewpoint=element(17,:);
stdata.visibility=element(18,:);
stdata.nowweather=element(19,:);
stdata.temp=element(20,:);
stdata.midcloudshape=element(21,:);
stdata.highcloudshape=element(22,:);
stdata.indicate1=element(23,:);
stdata.indicate2=element(24,:);
stdata.varytemp24h=element(25,:);
stdata.varypressure24h=element(26,:);