




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、GPS高程拟合的转换方法(河北理工大学)摘要:本文从 GPS定位系统的组成和介绍开始,分析研究了GPS高程的使用意义和不足,说明大地高、正常高和正高的概念及转换关系,阐述了GPS高程拟合的原理。讨论了绘等值线图法,解析内插法中的多项式曲线拟合方法、三次样条曲线拟合方法和Akima曲线拟合方法,曲面拟合中的多项式曲面拟合法、多面函数曲面拟合法、 曲面样条拟合法和移动曲面拟合方法。研究并分析了GPS水准拟合的精度评定。关键词:大地高,正常高,高程异常,高程拟合1引言GPS系统由GPS卫星星座(空间部分)、地面监控系统(地面控制部分)和GPSB号接收机图1.1全球定位系统(GPS)构成示意图GPS(
2、 Global Positio ning System/全球定位系统)系统是一种以空间卫星为基础的无线电导航与定位系统,是一种被动式卫星导航定位系统,能为世界上任何地方,包括空中、 陆地、海洋甚至于外层空间的用户,全天候、全时间、连续地提供精确的三维位置、三维速 度及时间信息,具有实时性的导航、定位和授时功能。GPS卫星星座21颗工作卫星和 3颗在轨备用卫星组成,这24颗卫星均匀分布在 6个轨道平面上。卫星轨道平面相对地球赤道平面的倾角约为55 ,各轨道平面升交点的赤经相差60,在相邻轨道上,卫星的升交距角相差300 ,轨道平均高度约为 20200km,卫星运行周期为11时58分(恒星时12h
3、)。这一分布方式,保证了地面上任何时间、任何地点至 少可同时观测到4颗卫星。GPS卫星的作用是接收和播发由地面监控系统提供的卫星星历。地面监控系统由主控站、 注入站和监测站三部分组成,它们主要负责编算 GPS星历将其发射到GPS卫星上,监测和控制 GPS卫星的“健康”状况,保持各颗卫星处于同一时间标准,即 处于GPS时间系统.。GPS信号接收机的主要任务是接收GPS卫星发射的信号,以获取必要的导航定位信息,并经数据处理而完成导航定位工作。当GPS卫星在用户视界时, 接收机能捕获到按一定卫星高度截止角所选择的待测卫星,并能跟踪这些卫星的运行:对所接收到的GPS信号具有变换、放大和处理的功能,以便
4、测量出GPS信号从卫星到接收机天线的传播时间,解译出GPS卫星所发射的导航电文,实时地计算出测站的三维坐标位置,甚至三维速度和时间。GPS水准高程拟合模型的研究是为了实现由GPS大地高求出水准高,用GPS水准代替几何水准,提高GPS水准测量的精度,发挥 GPS技术提供三维坐标的优越性;结合实际工程, 将各种模型应用到不同的地形条件下,经过模型优选,分析比较和精度评定,得出具体的结论,指导生产和工程实践。数据的处理,采用五种拟合方法编写程序实现似大地水准面拟合 的模型程序化,在一定条件下,拟合所得到的结果可以满足四等水准测量的精度要求,具有一定的使用价值。2 GPS水准高程基本概念图2.1大地高
5、、正高、正常高关系 (不考虑垂线偏差)大地高程系统是以参考椭球面为基准面的高程系统,它是一个几何量,不具有物理意义。大地高的概念:由地面点沿通过该点的椭球面法线,到参考椭球面的距离, 通常以H表示。GPS定位测量获得的是 WGS-84椭球空间直角坐标系中的成果,其中的高程值是地面点相 对于WGS-84椭球的大地高H。P的正高Hg定义为由地面点(水准面)中的一个,由于实际上是一条曲线, 正高正高系统是以大地水准面为基准面的高程系统,地面某点P沿垂线方向至大地水准面的距离。大地水准面是一族重力等位面水准面之间不平行, 所以,过一点并与水准面相垂直的铅垂直线, 的计算公式为:1Hg5gd H(2.1
6、)9m H g式中, g d H是地面点至大地水准面之间的位能差;gm为由地面点沿垂线方向至大地水准面的平均重力加速度。由于gm无法直接测定,所以严格的讲,正高是不能精确确定的。由于正高是以大地水准面为基准面,具有非常重要的物理意义,所以它在水利建设、管道和隧道建设等精密工程技术方面有着广泛的应用。若以hg表示大地水准面和椭球面之间的差距,则正高与大地高的关系按如下公式可得:=:g - hg(2.2)由于gm无法直接测定,导致正高无法严格确定。为了方便使用,根据前苏联大地测量学学 者莫洛金斯基的理论,建立了正常高系统。其公式为:任意点处的大地水准面与椭球面的差值称为高程异常,正常高与大地高的转
7、换关系为::.-(2.3)其中 为似大地水准面的高程异常。1H rgd H( 2.4)m3 GPS测量高程原理GPS是一种三维观测系统,通过相位观测值可求出网中每两点间的地心WGS-84坐标系中的坐标差 X= C:Xo,厶丫0,Zo)T,提供了地面点间位置和高程信息。如何求出地面点的 高程(正常高)需要经过一些中间步骤,现介绍其基本过程。GPS测得的基线向量 厶X,以坐标为未知参数进行自由网平差,求出该网点三维地心坐标。取网中至少三个己知地面控制点,其点位的大地坐标经度L、纬度B和大地高H为已知,将这些点上的己知数转换到相应椭球的三维直角坐标系中,求出坐标值。转换公式为:2一(N (B) +
8、H ) cos B cos LX =y(N (B) + H ) cos B sin Lz_i i2(1 _e ) N (B) + H ) sin B式中= :- / 1 -e2sin匕,:- f是相应的椭球参数。以已知点上大地直角三维坐标X为控制,采用七参数法,将 WGS-84系直角坐标转换至与控制点相应的直角坐标系中,公式为:Ay710 0Xs0-zsysyys+010yszs0xszs 一001zs-ysxs0 J占z;y式中 X 冷、工为平移参数,m为尺度比,y、z为旋转参数。由此求得 GPS的直角坐标,再经下式变换,即得与已知点相同椭球的经度、纬度和大地高:22B = arctge2N
9、 (B)N(B) + H L = arctg( y / x)2 2x y(正常高),两个基准面cos B N (B)众所周知,大地高是地面点至椭球面的高程,我们需要的是海拔高程 之差为该点的高程异常,即椭球面至似大地水准面之间的高差,表达式为:H y = H _ Q式中Hr,H分别为正常高、大地高,为高程异常。另一种过程是求出正常高高差,现简述其原理如下。由GPS获得的基线向量(Axs, Ays, Azs),通过网中至少三个已知点,经 (2.8)式变换到三维直角坐标系中,采用四参数法,公式为:0zs-ys一 zs0xsys22(2.10)式,变换到椭球面求得GPS测点从属于地面控制坐标系的三维
10、坐标差。有类似的 上,求得基线向量的精度差、纬度差和大地高差。再由下式求得正常高高差:占H 丫=-2由于坐标系转换时采用了七参数或四参数法,GPS网点经转换后与相应地面测点仍有间隙,亦即转换后GPS网点坐标与地面网系统并不兼容,这对于城市网和工程控制网而言,不能说不是一个问题。因此,在我国自行制作的软件中采取了一系列措施,使GPS网点坐标和地面网系统兼容。在地形起伏较大的地区及位于较高建筑物上的未知点,用传统水准测量方法测定这些点的高程较为困难,通常采用三角高程测量的方法。三角高程测量是以水平面为基准面和视线成直线为前提的。因此,只有当A、B两点间的距离很短时,才比较准确。当A、B两点距离较远
11、时,还应该考虑地球弯曲和大气折光的影响了,另外三角高程加测天顶距, 计算边长,则会不可避免的带来误差传播,并增加工程量。对于一般工程单位而言, 无法获得必要的重力数据和天文数据,因此重力方法在一般的工程单位很难施行。相比于前面的方法,GPS测高只需摆好接收机,开机等待即可,测站距离可以达几 km,劳动强度小,不受气候影响, 工作效率高,内业较为简单。因此在点多量大的水利工程测量中,如果能用GPS测高代替四等水准,其优越性是非常明显的,特别是在山区更为明显。4 GPS高程转换方法GPS水准高程,是目前 GPS作业中最常用的一种方法。目前,国内外用于 GPS水准计算的各种方法主要有:绘等值线图法;
12、解析内插法;曲 面拟合法等。下面重点介绍解析内插法和曲面拟合法的多项式曲面拟合法和多面函数拟合法。(1) 线形解析内插法1, Akima曲线拟合法Akima法的原理是:Akima插值法规定在两个实测点之间进行内插,除需要用到这两个实测值外,还要用这两个点相近邻的四个实测点上的观测值。也就是说,在两个实测点之间进行内插,共需六个实测点。设已知数据点为(Xj, y(i =1,2,3,n),现在要找一条光滑曲 线y = f(x),使得满足i = fg)。所谓“光滑”是指 f(x)具有连续的一阶导数。同时还约定, 在任意两相邻的数据点之间,用三次多项式来逼近。设给定n个不等距GPS测点为x0:乂勺:x
13、2xnj(或y0:y: y2”:”: yn = )。相应的 高程异常值为 i(i =0,1,2,n -1)若在子区间上Xk,Xk.J (k =0,1,2,,n 2)上两个端点处满足以下4个条件:二k = f (xk)-k-1 = f (Xk +)k(3.1)-k =g k严k 卡=g k 十式中gk, gk+i可由从Akima条件唯一确定。则在xk, xk “ (k = 0,1,2,n 一 2)区间上确定 一个三次多项式:(3.2)23s(x)二 S。 Sjx Xk) S2(X Xk) S3(X Xk)由上式可计算该子区间插值点t处的高程异常值-。其中:s0 = (3.3)51 =tk52 =
14、3(:_k 卡 一 rk)/(Xk* Xk) 2tk tk*/(Xk 卡一Xk)严严253 -tk k1 _2( k 1 - k)/(Xk:1 - xk)/(xk:1 - XQ上式中的tk,tk+1为k号和k+1号点实测要素的斜率,tk用k-2、k- 1、k、k+1、k+2已知点 计算,tk+1用k- 1、k、k+1、k+2、k+3已知点计算,一般计算公式为:ti =(同十-mm-+|m一mi口)/(血十一m+|mm)(3.4)式中 i 二k,k 1 mi =( )/(Xi .1 Xi)当上式分母为零时,ti =1/2(m亠mJ或L =mi2三次样条曲线拟合法当测线长,已知点多,变化大是,在R
15、二min条件下解各ai误差会增大,故通常总采取分 段计算。这样使曲线在分段点上不连续,也影响拟合经度。为此,采用三次样条法来拟合。三次样条曲线实际上是由一段一段的三次多项式曲线拼接而成的连续曲线。在连接点处,不仅函数自身是连续的,而且其一阶导数和二阶导数也是连续的。这样既保留了多项式在表达上的简便性,又克服了单个多项式不灵活、不稳定的缺点,且计算较简单,故在长测线似大地水准面拟合中得到了应用。设过n个已知点,和(或yi或拟合坐标)在区间xi ,xi J (i =0,1,2,n -1)上有三次样条 函数关系(X)二(Xi) (x Xi) (Xi , Xi 1) (x - Xi )(x - Xi
16、J (x,Xi , Xi 1)(3.5)式中X为待定点坐标,Xi , Xi 1为待求点两端已知点的坐标。其中(Xi,Xir)为一阶差商。-(X i , x i 1 ) = ( - i 1 - i ) /( X i 1 - xi );-. (X , x i , xi 1) 为二阶 差 商 , (x,Xj, Xi 1)= 1 (Xi ) (x) (Xi .1),而 (xj (i =1,2,,n 一1),满足系数矩阵 6为对称三角阵的线性方程组(Xi Xi) (Xi J 2 (Xi 1.Xi丄厂(Xi) (Xi 1xj (Xi 1)=6 (Xi,Xj 1 ) -:.(Xi,Xi)(3.6)(x0)=
17、(xn )=0用追赶法解(3.10)方程组,可求出(xi )和(xi,xi ,),而(3.7)(x)(Xi ) (X - Xi ) (Xi , Xi .1)事实上,我们可查询有关数值分析和其论文,一般采用下式来拟合:3r.(X i - X)(X-Xi)i (x)-Mi亠 M i -6hi6hiMi ; yi6 IXj X(i =1,2, n 1)(3.8)hi=XB ;(3.12)=XB ;(3.12)式中i(Xi), hjMXiXi, y i (xi)。待定参数 Mi同样可由对称三角阵的线性方程组解得。(2)曲面拟合法当GPS点布设成一定区域面时,在GPS网中联测一些水准点, 然后利用这些点
18、上的正高和大地高求出它们的高程异常值,再根据这些点上的高程异常值与坐标的关系,用最小二乘的方法拟合出测区的似大地水准面,利用拟合出的似大地水准面,内插出其他GPS点的咼程异常,从而求出各个未知点的正常咼。1多项式曲面拟合设点的与平面坐标x、y (或大地坐标B、L或二B、二L),有以下关系:=f(x,y) + ;(3.9)或厶=f ( :x, y) + :;(3.10)nJ其中 f(X, y)为中趋势值,:为误差、Ax =X-Xo、Ay 二y-yo、x Xi) / n ,t=0n丄y = ( yi) / nt zfi设 f (x, y)二 a0 .a2ya3x2a4xy a5y2a6x3a7x2
19、y a8y2x a9y3 -(3.11)写成矩阵形式有:=XB ;(3.12)其中:对于每一个已知点,均可以列出上式方程,在-1X。y。X:1X =1a x12X1亠1Xn丄yn27xnZ2Z=min条件下,可求解系数阵T1T “B = (X X ) -(X )(3.13)再由已知高程异常的权阵情况下,上式可改写为:1 .B =(X PX ) (X P )(3.14)系数求出后,再按(3.12)式求出待求点的Z,从而求出Hr。实践表明,在地势较为平坦的地区,当已知高程异常的点,密度适当,分布比较均匀时,该法计算高程异常的精度,可达厘米级。在难以达到代替三、四等水准测量要求的大地水准面的起伏较大
20、的山区,一般可采用根据测区地形情况,适当增加已知高程异常点的密度,并改善其分布;综合利用测区的重力观测资料,改善模型对高程异常的分辨率;考虑到高程异常与测区地形的密切相关性,在模型中可引入地形影响的改正项,以提高确定高程异常的精度等方法加以改善。下面列出几个常用的多项式内插公式:平面模型:f (x, y) = a0 a1x - a2y平面相关模型:2f (X, y)二 a。a1x - a2 y a3xa4xy二次曲面模型:2 2f (x,y) =a。 ax a?y asXaqXy a三次曲面模型:f (x, y)二a。 ax a? y a3x2 a4xy -a5y2a6x3a?x2ya&y2x
21、 agy3在本文论述的其它拟合方法中,其中移动曲面的基本数学模型采用的是多项式形式,由此可见对它的模型的研究和探讨有着重要的意义。其次,要说明的是,最小二乘原理和参数估计始终贯穿所有数学模型的计算,同时也是测量平差的基础。182多面函数曲面拟合多面函数拟合法,1971年由美国哈笛(Hardy)提出。1976年将此法应用于美国大地测量、 拟合重力异常、大地水准面差距、垂线偏差等,1978年将此法用于地壳形变。它的基本思想是任何一个规则或不规则的连续曲面均可以由若干简单面(或称单值数学面)来叠加逼近。具体做法是在每个数据点上建立一个曲面,然后在方向上将各个旋转曲面按一定比例叠加成一张整体的连续曲面
22、使之严格地通过各个数据点。多面叠加的数学表达式为:m二 f (x, y)aiQ(x,y,Xi, yj(3.15)i =1式中ai为待定系数;Q(x, y,Xi,yi)为x和y的核函数,其中心在(Xi,y)处,简称多面函数。常 用的简单核函数有以下几种。即:锥面Q(x,y,Xi, y二 C ; (x - Xi r (y 一 yi佝双曲面: 2 22Q(x, y,xi, yi) = (x - Xi)(y - y:)、(3.17)倒曲面i1222-Q(x, y,Xi, y=1/ . (x - x(y-y、(3.18)这里为非零参数。上式表示一段双曲面绕竖轴旋转而成的曲面,当时,此曲面就退化成圆锥面。
23、三次曲面3 Q (x, y, Xi , yi C (XXi)2 (yyj2、22(3.19)上式是母线为三次曲线的旋转面。在上式各式中,(x Xi)2 (y -yj2为内插点到参考点之间的水平距离的平方。旋转面Qa式中a为参数。母线旋转面Q = C0 exp_a2Dj2上式是以高斯曲线为母线的旋转面,C。和a为两个参数,其中:为光滑系数。设已知参考点有 m个,选择其中n个点作为节点(j =1,2,n_m),记B =22,an)T,则同样有(3.佝误差方程,式中X为核函数矩阵。当已知数m大于节点数n时,利用最小二乘法原理,可以根据(3.17)计算方程系数,在已知核函数权阵的情况下,可以根据(3.
24、18)计算方程系数。系数求出后,同样按(3.15)式求出待求点的从而求出Hr。如果已知数m等于节点数n时,方程唯一解,可以按线性方程求解,如下式:根据参与拟合计算已知点高程异常值i与计算后得到的高程异常值:,用Vi= 7 - i求拟这时待定点的解为:1 .p = X pB = X pX 一(3.21)这里补充说明,对于已知参考点个数等于误差方程系数的情况下,方程存在唯一解,只可以线性求解。m个已知点,要求是高程异常显著点, 即能很好描述该区域内高程异常分布 的特征点,最好位于最高、最低及坡度变化处。这种拟合法,虽然理论上严密,但使用上有些困难,特别是、:和核函数的选取,需要不断试验改进选取。下
25、面就其中三个问题进行深入分析。1. 核函数的选择多面函数是属于纯数学逼近方法,其差值精度与选用的核函数有密切的关系,高程变化大的地区内插的精度往往较低,所以我们要选择一个一成不变的核函数对各种地形进行最小 二乘推估。只要高程变化大的地区保证了精度,对于平坦地区的推估,其核函数影响是不大的。对于待求点的估值实质上就是数据点的加权平均值。因而选择优良核函数的准则是:由它计算的权系数只要相对集中,并且数值稳定性好。在对地壳垂直运动中速率面的拟合采用 正双曲线函数;2. 光滑系数:的确定优化的/选取对其高拟合效果有作用,但比较困难。试验表明越大,内插的曲面越平滑,但当大于某一值时将使结果出现很大的偏差
26、。/对于倒双曲面函数,:必须大于零,否则无法计算,同时、的取值也有正双曲面的性质,/越大,内插的曲面越平滑。 但当大于某一值时也 将使结果出现很大的偏差。3. 核函数节点的选取在大范围有较多的 GPS水准点,即已知高程异常的点较多,可选其中部分作为节点,其余作用拟合高程异常的检核。此时,拟合结果就与所选节点不同而异,这是需要试验研究的问题。对于局部 GPS网,GPS水准点不会很多,就可全部用来作为节点,进行拟合。5拟合精度评定为了能客观地评定 GPS水准计算的精度,采用内符合精度和外符合精度检核,内符合 精度表明的是结点与所选模型的符合程度,而外符合精度表示的是检核点与所选模型的符合程度。1内符合精度合残差V,按下式计算GPS水准拟合的内符合精度 U : VV /(n 1)(5.1)式中n为参参与计算的已知点个数。2外符合精度根据参与检核点的高程异常值与计算后得到的高程异常值:,用m= - i求拟合残差Vi,按下式计算GPS水准的外符合精度 M :M 二::;;W /(n 1)(5.2)式中n为参与检核点的个数3 GPS水准精度评定内
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年云计算服务模式演变下的云计算服务市场研究报告
- 2025年元宇宙虚拟艺术品市场交易活跃度分析与未来趋势研究报告
- 数字化金融生态2025年开放银行构建与合作模式创新趋势研究报告
- 2025年医药行业CRO模式下的临床试验方案设计与优化报告
- 新一代大学英语(第二版)综合教程1-U1-教师用书 Unit 1 A new journey in life
- 2025年医药企业研发外包(CRO)服务标准化与行业规范化报告
- 线下演出市场复苏中的市场潜力分析与竞争格局报告
- 2025年船舶制造行业订单分布与节能环保造船技术研究报告
- 工业互联网平台SDN网络架构优化与工业互联网平台可持续发展报告
- 北京安全监理试题及答案
- 2025年中国PCR仪市场全面调研及行业投资潜力预测报告
- 2025年医院信息科数据安全管理计划
- 智能工厂整体解决方案
- 2025年四川大学自主招生个人陈述的自我定位
- 2025年福建省建工集团及下属集团招聘235人高频重点提升(共500题)附带答案详解
- 紧急情况的处理措施、预案以及抵抗风险的措施
- 2025中智集团招聘重要岗位高频重点提升(共500题)附带答案详解
- 水暖维修培训课件
- 大学生心理健康教育知到智慧树章节测试课后答案2024年秋宁波大学
- DB33T 1152-2018 建筑工程建筑面积计算和竣工综合测量技术规程
- 临床路径变异分析
评论
0/150
提交评论