GPS高程测量原理及方法探讨.docVIP

广东科技2010.2总第231期 GPS 高程测量原理及方法探讨 谢劲松 (广东省广州市510000 1引言 GPS (欲了解更多?请见本期【科技“生词”解释】技术的出现,为确定大地水准面高提供了新的途径,提高了作业的效率。然而我们的实用高程采用的是以似大地水准面为基准的正常高。因此,我们必须要实现GPS 大地高向正常高的转换,从理论上讲,实现GPS 大地高向正常高转换最好的方法是综合利用GPS 测量数据、重力测量数据和地球重力场模型进行转换。然而,对于一般工程单位来说,考虑到作业成本的问题,人们不可能花那么多的经费去获取昂贵的重力资料。本文就是探讨如何结合工程实际 情况,利用较少的经费获取较高精度的GPS 高程, 从而实现低成本、高效率、高质量的测量成果。 2常用高程系统的基本定义 大地高:大地高系统是以参考椭球面为基准面的高程系 统。某点的大地高是该点到通过该点的参考椭球的法线与参考椭球面的交点间的距离。一般用H 表示。 正高:正高系统是以大地水准面为基准面的高程系统。某点的正高是该点到通过该点的铅垂线与大地水准面的交点之间的距离。一般用H 正高表示。 正常高:正常高系统是以似大地水准面为基准的高程系统。某点的正常高是该点到通过该点的铅垂线与似大地水准面的交点之间的距离。一般用h 表示。 高程异常:高程异常是指似大地水准面到参考椭球面的距离。一般用ζ表示。 大地高与正常高之间的关系可以表示为: ζ=H-h (13GPS 高程测量原理 利用GPS 可以测量出高精度的WGS84三维坐标(即大地经纬度和大地高,但是我国大多数用户需要的是正常高(例如 1985国家高程基准等 ,那么就需要将大地高转换成为我们需要的正常高,由图1可以知道,h=H-ζ, 由此可以知道,将GPS 大地高转换成正常高的关键就是求出GPS 点上的高程异常值 ζ。因此, 在利用GPS 确定了高精度的大地高后,求正常高的过程实际上就是求高程异常的过程。 4常用GPS 水准高程计算方法介绍 通过前面的介绍我们知道,由大地高转换到正常高的过程 就是求高程异常值ζ的过程。目前,求高程异常值的方法有很 多,比如GPS 三角高程、曲面拟合法、GPS 重力高程、 绘等值线图法、解析内插法等等。在实际工程应用中考虑到生产成本的因素,我们要根据实际情况选择既能提高工作效率,又能提高成果精度的方法。 下面简单介绍一下几种常用的确定高程异常的方法: 4.1绘等值线图法 绘等值线图法是较早使用的GPS 水准方法。其原理是:假 设在某一测区,有m 个GPS 点, 用传统的几何水准方法联测其中n 个点的正常高(即当作已知点,根据GPS 观测获得这些点的大地高,按式(1求出这n 个已知点的高程异常。然后,选定 适合的比例尺,将这n 个已知点的平面坐标 (平面坐标经GPS 网平差后获得展绘在图面上,并标注上相应的高程异常,再根据测区测量精度要求选取1~5cm 的等高距,绘出测区的高程异常图。在图上展绘其他待求点的平面位置,通过内插法计算出高程异常,从而求出这些待求点的正常高。 4.2曲面拟合法 当GPS 点布设成网状时,可用一个多项式曲面来拟合高程 异常。其原理是,根据测区中已知点的平面坐标x 、y 和ζ值, 用数值拟合法,拟合出测区似大地水准面,再内插出待求点的高 程异常ζ, 从而求出待求点的正常高。设测站点的高程异常ζ与坐标x,y 之间存在以下函数关系:ζi =f (x i ,y i +v i 式中:f (x,y 为的趋势值,v 为误差。我们可用以下空间曲面表达式来拟合: f (x i ,y i =a 0+a 1x i +a 2y i +a 3x i 2+a 4x i y i +a 5y i 2+a 6x i 3+a 7x i 2y i +a 8x i y i 2+a 9y i 3 (2式中:a i 为待定系数,于是,如果有n 个点,则可以得到下面的矩阵形式。 ζ=AX+v 式中:ζ=[ζ1,ζ2,ζ3,……,ζn ]T v=[v 1,v 2,v 3,……,v n ]T X=[a 1,a 2,a 3,……,a n ]T A=1 x 1y 1x 1y 11x 2y 2x 2y 2 1x 3y 3x 3y 3 ………… 1x n y n x n y n 对于每个已知点,都可列出以上方程,在∑v 2=min 条件下, 摘要:在控制测量过程中,GPS 测量由于其高效性、 经济性、实时性等特点已经在平面控制测量中得到了广泛的应用。但是在高程控制测量方面却没有得到更好的推广。长期以来, 水准测量是提供正常高的主要技术手段,它具有原理简单、误差易于检验的优点,但是其长距离作业时效率低、 作业强度大等缺点却也是不可避免的。因此,人们迫切期望能够用GPS 高程测量代替传统水准测量。本文根据近年来以美国