应用多面函数最小二乘推估法计算GPS高程.docVIP

应用多面函数最小二乘推估法计算GPS高程 　　摘要：本文阐述了多面函数法计算GPS高程的原理和方法，结合生产实例，采用此方法计算的数据与已知常规测量方法、GPS随机软件计算方法即地球模型法数据进行了对比，从数据对比结果分析，不同地区，不同GPS高程计算方法与常规测量方法数据差异不同，总有一种方法与常规测量方法差异最小化，达到长输管道测量规范中的高程精度的要求。因此针对长输管道分布不同地区，要求计算GPS高程方法多样化，提出一种适应地区变化且精度较高GPS高程计算方法。 　　关键词：GPS 拟合；大地水准面差距；多面函数；误差 　　中图分类号：P228.4文献标识码：A 　　一、 引言 　　众所周知, GPS技术在测绘领域中应用取得了巨大的成功, GPS高程应用也得到了一定程度发展，由于GPS施测的WGS-84椭球面上的大地高, 是一个几何分量, 而生产中需要的是基于测区的似大地水准面上的正常高, 两者存在着一个差数即高程异常, 为一个变量, 此差数因地区的不同而不同, GPS高程应用关健是应用各种数学模型求取测区高程异常值，也可以理解为将基于WGS-84椭球面的GPS施测数字高程模型转换为基于测区水准面的数字高程模型DEM，转换方法多种多样。但并不是每一种方法计算GPS高程在全球所有地区都适合，而GPS厂家随机软件大都带有一种计算GPS高程的方法，这在生产中应用很不方便。 专家、学者对计算GPS正常高方法进行了探索，多面函数最小二乘推估法就是其中的一种。 　　二、 多面函数最小二乘推估法的原理 　　多面函数最小二乘推估法是美国HARDY教授于1997年推出的一种DEM求解的方法，基本指导思想：任何一个光滑的数学表面，总可以用一系列的有规则的数学表面的总和以足够的精度逼近。地形表面总是可以划分成不同的区域，这些区域可以认为是局部光滑、连续的。可以选择适当的数学曲面加以描述并逼近，以求得其表面DEM。结合实际将此方法应用到计算GPS正常高方面，提高将GPS施测的大地高转换为测区正常高的精度，探索出一种新计算GPS正常高有效途经 　　依据以上原理，设（X，Y）为地形表面上的任意一点，其高程可以表示为多个数学曲面的总和①，即 　　n 　　Z＝#402;（X，Y）＝Σαj q(X,Y,Xj,Yj) 　　j=1 　　 =α1q(X,Y,X1,Y1)+ α2q(X,Y,X2,Y2)+…+αNq(X,Y,XN,YN) 　　(1) 　　上式中q(X, Y, Xj, Yj)称为拟合核函数，（Xj,Yj)为中心点，核函数根据需要选用，在工程界常用的用于拟合GPS高程的二次曲面函数也是核函数的一种，此次还有如垂线偏差法、平面转换法、地球模型水准面法都属于核函数。这里选用对称的双曲面和圆锥面作为核数： 　　 　　q(X,Y,Xj,Yj)＝［（X－Xj ）2＋（Y－Yj）2＋δ］1/2 a 　　 　　q(X,Y,Xj,Yj)＝［（X－Xj ）2＋（Y－Yj）2］1/2 b 　　 (2) 　　 　　其中第一式a是一个对称的双曲面核函数，第二式b是一对称的圆锥面核函数，后者较间单，前者在数据点处能够保证更好光滑性，也就是计算GPS 高程外符合性很好。 　　 　　设数据点总数为m，其中有n个作为核函数的中心点（m≥n），则其余各数据点（Xi ,Yi）点应当满足下式： 　　n 　　Z = Σ αj qij(i=1,2,…,m)(3) 　　 j=1 　　其中 　　qij＝q(Xi ,Y i , Xj , Y j )(4) 　　 　　由式（3）可列出m个误差方程，矩阵形式为 　　 　　 V=Qα-Z 　　m,1 m,n n,1 m,1 　　 (5) 　　将上式法化，并得到多面函数的系数解为 　　α ＝（QT Q ）－1 QTZ 　　 m,1n,m m,n n,m m,1 　　 (6) 　　求得系数后，就意味着多面函计算的GPS 与已知数据点实现了最佳拟合，任意一点k的高程由（5）式中得到： 　　 　　ZK = qkα=qk(Q TQ)-1QTZ 　　1,n n,1n,mm,nn,mm,1 　　 (7) 　　上式中 　　 qk = [qk1 , qk2 ,…,qkm ] 　　 (8) 　　 　　qkj= q ( Xk ,Yk , Xj , Yj ) (9) 　　此外，常用的GPS高程计算的核函数有以下几种： 　　 　　 -Kdj 2 　　 q(dj)= e 　　 (10) 　　 　　 q(dj)=1/(1+( dj /K)2(11) 　　以上两式K为选定的常数。 　　 　　q(dj)=1+ dj2 （12） 　　3 　　q(dj)= ∑bkd