分析小地区GPS高程拟合方法.docVIP

分析小地区GPS高程拟合方法 　　[摘要]本文对小区域GPS高程拟合方法精度进行了分析。首先从GPS高程测量原理出发，探讨了多项式函数拟合法的数学模型，对高程异常进行估算。 　　[关键词]GPS 似大地水准面 正常高 精化模型 　　[中图分类号] P228.4 [文献码] B [文章编号] 1000-405X（2013）-11-119-2 　　随着GPS定位技术的广泛应用，人们已经能够在10-6至10-9的精度量级上简捷而经济地获得所测点位的平面精度，但却一直未能以相应的精度解求点的高程，原因是虽然GPS能给出高精度的大地高，却由于没有一个具有相应精度高分辨率的似大地水准面模型，致使GPS大地高到GPS海拔高的转换过程中精度严重丢失。 　　1GPS高程测量原理 　　全球定位系统可同时获得相对精度较高的测站点的三维坐标，即大地经度L、大地纬度B和大地高H。对大地经度L、大地纬度B可以采用严密的数学公式将其转化成高斯平面坐标（x，y），其精度已被人们认识和接受。但GPS测高数据却没有像平面成果那样被广泛应用，这主要是因为GPS测量是在WGS-84地心坐标系上进行的，它所测得的高程是测站点相对与WGS-84椭球面的大地高，而在普通测量和城市测量系统中，地面点的高程采用正常高系统，即地面点沿铅垂线到似大地水准面的距离。由高程系统理论可知，测站点的大地高H与正常高h之间有如下关系： 　　h=H-ζ..........（1） 　　式中，F称为高程异常。 　　由式（1）可以看出，若能求出GPS点的高程异常，就可确定GPS点的正常高。因此，GPS高程转换的关键在于高程异常的精确求得。 　　2多项式函数拟合法的数学模型 　　多项式函数拟合法的基本思想：在小区域GPS网内，将似大地水准面看成曲面或平面，将高程异常F表示为平面坐标（x，y）的函数，通过网中起算点（既进行了GPS测量又进行了几何水准联测的点）已知的高程异常确定测区的似大地水准面形状，求出其余各点的高程异常，然后根据式（1）求出其他点的正常高，其数学模型为： 　　ζ=f（x，y）+ε..........（2） 　　式中，f（x，y）是拟合的似大地水准面；E是拟合误差，而： 　　f（x，y）=a0+a1x+a2y+a3x2+a4xy+a5y2+...（3） 　　式中，a0、a1、a2、a3、a4、a5，为拟合待定参数；x，y为各GPS点的平面坐标。 　　取式（3）中的一、二次项，合并（2）、（3）式后即得二次曲面拟合模型： 　　ζ=[a0 a1 a2 a3 a4 a5]×[1 x y x2 xy y2]T+ε。。。。。。。。。。（4） 　　取式（3）中的一次项，合并（2）、（3）式后即得二次平面拟合模型： 　　ζ=[a0 a1 a2]×[1 x y]T+ε。。。。。。。。。。（5） 　　每个起算点可组成一个式（4）或式（5），在[E2=min]条件下，解算出ai即可求出网中其余点的高程异常，并利用式（1）求出各待定点的正常高h。 　　3高程异常的估算 　　3.1利用地球重力场模型计算高程异常 　　大地水准面可以近似定义为通过平均海水面的重力等位面，并可用重力场模型来表示，地球重力场模型通常是指地球扰动位的球函数级数展开的系数集。这些位系数一般由下列数据组合确定：卫星轨道分析（用于确定大地水准面长波特性）和地面重力（中、短波长特性）以及卫星测高得到的反映区域中、短波长特性的值。 　　根据给定的位系数，用下式计算各点的高程异常，并按模型所具有的精度逼近高程异常场。 　　两测站间的高程异常差可表示为： 　　Δζ=ζj+1-ζj..........（7） 　　由于全球重力场模型在局部地区的不适应性，近几年来发展了地区性的密合重力场模型。它是利用区域重力数据拟合修正全球重力场模型的高阶位系数，得到与本区域更为吻合的模型，一般称之为/密合模型0。目前，密合模型已逐渐受到人们关注，在一些区域重力场研究中得到了应用，如西德和欧洲似大地水准面确定中采用的IFE88E2（360阶）在欧洲的60个GPS/水准点上的比较精度为±0.32 m。类似的研究工作，在国内也已进行，如由中科院测地所研制的中国IGG93C（360阶），总参测绘研究所研制的DQM94A（360阶）等模型，都属于符合中国地区的密合模型，模拟中国区域似大地水准面精度得到了较大地提高，其中IGG93C（360阶）模型模拟中国东部地区似大地水准面精度为±0.74m，DQM94A（360阶）模型模拟中国区域的似大地水准面精度为±0.88m。利用地球重力场模型计算高程异常具有较强的整体概念，应用方便，适应性广，但受模型精度影响较大。随着地面重力资料的不断增加以及模型解算方法的改进，地球重力场模型的精度将会不断提高，其表示的（似）大地水准面精度也将得到