广西AB级网及似大地水准面精度检测分析
廖超明1,2 罗满建2
(1.武汉大学GPS工程技术研究中心 湖北武汉 430079;
2.广西壮族自治区测绘局 广西南宁 530023)
【摘 要】 本文作者针对广西A、B级GPS控制网和广西2′30″×2′30″格网高分辨率似大地水准面成果外符合精度检测方法进行了深入研究,推导了外符合精度检测公式,通过高精度数据处理与统计分析,求得被检成果外符合精度,对被检项目所达到的精度等级进行了客观、准确的评价。
【关键词】 GPS 精密数据处理 外符合精度 检测
0 引言
为客观、准确地评价广西A、B级GPS控制网和广西2′30″×2′30″格网高分辨率似大地水准面成果所达到的精度,本文作者对被检测成果外符合精度检测方法进行了研究,推导出外符合精度检测公式,并通过GPS精密数据处理和统计分析,得到了被检项目本次检测的外符合精度。
1 外业数据采集
为获得高质量的GPS检测数据,外业统一采用双频大地型GPS接收机进行数据采集。广西A、B级GPS控制网外业检测数据采集参照《全球定位系统(GPS)测量规范》(GB/T 18314-2001)相应等级技术参数执行,详见表1。
表1 GPS数据采集技术参数表 单位:小时
检测区域 |
检测时段长(h) |
时段数 |
连接方式 |
地理位置 |
检测区域 |
检测时段长(h) |
时段数 |
连接方式 |
地理位置 |
平果B级子网 |
12 |
2 |
- |
桂西 |
林宾D级网 |
4 |
2 |
边连接 |
桂中 |
北海A级子网 |
24 |
3 |
- |
桂南 |
梧州D级网 |
8 |
1 |
网连接 |
桂东 |
河池D级网 |
4 |
2 |
边连接 |
桂西北 |
钦州D级网 |
4 |
2 |
边连接 |
桂南 |
广西似大地水准面外符合精度检测所选的四个D级GPS控制网分布于广西的桂西北、桂中、桂东和桂南,基本能够代表了广西的山地、丘陵、平地三种地形地貌。因此,这四个区域的外符合精度检测数据能够客观、真实地反映广西高精度、高分辨率似大地水准面成果的整体精度。
2 GPS精密数据处理技术研究
2.1 测量基准
高程基准:1985国家高程基准。
地心基准:ITRF97,参考历元2000.0。
2.2 GPS精密数据处理方案
为保证检测成果的客观性和可靠性,作者查阅了有关GPS数据处理资料,研究制定了如下的GPS精密数据处理方案。
(1)精密数据处理软件
GPS基线处理采用美国麻省理工学院和Scripps研究所共同研制的GAMIT(10.21V)软件,该软件是国际公认最好的GPS基线数据处理软件之一。
(2)精密星历
卫星轨道精度是影响GPS基线解算精度的重要因素之一。因此,GPS基线数据处理时采用IGS精密星历(ITRF框架)SP3格式。
(3)坐标框架与历元
GPS相对定位基准是由卫星星历和基准站坐标共同给出的。因此,GPS控制网基线数据处理所采用的坐标框架与历元为观测期间IGS 事后精密星历所对应的框架和瞬时历元;GPS网平差再统一到ITRF97、参考历元2000.0上。
(4)起算坐标
为确保基线解算的精度和稳定度,获得与广西A、B级GPS网一致的高精度ITRF坐标,在GPS基线处理中,直接引入联测的一个A、B级网点ITRF坐标加以控制。
(5)基线数据预处理
GPS基线数据预处理主要包括了四步:生成RINEX数据格式、编制接收机天线高度表、编制已知坐标和概略坐标表、IGS精密星历及其它信息文件的下载和编辑。
(6)基线解算主要模型和参数
GPS基线数据处理主要考虑以下五个模型改正:卫星钟差模型改正、接收机钟差模型改正、电离层折射模型改正、对流层折射模型改正、卫星和接收机天线相位中心模型改正。
3 外符合精度检测方法
3.1 广西A、B级GPS控制网
广西A、B级GPS控制网外符合精度检测目的是要求定广西A、B级GPS控制网三维坐标分量测量中误差。由于检测子网所求解的三维坐标分量不是真值,不能直接由公式(1)、(2)、(3)求A、B级GPS控制网的三维坐标分量中误差。以下作者根据误差传播定律详细推导了公式。
按照GPS精密数据处理方法分别对平果B级子网、北海A级子网的检测数据进行基线处理和网平差,得到检测点的空间直角坐标。将检测点的检测成果(Xj,Yj,Zj)和同名点的广西A、B级GPS控制网成果(X,Y,Z)进行三维坐标分量求差:
△X=Xj-X,
△Y=Yj-Y, ……(1)
△Z=Zj-Z,
通过公式(2): ……(2)
可计算每个检测点Pi的点位不符值。
通过公式(3): ……(3)
求得点位不符值中误差(由于检测样本个数较少,用n代替n-1)。
由于广西A、B级GPS控制网外业观测技术参数与两个检测子网外业观测技术参数相同,统一采用GPS精密数据处理方法进行基线处理和网平差,可认为是两次同精度(等权)的独立观测值。为简化计算,设检测点三维坐标分量(Xj,Yj,Zj)和同名点(X,Y,Z)的中误差分别为 和 ,且mj=m。对式(1)、(2)进行全微分线性化:
…… (4)
根据误差传播定律得:
…… (5)
则 …… (6)
其中:mp为检测点点位不符值中误差。
m为广西A、B级GPS控制网三维坐标分量测量中误差。
3.2 广西似大地水准面
由于广西似大地水准面外符合精度检测的点位同时具有精密一、二等水准测量成果hs,可认为hs是检测点正常高的真值,即hs的测量中误差为mhs=0。通过精密数据处理,求得检测点大地经、纬度和大地高(B、L、H),然后通过广西似大地水准面计算程序求解检测点正常高hg。由下式(7)计算检测点hg与hs的正常高之差:
v = hg - hs …… (7)
根据误差传播定律得:

