CN111060943B - 一种符合实际的gnss观测值定权方法 - Google Patents
一种符合实际的gnss观测值定权方法 Download PDFInfo
- Publication number
- CN111060943B CN111060943B CN201911338241.6A CN201911338241A CN111060943B CN 111060943 B CN111060943 B CN 111060943B CN 201911338241 A CN201911338241 A CN 201911338241A CN 111060943 B CN111060943 B CN 111060943B
- Authority
- CN
- China
- Prior art keywords
- observation
- gnss
- satellite
- constructing
- value
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000008569 process Effects 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种符合实际的GNSS观测值定权方法,方法包括获取GNSS观测值;确定各观测值对应的卫星高度角,卫星方位角和载噪比信息;构建顾及卫星高度角,卫星方位角和载噪比的定权函数;将定权函数用于构建随机模型;实现GNSS导航定位。本发明顾及了传统观测值定权方法中忽略掉的一些实际环境因素引起的观测值质量变化,由于这些测站环境因素往往与卫星高度角,卫星方位角和载噪比密切相关,本发明充分考虑了这些因素,并建立了相应的观测值定权方法。该方法与传统定权方法相比,可以提高GNSS导航定位精度和可靠性。
Description
技术领域
本发明涉及卫星导航定位技术领域,具体涉及一种符合实际的GNSS观测值定权方法。
背景技术:
在全球导航卫星系统(Global Navigation Satellite System,GNSS)中,准确的评估观测值的精度,从而构建相应准确的定权函数,是实现高精度高可靠性GNSS导航定位的前提。
之前,通常认为GNSS观测值具有同方差性,即等权模型,认为所有的观测值精度相等,但是这种方法在大多观测环境条件下不适用。因此为了克服观测值同方差性方面的缺陷,往往采用高度角或载噪比定权函数来实现观测值异方差性。高度角定权函数模型认为低高度角的卫星受到更严重的大气影响,还会使接收机天线增益减少,从而导致信号质量下降,因此认为高度角越高的卫星权重越大。载噪比定权函数模型认为载噪比和GNSS观测值由同一个跟踪环路记录,因此载噪比能直接评价信号质量,即载噪比越大,数据质量越好。
事实上,载噪比定权函数模型对测站相关的误差较为敏感。具体来说,当载噪比值较低时往往意味着观测值受到了多路径甚至是非视距传播的影响,因此观测值精度降低。此外,当信号扭曲严重的时候,其观测值会偏离标称值。载噪比这一重要特点目前并未充分考虑,造成之前的载噪比定权函数模型仅仅是间接和部分的反应观测值精度。此外,不同的方位角,接收机收到的观测值信号会受到测站环境周围植物以及建筑物等的影响,观测值的精度也会随之受到影响。然而,目前的定权函数模型中并未完全考虑这一特性。
发明内容
本发明的目的在于提供一种符合实际的GNSS观测值定权方法,以解决现有技术中导致的GNSS导航定位误差大不精确的缺陷。
一种符合实际的GNSS观测值定权方法,所述方法包括如下步骤:
获取GNSS观测值;
确定各观测值对应的参数数据;
根据观测数据构建观测值精度估计公式;
根据观测值精度估计公式构建定权函数;
通过定权函数实现GNSS导航定位。
进一步的,所述参数数据包括卫星高度角,卫星方位角以及载噪比中的一种或多种。
进一步的,所述观测数据的确定方法包括如下步骤:
获取观测文件和星历文件;
根据观测文件得到某个时刻某个观测值对应的载噪比数据;
通过星历文件计算该时刻的卫星位置;
根据观测文件求出GNSS接收设备的位置范围;
根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角。
进一步的,所述GNSS观测值包括不同频率的相位观测值和伪距观测值。
进一步的,根据参数数据构建定权函数的方法包括如下步骤:
确定观测值精度,计算公式为:
式中,σ2为GNSS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|.|代表求绝对值;为观测值对应的方位角;/>为方位角/>所对应的由用户定义的惩罚函数。
进一步的,根据观测值精度估计公式构建定权函数的方法包括如下步骤:
将观测值精度估计公式用于构建定权函数,计算公式为:
式中,为第i个观测值的方差;Pi为第i个观测值的权。
进一步的,通过定权函数实现GNSS导航定位的方法包括如下步骤:
构建观测数据的函数模型;
通过定权函数构建观测数据的随机模型;
通过函数模型与随机模型对位置参数进行求解,实现GNSS导航定位。
本发明的优点在于:该种符合实际的GNSS观测值定权方法,
1、顾及了传统观测值定权方法中忽略掉的一些测站环境因素引起的观测值质量变化,该方法与传统定权方法相比,更符合实际情况;
2、通过充分考虑测站环境因素对卫星高度角,卫星方位角和载噪比的影响,建立了相应的观测值定权方法,该方案操作简便,可执行性强;
3、利用该定权函数构建的随机模型,由于更能准确可靠地反应观测值精度,因此可以提高GNSS导航定位精度和可靠性。
附图说明
图1为本发明的流程示意图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
如图1所示,一种符合实际的GNSS观测值定权方法,所述方法包括如下步骤:
步骤一:获取GNSS观测值:
获取不同频率的全球定位系统(GPS)相位观测值和伪距观测值;
步骤二:确定各观测值对应的参数数据,参数数据包括卫星高度角,卫星方位角以及载噪比:
所述确定各观测值对应的高度角,方位角和载噪比信息,其中首先获取观测文件和星历文件,根据观测文件得到某个时刻某个观测值对应的载噪比数据;接着通过星历文件计算该时刻的卫星位置,根据观测文件求出GPS接收设备的位置范围;最后根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角。
步骤三:根据观测数据构建观测值精度估计公式:
确定观测值精度,计算公式为:
式中,σ2为GPS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|.|代表求绝对值;为观测值对应的方位角;/>为方位角/>所对应的由用户定义的惩罚函数。此时,由于是开阔环境,因此该惩罚函数为
步骤四:根据观测值精度估计公式构建定权函数:
将观测值精度估计公式用于构建定权函数,计算公式为:
式中,为第i个观测值的方差;Pi为第i个观测值的权;
步骤五:通过定权函数实现GNSS导航定位:
首先,构建观测数据的函数模型;接着通过定权函数构建观测数据的随机模型;最后,对位置参数进行求解,实现GPS导航定位;
获取GNSS观测数据并进行粗差探测与剔除,数据中断与修复,周跳探测与修复等预处理,由观测值构建误差方程,如需要固定相位观测值的模糊度,需要求解其中相位观测值的模糊度浮点解并进行固定,再次构建误差方程,即最终形成观测数据的函数模型。
由本发明的定权方法构建定权函数,从而确定观测数据的随机模型。
通过一定的准则如最小二乘准则,联立上述函数模型与随机模型,最终对含位置在内等参数进行求解,实现GNSS导航定位。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
Claims (1)
1.一种符合实际的GNSS观测值定权方法,其特征在于,所述方法包括如下步骤:
获取GNSS观测值;
确定各观测值对应的参数数据;
根据参数数据构建观测值精度估计公式;
根据观测值精度估计公式构建定权函数;
通过定权函数实现GNSS导航定位;
其中,所述参数数据包括卫星高度角,卫星方位角以及载噪比;
所述参数数据的确定方法包括如下步骤:
获取观测文件和星历文件;
根据观测文件得到某个时刻某个观测值对应的载噪比数据;
通过星历文件计算该时刻的卫星位置;
根据观测文件求出GNSS接收设备的位置范围;
根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角;
所述GNSS观测值包括不同频率的相位观测值和伪距观测值;
根据参数数据构建观测值精度估计公式的方法包括如下步骤:
确定观测值精度,计算公式为:
式中,σ2为GNSS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|.|代表求绝对值;为观测值对应的方位角;/>为方位角/>所对应的由用户定义的惩罚函数;
根据观测值精度估计公式构建定权函数的方法包括如下步骤:
将观测值精度估计公式用于构建定权函数,计算公式为:
式中,为第i个观测值的方差;Pi为第i个观测值的权;
通过定权函数实现GNSS导航定位的方法包括如下步骤:
构建观测数据的函数模型;
通过定权函数构建观测数据的随机模型;
通过函数模型与随机模型对位置参数进行求解,实现GNSS导航定位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338241.6A CN111060943B (zh) | 2019-12-23 | 2019-12-23 | 一种符合实际的gnss观测值定权方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338241.6A CN111060943B (zh) | 2019-12-23 | 2019-12-23 | 一种符合实际的gnss观测值定权方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111060943A CN111060943A (zh) | 2020-04-24 |
CN111060943B true CN111060943B (zh) | 2023-12-01 |
Family
ID=70300797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911338241.6A Active CN111060943B (zh) | 2019-12-23 | 2019-12-23 | 一种符合实际的gnss观测值定权方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111060943B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112198535A (zh) * | 2020-12-07 | 2021-01-08 | 湖南联智科技股份有限公司 | 一种gnss监测设备及监测点数据质量分析方法 |
CN112799106B (zh) * | 2021-03-30 | 2021-06-29 | 腾讯科技(深圳)有限公司 | 卫星观测信息处理方法、移动终端及存储介质 |
CN113281796B (zh) * | 2021-07-23 | 2021-10-15 | 腾讯科技(深圳)有限公司 | 位置确定方法、速度确定方法、装置、设备和存储介质 |
CN114217328B (zh) * | 2022-02-21 | 2022-05-17 | 长沙金维信息技术有限公司 | 导航电文半周跳变检测方法 |
CN114485655B (zh) * | 2022-04-06 | 2022-07-29 | 深圳华大北斗科技股份有限公司 | Gnss/ins组合导航数据质量控制方法 |
CN117647830B (zh) * | 2024-01-29 | 2024-05-07 | 中国测绘科学研究院 | 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799525A (zh) * | 2009-07-10 | 2010-08-11 | 中国测绘科学研究院 | 全球导航卫星系统增强系统监测站的自主完备性监测方法 |
CN104102822A (zh) * | 2014-07-01 | 2014-10-15 | 同济大学 | 一种多频gnss观测值随机特性建模方法 |
CN106814382A (zh) * | 2017-01-11 | 2017-06-09 | 武汉大学 | 联合高度角和用户测距误差的gnss卫星观测值定权方法 |
CN110018507A (zh) * | 2019-05-08 | 2019-07-16 | 中国科学院国家授时中心 | 一种基于星座间作差的组合精密单点定位方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2260320B1 (en) * | 2008-03-20 | 2016-01-13 | Telespazio S.p.A. | Position estimation enhancement for a global navigation satellite system receiver |
-
2019
- 2019-12-23 CN CN201911338241.6A patent/CN111060943B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799525A (zh) * | 2009-07-10 | 2010-08-11 | 中国测绘科学研究院 | 全球导航卫星系统增强系统监测站的自主完备性监测方法 |
CN104102822A (zh) * | 2014-07-01 | 2014-10-15 | 同济大学 | 一种多频gnss观测值随机特性建模方法 |
CN106814382A (zh) * | 2017-01-11 | 2017-06-09 | 武汉大学 | 联合高度角和用户测距误差的gnss卫星观测值定权方法 |
CN110018507A (zh) * | 2019-05-08 | 2019-07-16 | 中国科学院国家授时中心 | 一种基于星座间作差的组合精密单点定位方法及系统 |
Non-Patent Citations (3)
Title |
---|
Zhetao Zhang 等.Site-Specific Unmodeled Error Mitigation for GNSS Positioning in Urban Environments Using a Real-Time Adaptive Weighting Model.《Remote Sens.》.2018,10(7),全文. * |
冯彪 等.一种改进的组合导航系统选星算法.测绘工程.2018,27(2),全文. * |
王敏 等.一种顾及GEO卫星多路径效应影响的BDS长距离相对定位方法.测绘科学技术学报.2016,33(2),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111060943A (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111060943B (zh) | 一种符合实际的gnss观测值定权方法 | |
US8412456B2 (en) | Loosely-coupled integration of global navigation satellite system and inertial navigation system: speed scale-factor and heading bias calibration | |
US8860609B2 (en) | Loosely-coupled integration of global navigation satellite system and inertial navigation system | |
US8717226B2 (en) | Method for processing signals of an airborne radar with correction of the error in the radar beam pointing angle and corresponding device | |
US20170026797A1 (en) | Systems and methods for using doppler measurements to estimate a position of a receiver | |
CN114488230B (zh) | 一种多普勒定位方法、装置、电子设备和存储介质 | |
CN112859133B (zh) | 一种基于雷达与北斗数据的船舶深度融合定位方法 | |
US20240159529A1 (en) | Systems and methods for extending the spatial coverage of a reference pressure network | |
CN117146830A (zh) | 一种自适应多信标航位推算和长基线的紧组合导航方法 | |
CN115616636A (zh) | 一种面向无人机的多星座低轨卫星融合定位方法及装置 | |
CN111060936B (zh) | 基于信噪比的bds/gnss多路径探测方法 | |
CN113267793A (zh) | 一种基于外部增强信息的gbas对流层参数生成方法 | |
Guan et al. | Signal quality analysis and quality check of BDS3 Precise Point Positioning in the Arctic Ocean | |
KR100264977B1 (ko) | 위치정보제공시스템의 사용자 위치 실시간 측정 방법 | |
US11714198B2 (en) | Single-epoch pseudo-range positioning under varying ionosphere delays | |
KR102677918B1 (ko) | 가변적인 전리층 지연 하에서의 단일-에포크 의사-거리 위치 확인 | |
Raghuvanshi et al. | Precise positioning of smartphones using a robust adaptive Kalman filter | |
Hwang et al. | The standard for the selection of the appropriate GPS in the outdoor environment & the analysis of the performance for the improvement of reception | |
US20240159914A1 (en) | Method for taking provided gnss-relevant route information into account in the gnss-based localization of vehicles | |
CN116299234A (zh) | 雷达标校误差的修正方法及系统 | |
Rinnan et al. | Including GNSS based heading in inertial aided GNSS DP reference system | |
Böder | HCU-HMSS: A Multi Sensor System in Hydrographic Applications | |
Elsobeiey et al. | Effect of second-order ionospheric delay on GPS orbit and precise point positioning | |
CN113589349A (zh) | 基于粒子滤波的gnss实时高精度单差测姿方法 | |
CN116907470A (zh) | 创建环境模型的方法 |
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 |