利用GPS拟合方法代替四等水准的可行性分析

梁先兵，李兆雄

（1.国家测绘地理信息局 重庆测绘院，重庆 400015）

利用GPS拟合方法代替四等水准的可行性分析

梁先兵1，李兆雄1

（1.国家测绘地理信息局 重庆测绘院，重庆 400015）

对GPS高程转换的3种方法(二次多项式拟合、三次多项式拟合、多面函数法)进行实际计算，通过综合比较和分析，给出了相同地形不同拟合方法的精度评价。结果表明，在一定条件下，利用GPS测得的高精度平面坐标和己知的大地高、正常高，采用多面函数拟合法所得到的结果可以满足四等水准测量的精度要求，在测量中可以取代传统几何水准测量。

大地测量定位；GPS高；高程拟合模型；内符合精度；外符合精度

由于测绘新技术的飞速发展，GPS测量技术在很大程度上代替了传统的导线或三角测量。然而，由于协议坐标系与大地坐标系的基准面不重合，GPS在获得高精度的高程成果时仍需借助传统的水准测量，这一直是GPS测量试图突破的瓶颈。

国外从20世纪80年代就开始探索用GPS测定正常高的理论与方法。早期研究表明，在地势相对较平缓的地区，用GPS水准方法能获得cm级精度的水准成果，可达到三、四等水准精度；在地形起伏较大的区域，若用GPS结合重力测量的方法也能获得cm级精度的水准成果[1]。1997年芬兰在10 000 km2测区内，用6台Ashtechz-12接收机，观测96个GPS点的2期数据，并利用OSU91A全球重力场模型和芬兰重力模型FIN95，获得点位的高程精度为15 mm[1]。武汉测绘科技大学在河北的实验结果[2]表明，联测1/5的几何水准，其GPS高程拟合的精度可达2 cm。青海石油局[3]在吐鲁番某区域测量近300个环，并用GPS拟合获得水准成果的精度为1.84 cm。郑州测绘学校[4]对西南山区某地GPS网的高程异常的拟合精度可达1.9～5.4 cm。该地区地形起伏较大，若加入地形改正，精度会显著改善。

1 GPS高程拟合原理和方法

大地高H（椭球高）是指地面点沿法线至WGS84坐标系参考椭球面的距离；正常高h是指地面点沿铅垂线方向至似大地水准面的距离；似大地水准面与参考椭球面的距离称为高程异常，见图1。其关系[5]为：ζ=H-h。只要有了大地高，再通过其他方法取得高程异常，就可以求得正常高。

在一定区域中，当测区中有一部分点已用GPS定位技术求得其大地高及用常规水准测量的方法求得其正常高，则可以计算出这些已测点的高程异常。若测区中已测点的数量足够多且分布较为均匀，就可根据测区内已测点上的高程异常值构造出某种曲面来逼近似大地水准面，进而推算出测区中未进行水准联测的GPS点的高程异常,从而获得未测水准点的正常高。

利用GPS高程进行拟合的方法[6]随着GPS技术的不断发展，取得了突飞猛进的发展，国内外的高程拟合法主要有：绘制等值线图法、多面函数拟合法、函数内插法、神经网络法、地球重力场模型法、数学模型抗差估计法和数学模型优化法等。

图1 大地高、正常高、高程异常关系略图

2 GPS拟合水准精度评定

为客观地评定GPS水准计算的精度，在布设几何水准联测点时，应适当多联测几个GPS点，起点应均匀地分布全网，以作外部检核用。

1）GPS水准精度评定。用GPS水准求出的GPS点间的正常高程差，在己知点间组成附合或闭合高程导线，再将计算出的闭合差与《国家三、四等水准测量规范》[7]中允许残差进行比较，以衡量GPS水准达到的精度。

3）外符合精度。根据检核点的ζi与拟合值之差，计算GPS水准拟合的外符合精度：其中n为参考点的个数。

3 项目实例

某项目测区地形较为复杂，部分为山区，高差较大，平面成果由该地区C级GPS观测数据经约束平差所得，1985国家高程基准成果由该测区二等水准测量平差所得，因此，成果数据精度较高满足本文拟合要求（见图2）。本文对曲面函数拟合法中的二次多项式、三次多项式拟合法和多面函数拟合法进行对比，利用已有水准路线观测成果，给出检核点到最近已知点的路线长，计算检核点拟合残差的限值，评定水准拟合的外符合精度和内符合精度。

由表1、2数据分析可知，二次多项式、三次多项式和多面函数拟合3种方法在内符合精度上差别不大，在外符合精度上多面函数拟合精度要明显优于二次多项式和三次多项式拟合模型。根据《国家三、四等水准测量规范》计算其限差均为超限。所以若多面函数拟合选择点足够多，且点位分布均匀，在拟合精度上能达到四等水准测量的精度。

表1 3种拟合模型拟合后项目控制点内符合精度统计表/m

表2 3种拟合模型拟合后项目控制点外符合精度统计表/m

P228

B

1672-4623（2016）03-0083-02

10.3969/j.issn.1672-4623.2016.03.026

2015-06-25。