CN110596737B - 一种gnss虚拟参考站自适应建站方法 - Google Patents

一种gnss虚拟参考站自适应建站方法 Download PDF

Info

Publication number
CN110596737B
CN110596737B CN201910786930.7A CN201910786930A CN110596737B CN 110596737 B CN110596737 B CN 110596737B CN 201910786930 A CN201910786930 A CN 201910786930A CN 110596737 B CN110596737 B CN 110596737B
Authority
CN
China
Prior art keywords
station
virtual reference
cors
rover
calculating
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.)
Active
Application number
CN201910786930.7A
Other languages
English (en)
Other versions
CN110596737A (zh
Inventor
冯威
黄丁发
程宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201910786930.7A priority Critical patent/CN110596737B/zh
Publication of CN110596737A publication Critical patent/CN110596737A/zh
Application granted granted Critical
Publication of CN110596737B publication Critical patent/CN110596737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain

Abstract

本发明公开了一种GNSS虚拟参考站自适应建站方法,设置一个满足流动站定位精度要求的精度指标;在CORS网络中,选取CORS站;标定各CORS站所在位置达到要求精度的所需流动站与虚拟参考站的最大距离;根据各CORS站所在位置对应的值,建立在整个CORS网络中上述最大距离的分布模型;根据流动站的位置,结合计算流动站所在位置的值,为流动站提供或建立虚拟参考站,该方法可以自适应化实现多个流动站同时共用一个虚拟参考站,在保障流动站定位精度的前提下,生成尽量少数的虚拟参考站,从而降低了数据处理中心的计算负担。方法能够根据流动站位置的分布信息,自动调整虚拟参考站的位置,实现CORS网络位置服务中虚拟参考站位置的优化确定。

Description

一种GNSS虚拟参考站自适应建站方法
技术领域
本发明涉及全球卫星导航技术领域,尤其涉及一种GNSS虚拟参考站自适应建站方法。
背景技术
全球卫星导航系统(GNSS)包括美国的GPS、中国的BDS、俄罗斯的GLONASS、欧洲的Galileo,以及相关的增强系统,如美国的WAAS(广域增强系统)、欧洲的EGNOS(欧洲静地导航重叠系统)和日本的MSAS(多功能运输卫星增强系统)等,还涵盖在建和以后要建设的其他卫星导航系统。早期的定位是基于伪距观测值的伪距单点定位,但因为各种误差的影响,其精度只能在10m量级。为提高定位精度,随之出现了差分技术,目前最广泛的是实时动态载波相位差分技术(RTK),精度能够达到厘米级(实时),但作用距离有限,且由于只有一个参考站提供差分信息,一旦该参考站出现故障或产生粗差,流动站往往不能正常定位,降低了定位的可靠性,鉴于此,NRTK应运而生。虚拟参考站(VRS)是目前使用最成熟和最广泛的NRTK技术,数据控制中心根据主参考站的观测值和各项误差改正数生成虚拟参考站处的虚拟观测值,流动站和虚拟参考站组成短基线进行差分解算。
传统的做法是:流动站向数据处理中心发送自己的概略位置信息,中心就在其概略坐标处建立一个VRS,并将VRS处的虚拟观测值播发给流动站,也就说一个流动站就要生成一个虚拟参考站,且为双向通信,该方法在流动站数量较少时并无不妥,但随着流动站数量的增加,将不堪负重,无法满足多流动站同时工作的需求。
在传统做法上提出了格网化VRS,采用某一参考站的格网间距经验值作为整网划分的标准,也就说整只采用一个格网间距,以偏概全,这样就可能导致有些测区的格网间距较大,达不到精度要求,而有些测区的格网过于密集。
发明内容
本发明的目的就在于为了解决上述问题而提供一种GNSS虚拟参考站自适应建站方法。
本发明通过以下技术方案来实现上述目的:
本发明包括以下步骤:
a设置一个满足流动站定位精度要求的精度指标σu
b选取CORS网中s(s>6)个CORS站并进行编号;
c将步骤b所述的CORS站作为流动站并在其附近建立虚拟参考站,计算该CORS站所在位置达到a要求精度的所需流动站与虚拟参考站的最大距离l;
d根据步骤c中s个CORS站所在位置对应的l值,利用二次曲面拟合数学模型建立l在整个CORS网中的分布模型;
e根据步骤d计算流动站所在位置的l值,为流动站提供或建立虚拟参考站。
进一步,所述步骤c具体包括以下步骤:
c1在以步骤b所述的某个CORS站为中心,半径为l0km的圆上建立4个虚拟参考站,组成正方形;
c2把步骤c1所述的某个CORS站作为流动站,采用实测数据,分析其固定解的定位精度;该流动分别和4个虚拟参考站的虚拟观测值进行差分解算,计算4组该流动站每个历元的三维坐标误差;
c3基于步骤c2得到4组三维坐标误差序列,计算每组坐标分量的中误差和点位中误差;
c4计算步骤c1所述的CORS站所在位置的l值;重复步骤c1、c2、c3、c4,计算步骤b所述的所有CORS站所在位置的l值。
进一步,步骤d还包括以下步骤:
d1确定二次曲面拟合数学模型为:
Figure BDA0002178351460000031
B,L分别表示任意点处的纬度、经度坐标,ai(i=1,2,3,4,5,6)为二次曲面拟合所确定的系数,
Figure BDA0002178351460000032
表示任意点处的l值;
d2根据最小二乘原理计算步骤d1中所述系数ai
d3基于步骤d1、d2得到l在整个CORS网中的分布模型。
进一步,所述步骤e具体包括以下步骤:
e1流动站根据步骤d3得到的l在整个CORS网中的分布模型,计算其所在位置的l值;
e2若在以流动站位置为中心,以l为半径的圆内,或者2l为边长的正方形范围内存在虚拟参考站,则将上述虚拟参考站中最近的提供给流动站;否则,根据流动站位置,新建一个虚拟参考站,并将其提供给流动站。
本发明的有益效果在于:
本发明是一种GNSS虚拟参考站自适应建站方法,与现有技术相比,本发明提出了GNSS虚拟参考站自适应建站方法。该方法综合考虑了电离层延迟和对流层延迟等空间相关误差的影响。可以自适应化实现多个流动站同时共用一个虚拟参考站,在保障流动站定位精度的前提下,生成尽量少数的虚拟参考站,从而降低了数据处理中心的计算负担。该方法能够根据流动站位置的分布信息,自动调整虚拟参考站的位置,实现CORS网络位置服务中虚拟参考站位置的优化确定,本发明适用于跨度较大、地势起伏较大地区的CORS网络的虚拟参考站的自适应建站。
附图说明
图1是本发明CORS站编号图;
图2是本发明CORS站建立图。
具体实施方式
下面结合附图对本发明作进一步说明:
设置一个满足流动站定位精度要求的精度指标σu
选取CORS网中9个CORS站并进行编号(填充的三角形),如图1所示;
将上述所述的CORS站作为流动站并在其附近建立虚拟参考站,计算该CORS站所在位置达到上述要求精度的所需流动站与虚拟参考站的最大距离l。
1)在以CORS站1为中心,半径为l0km的圆上建立4个虚拟参考站,组成正方形,如图2所示。
2)把CORS站1作为流动站,采用实测数据,分析其固定解的定位精度。流动站1分别和4个虚拟参考站的虚拟观测值进行差分解算,计算4组流动站1每个历元的三维坐标误差:
Figure BDA0002178351460000041
其中,j=1,2,3,4;
Figure BDA0002178351460000042
为流动站1和第j个虚拟参考站的虚拟观测值差分解算得到的每一历元的流动站在对应坐标系下的三维直角坐标;
(x1,y1,z1)T为CORS站1在流动站对应坐标系下的三维直角坐标的精确值;
Figure BDA0002178351460000043
为每个历元的三维坐标误差。
3)基于步骤2)得到4组三维坐标误差序列,计算每组坐标分量的中误差和点位中误差:
Figure BDA0002178351460000051
Figure BDA0002178351460000052
n1j为第j组观测历元个数,
Figure BDA0002178351460000057
分别为第j组三维坐标误差计算的X,Y,Z方向上的中误差,M1j为第j组点位中误差,取4个点位中误差的平均值:
Figure BDA0002178351460000053
4)满足流动站定位精度要求的精度指标σu,则CORS站1所在位置的l值为:
Figure BDA0002178351460000054
重复以上步骤2、3、4、5,便可得到其余8个CORS站的l值:l2,l3,l4,l5,l6,l7,l8,l9
根据上述内容中9个CORS站所在位置对应的l值,利用二次曲面拟合数学模型建立l在整个CORS网中的分布模型。
1)二次曲面拟合数学模型为:
Figure BDA0002178351460000055
B,L分别表示任意点处的纬度、经度坐标,ai(i=1,2,3,4,5,6)为二次曲面拟合所确定的系数,
Figure BDA0002178351460000056
表示任意点处的l值。
2)根据最小二乘原理计算系数ai
A)列出误差方程:
Figure BDA0002178351460000061
B)表示为总误差方程:
V=CX-L  (3-1)
V=[v1 v2 … v9]T
Figure BDA0002178351460000062
X=[a1 a2 a3 a4 a5 a6]T
L=[l1 l2 … l9]T
其中,vj表示lj的误差,Bj,Lj分别表示CORS站j所在位置的纬度、经度坐标,(j=1,2…9)。
通过最小二乘得(3-1)的解为:
X=(CTC)-1CTL
从而得到l在整个CORS网中的分布模型:
Figure BDA0002178351460000063
根据上述计算流动站所在位置的l值,按以下流程为流动站提供或建立虚拟参考站。
1)流动站根据2.4计算其所在位置的l值:
Figure BDA0002178351460000064
其中,Bu,Lu分别表示流动站所在位置的纬度、经度坐标,
Figure BDA0002178351460000065
良示流动站所在位置的l值。
2)若在以流动站位置为中心,以
Figure BDA0002178351460000066
为半径的圆内,或者以2
Figure BDA0002178351460000067
为边长的正方形范围内存在虚拟参考站,则将上述虚拟参考站中最近的提供给流动站;否则,根据流动站位置,新建一个虚拟参考站,并将其提供给流动站。
结论
本发明提供了一种GNSS虚拟参考站自适应建站方法:(1)设置一个满足流动站定位精度要求的精度指标;(2)在CORS网络中,选取s(s>6)个CORS站;(3)标定各CORS站所在位置达到(1)要求精度的所需流动站与虚拟参考站的最大距离l;(4)根据(3)中各CORS站所在位置对应的l值,建立l在整个CORS网络中的分布模型;(5)根据流动站的位置,结合(4)计算流动站所在位置的l值,按照以下过程为流动站提供或建立虚拟参考站:若在以流动站位置为中心,以l为半径的圆,或者以2l为边长的正方形范围内存在虚拟参考站,则将上述虚拟参考站中最近的提供给流动站;否则,根据流动站位置,新建一个虚拟参考站,并将其提供给流动站。GNSS虚拟参考站自适应建站方法综合考虑了电离层延迟和对流层延迟等空间相关误差的影响。该方法可以自适应化实现多个流动站同时共用一个虚拟参考站,在保障流动站定位精度的前提下,生成尽量少数的虚拟参考站,从而降低了数据处理中心的计算负担。该方法能够根据流动站位置的分布信息,自动调整虚拟参考站的位置,实现CORS网络位置服务中虚拟参考站位置的优化确定,本发明适用于跨度较大、地势起伏较大地区的CORS网络的虚拟参考站的自适应建站。
以上显示和描述了本发明的基本原理和主要特征及本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (1)