即:mv = mhg
由公式(8)检测点正常高之差的中误差来评定广西似大地水准面成果检测精度:
…… (8)
4 外符合精度检测结果
4.1 广西A、B级GPS控制网三维坐标精度检测结果
广西A级网共12个点,相邻点平均边长达到了217公里。本次外符合精度检测了4个点,抽检率达到了33.3%。以其中118号A级点为固定点,得到广西A级网检测点点位不符值中误差: =±0.029m;三维坐标分量中误差: =±0.021m。检测结果详见表2。
表2 广西A、B级GPS控制网外符合精度检核结果 单位:m
序号 |
点号 |
分量 |
检核网平差成果 |
广西A、B级GPS控制网成果 |
检核成果-AB级成果 |
点位不符值v |
备注 |
1 |
49 |
X |
-.232 |
-.232 |
0.000 |
0 |
固定点 |
Y |
.227 |
.227 |
0.000 |
Z |
.460 |
.460 |
0.000 |
2 |
55 |
X |
-.160 |
-.179 |
0.018 |
0.033 |
B级 |
Y |
.189 |
.217 |
-0.028 |
Z |
.876 |
.878 |
-0.002 |
3 |
63 |
X |
-.383 |
-.379 |
-0.004 |
0.018 |
B级 |
Y |
.241 |
.257 |
-0.016 |
Z |
.832 |
.840 |
-0.008 |
4 |
74 |
X |
-.329 |
-.325 |
-0.004 |
0.032 |
B级 |
Y |
.643 |
.615 |
0.028 |
Z |
.074 |
.059 |
0.015 |
5 |
76 |
X |
-.465 |
-.468 |
0.003 |
0.014 |
B级 |
Y |
.429 |
.419 |
0.010 |
Z |
.103 |
.094 |
0.009 |
6 |
98 |
X |
-.897 |
-.909 |
0.012 |
0.042 |
B级 |
Y |
.494 |
.526 |
-0.033 |
Z |
.049 |
.073 |
-0.024 |
7 |
101 |
X |
-.124 |
-.135 |
0.011 |
0.069 |
B级 |
Y |
.515 |
.573 |
-0.058 |
Z |
.450 |
.486 |
-0.036 |
8 |
106 |
X |
-.743 |
-.773 |
0.030 |
0.086 |
B级 |
Y |
.826 |
.900 |
-0.075 |
Z |
.630 |
.661 |
-0.030 |
9 |
110 |
X |
-.06 |
-.059 |
-0.001 |
0.019 |
A级 |
Y |
.023 |
.007 |
0.016 |
Z |
.239 |
.230 |
0.009 |
10 |
116 |
X |
-.130 |
-.135 |
0.005 |
0.030 |
A级 |
Y |
.469 |
.494 |
-0.026 |
Z |
.472 |
.488 |
-0.015 |
11 |
117 |
X |
-.503 |
-.483 |
-0.020 |
0.035 |
A级 |
Y |
.381 |
.353 |
0.028 |
Z |
.569 |
.561 |
0.008 |
12 |
118 |
X |
-.098 |
-.098 |
0.000 |
0 |
固定点 |
Y |
.504 |
.504 |
0.000 |
Z |
.570 |
.570 |
0.000 |
备注:三维坐标成果已省略了小数点前的整数位。
广西B级网共106个点,相邻点平均边长达到了65公里。本次外符合精度检测了8个点,抽检率达到了7.5%。以其中49号B级点为固定点,得到广西B级网检测点点位不符值中误差: =±0.049m;三维坐标分量中误差: =±0.035m。检测结果详见表2。
4.2 广西似大地水准面成果精度检测结果
广西似大地水准面外符合精度检测总共抽检了四个区域共65个C、D级控制点,1985国家高程基准下的一期成果外符合精度检测结果:mhg1=±0.034m;二期成果外符合精度检测结果:mhg2=±0.033m。1985国家高程基准下一期检测结果详见表3。
表3 广西似大地水准面(一期)外符合精度检测结果 单位:m
序号 |
点号 |
水准(Sz86) |
Gps86 |
v |
序号 |
点号 |
水准(Sz86) |
Gps86 |
v |
1 |
C164 |
.682 |
.735 |
0.053 |
2 |
C163 |
.123 |
.156 |
0.033 |
3 |
C160 |
.654 |
.699 |
0.045 |
4 |
D003 |
.174 |
.132 |
-0.042 |
5 |
D004 |
.263 |
.281 |
0.018 |
6 |
D005 |
.444 |
.452 |
0.008 |
7 |
D006 |
.386 |
.399 |
0.013 |
8 |
D008 |
.820 |
.844 |
0.024 |
9 |
D009 |
.463 |
.476 |
0.013 |
10 |
D010 |
.727 |
.756 |
0.029 |
11 |
D011 |
.729 |
.742 |
0.013 |
12 |
D012 |
.700 |
.715 |
0.015 |
13 |
D013 |
.557 |
.555 |
-0.002 |
14 |
D014 |
.837 |
.876 |
0.039 |
15 |
D015 |
.407 |
.442 |
0.035 |
16 |
D016 |
.105 |
.085 |
-0.02 |
17 |
D017 |
.994 |
.993 |
-0.001 |
18 |
D018 |
.710 |
.711 |
0.001 |
19 |
D019 |
.407 |
.383 |
-0.024 |
20 |
D020 |
.630 |
.611 |
-0.019 |
21 |
D004 |
.670 |
.675 |
0.005 |
22 |
D008 |
.034 |
.003 |
-0.031 |
23 |
D009 |
.570 |
.543 |
-0.027 |
24 |
D014 |
.832 |
.837 |
0.005 |
25 |
D015 |
.340 |
.357 |
0.017 |
26 |
D016 |
.804 |
.843 |
0.039 |
27 |
D017 |
.005 |
.948 |
-0.057 |
28 |
D018 |
.545 |
.516 |
-0.029 |
29 |
D025 |
.909 |
.907 |
-0.002 |
30 |
D027 |
.475 |
.486 |
0.011 |
31 |
D028 |
.504 |
.477 |
-0.027 |
32 |
D029 |
.834 |
.777 |
-0.057 |
33 |
D030 |
.505 |
.462 |
-0.043 |
34 |
D036 |
.440 |
.421 |
-0.019 |
35 |
D044 |
.141 |
.076 |
-0.065 |
36 |
D046 |
.820 |
.778 |
-0.042 |
37 |
D049 |
.830 |
.757 |
-0.073 |
38 |
D055 |
.242 |
.187 |
-0.055 |
39 |
D056 |
.618 |
.541 |
-0.077 |
40 |
D057 |
.688 |
.617 |
-0.071 |
41 |
C1033 |
.791 |
.718 |
-0.073 |
42 |
C1036 |
.371 |
.403 |
0.032 |
43 |
D001 |
.206 |
.227 |
0.021 |
44 |
D003 |
.374 |
.420 |
0.045 |
45 |
D004 |
.643 |
.650 |
0.003 |
46 |
D006 |
.088 |
.136 |
0.048 |
47 |
D007 |
.762 |
.750 |
-0.010 |
48 |
D009 |
.298 |
.280 |
-0.018 |
49 |
D010 |
.238 |
.243 |
0.005 |
50 |
D011 |
.160 |
.160 |
0.001 |
51 |
D012 |
.029 |
.022 |
-0.007 |
52 |
D013 |
.002 |
.012 |
0.010 |
53 |
D014 |
.235 |
.250 |
0.015 |
54 |
D015 |
.261 |
.240 |
-0.020 |
55 |
D019 |
.982 |
.020 |
0.039 |
56 |
D020 |
.473 |
.499 |
0.026 |
57 |
D021 |
.084 |
.105 |
0.021 |
58 |
D022 |
.002 |
.99 |
-0.014 |
59 |
D023 |
.710 |
.710 |
0.000 |
60 |
D024 |
.817 |
.842 |
0.025 |
61 |
D026 |
.509 |
.556 |
0.047 |
62 |
D027 |
.027 |
.068 |
0.041 |
63 |
D028 |
.581 |
.600 |
0.019 |
64 |
D029 |
.201 |
.218 |
0.017 |
65 |
D031 |
.810 |
.800 |
-0.014 |
|
|
|
|
|
备注:正常高值已省略了小数点前的整数位。
5 结论
由以上分析得到的广西似大地水准面外符合精度与内符合精度是一致的,检测结果能够客观真实地反映了广西似大地水准面所达到的精度指标。内符合精度与外符合精度比较详见表4。
表4 广西似大地水准面内外符合精度比对 单位:m
高程基准 |
MAX |
内符合精度MIN |
STD |
MAX |
外符合精度MIN |
STD |
1985国家高程基准一期 |
.088 |
-.084 |
.033 |
.053 |
-.077 |
±0.034 |
1985国家高程基准二期 |
.089 |
-.087 |
.034 |
.048 |
-.073 |
±0.033 |
广西
A、B级GPS控制网外符合精度检测公式推导过程中,作者假设三维坐标分量(X,Y,Z)测量中误差是相等的,通过用X、Z轴吸收了Y轴(高程方向)的部分误差,简化了公式推导。另外,由于A、B级网外符合精度检测抽样较小,其检测精度还不能完全代表了广西A、B级网的测量精度。
【参考文献】
[1] 刘基余,李征航等.全球定位系统原理及其应用[M].北京:测绘出版社,1995
[2] 廖超明.GPS接收机标准检定场建设及相关问题研究.武汉大学硕士论文,2004.4
[3] 廖超明,谭庆林等.GPS接收机检定原理及检定数据自动化处理思路探讨. 测绘通报(已收录)
[4] 广西现代空间定位基准的建立及似大地水准面的确定技术设计书.广西壮族自治区测绘局武汉大学,2003.12
[5] 广西现代空间定位基准的建立及似大地水准面的确定外业检测报告.广西壮族自治区测绘局武汉大学,2006.11
|