TEQC功能简介

本章首先介绍了TEQC软件的基本功能，包括格式转换、数据编辑、质量检查、单点定位。接着列出了一些TEQC软件的常用命令，并作了详细的注释。

1. TEQC功能简介

TEQC（Translate/Edit/QualityCheck/Coordinate）是由UNAVCOFaclity 研制的为地学研究GPS监测站数据管理服务的公开免费软件，可用于检查双频GPS接收机的动态和静态数据质量。它利用伪距观测值和载波相位观测值的线形组合来进行GPS数据中的误差估计，在快速评定GPS数据质量方面有非常大的优势，一方面速度快，没有繁琐的操作步骤，只用几条简单的命令即可；另一方面能对GPS观测数据进行多角度全方位的质量分析，能分别从卫星高度角、方位角、多路径效应、电离层延迟误差、电离层延迟变化率、信噪比等方面在Qcview配合用图形的形式直观的反映GPS观测数据的质量。

TEQC主要有格式转换（Translate）、数据编辑（Edit）、质量检查（Quality Check）单点定位（Coordinate）四方面的功能。其中，格式转换可将许多不同厂家的GPS接收机观测文件转换成标准格式RINEX文件；编辑功能可用于对RINEX文件的字头块部分，也可以进行数据文件的任意切割与合并、卫星系统的选择和删减、接收机通道的选择和卫星高度角的设置等；质量检查可以反映出GPS数据的电离层延迟、多路径效应、周跳、信噪比等信息；单点定位则可以粗略计算出点位在WGS-84坐标系中的坐标和在大地坐标系中的坐标。它的四个模块是相互独立、互不影响的，既可以单独使用其中的一个模块，也可以组合使用。它不仅可以处理GPS系统，也可以处理GLONASS系统和GNSS系统。但它也有明显的不足之处：基于DOS界面，对它的操作是建立在命令的基础上，因此要想用好它必须熟悉许多繁杂的命令；能很好的反映观测数据的质量，但在对观测数据的修复改正功能还不是很完善。

2. TEQC命令介绍

下面简单介绍TEQC软件的操作格式。TEQC软件的命令格式非常有规律，其基本格式为：teqc {options} [File1 File2…] > File，其中teqc是可执行程序的名字（如teqc.exe），options控制参数，teqc软件包含了三百种左右的参数，可以控制完成各种功能，如格式转换、数据编辑、质量检查、点位坐标计算和帮助等。在参数前总会包含有“＋”或“－”符号，“＋”表示打开某一参数功能，也可以表示输出数据到屏幕（或文件）。“－”表示关闭某一参数功能或输入数据到文件。File1、Fiel2为待处理的数据文件。>File表示把处理结果保存到File中，若File所在目录已有名为File的文件，则覆盖之，否则自动创建名为File的文件。

2.1 RINEX格式转换（Translate）

不同的GPS接收机采用不同的数据格式，这给GPS数据处理软件造成了一定的困难，为了解决后处理软件数据的交换和输入格式问题，伯尔尼大学天文研究所提出了独立于接收机类型的数据交换格式RINEX，这一数据格式在科学研究领域得到了广泛的应用。大部分GPS接收机生产厂家都提供了相应的转换软件，但并非完全标准，在实际使用中，经常遇到格式不兼容的问题。TEQC软件的格式转换具有很好的通用性和较强的功能，适用于目前较常用的多数GPS接收机，而且还在不断的扩展。它主要是通过读入数据文件的开头部分来自动识别接收机的类型。

目前，TEQC可以转换的接收机数据文件大多是国外比较有名的GPS接收机，比如：美国的AOA、ASHTECH        、MOTOROLA、TRIMBLE、TECOM INDUSRTIES、ROGUE,瑞士的LEICA,加拿大的CMC,瑞士的WILD等，但是却不能识别中国的HD系列接收机和南方的北极星系列接收机。

以TRIMBLE 4700为例，介绍把TRIMBLE 4700接收机所接收到的数据文件转换成RINEX观测数据文件、导航数据文件和气象文件，其他类型的接收机转换与之类似。命令格式如下：（注意TEQC区分大小写，所以我们在写命令时一定要注意。）

         teqc -tr do -week 1323 +nav trimbel.05ntrimble.dat > trimble. 05o