1.一种GNSS虚拟参考站自适应建站方法,其特征在于,包括以下步骤:
a设置一个满足流动站定位精度要求的精度指标σu
b选取CORS网中s个CORS站并进行编号;
c将步骤b所述的CORS站作为流动站并在其附近建立虚拟参考站,计算该CORS站所在位置达到a要求精度的所需流动站与虚拟参考站的最大距离l;
c1在以步骤b所述的某个CORS站为中心,半径为l0km的圆上建立4个虚拟参考站,组成正方形;
c2把步骤c1所述的某个CORS站作为流动站,采用实测数据,分析其固定解的定位精度;该流动分别和4个虚拟参考站的虚拟观测值进行差分解算,计算4组该流动站每个历元的三维坐标误差;
c3基于步骤c2得到4组三维坐标误差序列,计算每组坐标分量的中误差和点位中误差;
c4计算步骤c1的l值;重复步骤c1、c2、c3、c4,计算s个CORS站的l值;
d根据步骤c中s个的l值,利用二次曲面拟合数学模型建立l在整个CORS网中的分布模型;
d1确定二次曲面拟合数学模型为:
Figure FDA0004048391120000011
B,L分别表示任意点处的纬度、经度坐标,ai为二次曲面拟合所确定的系数,其中i=1,2,3,4,5,6,
Figure FDA0004048391120000012
表示任意点处的l值;
d2根据最小二乘原理计算步骤d1中所述系数ai
d3基于步骤d1、d2得到l在整个CORS网中的分布模型;
e根据步骤d计算流动站所在位置的l值,为流动站提供或建立虚拟参考站;
e1流动站根据步骤d3得到的l在整个CORS网中的分布模型,计算其所在位置的l值;
e2若在以流动站位置为中心,以l为半径的圆内,或者2l为边长的正方形范围内存在虚拟参考站,则将上述虚拟参考站中最近的提供给流动站;否则,根据流动站位置,新建一个虚拟参考站,并将其提供给流动站。
CN201910786930.7A 2019-08-24 2019-08-24 一种gnss虚拟参考站自适应建站方法 Active CN110596737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910786930.7A CN110596737B (zh) 2019-08-24 2019-08-24 一种gnss虚拟参考站自适应建站方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910786930.7A CN110596737B (zh) 2019-08-24 2019-08-24 一种gnss虚拟参考站自适应建站方法

