RTKLIB开发文档

1、软件开发文档软件名称RTKPOST项目组成员王会、刘琦起止日期2017年6月20日至2017年8月15日1.项目论证1.1 目标任务-2-1.2 项目背景及必要性1.3计划进展-2.需求分析 - *4 -3.软件设计3.1软件体系结构软件解算流程 软件设计模块 软件功能描述全局数据结构说明软件操作说明4.1 界面介绍 4.2软件安装4.3软件操作软件测试5.1 测试计划5.2 测试方法及工具 5.3测试用例- 9 - -11 -*15 -18 -18 - -19 - -19 -24 -24 - - *24 - -25 - 2.1 软件描述 2.2运行环境及外部接

2、口需求-25 - -27 -5.4测试结果 参考文献-1. 项目论证1.1目标任务根据比赛要求，经过小组讨论，我们选择了RTKLIB进行二次开发。 RKTLIB是全球导航卫星系统GNSS的标准&精密定位开源程序包，RKTLIB由日本东京海洋大学的高须知二开发。RTKLIB主要功能有:GPS, GLONASS，Beidou，(1 )支持多个GNSS系统的标准和精密定位算法，包括Galileo，QZSS 和 SBAS。(2 )支持多种GNSS实时和后处理定位模式：单点定位、DGP S/DGNSS,动态RTK、静态RTK、移动基站、PPP。(3)支持多种 GNSS 标准格式和协议：RINEX

3、2.10、RINEX2.11、RINEX2.12、RINEX3.00、RINEX3.01、RINEX3.02、RTCM2.3、RTCM3.1 、RTCM3.2、BINEX、NTRIP、NMEA0183SP3、ANTEX1.4、 IONEX1.0、NGS PCV、EMS 2.0。(4)支持多种GNSS接收机专有数据协议格式：NovAtel:OEM4/V/6，OEM3, OEMStar 、SuperstarII、Hemisphere 、Crescent、u -blox:LEA-4T/5T/6TSkyTraq、JAVADGW10-II/III 和 NVS。(5 )支持外部通信：Serial、TCP/

4、IP、NTRIP、本地日志文件(记录和播发)和FTP/HTTP。(6)提供许多函数库和 API (appi icati on p rogram in terfaces):卫星和导航系统函数、矩阵和向量函数，时间和字符串函数、坐标的转换，输入和输出函数、调试跟踪函数、平台依赖函数、定位模型、大气模型、天线模型、地球潮汐模型、大地水准面模型、基准转换、RINEX函数、星历和时钟函数、精密星历和时钟、接收机原始数据函数、函数、流数据输入和输出函数、整周模糊度解算、标准定位、精密定位、后处理定位(解算)RTKLIB同样也有一些缺陷，RTKLIB在数据预处理（粗差探测），周跳探测方面的算法有些问题，RT

5、KLIB的源代码并不支持北斗系统的解算。同时,RTKLIB的框架过于过时，它的GUI是用的早期Borland C+ 设计，无法满足跨平台需求，基于这些问题，我们小组经过讨论，决定对这些问题进行修复，改进单频解算算法，同时基于Qt5对GUI层进行重新设计。1.2项目背景及必要性中国北斗卫星导航系统是中国自行研制的全球卫星导航系统。是继美国GPS、俄罗斯GLONASS之后的第三个成熟的卫星导航系统。目前北斗系统已具备覆盖亚太地区的定位、导航和授时以及短报文通信服务能力，预计2020年建成覆盖全球的北斗导航系统。目前多数接收机并不存在单纯的使用GPS进行卫星定位，而是在支持GPS的基础上还加入了 G