式中, -tr 指明接收机的类型为trimble;do 指明输入文件为Dat 文件、输出为RINEX 观测数据文件(o文件) ; -week 1323 (可选) 指明观测日期对应的GPS 周数,或以年/ 月/ 日方式表示(即上述命令也可以表示成teqc  -tr do -week 2005／05／18 +nav trimbel.05ntrimbel.dat > trimbel. 05o); +nav 指明同时输出RINEX 导航数据文件;转换结果文件为观测数据文件trimbel.05o和导航文件trimbel.05n。

其中，上述命令的控制参数设置项根据GPS种类的不同可自行设置，常用的GPS类型对应的控制参数如下：

表3－1 常用GPS对应的设置项

控制参数	对应的GPS类型	控制参数	对应的GPS类型
-aoa	AOA(JPL)	-rock	Rockwell
-ash	Ashtech	-tr	Trimble
-cmc	Canadian Marconi	-ti	Tecom Indusrties
-lei	Leica	-rogue	Rogue
-mot	Motorola	-wild	Wild

2.2 数据编辑

TEQC具有强大的数据编辑功能。下面将把其主要的编辑功能命令逐一介绍：

2.2.1       头文件的编辑

RINEX的观测数据文件（*.??o）、导航数据文件(*.??n)和气象数据文件(*.??m)的文件头部分均可利用TEQC 软件进行设置和更改,还可添加新的注释行且原有的注释行保持不变。利用TEQC 软件对RINEX观测数据文件（*.??o）的字头块部分进行设置和更改的格式为:

teqc -O.* “ ” input file > output file

teqc -O.mo good trimble.05o > temptrimble.05o

式中,-O.mo good 指明将文件的“测站名”记录更改为“good”,temptrimble.05o 为更改后输出的RINEX 观测数据文件。在这里必须重新设置输出文件名,即不能与被处理文件trimble.05o 同名,否则输出文件为空文件。

除了可更改测站名外，还可以修改其他信息。其他常用的设置项如表3－2所示。

表3-2 o5文件常用设置项

指令	含义	指令	含义
-O.mn	设置测站点编号	-O.at	设置天线类型
-O.rn	设置接收机编号	-O.pe	设置天线高偏心改正
-O.rt	设置接收机类型	-O.int	设置采样间隔（s）
-O.an	设置天线编号	-O.c	追加注释

RINEX导航数据文件和气象数据文件的编辑工作与观测数据文件类似,只不过导航数据文件的设置项以-N. 开头,气象数据文件以-M. 开头,而观测文件是以-O. 开头。

2.2.2       RINEX文件的切割

在我们进行GPS数据观测时，一般情况下刚开始一段时间的观测数据精度较差，所以对观测数据文件中时间的选取就势在必行。在TEQC中，利用时间窗可对RINEX文件进行任意的切割，使得对RINEX文件的提取相当容易。时间窗的设置常采用一下几种格式：

1.[start] +d[ Y,M ,d ,h ,m ,s ]  指定以RINEX文件开始观测时间为上限的时间间隔+d[ Y,M ,d ,h ,m ,s ]。若提取的时间段以年（月/日/小时/分钟/秒）为单位则控制参数应为dY（dM/dd/dh/dm/ds),例如提取trimble中前10个小时的数据命令格式如下：

    teqc +dh 10 timble.05o>trimble10.05o

2.-d[ Y,M ,d ,h ,m ,s ] [ end ] 指定以RINEX 文件结束观测时间为下限的时间间隔-d[Y,M ,d ,h ,m ,s ]的数据。

3.-st [YYMMddhhmmss[ . sss ⋯] ] [end] 指定时间窗上限-st [ YYMMddhhmmss[ . sss ⋯] ] ,默认以RINEX文件结束观测时间为下限。

4.[start]-e[YYMMddhhmmss[. sss ⋯] ] 指定时间窗的下限- e[] ,默认以RINEX文件开始观测时间为上限。

