读取GPS观测数据O文件的matlab编程.pdfVIP

functionHeadO=ReadObsHead

[fname,fpath]=uigetfile(*.*O,选择一个O文件);

HeadO.O_filename=strcat(fpath,fname);

fid=fopen(HeadO.O_filename,rt);

if(fid==-1)

msgbox(文件不存在,警告,warn);

return;

end

k=0;

while(1)

s=fgets(fid);

if(strfind(s,RINEXVERSION))

HeadO.Rtype=str2double(s(6:10));

elseif(strfind(s,APPROXPOSITIONXYZ))

HeadO.XYZ(1)=str2double(s(1:14));

HeadO.XYZ(2)=str2double(s(15:28));

HeadO.XYZ(3)=str2double(s(29:42));

elseif(strfind(s,MARKERNAME))

HeadO.staname=s(1:4);

elseif(strfind(s,MARKERNUMBER))

HeadO.stanum=s(1:4);

elseif(strfind(s,#/TYPESOFOBSERV))

HeadO.PRNtype=G;

HeadO.nObstype=str2double(s(6));

fori=1:HeadO.nObstype

HeadO.Obs{i}={s(i*6+5:i*6+6)};

end

elseif(strfind(s,SYS/#/OBSTYPES))

k=k+1;

HeadO.PRNtype(k)=char(s(1));

HeadO.nObstype(k)=str2double(s(5:6));

if(HeadO.nObstype(k)13)

fori=1:13

HeadO.Obs{k,i}=s(i*4+4:i*4+6);

end

s=fgets(fid);

fori=1:HeadO.nObstype(k)-13

HeadO.Obs{k,i+13}=s(i*4+4:i*4+6);

end

else

fori=1:HeadO.nObstype(k)

HeadO.Obs{k,i}=s(i*4+4:i*4+6);

end

end

elseif(strfind(s,ENDOFHEADER))

break;

end

end

下面为读取数据程序，供参考

%ReadRINEXObservationData

nepoch=0;

nobs=0;

nlli=0;

rnx.epoch=[];

rnx.obs=[];

rnx.lli=[];

while1

iffeof(fid)==1,break,end

line=fgetl(fid);

ifisempty(deblank(line)),break,end

ifversio