6、LONASS以及北斗系统。接收机在定位时可以同时使用这两大卫星定位系统，这样就GPS进行定位的可以在某个卫星定位发生问题时，使用另外一个系统进行定位。或者在使用 同时，再使用 GLONASS或者北斗来进行辅助，以此来增强定位的速度及准确率。另外，众 多的个人或行业应用中，定位都成为一种必不可少的基础功能，个人用户的社交信息、行程足迹；行业中的车辆管理、物流监控等无不需要卫星定位功能。RTK定位是北斗导航系统应用的一个重要方面，基于GPS的RTK算法已经十分成熟,各种商业软件的推出也推动了GPS RTK算法的完善和发展。此外由于北斗导航系统现在 已经有十几颗可用卫星，基于北斗系统的应用和研究也已

7、经在全国乃至世界范围内展开, 所以在研究 RTK的算法时，有必要考虑将 BDS加入算法中去，基于 BDS /GPS的RTK定 位算法研究及软件编制显得十分必要。1.3计划进展表1.1计划进展日期说明6月15日明确任务以及了解整个比赛进度安排6月16日至30日在RTKLIB上加入 DCB改正7月01日至20日在RTKLIB上加入单频周跳探测7月21日至8月01日在RTKLIB上加入数据预处理，修正周跳探测、添加载波平滑伪据、粗差探测7月02日至8月05日在RTKLIB上加入北斗系统，修正RTKLIB 一些现存问题8月06日至10日对RTKLIB的GUI进行重写，并在 Windows 和Ubunt

8、u 上进行测试8月11日至15日数据测试，编写数据精度报告2. 需求分析1)2.1软件描述序号功能名称功能需求标识简要描述01RTK解算rtkpost单频rtk解算02RTK成图rtkp lot对解算结果成图03RTK查看rtkview查看解算结果04RTK平滑rtksmooth对基准站进行平滑表2.1软件模块2.2运行环境及外部接口需求(1 )硬件要求程序硬件要求:1)具有In ter 13 处理器且满足以下要求的计算机:2)最低512 MB 内存、最小 1GB 硬盘、鼠标、键盘(2 )软件要求操作系统:Win dows 7 Win dows 8 Win dows 10 Ubu ntu 16

9、.04(3 )用户接口采用Qt通用图形界面，用户友好。2)界面具有一致性，界面规范遵循软件界面的规范。3)提供错误信息。4)遵循国家关于计算机方面词汇的标准，用词正确、准确、无歧义。5）本系统的用户一般需要通过中断进行操作，进行主界面后点击相应的窗口，分别进入相对应的界面（如，输入界面、输出界面）。3. 软件设计3.1软件体系结构RNEX读取一粗差探测解算DCB0E1EJ模糊度固定J卡尔曼谑波r波动图轨迹图绘图查看图3.1软件体系结构3.2软件解算流程1流程图中*意味着对RTKLIB原有的模块进行了改进或为新增模块图3.2结算流程图3.3软件设计模块RTK定位解算、结果软件主要包括四个模块：读

10、取文件、伪距单点定位、数据预处理、 输出。(1 )读取文件其中RINEX观测文件RTKLIB已经有现成模块，但是RTKLIB 2.3.4源码无法读取北斗观 测文件，我们在 RTKLIB基础上通过修改频率优先级增加了读取北斗观测文件功能。同时，基 于比赛广播星历数据，我们特地针对北斗星历增加了读取北斗广播星历函数。(2)伪据单点定位先进行伪据单点定位，求取平均值作为基站近似坐标，同时得到卫星坐标以及钟差信息, 以便组成双差观测方程。(3)数据预处理数据预处理阶段，我们着重做了一下三个工作：粗差探测、周跳探测、载波平滑伪距。经过我们测试，RTKLIB中的周跳探测并不完善，因此我们在解算前再加入一次

