GPS高程异常拟合方法

浅谈GPS高程异常拟合方法

摘要：在gps定位技术中,由于其测量定位技术的物理机制,其平面位置的精度可以达到较高水准,已被人们所认识和接受,而其高程精度较其平面精度约低2～5倍。尤其是在wgs-84坐标系向地方坐标系的转换过程中,由于wgs-84的大地高仅有几何意义而无物理内涵,而高程系统的正常高既有几何意义,又有地球内部质量密度分布不均匀这样一个物理现象。本文重点对gps高程水准拟合模型及其精度进行了分析探讨。

关键词：gps；高程异常；测量；定位技术

引言

gps定位技术因其优点突出,因而在测绘领域得到了广泛的应用。采用相对定位技术,通过gps网平差,可以得到高精度的平面坐标(或大地坐标)和大地高差；如果网中有1点或多点具有精确的wgs-84坐标系的大地高程,则可求得各gps点的大地高程。

gps 测量得到是wgs-84 椭球的大地高，而我国采用的是正常高系统，它是以似大地水准面作为参考面的，因此，精确计算gps 点的正常高，就必须作一些相应的转换。目前求定地面点的正常高的方法主要有gps 水准高程、gps 重力高程、gps三角高程、转换参数、整体平差和神经网络法等方法。

重力法是根据点位信息，可直接求得该点的高程异常值。在一定区域内，只要有足够数量的重力测量数据，就可以比较精确地求

定该区域的高程异常值。对于实施水准测量比较困难的丘陵和山区，利用重力测量方法是比较实用且可靠的方法。但此法的缺点是需要足够多且精度足够高的重力测量资料。从目前我国实际情况来看，gps 重力高程的精度低于gps 水准高程。

三角高程是在gps 点上加测各gps 点间的高度角(或天顶距)，利用gps 求出的边长，按三角高程测量公式计算gps 点间的高差，从而求出gps 点的正常高的一种方法。联合平差法是当测区内具有天文、大地、重力测量、水准测量及gps 测量等多种观测数据时，我们即可用整体平差模型将这些观测数据进行联合平差，最终可求得地面点的平面坐标及高程的最优无偏估值。

求转换参数法的原理是：当一测区内，有一定数量点平面坐标和高程己知，按坐标转换原理，求出参考椭球面与似大地水准面(或大地水准面)之间的平移和旋转参数，把这些参数加入gps 网的平差，在己知点高程约束下，通过平差，在求出各gps点平面坐标的同时，求出点的正常高。

人工神经网络是一门新兴交叉科学，它是生物神经系统的一种高度简化后的近似，是处理非线性映射问题的有效工具。基于神经网络来转换gps 高程是一种自适应的映射方法，没作假设，理论上比较合理，能避开未知因素的影响，减少人为构造的数学模型误差。

对于一般工程单位而言，无法获得必要的重力数据和天文数据，故一些常规方法和重力方法难于普及。三角高程加测天顶距，计算

边长，带来误差传播，并增加工程量，也是不可选的方法。gps 水准高程是综合利用gps 与水准测量成果求出似大地水准面，计算出未联测几何水准gps点的高程异常，从而求出这些gps 点的正常高的方法为gps 水准。

重力法是根据点位信息，可直接求得该点的高程异常值。在一定区域内，只要有足够数量的重力测量数据，就可以比较精确地求定该区域的高程异常值。此法的缺点是需要足够多且精度足够高的重力测量资料。

三角高程是在 gps 点上加测各gps 点间的高度角(或天顶距)，利用gps 求出的边长，按三角高程测量公式计算gps 点间的高差，从而求出gps点的正常高的一种方法。联合平差法是当测区内具有天文、大地、重力测量、水准测量及gps 测量等多种观测数据时，我们即可用整体平差模型将这些观测数据进行联合平差，最终可求得地面点的平面坐标及高程的最优无偏估值。

求转换参数法的原理是：当一测区内有一定数量点平面坐标和高程己知，按坐标转换原理，求出参考椭球面与似大地水准面(或大地水准面)之间的平移和旋转参数，把这些参数加入gps 网的平差，在己知点高程约束下，在求出各gps 点平面坐标的同时，求出点的正常高。

基于神经网络来转换 gps 高程是一种自适应的映射方法，没作假设，理论上比较合理，能避开未知因素的影响，减少人为构造的

数学模型误差。对于一般工程单位而言，无法获得必要的重力数据和天文数据，故一些常规方法和重力方法难于普及。三角高程加测天顶距，计算边长，带来误差传播，并增加工程量，也是不可选的方法。

gps 水准是从几何解析的角度出发，在gps网中联测一些水准点，再利用这些点上的正常高和大地高求出它们的高程异常值，根据这些点上的高程异常值与坐标的关系，拟合出测区的似大地水准面，内插入其他gps 点的高程异常，从而求出各个未知点的正常高。

对于线路工程，应用最广泛的是解析内插法，对于区域工程应用最广泛的是数学曲面拟合法。

1 gps 高程数学曲面拟合方法

当gps 点布设成一定区域面时，可以应用数学曲面拟合法进行拟合，求得待求点的正常高。

1.1 平面拟合法

在小区域且较为平坦的范围内，可以考虑用平面逼近局部似大地水准面。

设某公共点的高程异常ζ与该点的平面坐标有关系式

式中，a1,a2,a3 为模型待定参数。

1.2 二次曲面拟合法

似大地水准面的拟合也可采用二次曲面拟合法，即对于公共点上的高程异常与平面坐标之间，假定存在如下数学模型

式中，a0,a1,a2,a3,a4,a5为模型待定参数。

1.3 多面函数法

多面函数拟合法的基本思想是任何一个规则或不规则的连续曲面均可以由若干简单面(或称单值数学面)来叠加逼近。

多面叠加的数学表达式为

这里，q(x,y,xi,yi)为参加插值计算的简单数学面，又称为多面函数的核函数； n为简单数学面的张数或多层叠加面的层数，它的值与分块扩充范围内参与点的个数相等；ki(i=1,2,3,……,n)为待定参数，它代表了第i 个核函数对多层叠加面的贡献。

1.4 曲面样条拟合法

曲面样条拟合法是基于无限大平板小挠度方程的数学模型，设点的ζ与点的x、y 坐标存在如下样条关系

式(4)中ri2=(xxi)2 +(yyi)2 +δ，其中δ为平滑因子，一般可取δ为1。

由式 (4)可见，当有m 个已知的gps 水准联测点时，则有m + 3 个未知参数。(a0,a1,a2,f1,f2,……fm)，而方程个数亦为m +3个，因此可解此方程组求得这m + 3个参数，当gps水准联测点有n 个时，n>m，可按最小二乘法求解出这m + 3个未知参数，进而可推求其余各插值点的高程异常值ζi。

1.5 移动曲面法

移动曲面法在计算时，通常采用契比雪夫多项式为移动多项式。