Publications (2)

Publication Number Publication Date
CN110596737A CN110596737A (zh) 2019-12-20
CN110596737B true CN110596737B (zh) 2023-04-14

Family

ID=68855400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910786930.7A Active CN110596737B (zh) 2019-08-24 2019-08-24 一种gnss虚拟参考站自适应建站方法

Country Status (1)

Country Link
CN (1) CN110596737B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111954226B (zh) * 2020-07-30 2021-04-20 宁波冶金勘察设计研究股份有限公司 基于Delaunay异构CORS系统的基准站三角形的构网方法
CN114666804B (zh) * 2022-03-28 2023-06-23 北京四维图新科技股份有限公司 一种基于不同环境场景选取基站架设坐标的方法、装置及设备

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101943749A (zh) * 2010-09-10 2011-01-12 东南大学 基于星型结构的虚拟参考站网络rtk定位方法
CN104729486A (zh) * 2015-03-16 2015-06-24 东南大学 一种基于似大地水准面精化的无验潮水深测量方法
WO2015191854A1 (en) * 2014-06-13 2015-12-17 Trimble Navigation Limited Global navigation satellite system receiver system with radio frequency hardware component
CN105738934A (zh) * 2016-02-06 2016-07-06 武汉大学 附加大气信息动态约束的urtk模糊度快速固定方法
CN106772483A (zh) * 2016-12-20 2017-05-31 广州吉欧电子科技有限公司 一种基于cors系统的数据后处理方法及装置
CN107678050A (zh) * 2017-09-05 2018-02-09 西南交通大学 基于粒子滤波的glonass相位频间偏差实时追踪和精密估计方法
CN109143297A (zh) * 2018-09-12 2019-01-04 东南大学 一种兼容似大地水准面模型的实时网络rtk定位方法
CN109669193A (zh) * 2019-02-28 2019-04-23 成都北斗易寻科技有限公司 Gnss虚拟参考站手动设站方法及系统
CN109901204A (zh) * 2019-03-27 2019-06-18 北京航空航天大学 一种基于伪距误差分布模型的gbas完好性性能评估方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102498415B (zh) * 2009-09-19 2014-04-16 天宝导航有限公司 具有漫游器模糊度固定的gnss信号处理

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101943749A (zh) * 2010-09-10 2011-01-12 东南大学 基于星型结构的虚拟参考站网络rtk定位方法
WO2015191854A1 (en) * 2014-06-13 2015-12-17 Trimble Navigation Limited Global navigation satellite system receiver system with radio frequency hardware component
CN104729486A (zh) * 2015-03-16 2015-06-24 东南大学 一种基于似大地水准面精化的无验潮水深测量方法
CN105738934A (zh) * 2016-02-06 2016-07-06 武汉大学 附加大气信息动态约束的urtk模糊度快速固定方法
CN106772483A (zh) * 2016-12-20 2017-05-31 广州吉欧电子科技有限公司 一种基于cors系统的数据后处理方法及装置
CN107678050A (zh) * 2017-09-05 2018-02-09 西南交通大学 基于粒子滤波的glonass相位频间偏差实时追踪和精密估计方法
CN109143297A (zh) * 2018-09-12 2019-01-04 东南大学 一种兼容似大地水准面模型的实时网络rtk定位方法
CN109669193A (zh) * 2019-02-28 2019-04-23 成都北斗易寻科技有限公司 Gnss虚拟参考站手动设站方法及系统
CN109901204A (zh) * 2019-03-27 2019-06-18 北京航空航天大学 一种基于伪距误差分布模型的gbas完好性性能评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Virtual Satellite Positioning Technology Based on Differential Pressure Altimetry and PDR";Zhongliang Deng 等;《2018 Ubiquitous Positioning, Indoor Navigation and Location-Based Services (UPINLBS)》;20181231;1-7 *
"附有空间误差梯度信息的虚拟基准站服务技术";张熙 等;《西南交通大学学报》;20180914;1005-1011 *

