CN114355404A - 一种自适应网格点vrs改正数生成方法 - Google Patents
一种自适应网格点vrs改正数生成方法 Download PDFInfo
- Publication number
- CN114355404A CN114355404A CN202110288014.8A CN202110288014A CN114355404A CN 114355404 A CN114355404 A CN 114355404A CN 202110288014 A CN202110288014 A CN 202110288014A CN 114355404 A CN114355404 A CN 114355404A
- Authority
- CN
- China
- Prior art keywords
- reference station
- grid
- vrs
- correction number
- generating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种自适应网格点VRS改正数生成方法,所述方法包括以下步骤:通过实时区域大气延迟模型和大气误差阈值计算大气格网点坐标,再根据大气格网点坐标内插生成单差电离层延迟改正数和单差对流层延迟改正数;根据事先设定的几何格网点坐标,依次生成各格网点的几何距离改正数;根据单差电离层延迟改正数、单差对流层延迟改正数以及几何距离改正数生成各格网点的VRS改正数。以解决现有技术定位的可靠性低、定位效果差和高动态用户在会面临频繁初始化的问题。
Description
技术领域
本发明涉及姿态估计技术领域,尤其涉及一种自适应网格点VRS改正数生成方法。
背景技术
现有的格网化VRS技术存在以下问题:
首先,虽然格网化VRS技术有效破解了用户数量对网络RTK技术约束,从而使得GNSS高精度定位服务进一步扩大应用领域成为可能。但是,也需要注意到格网化VRS技术本身的问题,即格网划分标准全部依赖于操作人员的经验判断或实测结果验证。考虑到大气环境在不同时间、空间环境下的多变性和复杂性,基于经验设定的格网划分标准实际上增加了定位效果的不可控性,降低了定位的可靠性。
其次,大气格网点与当前大气延迟直接相关,又由于大气延迟随时间不断变化,故大气格网点的密度和点位也在不断变化。结合上述理论,VRS站点变更会导致用户需要重新收敛固定模糊度,频繁的变更站点会严重降低用户的定位效果。
最后,现有VRS改正数生成方法中,电离层延迟改正数、对流层延迟改正数、轨道误差改正数以及几何距离改正数均使用同一组虚拟基准站坐标,这导致现有VRS技术在面对高动态用户时,受限于自身技术方案,需要频繁更换虚拟基准站位置,以保证与用户之间距离,因此高动态用户在使用网络RTK服务时也面临虚拟基准站不断迁站导致的频繁初始化。
发明内容
为了解决以上现有技术的缺点和不足之处,本发明的目的是提供一种自适应网格点VRS改正数生成方法。
本发明的技术方案是:一种自适应网格点VRS改正数生成方法,所述方法包括以下步骤:
通过CORS网各基线实时大气信息基于线性内插或反距离内插生成实时大气误差模型,并基于实时大气误差模型和大气误差阈值计算大气格网点坐标,再根据大气格网点坐标内插生成单差电离层延迟改正数和单差对流层延迟改正数;
根据事先设定的几何格网点坐标,依次生成各格网点的几何距离改正数;
根据单差电离层延迟改正数、单差对流层延迟改正数以及几何距离改正数生成各格网点的VRS改正数。
进一步地,所述大气格网点坐标的计算方法包括以下步骤:
S01、计算大气格网步长;
S02、划分大气格网。
进一步地,所述大气格网步长的计算方法如下:
S01-1、计算本区域电离层延迟和对流层延迟在经线和纬线上的拟合因子;
S01-2、通过拟合因子和大气误差阈值计算格网在经线和纬线上的划分步长。
进一步地,所述拟合因子的计算方法为:
Coru=a1·ΔXu+a2·ΔYu;
式中,Err表示各基线的大气误差,ΔX表示基线在经线方向的映射长度,ΔY表示,Coru表示基线在纬线方向的映射长度,a1表示本区域电离层延迟和对流层延迟在经线上的内插系数,a2表示本区域电离层延迟和对流层延迟在纬线上的内插系数。
进一步地,所述格网在经线和纬线上的划分步长的计算方法为:
式中Δx表示格网在经线上的划分步长,Δy表示格网在经线上的划分步长。
进一步地,所述VRS改正数的生成方法如下:
式中,下标VRS代表虚拟基准站,下标MAS代表主基准站,其中主基准站选择离用户概略坐标最近的GNSS基准站,下标1代表卫星L1载波信号,下标2代表卫星L2载波信号;λ为载波信号波长,为载波相位观测值,ΔN为单差模糊度,P为伪距观测值,为虚拟基准站和主基准站关于卫星i的单差星地距改正数;为虚拟基准站和主基准站关于卫星i的单差L1电离层延迟改正数;为虚拟基准站和主基准站关于卫星i的单差L2电离层延迟改正数;为虚拟基准站和主基准站关于卫星i的单差对流层延迟改正数;为主基准站关于卫星i的残余误差改正数;为虚拟基准站和主基准站关于卫星i的单差轨道误差改正数。
进一步地,几何格网点的密度为每个基准站一个几何格网点。
本发明的有益效果是:与现有技术相比,
1)本发明在格网化VRS技术的基础上,提出了基于星-站大气延迟信息的自适应格网化VRS改正数的生成方法,该方法通过对CORS网的实时大气延迟模型进行精度分析,依据大气环境的变化情况自主调节格网划分步长,使得用户在CORS服务范围内的任意地点与改正数之间的大气延迟误差小于规定值,从而确保定位服务的精度可控和可靠性;
2)本发明对几何格网和大气格网进行分离,大气格网随大气信息实时变更坐标位置以保证用户的定位精度;
3)传统模式移动距离过长会重新初始化,而本发明几何格网点位置固定,保证用户了测量的连贯性,可将几何格网设置较为稀疏,通过两套格网的组合,相较于传统VRS技术和第一代的格网化VRS改正数技术需要频繁初始化,本发明长距离移动无需初始化,可以实时保障用户定位精度,同时提高对高动态用户的服务质量。
附图说明
图1为本发明实施实例1的流程图。
具体实施方式
下面结合附图及具体的实施例对发明进行进一步介绍:
实施实例1:参考图1,一种自适应网格点VRS改正数生成方法,所述方法包括以下步骤:
通过CORS网各基线实时大气信息基于线性内插或反距离内插生成实时大气误差模型,并基于实时大气误差模型和大气误差阈值计算大气格网点坐标,再根据大气格网点坐标内插生成单差电离层延迟改正数和单差对流层延迟改正数;根据事先设定的几何格网点坐标,依次生成各格网点的几何距离改正数;根据单差电离层延迟改正数、单差对流层延迟改正数以及几何距离改正数生成各格网点的VRS改正数。
用户在使用VRS改正数进行差分定位时,假设对流层延迟改正数、电离层延迟改正数、轨道误差改正数这些与距离相关的误差改正数被完全消除,则观测方差可以简化为:
式中,下标VRS代表虚拟基准站,下标rover代表流动站,VRS-rover代表虚拟站到流动站,下标1代表卫星L1载波信号,下标2代表卫星L2载波信号,λ为载波信号波长,为载波相位观测值,为双差模糊度,双差载波相位观测值,为双差伪距观测值,上标i和j表示卫星编号。
从上式可知,几何距离改正数和其他距离相关误差改正数之间是相互独立的。因此可以对几何距离改正数和其他误差改正数分开处理。又因为电离层延迟改正数、对流层延迟改正数、轨道误差改正数等由于其均为距离相关误差改正数,所以需要使用与用户较近的坐标生成改正数。但几何距离改正数由于其和观测值中的任何误差均不相关,其不影响用户最终的定位精度,因此几何距离改正数实际可以为任意值。传统VRS改正数中使用与大气延迟改正数相同的坐标计算几何改正数更多是希望提升VRS观测值的真实性,对定位结果并没有影响。但这一计算方法却制约了高动态用户的定位效果。在本发明中,对计算几何改正数和大气误差改正数使用的坐标进行区分,大气误差改正数使用大气格网点坐标,几何改正数所需坐标则可摆脱与用户间距离的限制。考虑到高动态用户频繁初始化这一问题,需要大幅降低用户在运动过程中VRS改正数中几何信息的变更频率,故计算几何改正数的几何格网密度要尽可能降低。但考虑到VRS观测值是在主站观测值的基础上生成的,所以最低需要为每个基准站设置一个几何格网点。由于几何格网点位置不影响用户的定位结果,故具体设置数量和位置可自行决定。
进一步地,所述大气格网点坐标的计算方法包括以下步骤:
S01、计算大气格网步长;
S02、划分大气格网。
进一步地,所述大气格网步长的计算方法如下:
S01-1、计算本区域电离层延迟和对流层延迟在经线和纬线上的拟合因子;
S01-2、通过拟合因子和大气误差阈值计算格网在经线和纬线上的划分步长。
进一步地,所述拟合因子的计算方法为:
Coru=a1·ΔXu+a2·ΔYu;
式中,Err表示各基线的大气误差,ΔX表示基线在经线方向的映射长度,ΔY表示,Coru表示基线在纬线方向的映射长度,a1表示本区域电离层延迟和对流层延迟在经线上的内插系数,a2表示本区域电离层延迟和对流层延迟在纬线上的内插系数。
进一步地,所述格网在经线和纬线上的划分步长的计算方法为:
式中Δx表示格网在经线上的划分步长,Δy表示格网在经线上的划分步长。
进一步地,所述VRS改正数的生成方法如下:
式中,下标VRS代表虚拟基准站,下标MAS代表主基准站,其中主基准站选择离用户概略坐标最近的GNSS基准站,下标1代表卫星L1载波信号,下标2代表卫星L2载波信号;λ为载波信号波长,为载波相位观测值,ΔN为单差模糊度,P为伪距观测值,为虚拟基准站和主基准站关于卫星i的单差星地距改正数;为虚拟基准站和主基准站关于卫星i的单差L1电离层延迟改正数;为虚拟基准站和主基准站关于卫星i的单差L2电离层延迟改正数;为虚拟基准站和主基准站关于卫星i的单差对流层延迟改正数;为主基准站关于卫星i的残余误差改正数;为虚拟基准站和主基准站关于卫星i的单差轨道误差改正数。
进一步地,几何格网点的密度为每个基准站一个几何格网点
最后,各格网点以大气格网点坐标进行区分,在根据用户概略坐标选取格网点VRS改正数时,也根据用户坐标与大气格网点坐标最小原则选择对应格网点VRS,但需要注意的是,播发改正数时,播发的VRS观测值坐标应为几何格网点坐标。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (7)
1.一种自适应网格点VRS改正数生成方法,其特征在于,所述方法包括以下步骤:
通过CORS网各基线实时大气信息基于线性内插或反距离内插生成实时大气误差模型,并基于实时大气误差模型和大气误差阈值计算大气格网点坐标,再根据大气格网点坐标内插生成单差电离层延迟改正数和单差对流层延迟改正数;
根据事先设定的几何格网点坐标,依次生成各格网点的几何距离改正数;
根据单差电离层延迟改正数、单差对流层延迟改正数以及几何距离改正数生成各格网点的VRS改正数。
2.根据权利要求1所述的自适应网格点VRS改正数生成方法,其特征在于,所述大气格网点坐标的计算方法包括以下步骤:
S01、计算大气格网步长;
S02、划分大气格网。
3.根据权利要求2所述的自适应网格点VRS改正数生成方法,其特征在于,所述大气格网步长的计算方法如下:
S01-1、计算本区域电离层延迟和对流层延迟在经线和纬线上的拟合因子;
S01-2、通过拟合因子和大气误差阈值计算格网在经线和纬线上的划分步长。
6.根据权利要求1所述的自适应网格点VRS改正数生成方法,其特征在于,所述VRS改正数的生成方法如下:
7.根据权利要求1所述的自适应网格点VRS改正数生成方法,其特征在于,几何格网点的密度为每个基准站一个几何格网点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110288014.8A CN114355404A (zh) | 2021-03-18 | 2021-03-18 | 一种自适应网格点vrs改正数生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110288014.8A CN114355404A (zh) | 2021-03-18 | 2021-03-18 | 一种自适应网格点vrs改正数生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114355404A true CN114355404A (zh) | 2022-04-15 |
Family
ID=81096007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110288014.8A Pending CN114355404A (zh) | 2021-03-18 | 2021-03-18 | 一种自适应网格点vrs改正数生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114355404A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114935768A (zh) * | 2022-07-13 | 2022-08-23 | 武汉大学 | 一种基于单基站的虚拟参考站的构建方法 |
CN116027357A (zh) * | 2023-03-30 | 2023-04-28 | 武汉大学 | 基于大气与地形的格网改正方法、装置、设备及存储介质 |
CN116125495A (zh) * | 2022-12-14 | 2023-05-16 | 北京六分科技有限公司 | 电离层改正数的确定方法、设备、存储介质及程序产品 |
CN116243351A (zh) * | 2023-05-11 | 2023-06-09 | 北京讯腾智慧科技股份有限公司 | 一种高精度gnss定位终端的图层化行为数据脱敏方法 |
CN116626725A (zh) * | 2023-07-24 | 2023-08-22 | 武汉大学 | 基于用户分布格网改正生成方法、装置、设备及存储介质 |
-
2021
- 2021-03-18 CN CN202110288014.8A patent/CN114355404A/zh active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114935768A (zh) * | 2022-07-13 | 2022-08-23 | 武汉大学 | 一种基于单基站的虚拟参考站的构建方法 |
CN114935768B (zh) * | 2022-07-13 | 2022-11-04 | 武汉大学 | 一种基于单基站的虚拟参考站的构建方法 |
CN116125495A (zh) * | 2022-12-14 | 2023-05-16 | 北京六分科技有限公司 | 电离层改正数的确定方法、设备、存储介质及程序产品 |
CN116125495B (zh) * | 2022-12-14 | 2024-04-16 | 北京六分科技有限公司 | 电离层改正数的确定方法、设备、存储介质及程序产品 |
CN116027357A (zh) * | 2023-03-30 | 2023-04-28 | 武汉大学 | 基于大气与地形的格网改正方法、装置、设备及存储介质 |
CN116243351A (zh) * | 2023-05-11 | 2023-06-09 | 北京讯腾智慧科技股份有限公司 | 一种高精度gnss定位终端的图层化行为数据脱敏方法 |
CN116626725A (zh) * | 2023-07-24 | 2023-08-22 | 武汉大学 | 基于用户分布格网改正生成方法、装置、设备及存储介质 |
CN116626725B (zh) * | 2023-07-24 | 2023-10-10 | 武汉大学 | 基于用户分布格网改正生成方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114355404A (zh) | 一种自适应网格点vrs改正数生成方法 | |
KR102547522B1 (ko) | 쾌속 정밀 측위 방법 및 시스템 | |
EP2191290B1 (en) | Method for accuracy estimation of network based corrections for a satellite-aided positioning system | |
CN108828626B (zh) | 基于实时网格的网络rtk电离层延迟内插方法及*** | |
CN110045407B (zh) | 一种分布式伪卫星/gnss优化定位方法 | |
KR20210006954A (ko) | 정밀 단독 측위 및 실시간 키네마틱 조합(ppp-rtk)의 측위 방법 및 설비 | |
CN112230252B (zh) | 终端定位方法、装置、计算机设备和存储介质 | |
CN105807300B (zh) | 一种用北斗双频接收机进行高精度动态单点定位的方法 | |
CN110261879B (zh) | 广域地基增强位置服务的格网化虚拟参考站方法 | |
CN108196284B (zh) | 一种进行星间单差模糊度固定的gnss网数据处理方法 | |
CN102269819B (zh) | 一种用于精密单点定位方法的估计方法 | |
WO2017070732A1 (en) | A method of analysing a signal transmitted between a global satellite navigation satellite system and a receiver | |
CN104680008A (zh) | 一种基于多参考站的网络rtk区域大气误差建模方法 | |
CN111983641B (zh) | 一种用于实时生成北斗星基增强***完好性参数的方法 | |
CN107966722B (zh) | 一种gnss钟差解算方法 | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN111694030A (zh) | 一种基于格网虚拟观测值的bds局域差分方法及*** | |
CN108345017B (zh) | 新的网络rtk大气内插方法 | |
CN116027357B (zh) | 基于大气与地形的格网改正方法、装置、设备及存储介质 | |
CN111381264A (zh) | 网络rtk中长基线模糊度固定方法和平台 | |
CN112902825A (zh) | 一种适用于高精度形变监测的北斗/gnss网络rtk算法 | |
Banville et al. | Improving real-time kinematic PPP with instantaneous cycle-slip correction | |
CN112596077A (zh) | 一种针对终端载体为低轨卫星的卫星导航信号仿真方法 | |
CN116626730A (zh) | 一种顾及nwp的海上区域cors增强ppp方法 | |
Charoenkalunyuta et al. | Performance of network-based RTK GPS in low-latitude region: A case study in Thailand |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |