GPS高程拟合方法.

GPS高程拟合方法

3.1等值线图示法

等值线图示法是最直接的求算高程异常的方法。这种方法的核心思想就是内插的思想，绘制高程异常的等值线图，然后采用内插法来确定未知点的高程异常值。具体操作十分的简单，在测区内制定分布均匀的GPS点，用水准测量的方法来测定这些点的水准高，根据公式ζ=H-Hr求出这些点的高程异常，选择适当的比例尺按照已知点的平面坐标展会在图纸内，对已知点标注出高程异常值，再确定等高距，绘制出高程异常值的等值线图。之后就可以内插出待测点的高程异常值，进而求出待测点的正常高。这种方法只适用地形相对平坦的地方，在此种测区内采用这种方法拟合的高程精度可达到厘米级。测区的地形相对复杂内插出的高程异常值就不准确，而且这种内插法的精度往往取决于两个方面，分别是测区内GPS点的分布密度和已知点大地高的精确度。首先GPS点的分布比较密集，那么内插精度就相对较高，如果比较稀疏这时候就要借助于此测区的重力测量资料，提高内插精度。且还要注意GPS点间高程异常的非线性变化。另外就是水准点的精度，联测时尽量选取高精度的正常高，尽可能使得出的高程异常值准确，进而才能内插出待测点高精度的高程异常值。这种方法虽然简单易操作，但是有其弱点，就是精度不高，只有当对拟合精度要求不高的时候才使用此种方法(注：等值线法不需构造数学模型)。

3.2狭长带状区域线性拟合

解析内插法作为拟合高程最常用的方法，主要思想是把似大地水准面用数学曲面近似拟合，建立所在测区内最为接近似大地水准面的数学模型，以此来计算测区内任意点的高程异常值，从而计算出正常高。这种方法计算出的高程异常值的精度是由所采用的数学模型和似大地水准面的拟合程度所决定的。

解析内插法在选择数学模型时，首先要考虑的就是GPS点的分布情况。GPS 点的分布情况可分为带状分布和面状分布。若GPS点是呈线状布设，而且是以

沿线似大地水准面为一条连续且光滑的曲线，这时就可以采用相对于狭长带状区域的解析内插法来内插出待定点的高程异常值，从而求出待定点的正常高。这种线状分布的内插原理是：测区内已知水准点，用GPS测出其GPS高程，计算出已知水准点高程异常值，根据已知点的平面坐标和计算得出的高程异常值，构造出一个插值函数，这个函数是用来拟合GPS分布线上的似大地水准面的。用这

个函数内插出位置点的高程异常值。下面是两种用来拟合线状分布的GPS高程

的内插法。

3.2.1 多项式曲线拟合法

多项式曲线拟合是线状分布拟合的主要方法。多项式拟合顾名思义其插值函数是一个m次的代数多项式，若高程控制点的高程异常为ζ，坐标为x l(或y i或d i 或拟合坐标或x i-x0或y i-y0)的函数关系为下式：

ζl=a0+a1x l1+a2x l2+a3x l3+……a m x l m3-1

各高程控制点的已知高程异常与其拟合值之差为下式所示：

r i=ζl(x)-ζi(i=0，1，2…n)3-2

上式我们称之为离差。(3-1)中x l是拟合点到参考点(x0，y0)的直线距离，x0，y0为设定的常数值。在一般情况下都认为，x0，y0就是测区内已知点坐标的均值。

多项式曲线拟合使用起来非常方便，但是它有自身的局限性，即是使用这种方法的时候，所测路线不能太长，要限制控制点到测点的距离不能太远，通常把距离控制在300米以内。这个要求是因为使用多项式曲线方法拟合似大地水准面，如果它拟合的范围太大，点位的高程异常变化就越复杂，削高补低的方法不能满足我们所要求的精度。随着多项式阶数的增大，也会使拟合出的曲线振荡的更厉害，从而造成拟合的误差增大。这些造成了多项式曲线拟合的缺陷，但是在路线较短的情况下，这种方法有足够的精度来拟合GPS点的正常高程。

在式(3-1)中用m次多项式拟合似大地水准面，这个m的值如何取定，一般

情况下如果测区不是很长，地形相对平坦，那么我们通常取m取为3。也就是说多项式为三次多项式。若测区比较长或者是测区地形比较复杂就要依情况而定，增加多项式的次数，提高拟合精度。依上述分析m的取值主要和测区长度以及测区的复杂程度有关。

3.2.2三次样条曲线拟合法

三次样条曲线拟合法针对测线长，已知点多的测区GPS高程拟合问题。由上述可以知道，当测线比较长已知点较多的时候，就需要构造高次的拟合多项式，当m值比较高的情况下，会出现不稳定的现象，对求解高程异常值会有比较大的影响，并且最小二乘法在求多项式系数中也会增大削高补低的误差，因此为了避免测线长、已知点多这种情况下所出现的问题，通常采用分段拟合的方法，采用三次样条函数拟合数学模型。这种方法很好的解决了因测线长而引起的问题。

三次样条曲线的实质就是一个拼接而成的连续函数，在把测线分为多段的情况下，每段设为三次多项式函数，然后将这些多项式函数组成三次样条函数。为了计算准确，应用中要求这种构造的曲线不仅在连接点处函数要连续。而且还要求这个函数的一级导数还有二级导数全部要是连续的，才能保证在分段之后构造的三次样条函数后期运算中能够计算出准确的高程异常值。

设过n个已知点，ζi和x i(或y i或拟合坐标)在区间[x i,x i+1](i=1,2,…,n-1)上有三次样条函数关系:

ζ(x)=ζ(x i)+(x-x i)ζ(x i,x i+1)+(x-x i)(x− x i+1)ζ(x,x i,x i+1)3-3

式中，x为待定点坐标,ζ(x i,x i+1)为一阶差商,ζ(x i,x i+1)=(ζi+1−ζi)/(x i+1-x i)；ζ(x，x i,x i+1)为二阶差商,ζ(x,x i,x i+1)=1

[ζ,,(x i)+ζ,,(x)+ζ,,(x i+1),而

6

ζ,,(x)(i=1,2,…,n-1)，满足系数矩阵为对称三角阵的线性方程组

(x i-x i+1)ζ,,(x i−1)+2(x i+1,x i−1)ζ,,(x i)+(x i+1,x i)ζ,,(x i+1)

=6[ζ,,(x,x i+1)-ζ(x,x i)]3-4

ζ(x0)=ζ,,(x n)=0