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

GPS高程拟合法的比较分析(机械工业勘察设计研究院测量公司)摘要：工程中需要把GPS高程测量的大地高转换为正常高。通常的做法是采用拟合法建立研究区域的似大地水准面。本文介绍了两种不同的拟合方法：二次曲面拟合法、多面函数拟合法。并结合某区域一定数量已知GPS高程异常点来内插和外推研究区域内的任一点的高程异常。通过比较发现多面函数拟合法拟合的精度要比二次曲面拟合的精度高。关键词：高程转换；二次曲面拟合法；多面函数拟合法TheelevationofGPSfittingtothecomparisonandanalysis（Machineryindustrysurveyanddesigninstituteofmeasuringcompany）Abstract:GPSheightmeasurementoftheearthshouldbeconvertedtonormalhighinengineering.Itisusuallytoestablishthequasi-geoidoftheresearchareabythefittingmethod.Thisarticleintroducestwodifferentfittingmethods:quadraticsurfacefittingandmultiple-surfacefunctionfitting.CombinedwithacertainnumberofaregionknownGPSelevationanomalypointstotheinterpolationandextrapolationoftheheightanomalyatanypointwithinthestudyarea.Bycomparison,themultiple-surfacefunctionfittingtotheprecisionishigherthanthequadraticsurfacefitting.Keywords：Elevationconversion;Quadraticsurfacefitting;Multiple-surfacefunctionfitting1.引言传统的几何水准测量虽然精度高，但耗时长、耗费多、工作效率低。GPS由于自身测量精度高、速度快、工作效率高等优点被广泛应用于高程测量。GPS测量的高程坐标是在WGS-84坐标系下的大地高[1]，大地高是地面一点沿参考椭球面的法线到参考椭球面的距离，用符号H表示。实际应用中需要把GPS测得的大地高转换为正常高，正常高是地面点到通过该点的铅垂线与似大地水准面的交点的距离，用符号rH表示。似大地水准面到参考椭球面之间的距离称为高程异常，用符号表示。因此大地高与正常高之间的关系为：rHH（1）由于我国采用的高程系统是相对于似大地水准面的正常高，因此如何进行GPS高程转换成为当前研究的热点问题。拟合法是GPS高程转换中比较常用的方法，主要的拟合模型有二次曲面法、多面函数拟合法。2.二次曲面拟合法假设测量区域内任一点的坐标为(,)Pxy，其高程异常为(,)xy。则二次曲面的拟合法的数学模型可以表示为[2]：22012345(,)xyaaxayaxayaxy（2）上式中有六个未知的拟合系数012345,,,,,aaaaaa，若已知测区内高程异常点的个数大于6个，则可根据最小二乘原理minTVPV求出以上未知数[3]，由此测区内任一未知点(,)Qxy的高程异常带入（2）式便可求得。再根据（1）式便可求得任意点的正常高。3.多面函数拟合法美国Hardy教授于1977年提出了多面函数拟合的方法[4]，多面函数拟合法的理论基础是，任何一个圆滑的数学表面总可以用一系列有规则的数学表面的总和以任意精度逼近。根据这一思想，假设测量区域内任一点的坐标为(,)Pxy，其高程异常为(,)xy。则多面函数拟合法的数学模型可以表示为：1(,)(,,,)ujjjjxyaQxyxy（4）上式中，ja为待定系数，(,,,)jjQxyxy为核函数，(,)xy为高程异常未知的点的坐标，(,)jjxy为高程异常已知的点的坐标。核函数一般可取[5-8]：22(,,,)()()kjjjjQxyxyxxyy（5）其中称为光滑因子。通常k取值为1/2，-1/2，3/2。分别对应于以下三种类型：正双曲面型1/222(,,,)()()jjjjQxyxyxxyy（6）倒双曲面型1/222(,,,)()()jjjjQxyxyxxyy（7）三次曲面型3/222(,,,)()()jjjjQxyxyxxyy（8）式（4）可表示为：Qa（9）待求参数a可根据已测点的垂直速率值按最小二乘法计算，可得：1()TTaQQQ（10）于是未测点p的高程异常可求得：1()TTpppQaQQQQV（11）再根据（1）式便可求得任意点的正常高。4.实例分析研究区域内布设有16个GPS控制点组成的控制网，记为G1-G16。进行平面观测时采用GPSC级标准，高程采用二等水准观测，利用高精度平差软件解算后获得各点的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高程拟合的方法，并结合实例进行比较分析，二次曲面拟合与多面函数拟合都能满足实际工程中的应用，但就精度而言，多面函数拟合的精度要高。参考文献：[1]李征航，黄劲松.GPS测量与数据处理[M].武汉：武汉大学出版社，2005[2]雷伟伟，郑红晓.二次曲面拟合法在区域似大地水准面精化中的应用[J].测绘与空间地理信息，2008,31（6）：38-42[3]武汉大学测绘学院测量平差学科组.误差理论与测量平差基础[M].武汉大学出版社，2003[4]Hardy.R.L.Theapplicationofmulti-quadricequationsandpointmassanomalymodelstocrustalmovementstudies[J],NOAATechnicalReportNGS76,NGS11,1978[5]刘万林，王利，赵超英.GPS水准的有限元法与多面函数法的加权综合模型[J].地球科学与环境学报,2004,26(3)：49-51[6]江在森.中国西部大地形变监测与地震预测[M].北京：地震出版社,2001,17-20[7]黄立人，陶本藻，赵承坤.多面函数拟合在地壳垂直运动研究中的应用[J].测绘学报,1993,22（1）：25-31[8]张菊清，刘平芝.抗差趋势面与正交多面函数结合拟合DEM数据[J].测绘学报,2008,37（4）：527-530