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

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

function HeadO=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,'RINEX VERSION'))

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

elseif(strfind(s,'APPROX POSITION XYZ'))

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

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

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

elseif (strfind(s,'MARKER NAME'))

HeadO.staname=s(1:4);

elseif (strfind(s,'MARKER NUMBER'))

HeadO.stanum=s(1:4);

elseif (strfind(s,'# / TYPES OF OBSERV'))

HeadO.PRNtype='G';

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

for i=1:HeadO.nObstype

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

end

elseif (strfind(s,'SYS / # / OBS TYPES'))

k=k+1;

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

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

if(HeadO.nObstype(k)>13)

for i=1:13

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

end

s=fgets(fid);

for i=1:HeadO.nObstype(k)-13

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

end

else

for i=1:HeadO.nObstype(k)

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

end

end

elseif (strfind(s,'END OF HEADER'))

break;

end

end

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

% Read RINEX Observation Data

nepoch=0;

nobs=0;

nlli=0;

rnx.epoch=[];

rnx.obs=[];

rnx.lli=[];

while 1

if feof(fid)==1,break,end

line=fgetl(fid);

if isempty(deblank(line)),break,end

if version < 3 % Version 2 RINEX

% get epoch flag and number of satellites

EpStr=line(1:26);

EpFlg=sscanf(line(27:29),'%d');

NumSat=sscanf(line(30:32),'%d');

if EpFlg>1

% epoch flags

% 0 - o.k.

% 1 - power failure between previous and current epoch

% event flags (epoch flag continued)

% 2 - start moving antenna

% 3 - new site occupation (end of kinematic data); at least a record

% with a new m