5.[ YYMMddhhmmss[ . sss ⋯] ] - e [ YYMMddhhmmss[ . sss ⋯] ]，提取时间段st [ YYMMddhhmmss[ . sss ⋯] ]到e[ YYMMddhhmmss[ . sss ⋯] ]的观测数据。如：要想提取trimble中10~13点的数据，命令如下：

    teqc -st 20040420100000 -e 20040420120000trimble.05o>trimble1012.05o

6.-st[ YYMMddhhmmss[ . sss ⋯] ] + d[ Y,M ,d ,h ,m ,s ],提取从st[ YYMMddhhmmss[ . sss ⋯] ]开始顺延d[Y,M ,d ,h ,m ,s ]时间的数据信息。

2.2.3       RINEX文件的合并

TEQC既然可以对RINEX进行切割，就可以对其进行合并。但要注意的是要合并的文件除了要求是RINEX格式外，还必须实在时间上连续的文件。合并的命令格式如下：

teqc > file1 file2… > myfile

表示把file1，file2…合并为myfile并输出。

2.2.4       卫星系统的选择和特定卫星的禁用

对于GPS/GLONASS双星接收机,可以使用TEQC 进行卫星的选用,如去掉GLONASS 卫星数据的指令是: teqc -R 输入文件 > 输出文件 ；禁用prn # 的GPS 卫星的观测数据指令是: teqc -G# 输入文件> 输出文件。其中prn ＃为卫星的编号。

2.2.5       设置卫星高度角

电离层延迟、多路径效应、接收机噪声是影响GPS数据质量的主要因素，在进行对流层和电离层延迟分析时,需要考虑低高度角卫星,用如下指令可以设置卫星高度角限值。命令格式为：teqc -set_maskx 输入文件> 输出文件。

如要把trimble.05o中卫星的高度角限值设为15度，命令为：teqc –set_mask 15 trimble.05o>trimble15.05o。

2.2.6       设置观测值类型

GPS/ GLONASS 的观测值类型一般用L 1 ( L 1载波相位) 、L 2 ( L 2载波相位) 、C1 ( L 1的C/ A 码伪距) 、P1 ( L 1的P 码伪距) 、P2 ( L 2的P 码伪距) 、D1 ( L 1的多普勒观测值) 、D2 ( L 2的多普勒观测值) 表示, 选取并按指定顺序形成RINEX文件的指令是: teqc -O. obs“观测值类型”输入文件> 输出文件,应用格式如下：

 teqc -O.obs“L1L2P1P2”trimble.05o > trimble. new

2.2.7       卫星范围的设定

在用GPS接收机的采集数据的时候，用户可以自行设定接收卫星的最大编号。命令格式为：teqc –n_GPS ＃ file1>file2。若接收数据中有大于设定的卫星编号，则此组数据将被删除。＃的范围是0<#<256,默认的卫星编号为32。

teqc –n_GLONASS # file1>file2  设定期望的GLONASS卫星的最大编号。若接收数据中有大于设定的卫星编号，则此组数据将被删除。其中：＃的范围是0<#<256,默认的编号为24。

teqc –n_gssp # file1>file2  设定期望的gssp卫星的最大编号，若接收数据中有大于设定的卫星的编号，此组数据同样也会被删除。＃的范围是0<#<256,默认的期望卫星最大编号为51。

2.2.8       接收机通道的设定

命令解析：teqc +ch file1>file2  用GPS接收机的所有的通道（默认）

               teqc -ch# file1>file2  禁用gps接收机的＃频道

2.2.9       多路径误差的有关设置

在TEQC中，可以对多路径效应进行相应的设置。主要有：

1．多路径效应的开启和关闭      

开启多路径效应的命令格式为：teqc ＋ma file1>file2

关闭多路径效应的命令格式为： teqc -mafile1>file2

2．多路径误差的均值设定

L1的期望的多路径误差均值设定命令：teqc –mp1_rms# file1>file2

L2的期望的多路径误差均值设定命令：teqc –mp2_rms# file1>file2

其中，＃是以厘米为单位的，L1的默认值为50cm，L2的默认值为65cm

3．平均移动点数

命令：teqc –mp_win #file1>file2.其中0<#<65536,默认为50