11、周跳探测。首先由L1、L2相位观测方程可得到 GF组合(消去了与频率无关的几何距离和钟差等)，由于相邻历元电离层变化很小，因此可忽略相邻历元差后的电离层残差。同时，在组成MW组合, 对MW 组合观测值做历元差分，以 3倍中误差为限，与 GF组合联合，通过求解方程组得到 精确周跳。RTKLIB中并没有单频周跳探测功能，于是我们增加了基于流动站的单频周跳探测功能。对每个历元k,卫星s，进行探测：若数据缺失超过阈值，标记周跳。若没有超过，则逐个计 算伪距与载波之差:d(s;k)md(s;kmd2(s;kR(s; k) (s;k)1 N1)- d(s；k 1)N i 11 N 21)- d2(s;k

12、1)N i 1Sd(s;k1) Jmd2(s;k 1) m：(s;k 1)若错误!未找到引用源。d(s；k) mj(s；k 1) nTSd(s；k 1)则标记周跳。GPS定位时，利用码观测在修复完周跳后，我们进行了粗差探测以及载波平滑伪据。在从而取出伪距粗差，提高伪值得到的结果精度较低，可以利用载波相位的变化来平滑伪据, 据观测定位精度。(4)RTK定位解算在解算过程中，我们基于RTKLIB添加了 DCB改正。GPS广播星历中卫星钟差参数是由L1P(丫)和L2P(Y)消电离层组合计算得到的，该参数中含有L1P(Y)和L2P(丫)在卫星端的硬件延迟影响。为扣除这一偏差对导航用户的影响,GPS引入

13、了TGD参数扣除。GPS卫星在发射前都会在地面对卫星 TGD参数进行提前标定。1999年4月后，GPS开始播发由JPL基于IGS基准站数据处理得到的TGD参数，且每四个月左右更新一次。2005年起发射的 GPSBlock IIR-M 卫星在L2频率上增加了 L2C民用信号，2010年起发射的 GPS Block IIF 卫星新增了 L5频率以及L5I5、L5Q5两种民用信号。与之对应，GPS在L1C/A、L2C、L5I5及L5Q5民用信号上相对于 L1 P(Y)分别定义了一个ISC 参数，即 ISCC/A、ISCL2C、ISCL5I 及ISCL5Q。ISC参数本质上也是 DCB参数，其中ISC

14、C/A是L1P(丫)与L1C/A 信号间的频内偏差参数,ISCL2C、ISCL5I 及 ISCL5Q 分别是L1P(Y)与L2C、L5I5、L5Q5 信号间的频间偏差参数。TGD也可以看作是 L2P(Y)信号上的ISC参数，不过要乘一个与频率有关的常数进行转换。考虑到广播的TGD及ISC参数与DCB参数之间的转换关系，GPS单双频用户也可以利用后处理的DCB参数改正硬件延迟误差后进行定位。clk L1P L2PflL1P i,x1 h'f2'2(5)结果输出基于RTKLIB增加了卫星系统输出，以及参考框架。最终输出还包括时间，状态，三维坐标，精度指标等内容。3.4软件功能描述表

15、3.1 POST模块功能项解算标识P OST功能描述对输入的观测文件进行解算使用角色用户主要事件用户操作系统响应事件流程单击Option按钮设置解算参数单击View按钮查看解算结果单击Plot按钮绘制波动图或定轨图单击Exec按钮开始解算单击Exit按钮退出前置条件用户需要设置 Op ti on参数，以及选择相应文件后置条件生成解算文件(*.pos)接口说明输入输出调用关系调用模块rtkp OS.C PP rin ex.c pp SatP ass.c pp DCP ass.c pp被调用模块rtkwork.c pp表3.2 View模块功能项查看标识View功能描述对输入的观测文件进行查看使用

16、角色用户系统响应事件主要事件用户操作流程单击Refresh按钮重新加载解算文件单击Read按钮重新选择解算文件前置条件单击Exit按钮退出后置条件输入接口说明调用关系输出调用模块pl otw in dow.c pp被调用模块post表3.3 Plot模块功能项绘图标识Plot功能描述对输入的观测文件进行成图使用角色用户主要事件流程前置条件后置条件接口说明调用关系用户操作单击Open按钮选择解算结果文件单击Refresh按钮单击Position按钮单击Gnd Trk按钮单击Save按钮单击Exit按钮系统响应事件重新绘制绘制波动图绘制轨迹图保存退出用户需要选择解算结果，并选择绘图类型输入输出调用

