《GPS数据处理》实验指导书资料(共20页).doc

《《GPS数据处理》实验指导书资料(共20页).doc》由会员分享，可在线阅读，更多相关《《GPS数据处理》实验指导书资料(共20页).doc（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上GPS数据处理实验指导书（供测绘工程专业用）（2013年）戴吾蛟 编二零一三年四月中南大学测绘与遥感科学系专心-专注-专业前言GPS数据处理实验是学习本门课程的重要环节之一，实验的设置是为了配合课堂教学，使学生熟练掌握和深入理解课堂教学内容，进一步强化GPS数据处理算法及流程。实验内容本门课程必做的实验有（在Visual C+.net 2008 环境下）：1. GPS卫星位置计算程序设计；2. 伪距单点定位程序设计；3. GPS数据质量分析；4. GPS工程控制网数据处理；实验1 GPS卫星位置计算程序设计一、实验目的 熟悉导航文件数据格式及参数意义； 掌握利用导航星

2、历数据进行卫星位置计算的详细步骤； 掌握卫星位置计算程序设计技巧；二、实验内容与要求设计一个卫星位置计算的函数；利用已有的RINEX格式导航电文数据读写函数(NavFileRead.h, NavFileRead.cpp)，从导航文件中读取数据，调用自己设计的卫星位置计算函数计算每颗卫星在星历发布时刻的卫星位置(double transtime; /电文发送时间(GPS星期秒)，并保存至一个文件中。计算的详细步骤：1.计算真近点角fs （1）计算平均角速度（2）计算归化观测时间tktk以参考时刻toe为基准的归化观测时间。t为观测时刻卫星的钟面时t经卫星钟差改正的观测时刻（3）计算观测时刻t的平

3、近点角Ms和偏近点角Es令Ms=Es进行迭代计算（4）计算观测时刻的真近点角2.计算升交距角及轨道摄动改正项参考历元的升交距角：u0=ws+fs 摄动改正项3.计算升交距角、卫星的地心距离及轨道倾角4.计算卫星在轨道坐标系中的坐标（x,y,z）5.计算观测时刻升交点的经度6.计算在协议地球系中的空间直角坐标实验2 伪距单点定位程序设计一、实验目的 熟悉RINEX观测文件的读写； 掌握卫星发射时刻位置的计算； 掌握伪距观测方程列立； 伪距观测方程系数计算及单点定位程序设计详细步骤。二、实验内容 设计一个伪距单点计算的函数，并用测试数据进行测试。，其中RINEX格式导航电文数据读写函数(NavFi

4、leRead.h, NavFileRead.cpp)、RINEX格式观测数据读写函数(ObsFileRead.h, ObsFileRead.cpp)已给出，卫星位置计算函数为实验一设计的函数，电离层和对流层改正暂略。单颗卫星观测数据结构体及一个历元观测数据类如下：/单颗卫星观测数据结构体typedef struct short satID; double vCA; double vY1; double vY2; double vL1; double vL2; double vD1; double vD2; float elevation; float azim; observation; /一

5、个历元的观测数据类class CEpochObs public:CEpochObs(void);CEpochObs(void);public:short nSatCount;/卫星数double GpsWeek;/观测时间GPS周double GpsSecond;/观测时间GPS秒 observation* Obs;bool SetSize(short nSat); /按照Rinex标准格式（2.0以上）读取一个历元的数据short ReadStandardEpochObs(FILE *pObsFile, CEpochObs *pEObs, char *ObsType, short& ObsT

6、ypeCount); /民用时至GPS时转换函数void DateTimeToWeekTOW(short year, short month, short day, short hour, short min, double sec, long *Week, double *TOW); 具体计算步骤：1. 将导航文件中所有卫星的星历全部读取，并存入ephemeris数组中 ；2. 调用相关函数读取一个历元的观测数据；3. 从ephemeris数组中找出当前历元第一颗观测卫星的导航数据，（注意一个导航文件中对一颗卫星来时可能有不同时间的星历数据，因此要查找当前历元最相近时间的星历数据）4. 计算

7、卫星的钟差 (dt)，调用卫星钟差改正函数Dt();5. 计算信号发射时刻 TTrans = TTag (C1/C + dt) 6.计算信号发射时刻卫星位置（调用实验一设计的函数）7.计算卫星至测站之间几何距离（考虑地球自转改正）do r0 = sqrt(pow(*rx,2.)+pow(*ry,2.)+pow(*rz,2.); theta = r0/C*W; /对距离分量进行地球自转改正 *rx = satx - recx+ theta*recy *ry = saty - recy- theta*recx; *rz = satz - recz; *r = sqrt(pow(*rx,2.)+po

8、w(*ry,2.)+pow(*rz,2.);while(fabs(r0-*r)1.e-3); /迭代求解8- 计算常数项（闭合差，不考虑对流层及电离层延迟） W = C1 (Ro2 - (c dt)9. 计算误差方程系数（设计矩阵）A1 = -Rox / Ro2A2 = -Roy / Ro2A3 = -Roz / Ro2A4 = c10.重复39步，计算全部所有卫星观测方程的系数阵A和常数项W11.组成法方程矩阵 N = AtA 与方程常数项U = ATW,权阵为单位阵，即 P=I；（使用已有的CMatrix类）12.测站点初始近似坐标与接收机钟差近似值设为0，进行迭代计算，求得测站点坐标及接

9、收机钟差；13.进行精度评定。先计算观测值残差，再计算单位权中误差；， u=4,n为观测历元数；14.保存伪距单点定位结果至文件中；15.循环读取下一历元观测数据，并重复314步完成该历元伪距单点定位计算；直至观测文件数据全部计算完毕。实验3 GPS数据质量分析一、实验目的 熟悉GPS数据质量分析的基本原理与方法 掌握TEQC软件的使用 熟悉TEQC软件分析结果的报告。二、实验内容1.利用DOS版TEQC软件剔除GPS+GLONASS双系统观测文件中GLONASS观测值；剔除高度角低于25度的观测数据；剔除某一颗卫星（任选一颗卫星）的观测值；提取单频GPS观测数据。 teqc R source

10、.*o GPS.*o teqc G source.*o GLONASS.*o禁用卫星。命令为：teqc G# source.*o result.*o在对对流层和电离层进行延迟分析或者为了保证数据的质量时，需要控制卫星的高度角。命令为：teqc set_mask # source.*o result.*o对于各种观测值类型，依据个人的需要进行观测值类型的定制。命令为：teqc O.obs P1+P2+L1+L2+ca source.*o result.*o2. 利用DOS版TEQC软件对一观测数据文件进行数据质量分析，查看数据质量分析文件；lite模式需要的文件是：O文件。命令为：teqc +q

11、c source.*o生成的文件有7个：*.*S, *.iod, *.ion, *.mp1, *.mp2, *.sn1, *.sn2。full模式需要的文件是：O文件和N文件。命令为：teqc +qc source.*o生成的文件有9个：*.*S, *.iod, *.ion, *.mp1, *.mp2, *.sn1, *.sn2, *.azi, *.ele。*.ion L1-L2电离层延迟；*.iod 电离层延迟的变化率；*.mp1/*.mp2 L1 上的C/A码或P码多路径观测误差和L2上的P码多路径观测误差；*.sn1/*.sn2 L1/L2观测值的信噪比；*.ele 卫星的高度角；*.a

12、zi 卫星的方位角。3.安装并使用WinTEQC软件，熟悉WinTEQC的各项功能，重点学习其数据质量分析功能。针对算例数据进行数据质量分析；4.利用WinTEQC的数据预处理功能处理一观测数据文件，然后再对预处理前后的数据文件进行质量分析，对比预处理前后数据的结果。具体实验步骤：详见TEQC使用说明、WinTEQC软件说明书实验4 GPS工程控制网数据处理一、实验目的1. 熟悉GPS基线解算基本步骤及参数设置含义；2. 掌握中海达GPS控制网数据处理程序操作；3. 掌握高精度基线解算质量控制策略及基本技巧；4. 掌握GPS控制网平差基本步骤；5. 掌握GPS高程拟合；6. 熟悉GPS控制网基

13、线解算及平差报告。二、实验内容利用中海达控制网数据处理软件，完成两套GPS控制网的基线解算、网平差及高程拟合的计算，并编制GPS控制网成果报告；具体包括：（1）GPS基线解算；（2）GPS三维网、二维网平差；（3）GPS 高程拟合实验步骤：六、实验步骤1. 新建项目2. 项目属性设置 选择控制网等级，基线解算与网平差将根据该设置值进行检核。3.导入文件 先把相应的观测文件拷贝到工程文件夹中。4.检查观测文件重点检查天线高与天线类型，检查是否有同一点但命名不同当所用天线在天线列表中不存在时，可用天线管理器进行编辑。5.坐标系统设置如果当前坐标系统中没有所需坐标系，则可以新建坐标系。新建坐标系的椭

14、球一般为BJ54或XA80，投影方式一般为高斯投影，关键是中央子午线与投影高程面。保存新建的坐标系6.静态基线解算6.1静态基线解算设置6.2处理基线 基线处理为自动过程，无需人工干预。6.3检查基线处理结果，处理解算失败的基线失败原因：Ratio值3; 固定解误差（RMS）大于相应等级的要求。处理方法：删除质量较差的观测值，根据残差图选定删除的数据。删除完后，重新处理选定基线，直到成功为止。6.4基线质量检核重复基线较差、异步环闭合差、同步环闭合差。主要是异步环闭合差，当异步环闭合差不满足条件时，要找出误差较大的基线，然后重新解算或者禁用该基线。6.5三维无约束平差主要用于检查质量较差的基线

15、X检验：先验精度与验后精度是否一致，通过调整比例因子。T检验：不通过时，说明该基线的误差还比较大，需进一步解算该基线7.网平差及高程拟合 网平差设置设置已知点二维平面已知点，输入XY, 高程拟合的水准点H, 平面与高程都已知的为XYH平差计算：当设置好后，处理为自动进行。平差后中误差超限的基线可以通过重新解基线提高精度，或删除基线。网平差平面点已知数据已知GPS点（54北京坐标系 高斯3度带投影 中央子午线114度）点号 X坐标 Y坐标等级备注GD01 .249.584D级可用0743 .975.474D级可用GE01.927.714E级可用高程拟合用已知水准高程点名H(m)备注GE0235.1041GE0333.4390GE0434.5178GE0545.1226GE0635.3699GE0738.9238GE0944.2127