具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图详予说明。
请参阅图1,一种基于北斗地基增强的网络RTK解算方法,包括:
接收北斗基准站发送的北斗卫星数据,所述北斗卫星数据包括北斗卫星星历数据和北斗卫星观测数据;
根据所述北斗卫星星历数据和北斗卫星观测数据,计算整周模糊度,并固定整周模糊度;
根据所述整周模糊度,计算北斗基准站上各基线的电离层延迟误差和对流层延迟误差,并分别根据所述电离层延迟误差和对流层延迟误差,构建电离层误差改正模型以及对流层误差改正模型;
建立虚拟参考站,并根据所述电离层误差改正模型以及对流层误差改正模型,获取所述虚拟参考站各可视卫星的误差改正信息。
从上述描述可知,本发明的有益效果在于:可全天候实时向用户终端提供厘米级导航定位。
进一步地,所述建立虚拟参考站,并根据所述电离层误差改正模型以及对流层误差改正模型,获取所述虚拟参考站各可视卫星的误差改正信息之后,进一步包括:
将所述误差改正信息实时播发至用户终端。
由上述描述可知,便于用户终端可根据误差改正信息解算自身精确位置。
进一步地,所述根据所述北斗卫星星历数据和北斗卫星观测数据,计算整周模糊度,并固定整周模糊度具体为:
根据所述北斗卫星星历数据,计算卫星坐标,并根据所述卫星坐标和所述北斗基准站的坐标,计算卫星和北斗基准站之间的距离,得到星站距离;
根据所述星站距离以及所述北斗卫星观测数据,构建单差观测方程,并根据所述单差观测方程,构建双差观测方程;
通过基于三差的探测周跳法,进行周跳探测;
根据所述双差观测方程,通过预设的宽巷组合,计算整周模糊度,并固定整周模糊度。
由上述描述可知,通过建立单差双差观测方程,可消除接收机钟差和卫星钟差对观测值的影响;采用基于三差的探测周跳法,不仅消除了卫星钟差和接收机钟差,同时还消除了整周模糊度。
进一步地,所述单差观测方程为λΔφ
a,b=Δρ
a,b-ΔI
a,b+ΔT
a,b+λΔN
a,b+Δt
a,b+ε,所述双差观测方程为
其中,Δ为单差算子,
为双差算子,a、b表示两个不同的北斗基准站,i、j表示两颗不同的北斗卫星,φ为载波相位观测量,λ为载波相位波长,ρ为卫星到北斗基准站的几何距离,I为电离层延迟误差,T为对流层延迟误差,N为整周模糊度,t为接收机钟差,ε为观测噪声。
进一步地,所述宽巷组合的公式为:
其中,
为双差电离层无关宽巷组合,L
1、L
2为载波观测量,
为双差载波,P
1、P
2为伪距观测量,f
1、f
2分别为L
1、L
2载波的频率。
由上述描述可知,上述宽巷组合可消除电离层延迟、卫星钟差、接收机钟差以及卫星至接收机之间的几何距离,仅受测量噪声和多路径误差的影响,并且这些影响可以通过多历元平滑减弱或消除。通过上述宽巷组合进行整周模糊度固定,可提高整周模糊度固定解的正确率和可靠性。
进一步地,所述根据所述整周模糊度,计算北斗基准站上各基线的电离层延迟误差和对流层延迟误差,并分别根据所述电离层延迟误差和对流层延迟误差,构建电离层误差改正模型以及对流层误差改正模型具体为:
根据所述整周模糊度和双差观测方程,计算北斗基准站上各基线的电离层延迟误差和对流层延迟误差;
根据所述的电离层延迟误差,建立电离层误差改正模型,所述电离层误差改正模型为vIa,b i=α1·Δxa,b+α2·Δya,b+α3,其中,vIa,b i为整周模糊度固定后单差天顶电离层延迟误差,Δxa,b和Δya,b分别为北斗基准站a与北斗基准站b之间的纬度差和经度差,α1、α2和α3为模型系数矢量;
根据所述对流层延迟误差,建立对流层误差改正模型,所述对流层误差改正模型为ZTWa,b(hb)=ZTWa,b-δZTWa,b,其中,ZTWa,b为整周模糊度固定后单差天顶对流层延迟湿分量,hb为北斗基准站b的高程,ZTWa,b(hb)为高程水平上的单差天顶对流层延迟湿分量的距离相关部分,δZTWa,b为通过预设的先验模型计算的对流层延迟高程方向相对偏差量。
进一步地,所述建立虚拟参考站,并根据所述电离层误差改正模型以及对流层误差改正模型,获取所述虚拟参考站各可视卫星的误差改正信息具体为:
接收用户终端发送的概略坐标,并根据所述概略坐标生成虚拟参考站;
选取距离所述虚拟参考站最近的北斗基准站作,作为参考基站;
根据所述虚拟参考站的坐标以及所述电离层误差改正模型和对流层误差改正模型,对所述参考基站的观测数据进行修正,生成虚拟参考站的观测数据,作为所述虚拟参考站各可视卫星的误差改正信息。
本发明还提出了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如上所述的方法。
实施例一
请参照图2,本发明的实施例一为:一种基于北斗地基增强的网络RTK解算方法,可应用于北斗卫星导航系统,适用于中长距离网络RTK解算,可以实现厘米级高精度定位结果。如图2所示,包括如下步骤:
S1:北斗解算中心接收北斗基准站发送的北斗卫星数据,所述北斗卫星数据包括北斗卫星星历数据和北斗卫星观测数据。
各个北斗基准站连续采集观测数据,并实时传输到北斗解算中心。北斗基准站发送的数据主要包括北斗卫星数据包括北斗卫星星历数据和北斗卫星观测数据两类,其中,北斗卫星星历数据用于计算卫星位置等参数,北斗卫星观测数据用于建立基线中的数学模型。利用两类数据结果可计算基线整周模糊度。
在实际应用场景中,北斗基准站可实时自动地向北斗解算中心提供加密的差分改正数据,并支持差分数据传输、支持标准NTRIP协议传输,其他厂家仪器可以接入系统。系统长时间无人值守稳定运行,高并发处理网络数据。
S2:根据所述北斗卫星星历数据,计算卫星坐标,并根据所述卫星坐标和所述北斗基准站的坐标,计算卫星和北斗基准站之间的距离,得到星站距离。
即使用星历数据计算卫星轨道位置,得到卫星坐标,由于北斗基准站的精确坐标已知,因此根据卫星坐标和北斗基准站的坐标,即可计算得到卫星和北斗基准站之间的精确距离。
S3:根据所述星站距离以及所述北斗卫星观测数据,构建单差观测方程,并根据所述单差观测方程,构建双差观测方程。
通过建立单差双差观测方程,可消除接收机钟差和卫星钟差对观测值的影响。
具体地,单差模型即将两个不同的北斗基准站上同步观测相同卫星所得的观测量求差。使用站间单差模型可以消除卫星钟差的影响,其观测方程如下:
λΔφa,b=Δρa,b-ΔIa,b+ΔTa,b+λΔNa,b+Δta,b+ε
其中,Δ为单差算子,a、b表示两个不同的北斗基准站,φ为载波相位观测量,λ为载波相位波长,ρ为卫星到北斗基准站的几何距离,I为电离层延迟误差,T为对流层延迟误差,N为双差整周模糊度,t为接收机钟差,ε为观测噪声。
双差模型即在站间单差的基础上,同步观测两颗不同的卫星,即所得单差模型再次求差。使用双差模型可以消除接收机钟差的影响,其观测方程如下:
其中,
为双差算子,a、b表示两个不同的北斗基准站,i、j表示两颗不同的北斗卫星,φ为载波相位观测量,λ为载波相位波长,ρ为卫星到北斗基准站的几何距离,I为电离层延迟误差,T为对流层延迟误差,N为双差整周模糊度,ε为观测噪声。
由于北斗基准站精准坐标已知,即北斗卫星与北斗基准站的真实距离已知,若想要建立电离层与对流层模型,则必须求得整周模糊度N的值。
S4:通过基于三差的探测周跳法,进行周跳探测。
在数据处理的过程中,常常伴有北斗卫星信号失锁或中断的情况,此时载波相位观测值是不连续的,会引起相位观测值的整周数发生跳变,这种现象称为周跳。因此,在计算最终的整周模糊度结果之前,必须对其中可能存在的周跳进行探测。
本实施例中,采用基于三差的探测周跳法,不仅消除了卫星钟差和接收机钟差,同时还消除了整周模糊度。其原理是在相邻历元间进行三差,如果三差观测值的粗差超过阈值,则判定出现周跳,对于出现周跳的观测值,不参与后续解算,或者进行周跳修复。
S5:根据所述双差观测方程,通过预设的宽巷组合,计算整周模糊度,并固定整周模糊度。
周跳探测完成后,可以利用双差观测方程求解整周模糊度。
具体地,先进行双差宽巷整周模糊度的求解。本实施例中,使用几何电离层无关宽巷组合求解双差宽巷整周模糊度,该组合的公式为:
其中,
为双差电离层无关宽巷组合,L
1、L
2为载波观测量,
为双差载波,P
1、P
2为伪距观测量;f
1、f
2分别为L
1、L
2载波的频率。
使用P(伪距观测量)进行计算得到的基准站坐标误差大,使用L(载波)观测量计算得到的基准站坐标误差小,但载波观测量计算复杂,因为载波需要求解整周模糊度,而伪距观测量计算简单,不需要求解整周模糊度,因此利用两个电磁波的特性,进行组合,即为电离层无关宽巷组合。
使用该宽巷组合求解整周模糊度的原因在于该组合消除了电离层延迟、卫星钟差、接收机钟差以及卫星至接收机之间的几何距离,仅受测量噪声和多路径误差的影响,并且这些影响可以通过多历元平滑减弱或消除,之所以称为宽巷组合是因为该组合的模糊度保持了整数特性,且波长较长,达到43cm,几乎可以直接确定宽巷模糊度。
固定宽巷模糊度后,可以利用传统求解整周模糊度方法计算每个波段的整周模糊度。具体地,首先利用kalman滤波器将每个波段无电离层浮点解求解出来,再结合之前得到的宽巷模糊度固定解,二者相结合,利用标准LAMBDA算法,解算每个波段的整周模糊度固定解。
S6:根据所述整周模糊度,计算北斗基准站上各基线的电离层延迟误差和对流层延迟误差,并分别根据所述电离层延迟误差和对流层延迟误差,构建电离层误差改正模型以及对流层误差改正模型。
具体地,计算得到整周模糊度的固定解后,由步骤S3中的双差观测方程,可求得北斗基准站网上各基线的电离层延迟误差和对流层延迟误差,以此建立电离层误差改正模型以及对流层误差改正模型。
对于建立电离层误差改正模型,可以建立内插模型,采用低阶曲面拟合法可逼近电离层延迟误差的空间相关性,模型如下:
其中,VIa,b i为整周模糊度固定后单差天顶电离层延迟误差,Δxa,b和Δya,b分别为北斗基准站a与北斗基准站b之间的纬度差和经度差,α1、α2和α3为模型系数矢量,使用最小二乘法解算模型系数矢量。
对于建立对流层误差改正模型,将Saastamoinen模型作为先验模型计算高程偏差,模型如下:
ZTWa,b(hb)=ZTWa,b-δZTWa,b
δZTWa,b=ZTWa(ha)-ZTWa(hb)
其中,ZTWa,b为整周模糊度固定后单差天顶对流层延迟湿分量,hb为北斗基准站b的高程,ZTWa,b(hb)为主参考站(即北斗基准站b)高程水平上的单差天顶对流层延迟湿分量的距离相关部分,δZTWa,b为对流层延迟高程方向相对偏差量,ZTWa(ha)为北斗基准站a高程水平上的天顶对流层延迟湿分量的距离相关部分,ZTWa(hb)为北斗基准站b高程水平上的天顶对流层延迟湿分量的距离相关部分。
S7:建立虚拟参考站,并根据所述电离层误差改正模型以及对流层误差改正模型,获取所述虚拟参考站各可视卫星的误差改正信息。
具体地,用户终端利用单点定位生成概略坐标并传给北斗解算中心作为虚拟参考站(VRS),北斗解算中心选取距离虚拟参考站最近的北斗基准站作为参考基站,根据虚拟参考站坐标以及解算出的电离层误差改正模型和对流层误差改正模型,对参考基站的观测数据进行修正,生成虚拟参考站的观测数据,即得到虚拟参考站(VRS)当前历元所有可视卫星的误差改正信息。
S8:将所述误差改正信息实时播发至用户终端,用户终端根据误差改正信息解算自身精确位置。
进一步地,在实际应用场景中,对流层延迟干分量可用Saastamoinen模型和NMF映射函数进行修正。接着,利用电离层误差改正模型和对流层误差改正模型进行修正,将虚拟参考站(VRS)带入两个模型中进行伪距和载波相位观测量的改正,将差分改正值播发给用户。
本实施例中,北斗解算中心在线解算北斗基准站网内各条基线的载波相位整周模糊度值;在北斗基准站网络间的整周模糊度确定后,利用网内相位观测值反算出每条基线上各种误差源的实际或综合误差影响值,并依此建立电离层、对流层等距离相关误差的空间参数模型;用户终端使用单点定位确定的用户概略坐标,并通过数据网络传送给北斗解算中心,中心根据该位置创建虚拟参考站(VRS),计算用户、虚拟参考站(VRS)和卫星的相对几何关系,通过内插得到VRS上各误差源影响的改正值,并按照RTCM格式发给流动用户,流动用户根据改正值解算自身精确位置。
本实施例能够全天候实时向用户终端提供厘米级导航定位,相比传统RTK,提高了作业距离和定位精度,并使用几何电离层无关宽巷组合进行双差整周模糊度固定,提高了整周模糊度固定解的正确率和可靠性。
实施例二
本实施例是对应上述实施例的一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如上述的基于北斗地基增强的网络RTK解算方法实施例中的各个过程,且能达到相同的技术效果,为避免重复,这里不再赘述。
综上所述,本发明提供的一种基于北斗地基增强的网络RTK解算方法及存储介质,北斗解算中心在线解算北斗基准站网内各条基线的载波相位整周模糊度值;在北斗基准站网络间的整周模糊度确定后,利用网内相位观测值反算出每条基线上各种误差源的实际或综合误差影响值,并依此建立电离层、对流层等距离相关误差的空间参数模型;北斗解算中心根据用户终端发送的概略坐标创建虚拟参考站,计算虚拟参考站各可视卫星的误差改正值,用户终端根据改正值解算自身精确位置。通过建立单差双差观测方程,可消除接收机钟差和卫星钟差对观测值的影响;采用基于三差的探测周跳法,不仅消除了卫星钟差和接收机钟差,同时还消除了整周模糊度;通过电离层无关宽巷组合,可消除电离层延迟、卫星钟差、接收机钟差以及卫星至接收机之间的几何距离,仅受测量噪声和多路径误差的影响,并且这些影响可以通过多历元平滑减弱或消除,且可提高整周模糊度固定解的正确率和可靠性。
本发明可应用于中长距离网络RTK解算,可以实现厘米级高精度定位结果,能够全天候实时向用户终端提供厘米级导航定位,相比传统RTK,提高了的作业距离和定位精度。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换,或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。