广西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/T18314-2001)相应等级技术参数执行,详见表1。表1GPS数据采集技术参数表单位:小时检测区域检测时段长(h)时段数连接方式地理位置检测区域检测时段长(h)时段数连接方式地理位置平果B级子网122-桂西林宾D级网42边连接桂中北海A级子网243-桂南梧州D级网81网连接桂东河池D级网42边连接桂西北钦州D级网42边连接桂南广西似大地水准面外符合精度检测所选的四个D级GPS控制网分布于广西的桂西北、桂中、桂东和桂南,基本能够代表了广西的山地、丘陵、平地三种地形地貌。因此,这四个区域的外符合精度检测数据能够客观、真实地反映广西高精度、高分辨率似大地水准面成果的整体精度。2GPS精密数据处理技术研究2.1测量基准高程基准:1985国家高程基准。地心基准:ITRF97,参考历元2000.0。2.2GPS精密数据处理方案为保证检测成果的客观性和可靠性,作者查阅了有关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备注149X-.232-.2320.0000固定点Y.227.2270.000Z.460.4600.000255X-.160-.1790.0180.033B级Y.189.217-0.028Z.876.878-0.002363X-.383-.379-0.0040.018B级Y.241.257-0.016Z.832.840-0.008474X-.329-.325-0.0040.032B级Y.643.6150.028Z.074.0590.015576X-.465-.4680.0030.014B级Y.429.4190.010Z.103.0940.009698X-.897-.9090.0120.042B级Y.494.526-0.033Z.049.073-0.0247101X-.124-.1350.0110.069B级Y.515.573-0.058Z.450.486-0.0368106X-.743-.7730.0300.086B级Y.826.900-0.075Z.630.661-0.0309110X-.06-.059-0.0010.019A级Y.023.0070.016Z.239.2300.00910116X-.130-.1350.0050.030A级Y.469.494-0.026Z.472.488-0.01511117X-.503-.483-0.0200.035A级Y.381.3530.028Z.569.5610.00812118X-.098-.0980.0000固定点Y.504.5040.000Z.570.5700.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)Gps86v序号点号水准(Sz86)Gps86v1C164.682.7350.0532C163.123.1560.0333C160.654.6990.0454D003.174.132-0.0425D004.263.2810.0186D005.444.4520.0087D006.386.3990.0138D008.820.8440.0249D009.463.4760.01310D010.727.7560.02911D011.729.7420.01312D012.700.7150.01513D013.557.555-0.00214D014.837.8760.03915D015.407.4420.03516D016.105.085-0.0217D017.994.993-0.00118D018.710.7110.00119D019.407.383-0.02420D020.630.611-0.01921D004.670.6750.00522D008.034.003-0.03123D009.570.543-0.02724D014.832.8370.00525D015.340.3570.01726D016.804.8430.03927D017.005.948-0.05728D018.545.516-0.02929D025.909.907-0.00230D027.475.4860.01131D028.504.477-0.02732D029.834.777-0.05733D030.505.462-0.04334D036.440.421-0.01935D044.141.076-0.06536D046.820.778-0.04237D049.830.757-0.07338D055.242.187-0.05539D056.618.541-0.07740D057.688.617-0.07141C1033.791.718-0.07342C1036.371.4030.03243D001.206.2270.02144D003.374.4200.04545D004.643.6500.00346D006.088.1360.04847D007.762.750-0.01048D009.298.280-0.01849D010.238.2430.00550D011.160.1600.00151D012.029.022-0.00752D013.002.0120.01053D014.235.2500.01554D015.261.240-0.02055D019.982.0200.03956D020.473.4990.02657D021.084.1050.02158D022.002.99-0.01459D023.710.7100.00060D024.817.8420.02561D026.509.5560.04762D027.027.0680.04163D028.581.6000.01964D029.201.2180.01765D031.810.800-0.014备注:正常高值已省略了小数点前的整数位。5结论由以上分析得到的广西似大地水准面外符合精度与内符合精度是一致的,检测结果能够客观真实地反映了广西似大地水准面所达到的精度指标。内符合精度与外符合精度比较详见表4。表4广西似大地水准面内外符合精度比对单位:m高程基准MAX内符合精度MINSTDMAX外符合精度MINSTD1985国家高程基准一期.088-.084.033.053-.077±0.0341985国家高程基准二期.089-.087.034.048-.073±0.033广西A、B级GPS控制网外符合精度检测公式推