2.2.10  信噪比的设定

载波L1上允许的最小信噪比设定： teqc –min_L1 # FILE1 >FILE2

载波L2上允许的最小信噪比设定： teqc –min_L2 # FILE1 >FILE2

其中：0<#<9,默认为0

2.2.11  电离层的有关设定

电离层设定有电离层延迟变化率和电离层延迟误差两个方面。命令分别如下：

1．电离层延迟。命令：teqc –ion_jump ＃ file1>file2 设置电离层延迟的最大值为＃cm，默认为1.00e＋035。

2．电离层变化率。命令：teqc –iod_jump# file1>file2 设置电离层延迟变化率为＃cm，默认为400cm/m。

2.2.12  GPS采样间隔的设定

在GPS测量中，一般观测间隔设定为5s或30s，设定格式如下：

teqc –O.int # file1>file2  其中＃>0

 

teqc –O.dec 30A.15O > B.15O  将采样间隔设置为30s，如果原始的采样间隔小于1s，比如0.5s，则需加上 –O.int 0.5，teqc –O.int 0.5 –O.dec30 A.15O > B.15O。如果文件头中已经写明internal，则不需加上 –O.int ** 选项

2.2.13  观测文件的信息查询

1.teqc +config file   通过此命令可以看到观测文件file 的各项设置。包括：GPS的频道设置，GPS、GLONASS、GLNASS卫星的编号设定，卫星周等信息。

2．teqc +meta file 通过此命令可显示观测文件的名称，大小，始末历元，采样率，测站名，站点编号，天线编号，天线类型，天线在大地坐标系中的坐标，天线高，接收机编号，接收机类型，坐标系等信息。

2.2.14  各种GPS接收机、天线类型的查询

我们可以通过命令 teqc +igs得到目前全世界范围内最流行的GPS接收机类型及天线类型。

2.2.15  TEQC 软件输出中的符号说明

在teqc程序中，我们常可以看到各种各样的符号。而通过命令teqc ++sym 可以查看各种符号代表的意义。

2.2.16  帮助信息

我们要想熟练使用任何软件，帮助信息必不可少，它可以帮助我们快速掌握我们所需的知识。在TEQC中弹出帮助信息的命令如下：

 teqc +help

通过此命令，我们不但可以看到teqc软件中的全部命令，而且他有英文注释，可以知道各具体命令的功能。

2.3 质量检核（Quality Check）

TEQC应用软件的数据质量检查功能可以处理静态或动态双频GPS和GLONASS导航定位系统的接收数据。只有单点数据并且包含广播星历信息才能进行数据质量检查，主要利用了伪距和载波相位观测值的线性组合方法。

根据是否利用导航文件信息，TEQC分为qc2lite和qc2full两种检核方式。

2.3.1       qc2lite 方式

如果输入文件只有RINEX 观测数据文件而没有导航数据文件,那么TEQC 将会在qc2lite 方式下运行。如运行teqc +qctrimble. 05o , TEQC 则对文件trimble. 05o 在qc2lite 方式下进行质量检核。通常在缺省状态下,质量检核的结果会生成报告文件trimble.05s 和数据文件trimble.ion(电离层延迟误差) 、trimble. iod(电离层延迟变化率) 、trimble.mp1 ( L1载波C/ A 码或P 码伪距的多路径影响) 、trimble.mp2 ( L2载波P 码伪距的多路径影响) 、trimble.sn1 ( L1载波的信噪比) 、trimble.sn2( L2载波的信噪比)。

2.3.2       qc2full 方式

如果输入文件为RINEX观测数据文件和导航数据文件,运行teqc +qc -nav trimble.05n  trimble.05o ;或者导航数据文件和观测数据文件在同一目录下,则TEQC 会自动搜索导航数据文件,而无需用-nav 指定,即运行teqc +qc trimble.05o,此时TEQC 则对文件trimble.05o 在qc2full 方式下进行质量检核。检核的结果除qc2lite方式下的报告文件和数据文件外,还增添了卫星和接收机天线的位置信息以及两个数据文件trimble. azi (方位角)和 trimble. ele(高度角)。

