发明内容
本发明针对目前VRS(虚拟参考站技术)的网络RTK中在高程差异较大地区对流层误差估计模型的缺陷,提出了一种GNSS网络差分定位***中的误差改正方法,它是一种顾及高程差异的三维内插方法,在线性组合模型的基础上,考虑了高程的影响,采用三维线性组合内插方法来提高网络RTK差分改正信息的准确性和可靠性。
本发明采用的技术方案如下:
一种GNSS网络差分定位***中的误差改正方法,差分定位***中,各个GNSS网络参考站连续采集观测数据,实时传输到数据处理与控制中心的数据库,进行网络计算,控制中心在线解算GPS参考站网内各独立基线的载波相位整周模糊度值,数据处理中心利用参考站网载波相位观测值计算每条基线上的双差综合误差,并据此建立距离相关误差(包括载波伪距误差、对流层延迟误差、电离层延迟误差、轨道误差及多路径误差)的空间参数模型;移动站将通过单点定位得到的NMEA格式的概略坐标发送给控制中心,控制中心在该坐标位置创建一个虚拟参考站VRS;控制中心根据参考站、移动站及GPS卫星的相对几何关系,通过内插计算模型得到移动站与参考站间的综合误差改正数,再根据虚拟观测值计算模型生成VRS处的虚拟观测值,控制中心把虚拟观测值作为网络差分改正信息发送给移动站;移动站接收网络差分信息与VRS构成短基线,通过常规RTK计算模型进行差分解算,确定移动站位置;其特征在于:通过内插计算模型得到移动站与参考站间的综合误差改正数时,考虑了高程的影响,采用三维线性组合内插方法来提高网络RTK差分改正信息的准确性和可靠性,利用参考站已解算的包括电离层、对流层误差的空间相关误差,根据移动站与周围参考站的空间位置关系,内插出移动站实时有效的差分改正信息,包括以下步骤:
(1)在CORS网络的三角形参考站网形中,A、B、C为参考站,其中以A站为主参考站,B,C站为辅参考站,u为移动站,构建虚拟位置VRS处虚拟双差观测值Φ:
同时系数ai满足下列约束条件,即
式中,
是移动站u与参考站i之间的双差观测值,X
i是各参考站位置及移动站位置的位置矢量,(2)式可整理得:
a2(XB-XA)+a3(XC-XA)=Xu-XA (3)
对Xi位置矢量采用三维平面坐标x、y、h进行求解,其中的x、y、h为:将WGS-84大地坐标转换为地方平面坐标,则各站点的平面坐标和高程值x、y、h可得到,对于(3)式可以转变成如下式的矩阵运算:
式中,X和Y为平面坐标,H是高程,则系数为:
α=(ATA)-1ATL (5)
其中,α2,α3可以作为相关系数内插出与距离相关的误差,而α1为主参考站相关系数;
由于轨道误差、电离层误差、对流层误差以及多路径误差影响通过观测值线性组合已经大削弱,因此组合双差观测方程可以简化为:
又由于关系式
同理,式(6)可化为:
把(7)式代入(8)中,可推出下式:
定义残差项
则移动站处u的综合误差改正数Lu,A为:
Lu,A=a2VB,A+a3VC,A
而对于多于3个参考站的CROS参考站网络,构成虚拟双差观测值:
式中,n为参考站的数量,r1,r2,r3,…rn为各个参考站,假设以r1为主参考站,u为移动站;
同时系数满足下列约束条件,即
式中,αi为组合的相关系数,Xu、Xri分别为移动站和第ri个参考站的平面坐标,n为参考站的数量,构成矩阵运算:
双差观测方程可转化为:
定义残差项
则移动站处u的综合误差改正数Lr1,u为:
Lr1,u=a2Vr1,r2+a3Vr1,r3+…+anVr1,rn (16)
综合误差改正数Lr1,u用于移动站位置虚拟观测值的生成,提高了虚拟位置虚拟观测值的精度,最终实现GNSS网络差分的高精度定位效果。
本发明可用于可应用于以下领域或行业:
(1)各级国土、测绘、勘查等部门建立GNSS连续运行卫星基准站网络用于实时定位获取位置信息。
(2)数字地面模型建立、数字地图综合等地学应用领域。
本发明的优点及有益效果:
(1)本发明提出的三维线性组合方法能准确修正高程方向偏差,精确估计移动站位置的对流层改正。相比较传统内插方法20-40cm的改正精度,本发明可提高改正精度至5-6cm,始终维持在厘米级的精度,且仅需3个参考站即可完成高精度内插,可靠性更高,适用于GPS/VRS参考站网络厘米级实时动态定位要求。
(2)本发明可提高省级和行业级乃至全国范围的大规模的网络RTK***中地形起伏较大地区的CORS网络建设,实现山区或高原等地形起伏大的地区的高精度内插改正数生成。
具体实施方式
参看图1,对于三个参考站的网络RTK三角网,A、B、C为参考站,其中以A为主参考站,B、C为辅参考站,u为移动站。
1)计算各参考站及移动站位置的三维位置矢量。
参考站及移动站位置的三维位置矢量计算方法:将WGS-84大地坐标以基准站当地子午线作为中央子午线进行高斯投影得到WGS-84平面坐标,再将WGS-84平面坐标通过7参数转换或4参数转换得到地方平面坐标,则各站点的平面坐标和高程值x、y、h可得到。
2)在CORS网络的三角形参考站网形中,构建虚拟位置VRS处虚拟双差观测值Φ。
同时系数ai满足下列约束条件,即
式中,
是移动站u与参考站i之间的双差观测值,X
i是各参考站位置及移动站位置的三维位置矢量。
转变成如下式的矩阵运算:
式中,X和Y为平面坐标,H是高程,则系数为:
α=(ATA)-1ATL
其中,α2,α3可以作为相关系数内插出与距离相关的误差,而α1为主参考站相关系数。
3)定义残差项
则移动站处u的综合误差改正数Lu,A为:
Lu,A=a2VB,A+a3VC,A
对于多于3个参考站的CROS参考站网络,首先构成虚拟双差观测值:
式中,n为参考站的数量,r1,r2,r3,…rn为各个参考站,假设以r1为主参考站,u为移动站,
同时系数满足下列约束条件,即
式中,αi为组合的相关系数,Xu、Xri分别为移动站和第ri个参考站的平面坐标,n为参考站的数量,构成矩阵运算:
同时,双差观测方程可转化为:
定义残差项
则移动站处u的综合误差改正数Lr1,u为:
Lr1,u=a2Vr1,r2+a3Vr1,r3+…+anVr1,rn
下面使用重庆CORS试验参考站作为算例,参看图2:
重庆CORS试验参考站网络包括长寿(CHSH)、丰都(FEDU)、南川(NACH)、彭水(PESH)和武隆(WULO)共5个连续运行参考站,试验采用2010年11月10日12:00至14:00的相位观测数据,采样间隔为15秒,试验网络平均基线长为88.325km,如图1-3所示。选用南川(NACH)作为主参考站,武隆(WULO)作为移动站,对本专利提出的误差改正方法的精度、实用性、可靠性等进行实地测试;
第一步,计算各参考站及移动站的三维位置矢量,即平面坐标X和Y及高程H,由于重庆参考站周边地区的地形起伏比较大,使得移动站与参考站之间的最大高程差异达到1182m,参看图3所示。
第二步,分别用常规内插方法及本发明三维线性组合内插方法对对流层延迟误差改正精度进行比较。
试验中,图3中的四条基线的双差对流层误差,是利用在参考站网络初始化之后,构造双差对流层延迟滤波器,实时计算参考站网络基线的对流层双差改正数得到真值。分别采用常规网络内插方法和三维线性组合方法对对流层延迟误差改正精度的进行比较。试验选取的参考卫星为PRN10,辅助卫星为PRN27、PRN26,这两颗卫星在2010年11月10日12:00至14:00这段时间内高度角的变化趋势分别为17.0度上升至55.3度、72.1度下降至13.2度,高度角的变化比较具有代表意义。其改正结果比较如图4a、图4b所示。
其中四种常规内插方法为:LIM为线性内插方法,LCM为线性组合方法,DIM为距离线性内插方法,LSM为低阶曲面方法。RAW为基线双差对流层延迟误差真值,NEW为本专利提出的三维线性组合内插方法。计算这五种内插方法的中误差,基线NACH-WULO对流层改正中误差如图5。
由图5可知,在重庆CORS的试验中,高程起伏的大小对对流层误差有很大的影响。传统的内插模型改正精度为20-40cm;对于三维线性组合模型的改正精度则为5-6cm,始终维持在厘米级的精度,使对流层改正精度从分米级提高到厘米级。综上所述,三维线性组合模型对高程差异比较大的区域改正精度比较明显,具有很高的优越性,对流层误差得到了很好消除,因此三维线性组合模型对流层改正模型精度大大高于常规模型结果,可以在实际中进行广泛的应用。