teqc进行GNSS数据质量分析时生成的s文件介绍

本文介绍使用teqc进行GNSS数据质量分析时生成的s文件。来源于n年之前写的一篇文章，重新整理出来放到这里，有不正确的地方欢迎指正。

第一行中“|”为时段分隔符，例子中一个小时的数据被分成六个时段，每个时段 10 分钟，“-“ 也是一个时段符号。

开始字符为 5到 12的行为对应 PRN 号卫星的状态

1	~	A/S 关闭（antispoofing：反电子欺骗关闭），只有 L1，C/A，L2，P2 数据
2	c	钟跳
3	I	电离层跳变
4	1	只有 P1 的多路径影响跳变
5	2	只有 P2 的多路径影响跳变
6	-	卫星高于截止角但没有数据
7	L	失锁指示
8	*	A/S 关闭（antispoofing：反电子欺骗关闭），只有 L1，P1，L2，P2 数据
9	_	卫星高于地平线低于截止角，无可用观测
10	+	卫星低于截止角，有可用观测值

Obs 开始的行对应的是每一时段的收星数目

Pos 定位状态

1	^	位置异常-有较大变动
2	X	码定位失败
3	C	定位比较散乱
4	H	水平定位不准
5	V	高程定位不准
6	T	总定位不准
7	>	动态测量正常
8	o	静态测量正常
9	O	观测量不足
10	E	星历不足
11	S	卫星数目不足

Time of start of window : 2013 May 16  23:00:05.000
Time of  end  of window : 2013 May 17  00:00:00.000
Time line window length : 59.92 minute(s), ticked every 10.0 minute(s)

以上三行分别表示,

• 观测文件开始时间
• 观测文件结束时间
• 观测时间跨度 .

注意下边这7行只有指定星历文件时才会输出(teqc +qc -nav nav_filename obs_filename)，只用观测文件进行质量分析的话是不会输出的。这部分信息是teqc解算出的天线位置。值得注意的是teqc进行质量分析并不需要知道天线的高精度位置，teqc使用天线位置只是为了得到“星图”而已，因此一个粗略的坐标就足够了。
这里仅解释一下，|qc - header| position，这里标识的是teqc解算出的坐标和rinex文件头的坐标之间的距离，通常这个值是小于40米的。

 mean antenna; # of pos : 79
  antenna WGS 84 (xyz)  : -2188528.2434 4325271.3127 4132189.6898 (m)
  antenna WGS 84 (geo)  : N  40 deg 38' 9.49" E 116 deg 50' 19.64"
  antenna WGS 84 (geo)  :   40.635970 deg   116.838788 deg
          WGS 84 height : 553.0102 m
|qc - header| position  :  15.2284 m
qc position offsets     :  14.9719 m vertical   2.7831 m horizontal

Observation interval    : 5.0000 seconds

采样间隔，即观测文件中的两个历元之间的间隔

Total satellites w/ obs : 9
NAVSTAR GPS SVs w/o OBS :  1   2   3   4   6   7   8  10  11  13  14  16 
                          17  19  20  23  25  27  28  29  30  31  32 

• 第一行表示卫星总数
• 第二行表示观测信息不全的GPS卫星 PRN 序列

Rx tracking capability  : unknown
Poss. # of obs epochs : 720
Epochs w/ observations  :    720

这三行表示,

• 接收机最多可同时追踪的卫星数
• 理论上可观测到的历元数（对于采样率为 5 秒的 1 个小时的观测文件，理论上应观测到 720个观测值）
• 拥有“完整观测”的历元数

完整观测是teqc中很重要的一个概念，teqc质量分析的算法是基于完整观测的。所谓“完整观测”是指某观测包括 L1、L2 上的载波观测值，伪距观测值和载噪比。

Epochs repeated         :      0  (0.00%)
Possible obs >   0.0 deg:   5869
Possible obs >  10.0 deg:   4910
Complete obs >  10.0 deg:   4694
 Deleted obs >  10.0 deg:      1
  Masked obs <  10.0 deg:    215
Obs w/ SV duplication   :      0  (within non-repeated epochs)

这七行表示,

• 观测文件中重复历元出现的次数，0次是正常的
• 0高度角以上理论观测数
• 10度高度角以上理论观测数
• 10度高度角以上完整观测数
• 10度高度角以上剔除观测数
• 完整观测中截止角小于10度的数目
• 同一历元中同一卫星有重复观测（多于一个）的个数

Moving average MP12     : 0.755136 m
Moving average MP21     : 1.151419 m

MP1 和 MP2 是多路径的移动平均的 RMS 值，单位：米。是伪距和载波相位的线性组合。

Points in MP moving avg : 50

使用移动平均方法计算时的窗口宽度（点数）.

No. of Rx clock offsets : 0
Total Rx clock drift    :  0.000000 ms
Rate of Rx clock drift  :  0.000 ms/hr
Avg time between resets : Inf minute(s)

这四行表示：

• 重启次数
• 接收机钟漂
• 接收机平均钟漂的估计值
• 平均重启间隔

Freq no. and timecode   : 2 15008 003000
Report gap > than       : 10.00 minute(s)
       but < than       : 90.00 minute(s)

• 第一行的第一个数字表示频段数目，例子中表示双频输出，第一行的剩余内容提供给在线UNAVCO 接口
• 第二行表示数据分析最短数据长度

IOD signifying a slip   : >400.0 cm/minute
IOD slips               :     29
IOD or MP slips         :     31

第一行是一个门限值，表示电离层最大变率，大于这个门限值的认为是跳变.下边这两行teqc 如此计算，teqc 会计算电离层延迟的导数，与门限比较大于门限的视作跳变，然后统计跳变次数,

first epoch	last epoch	hrs	dt	#expt	#have	%	mp1	mp2	o/slps
13 5 16 23:00	13 5 17 00:00	1.000	5	4910	4694	96	0.29	0.24	4694

• hrs 采样时段长
• dt 采样间隔
• #expt 表示大于截止角的所有可能观测的数目
• #have 表示完整观测的实际数目
• % 表示#have/#expt
• o/slps 为完整观测数/跳变数，此值一般为几百，当然越大越好