RINEX格式导航文件读取

读书报告《GPS原理及应用》题目RINEX格式导航文件读取学院土木建筑学院专业测绘科学与技术学生学号RINEX格式导航文件读取一、RINEX文件格式简介RINEX（ReceiverINdependentExchange）格式是与接收机无关的数据交换格式，该格式采用文本文件存储数据，数据记录格式与接收机的制造厂商和具体型号无关。

现在，RINEX格式已经成为了GPS测量应用等的标准数据格式，几乎所有测量型GPS接收机厂商都提供将其格式文件转换为RINEX格式文件的工具，而且几乎所有的数据分析处理软件都能够直接读取RINEX格式的数据。这意味着在实际观测作业中可以采用不同厂商、不同型号的接收机进行混合编队，而数据处理则可采用某一特定软件进行。RINEX格式文件包括6种文件类型：观测数据文件：GPS观测值导航电文文件：GPS卫星导航电文

气象数据文件：在测站处所测定的气象数据

GLONASS导航电文文件：GLONASS卫星导航电文

GEO导航电文文件：增强系统中搭载有类GPS信号发生器的地球同步卫星（GEO）的导航电文

卫星和接收机钟文件：包含卫星和接收机时钟信息

其中用的最多的是O文件、N文件，观测值文件的文件头存放有文件的创建日期、单位名、测站名、天线信息、测站近似坐标、观测值数量及类型、观测历元间隔等信息。导航电文的文件头存放有文件创建日期、单位名及其他一些相关信息,另外,还有可能会包含电离层模型的参数以及说明GPS时与UTC间关系的参数和跳秒等。气象数据文件的文件头则存放有文件创建日期、观测值类型、传感器信息和气象传感器的近似位置及其他一些相关信息。

RINEX格式文件的记录数据紧跟在文件头的后面,随文件类型的不同,所存放数据的内容和具体格式也不相同。在观测值文件中存放的是观测过程中每一观测历元所观测到的卫星及载波相位、伪距和多普勒等类型的观测值数据等,所包含的实际观测值类型与接收机所记录的类型及格式转换时的参数设置有关。在导航电文文件中存放的是所观测卫星的钟差改正模型参数及卫星轨道数据等。在气象数据文件中存放的是观测过程中每隔一段时间在测站天线附近所测定的干温、相对湿度和气压等数据。二、RINEX格式导航电文头文件如图一所示，导航电文头文件说明如下：RINEXVERSION/TYPE：RINEX版本号和观测值类型；PGM/RUNBY/DATE：创建本数据文件所采用程序的名称，创建本数据文件单位的名称，

创建本数据文件的日期。COMMENT：注释行。IONALPHA：历书中的电离层参数A0～A3；IONBETA：历书中的电离层参数B0～B3；DELTA-UTC:A0,A1,T,W：用于算UTC时间的历书参数。A0,A1:多项式系数；

T:UTC数据的参考时刻

W:UTC参考周数,为连续计数,不是1024的余数LEAPSECONDS：由于跳秒而造成的时间差ENDOFHEADER:头文件结束标签

图一RINEX头文件样例三、RINEX记录行格式

RINEX格式GPS导航电文文件数据记录节中的内容为按卫星和参考时刻存放的各颗卫星的时钟和轨道数据。每颗卫星一个参考时刻的数据占8行,第1行为卫星的PRN号和该卫星时钟的参考时刻及其改正模型参数,第2～8行为该卫星的广播轨道数据。由于导航电文通常每2h就更新一次,因此,某些卫星可能会有多个不同参考时刻的数据。PRN号/历元/卫星钟：

-卫星的PRN号

-历元:OC(卫星钟的参考时刻)

年(2个数字,如果需要可补0)

月,日,时,分,秒

-卫星钟的偏差(s)

-卫星钟的漂移(s/s)

-卫星钟的漂移速度(s/s2)

广播轨道—1：

广播轨道—2：

广播轨道—3：

广播轨道—4：

广播轨道—5：

广播轨道—6：

-卫星精度(m)

