GPS高程拟合方法的比较分析.docx

GPS高程拟合法的比较分析（机械工业勘察设计研究院测量公司） 摘要：工程中需要把 GPS 高程测量的大地高转换为正常高通常的做法是采用拟合法建立 研究区域的似大地水准面本文介绍了两种不同的拟合方法：二次曲面拟合法、多面函数拟 合法并结合某区域一定数量已知GPS高程异常点来内插和外推研究区域内的任一点的高 程异常通过比较发现多面函数拟合法拟合的精度要比二次曲面拟合的精度高关键词：高程转换；二次曲面拟合法；多面函数拟合法The elevation of GPS fitting to the comparison and analysis（Machinery industry survey and design institute of measuring company）Abstract: GPS height measurement of the earth should be converted to normal high in engineering. It is usually to establish the quasi-geoid of the research area by the fitting method. This article introduces two different fitting methods: quadratic surface fitting and multiple-surface function fitting. Combined with a certain number of a region known GPS elevation anomaly points to the interpolation and extrapolation of the height anomaly at any point within the study area. By comparison, the multiple-surface function fitting to the precision is higher than the quadratic surface fitting.Key words： Elevation conversion; Quadratic surface fitting; Multiple-surface function fitting1.引言传统的几何水准测量虽然精度高，但耗时长、耗费多、工作效率低。

GPS 由于自身测 量精度高、速度快、工作效率高等优点被广泛应用于高程测量 GPS 测量的高程坐标是在 WGS-84坐标系下的大地高［1］大地高是地面一点沿参考椭球面的法线到参考椭球面的距离, 用符号H表示实际应用中需要把GPS测得的大地高转换为正常高，正常高是地面点到通 过该点的铅垂线与似大地水准面的交点的距离，用符号H表示似大地水准面到参考椭球 r面之间的距离称为高程异常，用符号匚表示因此大地高与正常高之间的关系为：:二 H - H （1）r由于我国采用的高程系统是相对于似大地水准面的正常高，因此如何进行GPS高程转换成为当前研究的热点问题拟合法是GPS高程转换中比较常用的方法，主要的拟合模型 有二次曲面法、多面函数拟合法2. 二次曲面拟合法假设测量区域内任一点的坐标为P(x, y)，其高程异常为:(x, y)则二次曲面的拟合法的数学模型可以表示为［2］：:(x, y)二 a + ax + ay + a x 2 + ay 2 + a xy (2)0 1 2 3 4 5上式中有六个未知的拟合系数a ,a ,a ,a ,a ,a，若已知测区内高程异常点的个数大 012345于6个，则可根据最小二乘原理V7PV二min求出以上未知数［引，由此测区内任一未知点Q(x,y)的高程异常带入(2)式便可求得。

再根据(1)式便可求得任意点的正常高3. 多面函数拟合法美国 Hardy 教授于 1977 年提出了多面函数拟合的方法［4］，多面函数拟合法的理论基础 是，任何一个圆滑的数学表面总可以用一系列有规则的数学表面的总和以任意精度逼近根 据这一思想，假设测量区域内任一点的坐标为P(x,y)，其高程异常为:(x, y)则多面函 数拟合法的数学模型可以表示为：:(x, y )二》aQ (x, y , x , y ) (4)j=1 j j j上式中，a为待定系数，Q(x,y,x ,y )为核函数，(x,y)为高程异常未知的点的坐标, j j j(x , y )为高程异常已知的点的坐标核函数一般可取［5-8］：jjkQ ( x , y , x . , y . ) = ( x - x . )2 + ( y - y . )2 +§ (5)j j j j其中8称为光滑因子通常k取值为1/2, -1/2, 3/2分别对应于以下三种类型： 正双曲面型1/ 2Q ( x , y , x . , y . ) = ( x- x )2 + ( y - y . )2 +8 (6)倒双j j j j曲面型Q ( x , y , x , y ) =「( x - x.)2 + ( y - y .)2 + 8 1 -1/2 (7)j ］ L J J 」三次曲面型Q (x , y , x , y ) = 8 + L (x - x )2 + ( y - y )2」3/ 2 (8)J J J J式( 4)可表示为：9）待求参数a可根据已测点的垂直速率值按最小二乘法计算，可得：a 二（QtQ）-1 Qt： （10）于是未测点 p 的高程异常可求得：:=Q a = Q （QtQ）-1 QtV （11）p p p再根据（1）式便可求得任意点的正常高。

4. 实例分析研究区域内布设有16个GPS控制点组成的控制网，记为G1-G16进行平面观测时采 用GPS C级标准，高程采用二等水准观测，利用高精度平差软件解算后获得各点的GPS大 地高，通过水准观测求得各个GPS点的正常高，进而求得各点的高程异常值均匀选取控制网中10个点作为已知点，分别建立二次曲面拟合模型和多面函数拟合模 型，本文中多面函数模型选取正双曲面型通过对剩余的 6 个点进行高程拟合，分别比较拟 合所得的高程与已知的高程，来验证两种方法的拟合精度表 1 两种拟合方法精度对比点名已知高程二次曲面拟合高程残差多面函数拟合高程残差G219.76819.76120.015719.76120.0068G520.43520.5376-0.134620.5376-0.1026G818.49218.24580.358618.24580.2462G1121.98122.2979-0.428122.2979-0.3169G1417.86518.1623-0.381618.1623-0.2973G1617.53917.35560.256117.35560.1834从上表可以看出，两种拟合方法都能满足要求，多面函数拟合的精度总体而言要比二次曲面拟合的精度要高。

5.结论由于GPS测量自身的优点，GPS高程拟合在实际工程中得到了越来越广泛的应用，本文介绍了两种 GPS 高程拟合的方法，并结合实例进行比较分析，二次曲面拟合与多面函数拟合都能满足实际工程中的应用，但就精度而言，多面函数拟合的精度要高参考文献：⑴李征航，黄劲松.GPS测量与数据处理[M].武汉：武汉大学出版社，2005[2] 雷伟伟，郑红晓•二次曲面拟合法在区域似大地水准面精化中的应用[J].测绘与空间地理信 息， 2008,31(6)： 38-42[3] 武汉大学测绘学院测量平差学科组.误差理论与测量平差基础[M].武汉大学出版社，2003[4] Hardy.R.L. The application of multi-quadric equations and point mass anomaly models to crustal movement studies[J],NOAA Technical Report NGS 76,NGS 11,1978[5]刘万林，王利，赵超英.GPS水准的有限元法与多面函数法的加权综合模型[J].地球科学 与环境学报,2004,26 (3)： 49-51⑹江在森.中国西部大地形变监测与地震预测[M].北京：地震出版社,2001,17-20[7] 黄立人，陶本藻，赵承坤•多面函数拟合在地壳垂直运动研究中的应用[J].测绘学报,1993,22 (1)： 25-31[8] 张菊清，刘平芝•抗差趋势面与正交多面函数结合拟合DEM数据[J].测绘学报,2008,37(4)： 527-530。