GPS 高程拟合的方法及实现 - 湘潭市勘测设计院.PDF

GPS 高程拟合的方法及实现 颜家逵 湘潭市勘测设计院，湖南湘潭 411100 摘要： GPS 高程拟合已经拥有众多的方法，对于这些方法拟合的精准度，一直 以来都是业内的所关注的一个热点，本文对这些拟合方法进行了综合分析。结合 具体的实例，以已知的原始数据为基础，在设立测量区域的似大地水准面的过程 中，组合使用了包括二次曲面拟合法、二次曲面加权法、二次曲面多面函数法、 二次曲面最小二乘法 在内各种拟合方法，进行综合分析，从比对各种方法获得 的数据结果来看，组合方法能够使检验点的拟合更加精准，亦能够使外检验点的 拟合更加精准，是十分有效的拟合模型。因此综合利用多种拟合法，比单一地使 用二次曲面拟合法， 该高程拟合方法具备更高水准。 关键词:GPS 高程拟合，二次曲面法，组合法，分析 1 引言 对于 GPS 高程拟合，目前已有多种方法，如何利用这些方法，使拟合结果更 加准确可靠呢？ 2 高程系统的分类和概念及其它们之间的转换关系 2.1 高程系统包括三种，分别是大地高系统，正高系统，及正常高系统。 在我国的测绘工作中，统一使用的高程系统是正常高系统。 大地高系统的概念为：以地球椭球面为基准面，地面一点沿参考椭球面的法 线到参考椭球面的距离，在公式中通常写作 H， 正高系统的概念为：以大地水准面为基准面，地面点沿通过该点的铅垂线至 大地水准面的距离，在公式中通常写作Hg。 正常高系统的概念为：以似大地水准面为基准面，地面点沿通过该点的垂线 方向到似大地水准面的距离，在公式中通常写作 Hr。 2.2 各高程系统之间的转换关系。 大地水准面和地球椭球面之间的距离，叫做大地水准面差距，通常写成 N， 那么大地高和正高之间的关系则是H Hg  N 似大地水准面和参考椭球面之间的距离，叫做高程异常，通常写成 ζ， 那么大地高和正常高之间的关系则是 : H Hr ζ 3 多项式曲面拟合法的组合法 3.1 组合方法 为了使高程异常更加精准可靠，使用地球重力场模型，平均重力异常及数字 高程模型。然后利用从各种模型中获得的数据，推求出最可靠的高程异常。推求 公式为: ζ ζ ζ ζ ，各符号所代表的内容如下：[1] GM g  地球重力场， 代表似大地水准面之长波特性，在波长 100 到 200 千米 ζ GM ζ GM 之中呈现均匀的变化， ζ g 重力异常，ζ g 代表局部似大地水准面的特性，通常波长在 20 到 100 千米之中。 ζ 地形起伏对高程异常的影响, ζ 代表 20 千米 以下的波长特性。 通常，小区域的 GPS 网中 ，于高程异常的中长波项而言，大地水准面几乎接近 似大地水准面，同时没有明显的变化状态，合适的曲面拟合函数往往就可以解决 问题。而短波项就有所不同，它的产生原因在于起伏的地形，地形起伏若是微小， 似大地水准面就愈趋平整，地形起伏若是十分明显，它就愈为粗糙。对于普通的 工程来说，不会有过广的测绘范围，依据前面的知识，可以将高程异常划为两个 方面的内容， 中长波项，记作 ε，短波项 ，记作ζ ，写成公式①如下： ζ (ζ ζ )ζ  ζ GM g   ε表示高程异常的中长波项，它是一个趋势值 （与地面点坐标有关的函数）， ξr 是高程异常的短波项，它与地面地形信息有关。通过公式①，得知每个 GPS 水准点的高程异常都由两部分组成，即由中长波项的趋势值和短波项。假若只是 单一地采用多项式曲面拟合法，仅仅可以解决中长波项的拟合问题，对于短波项， 就必须进一步进行分析和研究。首先对观测方程和误差方程进行变换，得到如下 公式②③： ζ BX ζ    V BX ζ  f  通常，③中的最小二乘解基于如下目标函数 : T T V PV   P  min 1  2  前边提到的 “综合极小”法，得到的 ξ常常比实际的高程异常要小，这是由于 在测量时， 