报告文件trimle.05s称之为质量汇总文件，包括有各个接受卫星的数据质量状态。其中的参数信息有：

一、信息总结。有各颗卫星的接收数据情况，观测时间，文件名，首历元、末历元，GPS天线在WGS-84中的坐标，在大地坐标系中的坐标，采样率，没捕捉到的卫星数目及编号，在RINEX中的期望观测值数，采集百分比，多路径平均误差，平均移动点数，始终漂移，探测到的周跳数，卫星两个高度角之间的观测值数目信息。

二、观测统计量。如平均多路径误差，总平均高度角，周跳数目等。

三、QC设置参数。如接收机的通道最大值，电离层延迟变化率的最大允许值，期望的平均多路径误差等。

数据文件又成为视图文件，在GPS数据质量检核中，除汇总文件外的所有文件，通称为视图文件。生成的视图文件如下表所示：

表3－3 GPS质量检查中生成的视图文件后缀及其含义

文件类型	含义	文件类型	含义
.mp1	L 1载波C/ A 码或P 码伪距的多路径影响	.iod	电离层延迟变化率
.mp2	L 2载波P 码伪距的多路径影响	.ion	电离层延迟误差
.sn1	L 1 载波的信噪比	.azi	卫星方位角
.sn2	L 2 载波的信噪比	.ele	卫星高度角

我们可以利用软件QCVIEW对数据质量检查结果进行分析，使得观测数据的质量在图像上得到直观的反映。

2.3.3       各历元伪距单点定位计算（coordinate）

虽然TEQC是一专业的GPS数据预处理软件，但是也可进行各历元伪距的单点定位。计算单点在空间坐标系中的坐标命令格式如下：

    teqc +qc +eepx trimble.05o>trimblexyz.xyz

计算单点在大地坐标系中的坐标命令格式如下：

teqc +qc +eepgtrimble.05o>trimblebl.xyz

注意：虽然在上述命令中没有用到导航文件（.N），但是在所进行操作的目录下必须有导航文件，否则会显示： The instruction at 0x0043fe01 referenced memory at 0x0000000c. The memory could not be read。错误信息，不能运行上述命令。

3. TEQC处理实例

在本例中，将对2004年4月14日使用TRIMBLE 4700类型的3台接收机观测1小时的数据cann105a.dat,irym105a.dat,melb105a.dat进行处理（在此我们仅对cann105a.dat的处理给出详细说明，其他两个观测文件的处理过程与此相同），过程如下：

3.1 格式转换

1.首先用命令：teqc –tr do –week1266 +nav cann105a.04ncann105a.dat >cann105a.04o 将cann105a.dat文件转换成标准文件RINEX的格式。通过本命令后，将会生成观测文件cann105a.04o和导航文件cann105a.04n。

3.2 质量检查

在转换生成RINEX文件的同时就可以编辑RINEX文件的数据，当然也可以在生成RINEX文件后在进行编辑。在此，由于TEQC的各个模块就独立的，我们就暂时跳过这一步，直接查看观测数据的质量，然后对数据进行适当的编辑。

使用命令teqc +qc cann105a.04o进行质量检核。运行此命令后，将会生成cann150a.mp1,cann150a.mp2,cann150a.sn1,cann150a.sn2,cann150a.ele,cann150a.azi,cann150a.iod,cann150a.ion,cann150a.04s等文件。

1．汇总文件cann150a.04s的有关信息。用word程序打开cann150a.04s，我们可以看到观测文件的名称为cann105a.04o，观测时间段为2004年4月14日00：00：00到00：59：55.001，接收机接收到的卫星编号及质量，天线在空间坐标系中的坐标x: -4346271.2673 y: 2595234.2754  z:3867368.9447 (m)以及在大地坐标系中的坐标-37.564855deg   149.157839 deg 158.5777 m，采样的时间间隔为5秒，接收到11颗卫星的信号，没有接收到的卫星编号分别为2   5  6   7   9 10  11  12 14  15  17 18  21  22 23  25  26 28  29  30  32 ，卫星高度角的限值为10o，观测值总数为7842个，其中高度角大于10o的有5338个，高度角大于10o但由于其他误差超限而被删除的观测值有194个，载波在L1上的多路径误差均值为0.071614 m，在L2上的均值为 0.113825 m，平均移动点数为50个，最大的电离层延迟率为400.00 cm/min，L1的期望多路径平均值为50.00cm，L2的期望多路径平均值为65.00cm，以5o为单位各个高度角范围内的观测值数目分别为：

