曲面拟合法求GPS网正常高的几点认识.docVIP

曲面拟合法求GPS网正常高的几点认识 吉渊明　赵水泉(浙江省测绘产品质检站　310012) 　　【摘　要】　本文分析了二阶多项式函数模型的几何特征。通过曲面拟合测区的似大地水准面时对高程拟合起算点的可靠性分析、整体拟合和分区拟合的选择、GPS高程拟合精度的分析等方法提出了自己的见解，并通过实例，介绍了上述方法的具体应用。　　一、前　言　　GPS在控制测量中的应用越来越受到测量工作者的青睐，特别是平面控制测量，利用GPS可以很灵活地布网，测量速度快效益好，精度高且均匀，减少了手工操作和内业计算工作。然而GPS正常高的确定受高程异常的影响，用曲面拟合时需要多次试算，选择不同的起算点，不同的分区均对GPS高程拟合精度有较大影响。笔者通过对曲面拟合模型(二阶多项式函数)的分析，提出利用高程异常图选择整体拟合和分区拟合，对高程拟合起算点的可靠性分析和GPS高程拟合精度的分析等方法，有助于提高GPS高程拟合精度。通过实例，介绍了上述方法的具体应用。　　二、大地高与正常高的关系　　国家高程系统采用的是正常高或正高系统，我国采用正常高系统。正常高(或正高)是由地球重力场定义的高程参考面，即似大地水准面(或大地水准面)起算的。　　GPS测定的是大地高差，获得的是大地高而不是正常高，这是用GPS测定高程的主要问题。　　经过GPS网的三维平差，可获得各点的大地高hi，如果在其中一些点上同时进行水准测量可得出该点的正常高Hi，则该点高程异常可按下式求出： ζi=hi-Hi　　　　　　　　　　　　　　　　　　　　　(1) 　　若由GPS和水准测量共测点所得的ζi为已知，内插其他GPS测点高程异常，从而获得各点的正常高，这是GPS测定高程控制的基本方法。　　三、用二阶多项式函数拟合高程异常(曲面拟合)　　曲面拟合法通常选取二阶多项式函数，同济大学研制的TGPPS软件也采用如下模型： ζi＝a0+a1xi+a2yi+a3x2i+a4y2i+a5xiyi+εi　　　　　　　　　　　(2) 式中，ζi表示高程异常，a0,…，a5为拟合系数，xi、yi为平面坐标，εi为残差。　　通过6个以上已知高程异常值的GPS测点，用最小二乘法(［εi2］=min)求拟合系数，然后利用式(2)(不考虑残差)求其他GPS点高程异常，最后获得正常高。　　下面来分析一下模型(2)的几何特征。　　对式(2)求一阶偏导： 　　　　　　　　　　　　　　　(3) 　　求二阶偏导： 　　　　　　　　　　　　　　　　　(4) 　　一旦式(2)确定，a0，…，a5均为一常数，模型(2)与某一切平面存在一个切点，说明式(2)几何表示为一“抛物单曲面”，仅有一个凹面或凸面。了解二阶多项式函数模型的几何特征非常重要，如果高程异常图像鞍部、水波浪形等存在多个凹凸面时，用一个二阶多项式函数就无法以数学形式表示，采用TGPPS软件计算时结果就出现较大残差；选择不同起算点，无论怎样试算也无法得以改善。　　四、对高程拟合起算点的可靠性分析　　似大地水准面的形状可以通过高程异常图来反映。从式(1)可知高程异常值与正常高和大地高有关，等值线的绘制与分布点的高程异常值有关。在高程异常图的绘制过程中，如果有某点的高程异常值与等值线趋势明显地不协调，我们应对该点的高程(正常高和大地高)进行分析和检测。对大地高的数据分析可以打开野外数据采集时仪器高的数据文件核查仪器高的输入与手簿记录是否一致，是否二次设站；对正常高的数据分析采用常规方法；必要时实地检测大地高或正常高。　　笔者最近对某市690 km2 GPS控制网的检验中通过高程异常图的绘制(图1)发现443号点的高程异常值与等值线趋势存在20 cm的偏离。从二次设站观测方法可以判定该点大地高数据可靠；正常高的分析采用实地检测方法，用三等水准附合路线检测该点高程，也得出原高程正确的结论。 图 1 　　从高程异常图可以看出443号点周围高程异常变化较大，为了反映高程异常细部，有必要在该点周围增加水准联测点密度，有利于提高等值线绘制精度和高程拟合精度。 　　五、整体拟合和分区拟合的选择　　二阶多项式函数模型的几何特征为一“抛物单曲面”，然而一个测区的似大地水准面并不会如此规整，为了使二阶多项式函数与似大地水准面拟合得更好，有必要为二阶多项式函数定义一个区域，即选择整体拟合和分区拟合。　　从二阶多项式函数模型的几何特征我们知道如果似大地水准面的形状仅存在一个凹面或凸面，选择整体拟合能够较好地获得表示似大地水准面形状的二阶多项式函数。在使用TGPPS软件时由于受内存的限制，三维平差一次最多平差70点左右，对于较大范围控制的测区可以分区三维平差，获得的大地高采用同一基准，然后利用高程计算中数据编辑功能把各分区的高程数据文件组成一个整体拟合的