17、模块qcusto mpl ot.c pp被调用模块rtkwork.c pp功能项选项功能描述设置解算策略使用角色用户主要事件流程前置条件后置条件接口说明调用关系用户操作单击OK按钮单击Cancel按钮表3.4 Option 模块标识Op ti on设置并保存解算参数退出用户需要设置 Op ti on参数，以及选择相应文件生成解算文件（*.pos）输入输出调用模块系统响应事件op ti on dlg.c pp controlpl us.c pp被调用模块post表3.5 Smooth 模块2我们将该功能单独剥离了出来，可以在 bin目录下直接找到。前置条件-in file <in file

18、> -outfile <outfile> -dt <dt>后置条件生成平滑后的观测文件（*.*o）接口说明输入输出调用关系调用模块smoothe.c pp被调用模块post3.5全局数据结构说明除了 RTKLIB原有文件外，我们新增加了5个头文件:表3.6新增文件序号头文件实现文件功能01MathBase. h常用排序函数，以及统计函数02Sat Pass.hSatP ass.c pP将观测数据按照卫星编号分类存储，可用于载波平滑伪据03Segme nt.h存储了 Sat Pass中卫星历元连续信息，用于探测周跳04DCP ass.hDCP ass.c pPM-W

19、探测周跳，修复周跳05Stats.h常用数学函数，最大值，最小值，方差等软件中全局数据结构如下表所示:表3.7全局宏定义序号宏定义所属文件说明01PIrtklib.h圆周率02CLIGHTrtklib.h光在真空中的速度03OMGErtklib.h地球自转角速度04FEQ1rtklib.hGPS L1频率05FEQ2rtklib.hGPS L2频率06FEQ1_C MPrtklib.h北斗B1频率07FEQ2_C MPrtklib.h北斗B2频率08FEQ3_C MPrtkiib.h北斗B3频率09D2Rrtkiib.h度转弧度10R2Drtklib.h弧度转度11WLNWINDOWS*DCP

20、 ass.c pPWL周跳探测滑动窗口宽度12MAXWINDOWSIZE*rtkpo s.c pp单频周跳探测最大滑动窗口宽度 *表明新增宏，开发过程中我们尽量不在添加新有宏定义，所以开发过程中我们尽量沿用了原有RTKLIB中的宏定义。同时为了增加代码兼容性和可读性，我们没有增添新的结构体，只对少部分的结构体做了改动。4. 软件操作说明4.1界面介绍图4.1为软件主界面图，其中从上到下分别是流动站文件，参考站文件，广播星历文件, 精密星历文件，精密种差等 IGS文件，DCBDCB 文件可以在 IGS 官网下载： ftp：/cddis gsfc nasa gev/pub/gps/prcducts

