浅谈GPS高程拟合技术.docVIP

浅谈GPS高程拟合技术 摘 要介绍了GPS技术在高程拟合方面的应用。 关键词 GPS技术 高程拟合 高程异常。 前言 GPS（Global Positioning System）即全球定位系统，是1973年美国国防部为了满足军事部门对海上、陆地和空中设施进行高精度导航和定位的要求而研究的新一代高精度卫星导航系统。GPS是以人造卫星为基础的无线电导航系统，它是利用天空中均匀分布的24颗GPS卫星轨道参数及其载波相位信号，通过地面接收设备接收其发射信息，实时地测定地面接收载体的三维位置。 我院从1999年开展了GPS技术在公路勘测中的应用研究。几年来的生产实践，我们认识到了GPS技术在平面控制测量和路线中桩、边桩放样方面具有传统测量工作不可比拟的优势，可以极大的降低劳动作业强度，提高作业效率，但GPS技术在高程测量方面的应用还一直处于研究状态。本文结合几年来的生产实践仅就GPS技术在高程拟合方面的应用谈谈自己的观点： 高程异常 GPS测得三维坐标高程为各GPS点在WGS—84坐标系中的大地高H，而公路勘测所用的地面高程是相对于似大地水准面的正常高H正，两者之间的差值称为高程异常，用公式可表示为： ζ=H—H正 式中：ζ—为高程异常 要将GPS所求的大地高转换成正常高，关键是求得精确的高程异常ζ。 目前通常采用二次曲面函数对高程异常进行曲面拟合，对于GPS水准联测点PK拟合模型可写为 ζK=a0+a1Δxk+a2Δyk+a3Δx2k+ a4Δy2k+ a5ΔxkΔyk—εk 式中Δxk=xk—x0 Δyk=yk—y0 x0 ，y0是参考点的坐标，一般取重心坐标；xk ，yk是Pk点的平面坐标，也可是大地纬度和大地经度；εk为拟合残差。按最小二乘法可求得拟合系数a为 a=（ATA）-1ATζ 式中a=[ a0 a1…an]T ζ==[ζ0 ζ1…ζn]T 1 Δx1 Δy1 Δx21 Δy21 Δx1Δy1 1 Δx2 Δy2 Δx22 Δy22 Δx2Δy2 A= … … … … … … … … … … … … 可以看到，在采用二次曲面拟合时，至少应有6人GPS水准联测点，当少于6个时，则应去掉二次项拟合系数σ3，σ4，σ5，即采用平面系数拟合，此时拟合模型为 ξk=σ0+σ1Δxk+σ2Δyk-εk 求得拟合系数后，任意点Pi的高程异常为 ξi=σ0+σ1Δxi+σ2Δyi-σ3Δxi2+σ4Δyi2+σ5ΔxiΔyi 或为ξi=σ0+σ1Δxi+σ2Δyi 当然还有基它一些方法：绘等高直线图法、解析内差法、滤波推估法，但这些方法在实际操作中计算量过大。为了更加准确、快速的消除高程异常ξ的影响，求的准确的大地水准高，我们可以借助GPS500的SKI-Pro软件解决。 生产实践 在2003年龙头至色家二级路（K12+280-K13+580段）勘测中，首先我们利用三台瑞士Leica公司生产的GPS500双频接收机对路线走向进行了静态观测，有效观测时间均大于60min，随后借助SKI-Pro专业软件对静态观测成果进行基线处理和网平差，求得满足要求的转换参数。在中桩放样中我们采用RPK技术，同时对BM13至BM15及其之间的中桩逐一进行了测量，随后利用Leica　GPS500的SKI-Pro专业软件对BM13单个水准点、BM13至BM14两个水准点、BM134、BM15三个水准点间的中桩地面高程分别进行拟合： （1）新建项目； （2）输入原始数据既测量数据，并将原始数据分配到项目中； （3）新建坐标集； （4）输入ASCII数据，并将ASCII数据分配到新建坐标集中，ASCII数据应以记事本形式以如下格式（点标识，东坐标（y），北坐标（x），高程（h）给出： BM13，18769．1838，200296．2304，109．408 BM14，18070．6723，200402．9440，101．894 其中坐标为地方坐标，且需与项目中对应点坐标保持一致，高程为水准基平测得的大地水准高。 （5）基准与投影，上面选中项目，下面选中坐标集，利用插值法（平面+高程（正高））进行匹配，结果显示残差为零，将此坐标系统保存； （6）回到项目，选择上面图标中LOCAL以及格网坐标，则显示所有点的当地坐标及拟合后的正高。 结果显示两个或三个水准点间地面高程拟合误差较小，下面是BM13至BM14间地面高程拟合对比： 表1GPS与中平测量高程对比 桩号 中平测值 GPS测值 差值 BM13 109.408 K12+240