下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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;endk=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.nObstypeHeadO.Obs{i}={s(i*6+5:i*6+6)};endelseif(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:13HeadO.Obs{k,i}=s(i*4+4:i*4+6);ends=fgets(fid);fori=1:HeadO.nObstype(k)-13HeadO.Obs{k,i+13}=s(i*4+4:i*4+6);endelsefori=1:HeadO.nObstype(k)HeadO.Obs{k,i}=s(i*4+4:i*4+6);endendelseif(strfind(s,'ENDOFHEADER'))break;endend下面为读取数据程序,供参考
%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
ifversion<3
%Version2RINEX
%getepochflagandnumberofsatellites
EpStr=line(1:26);
EpFlg=sscanf(line(27:29),'%d');
NumSat=sscanf(line(30:32),'%d');
ifEpFlg>1
%epochflags
%
0-o.k.
%
1-powerfailurebetweenpreviousandcurrentepoch
%eventflags(epochflagcontinued)
%
2-startmovingantenna
%
3-newsiteoccupation(endofkinematicdata);atleastarecord
%
withanewmarkernamemustfollow
%
4-newheaderinfofollows
%
5-externalevent(epochissignificant)
%
6-cyclesliprecordsfollow(sameformatasobservations)
%
%thenumberofrecordswhichfollowisgivenbyNumSat
fori=1:NumSat
line=fgetl(fid);
end
continue
end
%getepochtime
t=sscanf(EpStr,'%f',6)';
ift(1)<80
t(1)=2000+t(1);
else
t(1)=1900+t(1);
end
%keyboard
%-adaptionrvbree25112008---
%
date
=datenum(t(1),t(2),t(3));%Theserialdatenumbersfrom1-Jan-0000tothedayofobservation
%
gpsweek
=fix((date-datenum(1980,1,6))/7);
%
secsinweek
=(date-datenum(1980,1,6)-gpsweek*7)*24*3600+t(4)*3600+t(5)*60+t(6);
%
%
EpTime
=secsinweek;
%
GPSW
=gpsweek;
%EpTime=datenum(t(1),t(2),t(3),t(4),t(5),t(6));
%----endadaption-----------------------------------------
%EpTime=datenum(t(1),t(2),t(3),t(4),t(5),t(6));
%readoptionalclockerror***wedonotneedtheclockerror***
%ClkOffset=sscanf(line(69:end),'%f');
%getsatelliteid's
%-sprncontainsthe3lettercodesforthesatellites
%-prncontainstheuniversalsatelliteid's
sprn=line(33:min(end,68));
forj=2:ceil(NumSat/12)
line2=fgetl(fid);
sprn=[sprnline2(33:min(end,68))];
end
prn=rxsatnum(sprn);
%gettheobservations,lliandsnrforallsatellites
obs=zeros(NumSat,NumRxObs);
fori=1:NumSat
line=fgetl(fid);line=[lineblanks(80-length(line))];
forj=2:ceil(NumRxObs/5)
line2=fgetl(fid);
line=[lineline2blanks(80-length(line2))];
end
%gettheobservations
p=0;
forj=1:NumRxObs
obsi=sscanf(line(p+1:p+14),'%f');
%testformissingobservations
ifisempty(obsi)|obsi==0,obsi=NaN;,end
obs(i,j)=obsi;
p=p+16;
end
%getthelli(lossoflockindicator)forcarrierphasedata
forjj=1:length(icarrier)
j=icarrier(jj);
flg=double(line((j-1)*16+15))-48;
ifrem(flg,2)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 电影行业安全生产工作总结
- 传统制造业技术职位展望
- 二零二五年度航空航天材料试验委托协议3篇
- 二零二五年度房屋收购合同环保验收与评估范本3篇
- 二零二五版养老院专业保洁及消毒服务合同2篇
- 二零二五版个人二手房购房合同与产权过户指导书
- 航空行业助理的职位介绍
- 汽车行业财务预测分析工作总结
- 二零二五年度产品责任纠纷民事答辩状范文3篇
- 二零二五年度木材市场树木买卖协议3篇
- 问题探究如何让城市不再看海(教学课件)高一地理
- 2024年人教版五年级数学(上册)模拟考卷及答案(各版本)
- 人教版八年级下册历史第1课 中华人民共和国成立 说课稿
- 《地球物理勘查》全册配套完整教学课件
- (正式版)JBT 5300-2024 工业用阀门材料 选用指南
- 九年级上册-备战2024年中考历史总复习核心考点与重难点练习(统部编版)
- 健康指南如何正确护理蚕豆病学会这些技巧保持身体健康
- 老客户的开发与技巧课件
- 26个英文字母书写(手写体)Word版
- GB/T 13813-2023煤矿用金属材料摩擦火花安全性试验方法和判定规则
- 日语专八分类词汇
评论
0/150
提交评论