【技术实现步骤摘要】
区域高程异常数据的多维粗差定位及定值方法
本专利技术涉及大地测量学领域,特别是涉及一种区域高程异常数据的多维粗差定位及定值方法。
技术介绍
区域为地球表面的一部分,高程异常是地球椭球面与地球似大地水准面之间的差异。确定区域高程异常是大地测量学及工程任务的重点工作之一。目前,确定区域高程异常的主要计算方法有以下两种:(1)重力似大地水准面的计算。按照莫洛金斯基理论求定高程异常,在计算过程中一般采用两次移去-恢复技术方法。但精度较低,难以满足工程需要。(2)利用GPS水准计算高程异常。该方法为在同一个点上获得该点的GPS大地高及精密水准正常高,两者的差值即为该点的高程异常,其精度比第一种方法要高。而且,如果在一个区域布设有足够多的GPS水准点,那么通过计算这些点的高程异常值,便可以建立该区域高程异常计算模型。方法(2)是目前确定区域(省市级)高程异常的常用方法。高程异常模型的精度,除了取决于GPS测定大地高的精度和水准测量正常高的精度外,还取决于高程异常值的数据质量,即高程异常值是否含有粗差。如果不能有效地对粗差进行定位并定值(用于对粗差改正),将会对参数的估计结果产 ...
【技术保护点】
1.区域高程异常数据的多维粗差定位及定值方法,其特征在于:包括以下步骤:S1:获取n个GPS点的高程异常值建立如下模型:
【技术特征摘要】
1.区域高程异常数据的多维粗差定位及定值方法,其特征在于:包括以下步骤:S1:获取n个GPS点的高程异常值建立如下模型:式(1)中,Li为第i个GPS点的高程异常值,Li=HGi-H0i,HGi为第i个GPS点的大地高、H0i为第i个GPS点的正常高,i=1,…,n;Xj,j=0,1,…,5,为二次曲面模型的待定参数;x为GPS点在高斯平面的纵坐标,y为GPS点在高斯平面的横坐标;将式(1)写成矩阵形式:L=AX+Δ(2)式(2)中,为高程异常值向量;为待定参数向量;为未知参数的系数矩阵,为随机误差向量;S2:根据式(3)得到X的最小二乘估计值根据式(4)得到单位权方差的最小二乘估计值其中,N=ATPA,P为权阵;V为随机误差向量的最小二乘估计,根据式(5)得到;式(5)中,H=A(ATPA)-1ATP;S3:根据式(6)得到ΔT(i):ΔT(i)=T(i)-Ta(6)式(6)中,T(i)为第i个GPS点的影响函数统计值,Ta为T(1),…,T(n)的算术平均值,根据式(7)得到,Ta根据式(8)得到;式(7)中,为删掉第i组数据后从剩余的n-1组数据回归X的最小二乘估计;通过以下方法判定强影响点:当ΔT(i)>0时,则第i点为强影响点;令所有强影响点个数之和为m;S4:步骤S4具体包括以下步骤:S4.1:根据式(9)计算检验统计量F:式(9)中,根据式(10...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。