21、/mgex/dcb文件以及解算结果文件。其中参考站文件，基准 站文件和广播星历文件为必选，其余的为可选项。针对这次的比赛数据，我们特地为北斗广 播文件增设了一个选择框，并且在源码中进行了识别。承 RTKLIBPOSTRINEX OES: Rover及 卓面单iRTEttBiJg2O16Og03-2tov6r. ISORINEX OES: easeE：、卓面单iRTEtBS!(tg20150903-2base. 16oRINEX MAV / IGS Product: SP3,CLK,lONEXE: 卓面、单iRTKttitl20160903-2rover. ISHE: '卓面、单iRTK

22、tt；l!tl20160903-2rover. ISCOCB FILE tlg2016942C£StG)(RAF_2(HE252(lD_CllD_lD_DCE.HSX.delationE:臺面、单频肮时匕務数壮PlotViewOptions ExecuteExit图4.1 RTKPOST 软件主界面图主要功能按钮介绍:Execute:执行解算gtlO":设置解算参数View:查看解算结果Plot:绘图4.2软件安装默认我们在bin文件夹下已经生成了 windows下的exe可执行程序，可直接使用,若在LinUX下则需要使用qmake对gui文件夹下的重新编译。

23、4.3软件操作以单频RTK比赛数据文件中的 20160909-2 下的数据为例:>> 宝31 * 单！RFilUtlSlP抵 2016Q«B-2LtP卜駅hUSf b日吕已16o2DlD/9d2 15=14163 TSrf牛"7004 KBn DSQMXRA己刃厲2TSTOQ0 OWXW pca esx1&5KB盍1 九1 1轧好幅re*1 rcnft-r.1-&N扣1 钠ZIP 1" 口IfiN忘性11 tJi=H片*rcJ¥tr.lli&0Ki： ) vli4U4 航图4.2演示用例1 ）首先点击Options按钮

24、，如图4.3所示，Settiing1选项卡中,In terval为采样10，lonosphere率，设置为1s，Elevaton Mask 为截止高度角（单位为度）这里设置为Correction 与Troposphere Correction为电离层与对流层模式，由于只是演示，这里 分别选择Broadcast和Saastamoinen ，卫星系统选择GPS和北斗。Setti ng2 选项卡中，Freqique ncies为频率选择，虽然流动站为单频，但是基准站为双频，所以还是选择L1+2，在解算过程中会程序会自动识别。Integer Ambiguity ResON。为模糊度固定方式，由于是连续

25、观测，这里选择Fix and Hold ，MinRatio to FixAmbiguity 为ratio值，默认为3 , Smooth为是否平滑以及粗差探测，这里选择对于基准站坐标，这里选择Average of Si ngle Position，即单点定位平均值，若已知基 准站坐标，则可以选择 X/Y/Z-ECEF(m)模式，从而输入基准站坐标。Out put中主要设置了输出选项。这里Solution Format输出坐标格式，选择E/N/U-Baseline，Time Format 为时间格式，选择 hh:mm:ss GPST ，Latitude/L on gitude Format为基准站

26、坐标格式，这里选择ddd.dddddd，其余的保持 默认即可。罔 Optiorc?XOutput1 line itarC tfST"圧fU.l/X:U UU Ulj«匚 T加 EiiJ（CT£T）1370/1/1:0 00-00J.rFositioi KcdaLEinaraticTInterval【：）1ZT _EliVtti JK20* _Tft7inipK»r*jflnTS-E.utuT oiitr.TEttdliti Iphntris/CloclfracistTt4 St4r>回阳s0日i©*11CttBLCel1衆 Option?

27、XStttinilEtttiftsSOutputFrei qucni czLL吃Intec*r口ity £«sFLk mH Ho：才Tin R-ti t t硏 riw Knti1?. 00：Srr*ocLUHrPas Sliti orI薦 QptiorScttincl卫农 O-c+jHatSolutioo F>rnfi.tTift* PcruiTLit 3 烈«4九宾日 trjJo ? arm4±U'ittnR；*i渔丰T1 ajL：±1 de 1 1,?XlA-tCEFwIhrn := CP STrUddcJd沁甜Pllipic

28、AlCFFT4t.UK I Cwc»l，r图4.3参数设置2 ）按照如图4.4所示选择解算文件。5图4.4解算文件选择3）上述步骤操作完成后，点击Execute，进行解算，解算完成后，点击View可以查看解算结果。5 DCB改正同时需要 CLK文件，CLK文件可以在IGS的FTP上下载：/JK.SMiUTKiaJO I O*3fr-2rMr.=.Di jib JI IV4诚» -珂 c-ViflitFiht PE帝曲C vtlvbi 4E. !<* F vil 如"伽* iLitOift I阳邛屮 1J9D対止 GK* L 尸

29、 > ull J 1V <hT tlTf = Swu tfiEvJAk； ryi fy-bp-t1 I ph piiTun VH : 4 -占八送耐冲埼not弼IS厲朋吟甘丁叶刁佃B '令血单:KIt|tS3aS'2«*lS«-i',b4».I ' *£ !p fl1trihl¥tttau».r»<r tu :£:总£'中电TZttFti后魄心DHRY Hku. tCCBowwue « s M 0 朋T <*Jd9n "如

30、Q u,i 应 Da oa E? Q OST (”亠帕：卫gk I b-PBati dlpuf-in B-d.ID I d列oEf也 h.Hr*tS-U Lt4r>>Jl14 Dr*dflwtIPI b.ilaD.£j± All 1 bilJ.a 0iiit|.IKS»9»D h£Lt a du (D.OQDQa血团亡b Mm D CKQ 0 KQMiL (?/fc<i-fcnHinL*CC5* l】 Im? Elbt.O.sLfa.-l dcpi6 xiih * TPI-cls 涎 *1“dMpSb HwyliAlrJ u吗

31、uJiah； 中 r>-54t±llitM 门<,1 D1駆備ZCJtjdP.'te D& £9 ZB.CIDI-45C1-aE3a=iT2-ED.SflJBD. riKJ-60MOft 汕 IS.OPU1上1-0買护ft4 »320 K岀11拠机iE/g/be U F5 ；9.KD1 44*1-aPTLB=499ED迂吕d口 4lEBD S36E=门氐B/n DS Z9.3D.CC-DJ-*他-aX.eV kBI?0. JbU4D.kkJ；.ifldiiosyit 04 PC M o1和H-0糾阳L.0A涯諒0势ft请2讹阻/OE 昭

32、 23 JZ.lZDDI.B3L-ai£DD冋4D.efiSE=址伽牺 K SP )1. UD.mi-0net*tr±LHo.iiaO.tJV25恥轲DP能出3弓.WD1祁西-0廖tc旧审0林0 6?5tMigju (A se ts.aa1 ELXtifl? Jl»0 ：4UD tlM2«6X咯诟 P& EB Ip.nu-uTW5>fl HTQUWQ创舊E(H64»/De K 四 PLffiQ1.E&0-0吨?eC ?CMo.aa?2o.jraU ：6 )1. »hk 啟ld«?t< dl)d0

33、A用Q 2倍血mfg D&.£3；i3.KPl.gET5-QPKfi£EC.15SD.P.UILisesaw汕 a4fT图4.5查看解算结果）在主界面点击Plot，可以查看绘图结果，如图 4.6如同肛®D.3jra-n. 51DCa DO0 2迥-fl. JH>Jftj 1也-n 79 BCa DOO.JFZ-0. _iWC DU0 1卿-0眷HCl an0 1坯L-D. idia0 DO0. kjui-aaptre DU0卄冷-D. LBd1D WD LlM-fl. nilt U0汁聊-q. WWG m0.沪旳-D. leei0 mn QfeH-n

34、 MUF mD.CM-D.衣L Wrlkt1-213I 2I 9 ML 豐血-QI-itbrtlDb I所示，右下角显示的是固定率。5. 软件测试5.1测试计划表5.1测试计划测试测试类型负责人时间1算法测试小组所有成员08月11日2界面测试小组所有成员08月13日3验收小组所有成员08月14日5.2测试方法及工具测试方法：在 Visual Studio 中编译改进后 RTKLIB源码，测试所有数据，同时在 Qt中测试界面逻辑。工具：Visual Studio 2013, Qt Creator 4.3 + Qt 5.95.3测试用例I '聞旦脑，空直b早领盯Ktt績说辱T应曰期<2011 &0逊-2£07/6/13 14:30丈肛201&CS09-1201 /e/10 14:3/文件衰J201讯9网-220-7/e/n i&：d7丈件突_ ?01 fXWU-lifi-i7-201GWH-220UZC/20 帖29VJoiMyiJU 1i5：jy201 M90a-21 201510252017/6/28 15:29丈件夹bin20151102的7 曲RT、;7q=(Dj2015122120打必，8 15：四20170220 02017/6/28四20170220-10