85 - 90      0    0  0.000000 

 80 - 85     0    0   0.000000 

 75 - 80     0    0   0.000000 

 70 - 75   115    0   0.013890 

 65 - 70   396    0   0.025700 

 60 - 65   369    0   0.018057 

 55 - 60   336    0   0.030205 

 50 - 55   294    0   0.034133 

 45 - 50   219    0   0.046327 

 40 - 45   441    0   0.041733 

 35 - 40   482    0   0.044777 

 30 - 35   387    0  0.064248 

 25 - 30   289    0   0.071027 

 20 - 25   617    0   0.085342 

 15 - 20   643    0   0.124734 

 10 - 15   742    1   0.144828 

  5 - 10   281    0   0.114023 

  0 - 5      0    0  0.000000 

3.3 视图文件查看

1．查看卫星高度角  命令：qcview cann105a.ele。图4－1是第19颗卫星的高度角曲线图，从图中我们可以看出它的高度角变化范围为0～19o之间，图4－2是第13颗卫星的高度角曲线图，从图中得知它的变化范围为50～74o之间。

2．查看卫星方位角  命令：qcview cann105a.azi。 从图4－3中我们可以看出第13颗卫星的方位角变化范围为－140～170o之间，而从图4－4中看出第19颗卫星的方位角变化范围为40～60o之间。

3．查看卫星的多路径误差  命令：qcview cann105a.mp1（查看载波L1上的多路径误差） 命令：qcview cann105a.mp2（查看载波L2上的多路径误差）。从图4－5我们可以看出第19颗卫星在L1上的多路径误差MP1的变化范围为-30～50cm，从图4－7可以看出第19颗卫星在L2 上的多路径误差变化范围为-2.5～1.0m，说明此颗卫星多路径误差影响较大。从图4－6可以得知第13颗卫星在L1上的多路径误差范围为-10～10cm。,而从图4－8中可得知第13颗卫星在L2上的多路径误差变化范围为-0.1～0.1m，说明此颗卫星的多路径误差比较小。

4.查看电离层延迟误差 命令：qcview cann105a.ion。从图4－9可以看出第19颗卫星的电离层延迟误差，其值在-0.3～3m之间，并且在第560历元处探测到一个“I”周跳。从图4－10中可以看到第13颗卫星的电离层延迟误差为-00.1～0.3。从中我们不难看到卫星19的电离层延迟误差大于卫星13的电离层误差。

5．查看电离层延迟变化率 命令：qcview cann105a.iod。从图4－11我们很容易看出第19颗卫星的电离层延迟变化率为-0.9～1.0之间，在第560个历元处同样也有一个“I”标记。从图14可以看出第13颗卫星的电离层延迟变化率为-0.1～0.1之间，变化幅度相当小。

6．查看信噪比  命令：qcview cann105a.sn1(sn2)。第19颗卫星的信噪比在1-6之间，第13颗卫星的信噪比在8－9之间。

4. TEQC结果可视化

1. 首先下载QCVIEW32.EXE以及QCVIEW.EXE，下载网址http://www.unavco.org/software/data-processing/preprocessing/qcview/qcview.html

2. 如果你是在windows下，可能没法正常显示画出来的图形，原因是QCVIEW是在DOS环境下运行的，现在普遍使用的Windows 7/8不支持DOS下的图形显示，所以可以在网上下载一个工具DosBox，在windows下模拟DOS窗口，直接在官网上下载
3.  在dosbox环境下使用qcview
(1). 双击进入dosbox
(2). 敲入mount C E:\   （E为你的qcview.exe所在的盘）
(3). 敲入C:   
(4). 跟在DOS中一样：进入你的qcview所在的目录，然后运行命令，如 qcview32 test.ion -mono，将之前teqc生成的质量文件用图形的形式显示出来