1. 初识TEQCGNSS数据预处理的“瑞士军刀”如果你刚开始接触GNSS数据处理面对一堆从接收机里导出来的、后缀名稀奇古怪的原始数据文件是不是有点无从下手别慌这几乎是每个GNSS工程师或研究生的必经之路。今天我要跟你分享的就是我在处理这些数据时最常用也最信赖的一个“老伙计”——TEQC工具。你可以把它想象成GNSS数据预处理领域的“瑞士军刀”虽然界面看起来就是个黑乎乎的终端窗口但功能却强大得惊人尤其是做RINEX格式转换和数据质量检核效率非常高。TEQC的全称是“Translation, Editing, and Quality Checking”翻译过来就是“翻译、编辑和质量检查”。这个名字非常直白地概括了它的三大核心功能。它最初是由UNAVCO一个地球科学领域的大学联盟开发并维护的在科研和工程领域应用了二十多年经久不衰。它的主要工作就是帮你把各家接收机厂商的私有数据格式比如天宝的.T02、徕卡的.M00、诺瓦泰的.obs等等转换成国际通用的RINEX格式。RINEX格式就像是GNSS数据界的“普通话”只有转换成它后续的各种高精度解算软件比如GAMIT/GLOBK、Bernese、RTKLIB才能顺利“听懂”并进行处理。除了格式转换TEQC更强大的地方在于它的质量检核QC功能。很多时候我们从野外辛辛苦苦采集回来的数据里面可能藏着各种“暗伤”比如卫星信号被树木或建筑物遮挡导致数据中断接收机偶尔“抽风”产生跳变或者电离层扰动引起周跳等等。如果带着这些问题数据直接去做解算结果很可能不可靠。TEQC的质量检核就像给数据做一次全面的“体检”它能生成详细的报告和图表直观地告诉你数据哪里好、哪里有问题比如卫星的可见性、信噪比变化、多路径误差等等。这对于评估数据是否可用、定位问题原因至关重要。接下来我就手把手带你从零开始搞定TEQC的安装和最基本的操作。2. 环境准备与第一步获取并运行TEQC2.1 获取TEQC可执行文件TEQC本身是开源免费的但官方并不直接提供编译好的Windows安装包。最可靠的获取方式是访问UNAVCO的官方网站。不过对于初学者更简单的方法是直接搜索“teqc executable download”通常可以找到一些大学或研究机构分享的已编译版本。你需要根据你的操作系统Windows、Linux或macOS下载对应的可执行文件。对于Windows用户下载下来的通常就是一个名为teqc.exe的文件。这里要特别强调一点TEQC不是一个有图形界面的软件它是一个纯粹的命令行工具。这意味着你别指望双击teqc.exe就能打开一个窗口点点鼠标。我刚开始用的时候也犯过这个傻双击后一个黑框一闪而过啥也没发生。正确的方式是需要在系统的命令行终端比如Windows的CMD或PowerShellLinux/macOS的Terminal里来运行它。你可以把teqc.exe文件放在任何一个你喜欢的文件夹比如D:\GNSS_Tools\。为了方便在任何目录下都能调用它我强烈建议你将这个目录添加到系统的环境变量PATH中。如果觉得配置环境变量麻烦也可以直接在命令行里切换到teqc.exe所在的目录进行操作。2.2 处理前的关键一步RINEX版本降级拿到TEQC你可能迫不及待地想处理手头的数据。但这里有一个常见的“坑”目前广泛使用的TEQC版本主要支持RINEX 2.xx格式的数据。而如今很多新型号的GNSS接收机默认输出的已经是RINEX 3.xx格式了。如果你直接拿一个3.02版本的观测文件后缀通常是.21o给TEQC处理它很可能会报错或者无法识别。所以预处理的第一步往往是把RINEX 3格式降级到RINEX 2格式。这就需要用到另一个小工具比如原始文章里提到的gfzrnx。这是德国地学研究中心GFZ开发的一个非常高效的RINEX处理工具转换速度很快。它的基本命令格式很简单.\gfzrnx_1.15-8044_win10_64.exe -finp 输入文件.21o -fout 输出文件.21o --version_out 2这条命令的意思就是读取输入的RINEX 3观测文件输出一个RINEX 2.11版本的观测文件。那个--version_out 2参数就是指定输出版本为2。原始文章命令里还有一个-no_prn C38,C39...参数这是用来在转换时排除指定的北斗卫星PRN号。这是因为在某些早期版本的解算软件中可能不支持所有北斗卫星排除它们可以避免兼容性问题。你可以根据自己后续处理软件的要求来决定是否使用这个参数。导航星历文件广播星历也需要进行同样的版本转换它的后缀通常是.21g混合系统或.21nGPS。命令是类似的.\gfzrnx -finp BRDC00IGS_R_20211550000_01D_MN.rnx -fout MA36-x33.21g --version_out 2 -satsys C这里的-satsys C是指定处理北斗系统的星历。完成这一步我们才得到了TEQC能“吃”得下去的RINEX 2格式数据接下来就可以请TEQC正式登场了。3. 核心实战一使用TEQC进行数据质量检核QC质量检核是TEQC的看家本领也是我们判断数据好坏的第一道关卡。运行一次完整的QCTEQC会生成一份文本报告和一系列用于绘图的数据文件让我们对数据质量有个全局的把控。3.1 执行基础质量检核最基本的QC命令格式如下.\teqc qc 输入文件.21o运行这个命令后TEQC会做以下几件事自动寻找导航文件它会在观测文件.21o所在的同一目录下自动寻找同名的导航文件.21n或.21g。如果找不到QC报告里关于卫星轨道和位置的信息就会缺失。在屏幕上滚动输出终端里会快速滚动显示检核的概要信息包括数据时长、采样间隔、观测到的卫星系统等。生成QC报告文件最重要的输出是一个后缀为S的文件例如输入文件.21o.S。这是一个纯文本文件里面包含了极其详细的统计信息比如每个卫星每个频段的观测值数量、信噪比SNR的平均值和标准差、周跳和失锁标记的数量、多路径误差的估计值等等。生成绘图数据文件同时还会生成一系列后缀为M多路径、N信噪比、O观测值、I电离层等的文件。这些文件是专门给另一个叫teqcplot的工具通常和TEQC一起发布来绘制质量检核图的。3.2 卫星系统与特定卫星的筛选在实际处理中我们可能只关心特定卫星系统的数据。比如我的项目只使用GPS和北斗那么我可以排除GLONASS、Galileo等系统的数据让报告更简洁处理速度也更快。TEQC提供了非常灵活的筛选选项-G排除所有GPS卫星-R排除所有GLONASS卫星-E排除所有Galileo卫星-C排除所有北斗BDS卫星-J排除所有QZSS卫星-S排除所有SBAS卫星这些选项可以组合使用。例如如果我只想分析GPS和北斗的数据命令可以这样写.\teqc qc -R -E -J -S 输入文件.21o这条命令排除了GLONASS(-R)、Galileo(-E)、QZSS(-J)和SBAS(-S)只留下了GPS和北斗的数据进行检核。更进一步我们还可以排除特定的卫星。比如我知道GPS的PRN 22号卫星在那天的健康状况不佳或者北斗的PRN C01卫星被建筑物长期遮挡我想在检核时忽略它们。命令如下.\teqc qc -G22 -CC01 输入文件.21o甚至你可以排除一个范围内的卫星比如排除GPS PRN 1到15号的所有卫星也许这个区间卫星的仰角一直很低.\teqc qc -G1-15 输入文件.21o3.3 解读QC报告的关键指标生成了.S报告文件后怎么读呢我挑几个最常用的指标说说mp1和mp2这是估算的L1和L2频段的多路径误差单位是米。一般来说这个值越小越好通常低于0.5米可以认为是质量很好的数据大于1米则说明多路径效应比较严重可能站点周围有反射面。sn1和sn2L1和L2的信噪比dB-Hz。信噪比高且稳定说明信号质量好。如果某颗卫星的信噪比曲线像锯齿一样剧烈波动很可能信号受到了干扰或遮挡。和这是数据完整性的一个粗略指标表示“观测值/期望值”的百分比。它会根据卫星的可见性理论值来估算你的接收机实际捕获到了多少数据。这个值越高说明数据中断越少。slips周跳的数量。周跳是载波相位观测值发生整周跳变的现象太多周跳会增加数据处理的难度。报告里会按卫星列出周跳数。看报告不能只看平均值一定要结合teqcplot生成的图表来看趋势。比如多路径误差在某个时间段突然增大很可能对应了车辆或行人从天线附近经过。信噪比在某个高度角区间系统性降低可能意味着有固定的遮挡物比如树梢。4. 核心实战二使用TEQC进行数据编辑与重采样QC告诉我们数据哪里有问题而TEQC的编辑功能则允许我们对数据进行一些简单的“手术”来优化它使其更适合后续的高精度解算。4.1 数据重采样降低数据量现代GNSS接收机的采样率可以很高比如1Hz、5Hz甚至10Hz。这意味着每秒钟就有很多组观测值。对于某些不需要高时间分辨率的研究比如静态大地测量或长时间的气候学研究如此高的采样率会导致数据文件非常庞大处理起来耗时耗力。这时我们就可以使用TEQC进行数据重采样降采样。假设我有一个1秒采样率的文件highrate.21o我想把它降到30秒采样率命令如下.\teqc -O.dec 30 highrate.21o lowrate.21o这里的-O.dec 30就是降采样参数表示每30秒保留一个历元第一秒的数据中间29秒的数据被丢弃。是命令行中的重定向符号意思是将TEQC处理后的结果输出保存到新的文件lowrate.21o中。非常重要的一点如果你的原始数据采样间隔小于或等于1秒在进行降采样时必须先用-O.int参数明确告诉TEQC原始的采样间隔是多少。例如原始数据是0.1秒10Hz采样的要降到60秒.\teqc -O.int 0.1 -O.dec 60 ultrahigh.21o ultrahigh_decim60.21o如果不指定-O.intTEQC可能会误判采样间隔导致降采样结果错误。4.2 简单的数据切割与拼接TEQC也可以用来按时间切割数据或者将多个同类型的观测文件首尾拼接起来。比如我有一整天的数据但只想要UTC时间06:00到18:00这12个小时的数据.\teqc -sta 我的站点 -day 2023 150 -dh 6 12 allday.21o daytime.21o这里的-day 2023 150表示年积日2023年第150天-dh 6 12表示从当天第6小时开始截取12小时的时长。拼接功能对于处理因接收机重启而产生的多个连续文件非常有用。假设我有site1.21o、site2.21o、site3.21o三个文件它们是同一站点连续观测的.\teqc site1.21o site2.21o site3.21o site_combined.21oTEQC会自动按照时间顺序将它们合并成一个文件。但务必注意这个功能对导航星历文件.21n要谨慎使用因为不同时间的星历文件可能包含重复的卫星星历直接拼接可能导致冲突。通常导航文件是直接取用时间跨度覆盖观测时段的那一个即可。5. 从接收机原始数据到RINEX格式转换实战虽然前面我们用gfzrnx处理了RINEX 3到RINEX 2的转换但TEQC最原始、最核心的功能其实是直接把各种接收机厂商的原始二进制数据转换成RINEX 2.11格式。这是数据预处理流水线的真正起点。5.1 基本转换命令当你从接收机比如Trimble、Leica、NovAtel等用配套软件下载到.dat、.t02、.m00等格式的文件后就可以使用TEQC进行转换。命令形式极其简单.\teqc 原始数据文件.dat 输出文件.21oTEQC会根据文件内容自动识别接收机类型和格式。同时它通常也会尝试在相同目录下寻找或生成对应的导航文件.21n。为了获得更完整的信息我习惯加上站点名、天线信息等元数据.\teqc -tr m -O.obs G L C -nav 输出文件.21n -sta SITE001 -ant AT504 -rcv TRIMBLE NETR9 原始文件.t02 SITE0010.21o我来解释一下这几个参数-tr m 指定接收机类型为“Trimble”具体字母代码需查TEQC手册。-O.obs G L C 指定输出文件中包含GPS(G)、GLONASS(L)、北斗(C)的观测值。-nav 输出文件.21n 指定输出的导航文件名。-sta SITE001 设置站点名这个会写入RINEX文件头。-ant AT504和-rcv TRIMBLE NETR9 设置天线和接收机型号这些信息对于高精度处理中的天线相位中心改正至关重要。5.2 转换中的常见问题与排查转换过程并非总是一帆风顺。如果你遇到TEQC报错或者转换出来的RINEX文件用其他软件打不开可以从以下几个方面排查文件格式识别错误TEQC可能无法自动识别你的特定接收机型号或固件版本输出的数据格式。这时需要显式指定格式例如对于某些天宝二进制格式可能需要尝试tfb或tr等选项。最好的方法是查阅TEQC官方文档中关于你所用接收机型号的说明。数据不完整或损坏如果原始数据在下载或传输过程中损坏TEQC会在转换过程中报错提示“unexpected end of file”之类的信息。解决方法是重新从接收机下载原始数据。时间系统问题有些接收机数据使用本地时间而RINEX标准要求是GPS时间或UTC。如果时间转换出错会导致所有观测值的时间标签错误。在命令中可以使用-week和-toh等参数手动指定GPS周和秒但这需要你知道准确的时间信息。内存不足处理超长时间如数周连续观测的高频数据时TEQC可能会因为内存不足而崩溃。这时可以先用前面提到的数据切割功能将大文件分成几个小段分别处理。我个人的经验是对于一款新型号的接收机第一次转换时最好先用一小段数据做测试并仔细检查生成的RINEX文件头信息是否正确再用QC功能快速浏览一下数据质量确认无误后再进行批量处理。磨刀不误砍柴工这个步骤能避免后续大量返工。6. 进阶技巧与脚本化批量处理当你熟练掌握了单个文件的操作后肯定会面临批量处理几十甚至上百个数据文件的任务。手动一条条输入命令是不可行的这时就需要借助脚本的力量。6.1 编写批处理脚本Windows Batch示例在Windows下我们可以写一个简单的.bat批处理文件。假设我有一个文件夹里面全是.t02格式的Trimble原始数据文件需要全部转换成RINEX并做质量检核。echo off set TEQC_PATHD:\GNSS_Tools\teqc.exe set GFZRNX_PATHD:\GNSS_Tools\gfzrnx.exe for %%f in (*.t02) do ( echo Processing %%~nf... %TEQC_PATH% -tr m -sta %%~nf -ant AT504 -rcv TRIMBLE NETR9 %%f %%~nf.21o if exist %%~nf.21o ( echo QC for %%~nf... %TEQC_PATH% qc %%~nf.21o ) ) echo All done! pause这个脚本会遍历当前目录下所有.t02文件对每个文件进行转换以文件名作为站点名并在转换成功后自动运行质量检核。你可以根据实际情况修改TEQC_PATH、天线型号等参数。6.2 在Linux/macOS下使用Shell脚本在Linux或macOS环境下使用Bash shell脚本会更加强大和灵活。下面是一个功能更全面的示例它包含了版本降级如果检测到RINEX 3、格式转换和质量检核并将所有QC报告汇总到一个HTML文件中需要借助其他工具如qc2html。#!/bin/bash # 定义工具路径 TEQC/usr/local/bin/teqc GFZRNX/usr/local/bin/gfzrnx for raw_file in ./raw_data/*.[Tt]02; do base_name$(basename $raw_file .t02) echo 正在处理: $base_name # 步骤1转换原始数据为RINEX 2.11 $TEQC -tr m -sta ${base_name:0:4} -ant AT504GG -rcv TRIMBLE NETR9 $raw_file ./rinex/$base_name.21o 2 ./log/convert.log # 步骤2检查并转换可能的RINEX 3导航文件 nav_file./nav/brdc${base_name: -3}0.21n # 假设导航文件命名规则 if [[ -f $nav_file ]]; then $GFZRNX -finp $nav_file -fout ./rinex/$base_name.21n --version_out 2 fi # 步骤3执行质量检核 if [[ -f ./rinex/$base_name.21o ]]; then $TEQC qc -R -E -J ./rinex/$base_name.21o ./qc/$base_name.qc.log 21 echo 已完成QC: $base_name fi done # 步骤4可选汇总QC结果 cat ./qc/*.qc.log | grep -A5 SUMMARY ./qc/all_stations_summary.txt这个脚本结构清晰加入了错误日志记录并且更健壮。你可以把它保存为process_gnss.sh然后在终端里用bash process_gnss.sh命令运行。通过脚本化你可以把整个预处理流程固定下来确保每次处理的方法和步骤都是一致的极大提高了重复工作的效率和可靠性。走到这一步你已经从一个需要手动敲命令的新手进化到可以搭建自动化预处理流水线的“玩家”了。TEQC工具虽然古老但其设计理念和功能在今天依然极具价值。它强迫你去理解数据处理的每一个步骤而不是当一个“黑箱”的点击工。这种理解对于你后续深入学习GNSS数据解算原理、诊断复杂问题有着莫大的好处。我至今还记得第一次看到自己处理的站点QC图发现下午信噪比周期性下降是因为旁边一棵树的影子时那种把抽象数据和现实世界对应起来的兴奋感。希望TEQC也能成为你探索GNSS世界的一把得力钥匙。