简介RTKLIB是开源GNSS数据处理软件库这份压缩包提供的是作者修改后的精密单点定位PPP功能代码面向GNSS定位技术学习者和从事北斗、GPS数据处理研究的开发者。原始RTKLIB对北斗BDS和GPS的多频PPP支持有限该代码重点补齐了三频、双频及组合PPP解算能力便于用户利用北斗和GPS观测数据进行高精度定位实验。包内共37个文件以C语言源码为主包括36个源文件和1个头文件整体大小仅496KB覆盖RTCM解析、星历计算、PPP解算等关键模块结构清晰适合直接阅读和二次编译。该资源已有4497人浏览学习受到不少RTKLIB二次开发者的关注。读者能从这套代码中获得可运行的RTKLIB改进实现理解北斗与GPS组合PPP的算法流程包括L1/L2/L5三频信号处理、双频电离层改正、卫星钟差和大气延迟误差修正等通过对比官方原版代码还能深入掌握多系统精密单点定位的工程实现细节适合用于课程设计、算法研究或工程项目入门。 做高精度定位的人手里基本都有一份RTKlib源码。这两年问得最多的一个问题不是“RTKlib怎么装”而是“RTKlib的PPP到底能不能直接用北斗”。我自己的答案是能用但默认版本对北斗的支持非常勉强。如果你想实现“精密单点定位PPP并且同时支持北斗和GPS”那你大概率需要自己动手改代码。这篇我把自己实际改过的思路、踩过的坑、跑通的配置全部写出来给正在折腾的朋友一个参考。RTKlib本身是开源里做PPP非常成熟的框架但“支持北斗”和“支持好北斗”完全是两回事。默认版本里北斗的观测值映射、信号频率定义、码间偏差处理都不完整尤其是BDS-3的B1C、B2a这类新信号经常解析出来直接被丢掉。所以真正能落地跑双系统PPP的版本基本都是改过的。这篇文章就是围绕“改了什么、为什么改、怎么跑通”来展开。1. 为什么默认RTKlib跑不好北斗PPP1.1 先说清楚RTKlib的PPP能力边界RTKlib目录下有个RTKPOST是后处理工具PPP功能主要落在pntpos.c和pppos.c这两个文件里。它本身支持GPS、GLONASS、Galileo、QZSS北斗也有对应代码也就是SYS_CMP。从架构上说RTKlib设计得相当灵活你只要给它观测值、精密星历、钟差这些输入它就能在伪距和载波相位上做滤波解算。但实际用下来RTKlib对北斗的支持更像“能用但没优化”。我碰到最直接的问题同一组RINEX数据GPS参与解算一切正常加入北斗后要么卫星数量异常要么残差明显偏大。排除了数据源质量问题之后基本可以断定是代码层面的映射和处理不到位。这里的核心矛盾在于RTKlib 2.4.3这套代码诞生的时间比较早当时BDS-3还没大规模组网。所以它的信号定义表里北斗部分其实是按照BDS-2的B1I、B2I、B3I这些老信号来写的。现在采集到的数据里B1C、B2a、B2b这些新信号已经很常见默认代码根本没处理。1.2 北斗支持真实的短板在哪里我总结下来默认RTKlib跑北斗PPP有三个明显的短板第一信号与频率映射不完整。RTKlib在rtklib.h里用NFREQ定义参与解算的最大频率数代码里观测值映射函数obs2code()和code2obs()只覆盖了北斗少数几个信号。B1C和B2a这种新信号如果出现在RINEX文件里解析时会走default分支最后变成“未知码”参与不了解算。第二DCB码偏差没有内置修正方案。PPP解算用的是精密钟差这个钟差通常基于无电离层组合或者特定频率定义。BDS-2和BDS-3之间、不同频率之间的码间偏差如果不处理伪距观测值会带一个系统性的硬件延迟直接反映到定位结果上就是位置偏差好几米甚至更大。RTKlib本身没有一个自动下载和应用DCB文件的机制需要你在数据准备阶段自己处理。第三时间系统差别容易被忽略。北斗的时间基准是北斗时BDT和GPS时差了整整14秒而且BDT不做闰秒调整。RTKlib里虽然提供了bdt2gpst()这个函数但实际读取RINEX时如果接收机端或者转换程序没处理好这个偏移量会导致钟差匹配对不上解算结果自然是一团糟。所以所谓“更改后代码”核心就是把这些缺失的环节补上。补完之后GPS和北斗才能真正同时在同一个PPP滤波里互相配合。2. 改代码前必须搞懂的三个关键点2.1 PPP主流程和数据结构先看整体流程。RTKlib处理PPP时大致是读取观测文件得到obs_t结构体存放每个历元所有卫星的伪距、载波、多普勒等观测值读取导航文件和精密星历填充nav_t结构体然后进入pntpos()做过粗定位、选星、构建观测方程最后进pppos()做扩展卡尔曼滤波估计接收机位置、钟差、对流层湿延迟、电离层参数和模糊度。你在改代码时主要打交道的就是obs_t和nav_t这两个结构体。obs_t里的每个obsd_t有一条卫星的观测值里面有卫星编号sat、观测码code、伪距P、载波L、多普勒D等。nav_t里则存放星历、精密星历、电离层参数、UTC参数等。如果要支持北斗你需要保证两个层面的贯通一是观测值层面RINEX里的BDS信号能被正确解析到obsd_t里并且被当成“可用观测值”二是导航星历层面精密星历文件里的BDS轨道和钟差能被读进nav_t的peph结构体里。这两个环节任何一处断了北斗卫星就进不了解算。2.2 卫星系统和频率编号的对应关系RTKlib里卫星系统用一组宏定义比如SYS_GPS、SYS_GLO、SYS_GAL、SYS_CMP等。北斗那个宏在代码里是SYS_CMP不是很多人以为的SYS_BDS。这个细节很容易引起困惑因为网上很多资料直接写SYS_BDS但你拿这个宏去编译根本过不了。频率编号也有讲究。rtklib.h里freq_t数组定义了不同系统的频率值比如GPS L1是1575.42MHzBDS B1I是1561.098MHzB3I是1268.52MHz。观测值代码里比如GPS L1对应的码是“1C”BDS B1I对应“2I”B3I对应“6I”。修改代码时一个常见的任务就是在obs2code()里把北斗新信号和已有频率编号对应起来比如B1C对应L1频点、B2a对应L5频点。改代码前先在rtklib.h里确认你手上版本的NFREQ和信号定义表再对照你RINEX数据里实际出现的北斗信号类型。别上来就改否则很容易出现“编译过了但观测值全是0”的诡异情况。2.3 精密星历和DCB怎么进解算PPP定位必须使用精密星历和精密钟差。RTKlib支持读取SP3格式的精密轨道文件.sp3和钟差文件.clk这两个文件通常从IGS或者MGEX分析中心下载。北斗的精密产品现在公开渠道已经很多比如武汉大学、德国地学中心等机构发布的multi-GNSS产品里面就包含BDS-2和BDS-3的轨道钟差。但即便你准备了精密星历还是绕不开一个实际问题代码在peph2pos()里对每颗卫星找精密星历的位置如果星历文件里该卫星的位置不在参考时刻附近或者钟差类型不匹配RTKlib会选择插值失败然后跳过这颗星。所以精密星历和观测数据的时间跨度要匹配我一般会下载覆盖观测时段前后至少半小时的产品。DCB文件的处理则更麻烦。RTKlib默认并不读取常见的DCB格式。你需要自己做一次转换把DCB改正值加到伪距观测值上或者写成RTKlib能识别的格式。实际操作中我见过有人直接在观测值上做减法/* 以BDS B1I/B3I双频组合为例 */ P1_corrected P_B1I dcb_B1I_B3I * 0.5; P3_corrected P_B3I - dcb_B1I_B3I * 0.5;这个公式的含义是把码间偏差按频率关系分配到两个伪距上。具体DCB数值可以从CAS或者DLR提供的月文件中读取。这个步骤不做北斗的伪距残差会明显偏大最终位置结果很不可靠。3. 代码实际改造步骤3.1 第一步补齐北斗信号和观测值映射我用的源码版本是RTKlib 2.4.3。打开rtklib.c找到obs2code()函数这里有一个很大的switch语句。你会看到GPS、GLONASS、Galileo都列出了很多信号北斗部分却只有几行。我的做法是在BDS分支里补上B1C、B2a、B2b这些新信号和对应频率的映射case SYS_CMP: if (freq1) return 2I; /* B1I */ if (freq2) return 7I; /* B2I */ if (freq3) return 6I; /* B3I */ if (freq5) return 1X; /* B1C */ if (freq7) return 5X; /* B2a */ break;注意这里的freq编号必须和rtklib.h里定义的频率数组下标一致。改完obs2code()之后还要同步改code2obs()它是反查函数负责把RINEX文件里的观测码转换成内部的频率和系统编号。这两个函数必须成对修改不然解析链路会中断。另一个容易遗漏的地方rtklib.h里NFREQ和NFREQG这两个宏控制参与解算的频率数量。如果你想让B1C、B2a参与双频组合很可能需要把NFREQ从3改成更大的值同时检查所有用NFREQ做循环的地方比如观测方程里对频率的遍历。3.2 第二步开启GPS和北斗双系统联合解算代码层面确认了北斗观测值能被解析接下来就是让解算器真正把北斗卫星纳入观测方程。RTKlib里处理选项prcopt_t结构体中有一个navsys字段用来控制参与定位的卫星系统。我一般在代码初始化时写成prcopt.navsys SYS_GPS | SYS_CMP; /* GPS 北斗 */如果你用的是RTKPOST图形界面这一步对应勾选“GPS”和“BDS”。但注意很多人死在这界面勾选了BDS代码里却没把SYS_CMP加进去因为界面的选项和内部prcopt_t不一定是直接映射。所以建议你直接在代码里写死后编译这样最可控。选星逻辑同样要检查。RTKlib在selstar()函数里按卫星系统筛选可用星。如果navsys里没包含北斗后面定位方程根本不会遍历到北斗卫星。另外pppos()里的遍历范围是MAXOBS只要北斗观测值能正确进入obs_t这部分通常问题不大。3.3 第三步处理BDT和GPST的时间基准这一步是我最想强调的。很多改完代码跑出来的位置偏差几十米第一时间怀疑是代码逻辑问题其实只是时间系统没对齐。RINEX 3.03之后北斗观测数据的时间系统默认是BDT而精密星历和钟差一般使用GPST或者对应的系统时间。RTKlib在读RINEX文件时对BDS观测值会调用bdt2gpst()做一次转换把时间统一到GPST。但如果你用的数据源是某些接收机厂商自己的格式转换之后的RINEX文件里时间系统标签可能会混乱。我的经验是在处理完RINEX之后先抽查几个历元extern gtime_t bdt2gpst(gtime_t t) { return timeadd(t, 14.0); /* BDT比GPST快14秒 */ }确认转换正确了再往下走。同时精密星历文件SP3里的第一行有时间系统标识如果你发现下载的精密星历是BDT而你的程序按GPST处理也需要在读取时做相应偏移。这条做不对后面的收敛和精度无从谈起。3.4 第四步编译和常见编译错误RTKlib支持Windows和Linux编译。Windows下我用Visual Studio打开sln工程直接编译RTKPOST和RTKLIB库Linux下用make编译整个lib。改完代码后最容易遇到两类编译错误。第一类是宏定义或结构体字段名拼写错误。比如SYS_CMP写成SYS_BDS会报未定义标识符。第二类是数组越界警告。因为修改了频率定义和信号映射后某些循环的边界条件没跟上比如obs2code()里对freq的默认值返回NULL上层没有做空指针判断运行到北斗新信号时就崩溃或者跳过。我的建议是每一步修改后先编译一次库再编译调用程序确认无误后再继续下一步。不要一次改动太多否则排查成本很高。编译通过之后先用一份短时长的静态观测数据做验证确认北斗卫星能被识别、参与解算再跑长时间动态数据。4. 跑通PPP需要准备的数据与参数4.1 RINEX观测数据怎么来改好代码后数据准备的质量直接决定定位结果。PPP最少需要双频伪距和载波相位观测值。对于GPS至少要有L1和L2对于北斗至少要有B1I和B3I或者B1C和B2a。市面上常见的测量型接收机比如中海达、华测、司南等都支持输出包含北斗信号的RINEX 3.02以上格式。如果你手头没有接收机可以到IGS数据中心下载一些MGEX测站的静态观测数据来测试。我经常用BKG或者CDDIS下载某个测站的RINEX 3文件。注意下载时要选同时包含GPS和BDS数据的文件很多老测站只保留了GPS观测值。4.2 精密星历、钟差和DCB文件PPP需要的第二个数据源是精密产品。IGS的最终产品里北斗覆盖不够全所以一般推荐下载MGEX分析中心的产品。武汉大学的WUM产品、德国地学中心的GBM产品都是公开可下载的。这些产品里包含SP3精密轨道文件和CLK钟差文件。下载时重点看文件的星历覆盖时间段一定要覆盖你的观测时段否则插值外推会导致很大的轨道误差。至于DCB文件我一般在CAS的ftp下载BDS DCB月文件然后自己写个小脚本把需要的那天数据提取出来应用到观测值上。4.3 RTKPOST里最关键的一组参数如果你用RTKPOST跑后处理最影响结果的参数有这些Positioning Mode选PPP卫星系统勾选GPS和BDS其他系统保持默认也问题不大但初期建议先只做GPSBDS方便排查双频组合方式选L1L2北斗端也要对应改成B1IB3I或B1CB2a这取决于你的数据高度截止角静态数据设10度动态数据可以适当提高对流层延迟选Estimate ZTD这是PPP的常规做法电离层延迟选IFLC无电离层组合PPP收敛会快很多还有一个容易被忽略的参数观测值采样率。如果你手里是30秒采样率的静态数据没问题如果是1秒采样率的动态数据RTKPOST运行时对内存和处理时间要求会明显提高建议先用30秒数据做验证。5. 实战GPSBDS联合PPP跑通全流程5.1 完整跑一次的操作流程我以一组静态测站数据为例。数据是某个MGEX测站某天的RINEX 3观测文件时长24小时采样率30秒同时包含GPS和BDS观测值。精密星历用GBM的最终产品DCB用CAS当天文件。操作步骤大致如下把观测文件、精密星历、精密钟差、DCB文件放到同一个目录先用自己写的脚本对RINEX观测值做DCB改正生成一个新的RINEX文件。这一步我会在伪距上应用DCB改正载波相位不动打开RTKPOST观测文件选择改正后的RINEX精密星历选择SP3文件钟差选择CLK文件设置PPP模式勾选GPSBDS选择无电离层组合对流层选Estimate ZTD点击Execute开始处理后台会跑出.pos文件整个处理过程我这边大概需要几分钟。如果你的数据量很大要注意电脑内存RTKlib处理长时段数据时内存占用会上升尤其是双系统观测值多的情况下。5.2 结果怎么看位置、残差与收敛跑完之后打开.pos文件里面记录了每个历元的XYZ坐标、Q标志、NS参与解算的卫星数、GDOP等。我一般会正常化到经纬度和高程然后画时间序列图。判断结果是否正常的几个标准NS值在双系统下通常比单GPS多出5到10颗星GDOP值明显下降这是观测量几何结构变好的直接证明位置时间序列在静态场景下初始阶段波动比较大应该在30到60分钟内收敛到稳定值北东高三个方向误差双系统收敛后一般能到厘米级至少比单GPS快残差信息在RTKLIB的输出选项里可以打开最好也看一下。北斗伪距残差如果始终带着固定偏移基本可以确定DCB没改好回头检查数据预处理环节。5.3 我踩过的三个坑第一个坑是RTKPOST界面勾选了北斗但核心代码没加SYS_CMP。表现就是定位结果只有GPS卫星参与NS和单GPS完全一样。我后来习惯先看NS有没有增加再判断北斗是否真的参与了解算。第二个坑是SP3精密星历文件时间覆盖不够。我有一回下载了整周的SP3文件但RTKlib在插值某个时刻的卫星位置时如果文件里对应卫星的位置在前后各2小时窗口内缺失插值会失败。表现是某颗北斗卫星一直不进解算但GPS正常。解决方式是确保精密星历覆盖观测时段且轨道数据尽量使用最终产品精度和连续性都更好。第三个坑也是最隐蔽的RINEX数据本身从某些转换工具生成后BDS观测值的“code”字段和RTKlib内部的obsid不匹配。比如同样的B1I信号RINEX里标成“2I”但有些工具标成“2X”。代码里obs2code()返回的是NULL观测值直接被跳过。解决方法是打印一份观测值解析日志检查被跳过的卫星针对性补全映射。6. 常见问题排查速查表这里整理了我改代码和跑数据过程中频繁遇到的问题按“现象、原因、解决方式”整理成了速查表。现象可能原因解决方式北斗卫星完全不参与解算代码里navsys没加SYS_CMP检查prcopt_t初始化补上SYS_CMPNS比单GPS没增加北斗观测值在RINEX解析阶段被丢弃打开观测值解析日志检查code2obs()映射北斗残差有固定大偏移DCB码偏差未处理在伪距观测值上应用DCB改正位置结果初始偏差几十米BDT和GPST时间系统未对齐检查bdt2gpst()调用和时间标签某颗星经常不收敛精密星历插值窗口不够更换覆盖更全的SP3文件编译不通过提示SYS_BDS未定义宏名写错统一使用SYS_CMPRTKPOST界面卡死或内存暴涨观测文件过大或采样率过高先降采样到30秒再处理静态结果持续漂移不收敛对流层或电离层处理策略不对确认开启ZTD估计和IFLC组合排查技巧上我个人的习惯是“从头卡数据链路”。先确认RINEX能被解析出北斗卫星再看精密星历能不能给北斗卫星插值出轨道和钟差最后看解算结果有没有北斗参与。每一步都用日志或输出量来验证不要凭感觉判断。在处理过程中还要特别注意输出文件的“Q”标志。Q为1表示浮点解Q为5表示固定解。PPP通常是浮点解如果结果里Q标志一直不达标先别急着怀疑代码回顾一下数据质量和参数配置。7. 后续还可以往哪些方向扩展改完“GPSBDS”双系统PPP之后后续可以继续扩展的方向不少。比如加入Galileo和QZSS多系统联合后卫星数更多城市峡谷环境下定位连续性会好很多。再比如把代码从后处理挪到实时处理对接NTRIP协议或者串口数据流那就是一套实时PPP的原型。还有一个方向是模糊度固定的实现。RTKlib 2.4.3对PPP模糊度固定支持比较有限如果你有研究需求可以基于现有的双系统PPP代码把宽巷和窄巷小数偏差FCB估计加进去这是一条很值得深入的路。我个人在实际操作中的体会是RTKlib这套代码虽然看着老但胜在结构清晰、可控性强非常适合做研究和二次开发。改代码支持北斗这件事最难的不是算法本身而是把数据链路的每个环节都打通。正因为这样跑通第一份GPSBDS联合PPP结果的时候那种成就感确实很实在。希望这篇内容能帮正在折腾RTKlib的朋友少走点弯路。本文还有配套的精品资源点击获取
RTKLIB实现GPS+BDS双系统PPP的代码改造实战指南
简介RTKLIB是开源GNSS数据处理软件库这份压缩包提供的是作者修改后的精密单点定位PPP功能代码面向GNSS定位技术学习者和从事北斗、GPS数据处理研究的开发者。原始RTKLIB对北斗BDS和GPS的多频PPP支持有限该代码重点补齐了三频、双频及组合PPP解算能力便于用户利用北斗和GPS观测数据进行高精度定位实验。包内共37个文件以C语言源码为主包括36个源文件和1个头文件整体大小仅496KB覆盖RTCM解析、星历计算、PPP解算等关键模块结构清晰适合直接阅读和二次编译。该资源已有4497人浏览学习受到不少RTKLIB二次开发者的关注。读者能从这套代码中获得可运行的RTKLIB改进实现理解北斗与GPS组合PPP的算法流程包括L1/L2/L5三频信号处理、双频电离层改正、卫星钟差和大气延迟误差修正等通过对比官方原版代码还能深入掌握多系统精密单点定位的工程实现细节适合用于课程设计、算法研究或工程项目入门。 做高精度定位的人手里基本都有一份RTKlib源码。这两年问得最多的一个问题不是“RTKlib怎么装”而是“RTKlib的PPP到底能不能直接用北斗”。我自己的答案是能用但默认版本对北斗的支持非常勉强。如果你想实现“精密单点定位PPP并且同时支持北斗和GPS”那你大概率需要自己动手改代码。这篇我把自己实际改过的思路、踩过的坑、跑通的配置全部写出来给正在折腾的朋友一个参考。RTKlib本身是开源里做PPP非常成熟的框架但“支持北斗”和“支持好北斗”完全是两回事。默认版本里北斗的观测值映射、信号频率定义、码间偏差处理都不完整尤其是BDS-3的B1C、B2a这类新信号经常解析出来直接被丢掉。所以真正能落地跑双系统PPP的版本基本都是改过的。这篇文章就是围绕“改了什么、为什么改、怎么跑通”来展开。1. 为什么默认RTKlib跑不好北斗PPP1.1 先说清楚RTKlib的PPP能力边界RTKlib目录下有个RTKPOST是后处理工具PPP功能主要落在pntpos.c和pppos.c这两个文件里。它本身支持GPS、GLONASS、Galileo、QZSS北斗也有对应代码也就是SYS_CMP。从架构上说RTKlib设计得相当灵活你只要给它观测值、精密星历、钟差这些输入它就能在伪距和载波相位上做滤波解算。但实际用下来RTKlib对北斗的支持更像“能用但没优化”。我碰到最直接的问题同一组RINEX数据GPS参与解算一切正常加入北斗后要么卫星数量异常要么残差明显偏大。排除了数据源质量问题之后基本可以断定是代码层面的映射和处理不到位。这里的核心矛盾在于RTKlib 2.4.3这套代码诞生的时间比较早当时BDS-3还没大规模组网。所以它的信号定义表里北斗部分其实是按照BDS-2的B1I、B2I、B3I这些老信号来写的。现在采集到的数据里B1C、B2a、B2b这些新信号已经很常见默认代码根本没处理。1.2 北斗支持真实的短板在哪里我总结下来默认RTKlib跑北斗PPP有三个明显的短板第一信号与频率映射不完整。RTKlib在rtklib.h里用NFREQ定义参与解算的最大频率数代码里观测值映射函数obs2code()和code2obs()只覆盖了北斗少数几个信号。B1C和B2a这种新信号如果出现在RINEX文件里解析时会走default分支最后变成“未知码”参与不了解算。第二DCB码偏差没有内置修正方案。PPP解算用的是精密钟差这个钟差通常基于无电离层组合或者特定频率定义。BDS-2和BDS-3之间、不同频率之间的码间偏差如果不处理伪距观测值会带一个系统性的硬件延迟直接反映到定位结果上就是位置偏差好几米甚至更大。RTKlib本身没有一个自动下载和应用DCB文件的机制需要你在数据准备阶段自己处理。第三时间系统差别容易被忽略。北斗的时间基准是北斗时BDT和GPS时差了整整14秒而且BDT不做闰秒调整。RTKlib里虽然提供了bdt2gpst()这个函数但实际读取RINEX时如果接收机端或者转换程序没处理好这个偏移量会导致钟差匹配对不上解算结果自然是一团糟。所以所谓“更改后代码”核心就是把这些缺失的环节补上。补完之后GPS和北斗才能真正同时在同一个PPP滤波里互相配合。2. 改代码前必须搞懂的三个关键点2.1 PPP主流程和数据结构先看整体流程。RTKlib处理PPP时大致是读取观测文件得到obs_t结构体存放每个历元所有卫星的伪距、载波、多普勒等观测值读取导航文件和精密星历填充nav_t结构体然后进入pntpos()做过粗定位、选星、构建观测方程最后进pppos()做扩展卡尔曼滤波估计接收机位置、钟差、对流层湿延迟、电离层参数和模糊度。你在改代码时主要打交道的就是obs_t和nav_t这两个结构体。obs_t里的每个obsd_t有一条卫星的观测值里面有卫星编号sat、观测码code、伪距P、载波L、多普勒D等。nav_t里则存放星历、精密星历、电离层参数、UTC参数等。如果要支持北斗你需要保证两个层面的贯通一是观测值层面RINEX里的BDS信号能被正确解析到obsd_t里并且被当成“可用观测值”二是导航星历层面精密星历文件里的BDS轨道和钟差能被读进nav_t的peph结构体里。这两个环节任何一处断了北斗卫星就进不了解算。2.2 卫星系统和频率编号的对应关系RTKlib里卫星系统用一组宏定义比如SYS_GPS、SYS_GLO、SYS_GAL、SYS_CMP等。北斗那个宏在代码里是SYS_CMP不是很多人以为的SYS_BDS。这个细节很容易引起困惑因为网上很多资料直接写SYS_BDS但你拿这个宏去编译根本过不了。频率编号也有讲究。rtklib.h里freq_t数组定义了不同系统的频率值比如GPS L1是1575.42MHzBDS B1I是1561.098MHzB3I是1268.52MHz。观测值代码里比如GPS L1对应的码是“1C”BDS B1I对应“2I”B3I对应“6I”。修改代码时一个常见的任务就是在obs2code()里把北斗新信号和已有频率编号对应起来比如B1C对应L1频点、B2a对应L5频点。改代码前先在rtklib.h里确认你手上版本的NFREQ和信号定义表再对照你RINEX数据里实际出现的北斗信号类型。别上来就改否则很容易出现“编译过了但观测值全是0”的诡异情况。2.3 精密星历和DCB怎么进解算PPP定位必须使用精密星历和精密钟差。RTKlib支持读取SP3格式的精密轨道文件.sp3和钟差文件.clk这两个文件通常从IGS或者MGEX分析中心下载。北斗的精密产品现在公开渠道已经很多比如武汉大学、德国地学中心等机构发布的multi-GNSS产品里面就包含BDS-2和BDS-3的轨道钟差。但即便你准备了精密星历还是绕不开一个实际问题代码在peph2pos()里对每颗卫星找精密星历的位置如果星历文件里该卫星的位置不在参考时刻附近或者钟差类型不匹配RTKlib会选择插值失败然后跳过这颗星。所以精密星历和观测数据的时间跨度要匹配我一般会下载覆盖观测时段前后至少半小时的产品。DCB文件的处理则更麻烦。RTKlib默认并不读取常见的DCB格式。你需要自己做一次转换把DCB改正值加到伪距观测值上或者写成RTKlib能识别的格式。实际操作中我见过有人直接在观测值上做减法/* 以BDS B1I/B3I双频组合为例 */ P1_corrected P_B1I dcb_B1I_B3I * 0.5; P3_corrected P_B3I - dcb_B1I_B3I * 0.5;这个公式的含义是把码间偏差按频率关系分配到两个伪距上。具体DCB数值可以从CAS或者DLR提供的月文件中读取。这个步骤不做北斗的伪距残差会明显偏大最终位置结果很不可靠。3. 代码实际改造步骤3.1 第一步补齐北斗信号和观测值映射我用的源码版本是RTKlib 2.4.3。打开rtklib.c找到obs2code()函数这里有一个很大的switch语句。你会看到GPS、GLONASS、Galileo都列出了很多信号北斗部分却只有几行。我的做法是在BDS分支里补上B1C、B2a、B2b这些新信号和对应频率的映射case SYS_CMP: if (freq1) return 2I; /* B1I */ if (freq2) return 7I; /* B2I */ if (freq3) return 6I; /* B3I */ if (freq5) return 1X; /* B1C */ if (freq7) return 5X; /* B2a */ break;注意这里的freq编号必须和rtklib.h里定义的频率数组下标一致。改完obs2code()之后还要同步改code2obs()它是反查函数负责把RINEX文件里的观测码转换成内部的频率和系统编号。这两个函数必须成对修改不然解析链路会中断。另一个容易遗漏的地方rtklib.h里NFREQ和NFREQG这两个宏控制参与解算的频率数量。如果你想让B1C、B2a参与双频组合很可能需要把NFREQ从3改成更大的值同时检查所有用NFREQ做循环的地方比如观测方程里对频率的遍历。3.2 第二步开启GPS和北斗双系统联合解算代码层面确认了北斗观测值能被解析接下来就是让解算器真正把北斗卫星纳入观测方程。RTKlib里处理选项prcopt_t结构体中有一个navsys字段用来控制参与定位的卫星系统。我一般在代码初始化时写成prcopt.navsys SYS_GPS | SYS_CMP; /* GPS 北斗 */如果你用的是RTKPOST图形界面这一步对应勾选“GPS”和“BDS”。但注意很多人死在这界面勾选了BDS代码里却没把SYS_CMP加进去因为界面的选项和内部prcopt_t不一定是直接映射。所以建议你直接在代码里写死后编译这样最可控。选星逻辑同样要检查。RTKlib在selstar()函数里按卫星系统筛选可用星。如果navsys里没包含北斗后面定位方程根本不会遍历到北斗卫星。另外pppos()里的遍历范围是MAXOBS只要北斗观测值能正确进入obs_t这部分通常问题不大。3.3 第三步处理BDT和GPST的时间基准这一步是我最想强调的。很多改完代码跑出来的位置偏差几十米第一时间怀疑是代码逻辑问题其实只是时间系统没对齐。RINEX 3.03之后北斗观测数据的时间系统默认是BDT而精密星历和钟差一般使用GPST或者对应的系统时间。RTKlib在读RINEX文件时对BDS观测值会调用bdt2gpst()做一次转换把时间统一到GPST。但如果你用的数据源是某些接收机厂商自己的格式转换之后的RINEX文件里时间系统标签可能会混乱。我的经验是在处理完RINEX之后先抽查几个历元extern gtime_t bdt2gpst(gtime_t t) { return timeadd(t, 14.0); /* BDT比GPST快14秒 */ }确认转换正确了再往下走。同时精密星历文件SP3里的第一行有时间系统标识如果你发现下载的精密星历是BDT而你的程序按GPST处理也需要在读取时做相应偏移。这条做不对后面的收敛和精度无从谈起。3.4 第四步编译和常见编译错误RTKlib支持Windows和Linux编译。Windows下我用Visual Studio打开sln工程直接编译RTKPOST和RTKLIB库Linux下用make编译整个lib。改完代码后最容易遇到两类编译错误。第一类是宏定义或结构体字段名拼写错误。比如SYS_CMP写成SYS_BDS会报未定义标识符。第二类是数组越界警告。因为修改了频率定义和信号映射后某些循环的边界条件没跟上比如obs2code()里对freq的默认值返回NULL上层没有做空指针判断运行到北斗新信号时就崩溃或者跳过。我的建议是每一步修改后先编译一次库再编译调用程序确认无误后再继续下一步。不要一次改动太多否则排查成本很高。编译通过之后先用一份短时长的静态观测数据做验证确认北斗卫星能被识别、参与解算再跑长时间动态数据。4. 跑通PPP需要准备的数据与参数4.1 RINEX观测数据怎么来改好代码后数据准备的质量直接决定定位结果。PPP最少需要双频伪距和载波相位观测值。对于GPS至少要有L1和L2对于北斗至少要有B1I和B3I或者B1C和B2a。市面上常见的测量型接收机比如中海达、华测、司南等都支持输出包含北斗信号的RINEX 3.02以上格式。如果你手头没有接收机可以到IGS数据中心下载一些MGEX测站的静态观测数据来测试。我经常用BKG或者CDDIS下载某个测站的RINEX 3文件。注意下载时要选同时包含GPS和BDS数据的文件很多老测站只保留了GPS观测值。4.2 精密星历、钟差和DCB文件PPP需要的第二个数据源是精密产品。IGS的最终产品里北斗覆盖不够全所以一般推荐下载MGEX分析中心的产品。武汉大学的WUM产品、德国地学中心的GBM产品都是公开可下载的。这些产品里包含SP3精密轨道文件和CLK钟差文件。下载时重点看文件的星历覆盖时间段一定要覆盖你的观测时段否则插值外推会导致很大的轨道误差。至于DCB文件我一般在CAS的ftp下载BDS DCB月文件然后自己写个小脚本把需要的那天数据提取出来应用到观测值上。4.3 RTKPOST里最关键的一组参数如果你用RTKPOST跑后处理最影响结果的参数有这些Positioning Mode选PPP卫星系统勾选GPS和BDS其他系统保持默认也问题不大但初期建议先只做GPSBDS方便排查双频组合方式选L1L2北斗端也要对应改成B1IB3I或B1CB2a这取决于你的数据高度截止角静态数据设10度动态数据可以适当提高对流层延迟选Estimate ZTD这是PPP的常规做法电离层延迟选IFLC无电离层组合PPP收敛会快很多还有一个容易被忽略的参数观测值采样率。如果你手里是30秒采样率的静态数据没问题如果是1秒采样率的动态数据RTKPOST运行时对内存和处理时间要求会明显提高建议先用30秒数据做验证。5. 实战GPSBDS联合PPP跑通全流程5.1 完整跑一次的操作流程我以一组静态测站数据为例。数据是某个MGEX测站某天的RINEX 3观测文件时长24小时采样率30秒同时包含GPS和BDS观测值。精密星历用GBM的最终产品DCB用CAS当天文件。操作步骤大致如下把观测文件、精密星历、精密钟差、DCB文件放到同一个目录先用自己写的脚本对RINEX观测值做DCB改正生成一个新的RINEX文件。这一步我会在伪距上应用DCB改正载波相位不动打开RTKPOST观测文件选择改正后的RINEX精密星历选择SP3文件钟差选择CLK文件设置PPP模式勾选GPSBDS选择无电离层组合对流层选Estimate ZTD点击Execute开始处理后台会跑出.pos文件整个处理过程我这边大概需要几分钟。如果你的数据量很大要注意电脑内存RTKlib处理长时段数据时内存占用会上升尤其是双系统观测值多的情况下。5.2 结果怎么看位置、残差与收敛跑完之后打开.pos文件里面记录了每个历元的XYZ坐标、Q标志、NS参与解算的卫星数、GDOP等。我一般会正常化到经纬度和高程然后画时间序列图。判断结果是否正常的几个标准NS值在双系统下通常比单GPS多出5到10颗星GDOP值明显下降这是观测量几何结构变好的直接证明位置时间序列在静态场景下初始阶段波动比较大应该在30到60分钟内收敛到稳定值北东高三个方向误差双系统收敛后一般能到厘米级至少比单GPS快残差信息在RTKLIB的输出选项里可以打开最好也看一下。北斗伪距残差如果始终带着固定偏移基本可以确定DCB没改好回头检查数据预处理环节。5.3 我踩过的三个坑第一个坑是RTKPOST界面勾选了北斗但核心代码没加SYS_CMP。表现就是定位结果只有GPS卫星参与NS和单GPS完全一样。我后来习惯先看NS有没有增加再判断北斗是否真的参与了解算。第二个坑是SP3精密星历文件时间覆盖不够。我有一回下载了整周的SP3文件但RTKlib在插值某个时刻的卫星位置时如果文件里对应卫星的位置在前后各2小时窗口内缺失插值会失败。表现是某颗北斗卫星一直不进解算但GPS正常。解决方式是确保精密星历覆盖观测时段且轨道数据尽量使用最终产品精度和连续性都更好。第三个坑也是最隐蔽的RINEX数据本身从某些转换工具生成后BDS观测值的“code”字段和RTKlib内部的obsid不匹配。比如同样的B1I信号RINEX里标成“2I”但有些工具标成“2X”。代码里obs2code()返回的是NULL观测值直接被跳过。解决方法是打印一份观测值解析日志检查被跳过的卫星针对性补全映射。6. 常见问题排查速查表这里整理了我改代码和跑数据过程中频繁遇到的问题按“现象、原因、解决方式”整理成了速查表。现象可能原因解决方式北斗卫星完全不参与解算代码里navsys没加SYS_CMP检查prcopt_t初始化补上SYS_CMPNS比单GPS没增加北斗观测值在RINEX解析阶段被丢弃打开观测值解析日志检查code2obs()映射北斗残差有固定大偏移DCB码偏差未处理在伪距观测值上应用DCB改正位置结果初始偏差几十米BDT和GPST时间系统未对齐检查bdt2gpst()调用和时间标签某颗星经常不收敛精密星历插值窗口不够更换覆盖更全的SP3文件编译不通过提示SYS_BDS未定义宏名写错统一使用SYS_CMPRTKPOST界面卡死或内存暴涨观测文件过大或采样率过高先降采样到30秒再处理静态结果持续漂移不收敛对流层或电离层处理策略不对确认开启ZTD估计和IFLC组合排查技巧上我个人的习惯是“从头卡数据链路”。先确认RINEX能被解析出北斗卫星再看精密星历能不能给北斗卫星插值出轨道和钟差最后看解算结果有没有北斗参与。每一步都用日志或输出量来验证不要凭感觉判断。在处理过程中还要特别注意输出文件的“Q”标志。Q为1表示浮点解Q为5表示固定解。PPP通常是浮点解如果结果里Q标志一直不达标先别急着怀疑代码回顾一下数据质量和参数配置。7. 后续还可以往哪些方向扩展改完“GPSBDS”双系统PPP之后后续可以继续扩展的方向不少。比如加入Galileo和QZSS多系统联合后卫星数更多城市峡谷环境下定位连续性会好很多。再比如把代码从后处理挪到实时处理对接NTRIP协议或者串口数据流那就是一套实时PPP的原型。还有一个方向是模糊度固定的实现。RTKlib 2.4.3对PPP模糊度固定支持比较有限如果你有研究需求可以基于现有的双系统PPP代码把宽巷和窄巷小数偏差FCB估计加进去这是一条很值得深入的路。我个人在实际操作中的体会是RTKlib这套代码虽然看着老但胜在结构清晰、可控性强非常适合做研究和二次开发。改代码支持北斗这件事最难的不是算法本身而是把数据链路的每个环节都打通。正因为这样跑通第一份GPSBDS联合PPP结果的时候那种成就感确实很实在。希望这篇内容能帮正在折腾RTKlib的朋友少走点弯路。本文还有配套的精品资源点击获取