Also Published As

Publication number Publication date
CN110596737A (zh) 2019-12-20

Similar Documents

Publication Publication Date Title
CN110261879B (zh) 广域地基增强位置服务的格网化虚拟参考站方法
CN106970404B (zh) 一种基于Delaunay三角网的多冗余网络RTK大气误差内插方法
US10690775B2 (en) Crowdsourcing atmospheric correction data
CN107390233B (zh) 一种低轨卫星导航增强电离层延迟改正参数方法
CN106569239B (zh) 一种广播式网络rtk定位技术
Matosevic et al. A comparison of accuracy using a GPS and a low-cost DGPS
US8072373B2 (en) Method for accuracy estimation of network based corrections for a satellite-aided positioning system
CN108490473A (zh) 一种融合gnss和uwb的无人机增强定位方法与系统
CN110297259B (zh) 基于格网化的基准站网定位增强信息可用性监测方法及系统
WO2017128871A1 (zh) 一种高精度实时卫星定位装置及其方法
CN104965207A (zh) 一种区域对流层天顶延迟的获取方法
CN105044747B (zh) 一种基于多星共视和滤波的时间同步装置及其方法
El-Mowafy Precise real-time positioning using Network RTK
CN112230252B (zh) 终端定位方法、装置、计算机设备和存储介质
CN110146908B (zh) 一种虚拟参考站观测数据的生成方法
CN110596737B (zh) 一种gnss虚拟参考站自适应建站方法
CN110426723A (zh) 一种卫星定位gga数据的获取与地图发布的方法
CN111314849B (zh) 定位方法、装置、移动终端及存储介质
CN111290005A (zh) 载波相位的差分定位方法、装置、电子设备及存储介质
CN110412634B (zh) 基于参考站网的伪距差分定位方法及装置
CN101140322A (zh) 用于布局卫星导航地基区域增强系统的装置及方法
WO2022046317A2 (en) System and method for providing gnss corrections
CN111983641A (zh) 一种用于实时生成北斗星基增强系统完好性参数的方法
Krypiak-Gregorczyk et al. Validation of approximation techniques for local total electron content mapping
Przestrzelski et al. Performance of real-time network code DGPS services of ASG-EUPOS in north-eastern Poland

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
GR01 Patent grant
GR01 Patent grant