-卫星健康状态(第1子帧第3字第17～22位)

-TGD(sec)

-IODC钟的数据龄期

广播轨道—7：

-电文发送时刻(单位为GPS周的秒,通过交接字(HOW)中的Z计数得出)

-拟合区间(h),如未知则为零

-备用

-备用四、导航电文提取算法过程1、打开文件2、按行提取每一行记录3、判断头文件结束标志“ｅｎｄｏｆｈｅａｄ”4、按行提取数据记录，按导航电文格式提取每一行数存到数组中5、关闭文件，文件提取结束输出内容读取数据记录存放到结构数组每一行从第60位开始取13位与endofhead比较，直至相等为止读到文件尾，关闭文件按行读取文件文件打开文件输出内容读取数据记录存放到结构数组每一行从第60位开始取13位与endofhead比较，直至相等为止读到文件尾，关闭文件按行读取文件文件打开文件图二简易流程图关键代码说明//*按指定位置及长度取字符串函数*//voidxstrmid(char*szDest,char*line,intstart,intlen){ char*src_index; char*dst_index=line;if(len>0&&start>=0) { src_index=szDest+start; while(len>0&&*src_index!='\0') { *dst_index=*src_index; src_index++; dst_index++; len--; } *dst_index='\0'; }}GMNGETGMNREC(GMNgmn,char*filename)//*定义读取文件记录函数*//{ intj=0; char*temp; charstr0[]="ENDOFHEADER";//*定义文件头结束标签*// charstr1[25]; charstr2[MAX_IN_ROW]; charstr3[21]; FILE*fp; fp=fopen(filename,"r");//*打开文件*// do { fgets(str2,MAX_IN_ROW,fp);//*按行读取文件每条记录*//xstrmid(str2,str1,60,13); }while(strcmp(str0,str1)!=0);//*判断头文件结束标志*// temp=fgets(str2,MAX_IN_ROW,fp); while(strlen(str2)>15&&temp&&j<MAX_Satellite) {inti; for(i=0;i<MAX_IN_ROW;i++) { if(str2[i]=='D')//*替换字符D为e,用于数值转换*// str2[i]='e'; }xstrmid(str2,str3,0,2);//*按RINEX文件格式逐个读取PRN号、时间、卫星钟偏差*//gmn.Record[j].PRN=atoi(str3);//*卫星钟漂移、卫星钟漂移速度*// xstrmid(str2,str3,3,2); gmn.Record[j].TOC.year=atoi(str3); if(gmn.Record[j].TOC.year>80)gmn.Record[j].TOC.year+=1900; elsegmn.Record[j].TOC.year+=2000;xstrmid(str2,str3,6,2); gmn.Record[j].TOC.month=atoi(str3);//*读取月份存储到记录数组*//xstrmid(str2,str3,9,2); gmn.Record[j].TOC.day=atoi(str3);xstrmid(str2,str3,12,2); gmn.Record[j].TOC.hour=atoi(str3);xstrmid(str2,str3,15,2); gmn.Record[j].TOC.minite=atoi(str3);xstrmid(str2,str3,17,5); gmn.Record[j].TOC.second=atoi(str3);xstrmid(str2,str3,22,19); gmn.Record[j].ClkBias=atof(str3);//*读取卫星钟偏差存储到记录数组*//xstrmid(str2,str3,41,19); gmn.Record[j].ClkDrift=atof(str3);//*读取卫星钟漂移存储到记录数组*//xstrmid(str2,str3,60,19); gmn.Record[j].ClkDriftRate=atof(str3);//*读取卫星钟漂移速度存储到记录数组*//temp=fgets(str2,MAX_IN_ROW,fp)//*读取星历发播时间IODE、、、*// for(i=0;i<MAX_IN_ROW;i++) { if(str2[i]=='D') str2[i]='e'; }xstrmid(str2,str3,3,19); gmn.Record[j].IODE=atof(str3);//*读取IODE存储到记录数组*//xstrmid(str2,str3,22,19); gmn.Record[j].Crs=atof(str3);xstrmid(str2,str3,41,19); g