CN107544081A - 顾及电离层约束的rtk定位方法 - Google Patents
顾及电离层约束的rtk定位方法 Download PDFInfo
- Publication number
- CN107544081A CN107544081A CN201710721898.5A CN201710721898A CN107544081A CN 107544081 A CN107544081 A CN 107544081A CN 201710721898 A CN201710721898 A CN 201710721898A CN 107544081 A CN107544081 A CN 107544081A
- Authority
- CN
- China
- Prior art keywords
- ionospheric
- double
- difference
- rtk
- constraints
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 11
- 239000005433 ionosphere Substances 0.000 title claims description 10
- 230000004807 localization Effects 0.000 title 1
- 239000013598 vector Substances 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 239000005436 troposphere Substances 0.000 claims description 2
- 230000003321 amplification Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种顾及电离层约束的RTK定位方法,基于GNSS非组合的双差观测值,通过附加电离层先验约束、空间域和时间域的约束信息,将电离层延迟误差当作未知参数进行实时估计,不仅高效解决了RTK定位中电离层误差影响,尤其是长基线条件下的电离层残差影响;而且相比常规的无电离层组合观测降低了观测噪声,使得长基线条件下RTK的定位精度更高更稳定。
Description
技术领域
本发明涉及一种RTK定位方法。
背景技术
载波相位差分实时动态定位(RTK)是精密定位中最主要的数据处理技术之一,基于高精度的载波相位差分观测值即可以实时获取厘米级定位精度,因此已经得到深入研究和广泛应用,但其适用性和可拓展性仍存在一些缺陷。首先,最为显著的问题是RTK定位精度随着基线距离增加逐渐衰减,当基线距离较长(大于10km)时,双差观测值难以准确消除电离层残差,导致定位精度降低,并且收敛时间加长,同时模糊度固定准确率下降;此外,对于较长基线的RTK,常采用无电离层组合观测值,将放大观测噪声并影响模糊度的快速固定。顾及电离层约束信息的RTK定位方法,将电离层延迟误差附加约束信息并进行参数化估计,较好的解决了长基线RTK定位的实际问题,大大提升了RTK的作业范围和作业领域。
发明内容
为了克服现有技术在长基线条件下定位精度衰减,收敛时间加长,模糊度固定准确率下降,且双差无电离层组合观测噪声放大等不足,本发明提供一种顾及电离层约束信息的RTK定位方法,基于GNSS非组合的双差观测值,通过附加电离层先验约束、空间域和时间域的约束信息,将电离层延迟误差当作未知参数进行实时估计,不仅高效解决了RTK定位中电离层误差影响,尤其是长基线条件下的电离层残差影响;而且相比常规的无电离层组合观测降低了观测噪声,使得长基线条件下RTK的定位精度更高更稳定。
本发明解决其技术问题所采用的技术方案包括建立多系统RTK双差观测方程模型和双差电离层残差约束模型两步。
第一步,基于多GNSS系统非组合的双差观测值,考虑电离层和对流层残余误差的影响,建立双差的RTK观测模型,实时估计用户和基准站间的基线向量以及模糊度、电离层、对流层参数;
第二步,确定双差电离层残差约束模型,将双差的电离层残差当作未知参数进行实时估计,并附加先验约束,空间域和时间域的信息约束,消除双差电离层残差影响。
本发明的有益效果是:
第一,解决了中长基线RTK定位精度衰减问题。本发明将双差电离层残差进行参数估计,高效解决了中长距离RTK中电离层残差难以准备处理的问题,有效的抵抗了常规RTK随基线长度增加定位精度的衰减。
第二,解决了中长基线RTK定位模糊度快速准确固定问题。对于中长基线RTK,常基于无电离层组合观测来消除电离层误差影响,但是存在噪声放大现象,导致模糊度的收敛较慢且准确性下降。本发明仍然采用非组合观测值,并基于附加约束信息的参数估计思想高效消除电离层误差,避免了噪声放大影响,解决了中长基线RTK定位模糊度快速准确固定问题。
第三,拓展了RTK作用范围和作业领域。本发明通过高效解决电离层误差,既保障了长距离RTK定位精度,又保障了模糊度的快速准确固定,将大大增加RTK的有效作用范围,同时作业领域也将得到延伸。
附图说明
图1是顾及电离层约束信息的RTK定位方法实施框图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明基于GNSS非组合的双差观测值,通过附加电离层先验约束、空间域和时间域的约束信息,将电离层延迟误差当作未知参数进行实时估计,不仅高效解决了RTK定位中电离层误差影响,尤其是长基线条件下的电离层残差影响;而且相比常规的无电离层组合观测降低了观测噪声,使得长基线条件下RTK的定位精度更高更稳定。基于该模型,长距离RTK用户仍可快速准确固定模糊度,实时获取厘米级动态位置信息,大大提升了RTK的作业范围和作业领域。
本发明包括建立多系统RTK双差观测方程模型和双差电离层残差约束模型两步。
第一步,建立多系统RTK双差观测方程模型
基于多GNSS系统非组合的双差观测值,通过精细考虑电离层和对流层残余误差的影响,建立双差的RTK观测模型,可实时估计用户和基准站间的基线向量以及模糊度、电离层、对流层参数。
第二步,确定双差电离层残差约束模型
通常,在短基线条件下,电离层误差可以通过双差观测进行有效消除;在长基线条件下,电离层误差可以通过组无电离层组合方程进行消除,但是它会放大观测噪声,影响模糊度的快速准确固定。本发明中,将双差的电离层残差当作未知参数进行实时估计,并附加先验约束,空间域和时间域的信息约束,高效解决了双差电离层残差影响。
本发明的技术方案主要包括两个核心技术:
(1)多系统双差RTK观测方程模型
GNSS伪距和载波相位的双差观测值线性化后可以表达为如下形式:
其中,角标i和j分别代表基准站和用户站,p和q表示卫星,并且p代表参考卫星。ij表示站间差分,pq表示星间差分。ρ是站星间几何距离,P,Φ分别为伪距和载波相位观测值,λ为载波波长。dX、dY、dZ为需要求解的基线向量,l、m和n代表卫星和接收机间的单位矢量。Amb代表载波相位模糊度,δI和δT代表电离层和对流层延迟。t表示信号接收时间,f1和f2为频率,ε为观测噪声,P1、P2、L1、L2分别表示两个频率上的伪距和相位观测。另外,双差变量按照如下方式定义:
这里δ为标准差,MT(t)为对流层延迟湿分量投影系数。
(2)双差电离层残差约束模型
本发明中,通过附加电离层先验约束,空间域和时间域的信息约束,将双差的电离层残差当作未知参数进行实时估计。
①先验信息约束
对于电离层先验信息约束,可以采用常规的Klobuchar模型或是Bent模型。本发明中,利用IGMAS分析中心提供的实时电离层模型进行先验信息约束,其精度为3-6TECU。电离层先验信息约束模型可以用以下方程表示:
这里,为电离层延迟的投影系数,它可以由单层模型计算。代表垂直方向上的总电子含量。
基于先验电离层信息约束,双差观测电离层残差的虚拟观测方程可以写成如下形式:
为先验信息噪声εprior的方差,可以设定为10-20cm2。
②空间约束
电离层VTEC的空间特性可以表达为与经差、纬差有关的平面函数,如公式(17)所示。
这里,m0,n0为平面模型的阶数,因对单站而言,电离层穿刺点的分布较近,可以定为2阶。λ为穿刺点的经纬度;λ0为测站的经纬度;Eij为平面模型的系数。当历元观测到的卫星个数多于4颗时,可以基于方程(18)求解平面模型系数,除首历元外当历元观测卫星个数小于4颗时,采用上一历元的系数代替本历元的系数。
基于空间信息约束,双差观测电离层残差的虚拟观测方程可以写成如下形式:
为空间约束εspace的噪声方差,可以设定为5-10cm2。
③时间域约束
根据电离层变化的特性,VTEC随时间变化缓慢,可以采用随机过程进行描述。本发明将双差的电离层残差用一个随机游走过程表达,如式(20)。
这里,代表双差电离层残差的历元变化,为时间约束εtemp的方差,可以设定为1-2(cm/s)2。
本发明的实施例具体包括以下步骤:
第一步,实时数据输入与预处理。输入基准站和用户站原始观测数据,广播星历产品和电离层先验信息。并对原始数据进行预处理分析,得到干净的数据。
第二步,RTK解算。首先,组成双差观测方程和电离层约束的虚拟观测方程,进行最小二乘解算,得到浮点解,并进行模糊度的固定和质量控制,得到最终的固定解结果。
第三步,结果输出。将RTK定位结果和精度信息通过通信手段实时输出给用户。
Claims (1)
1.一种顾及电离层约束的RTK定位方法,其特征在于包括下述步骤:
第一步,基于多GNSS系统非组合的双差观测值,考虑电离层和对流层残余误差的影响,建立双差的RTK观测模型,实时估计用户和基准站间的基线向量以及模糊度、电离层、对流层参数;
第二步,确定双差电离层残差约束模型,将双差的电离层残差当作未知参数进行实时估计,并附加先验约束,空间域和时间域的信息约束,消除双差电离层残差影响。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710721898.5A CN107544081B (zh) | 2017-08-22 | 2017-08-22 | 顾及电离层约束的rtk定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710721898.5A CN107544081B (zh) | 2017-08-22 | 2017-08-22 | 顾及电离层约束的rtk定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107544081A true CN107544081A (zh) | 2018-01-05 |
CN107544081B CN107544081B (zh) | 2022-12-02 |
Family
ID=60958480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710721898.5A Active CN107544081B (zh) | 2017-08-22 | 2017-08-22 | 顾及电离层约束的rtk定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107544081B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919305A (zh) * | 2018-08-07 | 2018-11-30 | 北斗导航位置服务(北京)有限公司 | 交通运输中北斗地基增强带状服务方法及系统 |
CN110007328A (zh) * | 2019-05-10 | 2019-07-12 | 国网浙江省电力有限公司信息通信分公司 | 基于北斗三号卫星四频观测值的非组合rtk定位方法 |
CN110346816A (zh) * | 2019-05-24 | 2019-10-18 | 武汉大学 | 短到长基线场景下电离层建模约束的实时差分定位方法 |
CN111123315A (zh) * | 2018-11-01 | 2020-05-08 | 千寻位置网络有限公司 | 非差非组合ppp模型的优化方法及装置、定位系统 |
CN111323748A (zh) * | 2018-12-13 | 2020-06-23 | 千寻位置网络有限公司 | 差分定位方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105204048A (zh) * | 2014-06-24 | 2015-12-30 | 东南大学 | 一种北斗兼容gps/glonass的网络rtk参考站间模糊度快速解算方法 |
CN105629279A (zh) * | 2015-12-18 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | 一种网络基准站间的宽巷模糊度固定方法 |
WO2016185500A1 (en) * | 2015-05-19 | 2016-11-24 | Spacearth Technology S.R.L. | Method for forecasting ionosphere total electron content and/or scintillation parameters |
CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
-
2017
- 2017-08-22 CN CN201710721898.5A patent/CN107544081B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105204048A (zh) * | 2014-06-24 | 2015-12-30 | 东南大学 | 一种北斗兼容gps/glonass的网络rtk参考站间模糊度快速解算方法 |
WO2016185500A1 (en) * | 2015-05-19 | 2016-11-24 | Spacearth Technology S.R.L. | Method for forecasting ionosphere total electron content and/or scintillation parameters |
CN105629279A (zh) * | 2015-12-18 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | 一种网络基准站间的宽巷模糊度固定方法 |
CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
Non-Patent Citations (2)
Title |
---|
ADENIYI IYIADE: "Ionospheric Effects on GPS Signals in Real Time Kinematic(RTK) Applications", 《2006 IEEE》 * |
王胜利等: "北斗系统高轨卫星高精度定位模型病态性分析及改进方法", 《导航定位学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919305A (zh) * | 2018-08-07 | 2018-11-30 | 北斗导航位置服务(北京)有限公司 | 交通运输中北斗地基增强带状服务方法及系统 |
CN111123315A (zh) * | 2018-11-01 | 2020-05-08 | 千寻位置网络有限公司 | 非差非组合ppp模型的优化方法及装置、定位系统 |
CN111323748A (zh) * | 2018-12-13 | 2020-06-23 | 千寻位置网络有限公司 | 差分定位方法及系统 |
CN111323748B (zh) * | 2018-12-13 | 2023-07-21 | 千寻位置网络有限公司 | 差分定位方法及系统 |
CN110007328A (zh) * | 2019-05-10 | 2019-07-12 | 国网浙江省电力有限公司信息通信分公司 | 基于北斗三号卫星四频观测值的非组合rtk定位方法 |
CN110346816A (zh) * | 2019-05-24 | 2019-10-18 | 武汉大学 | 短到长基线场景下电离层建模约束的实时差分定位方法 |
CN110346816B (zh) * | 2019-05-24 | 2023-08-15 | 武汉大学 | 短到长基线场景下电离层建模约束的实时差分定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107544081B (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544081B (zh) | 顾及电离层约束的rtk定位方法 | |
CN109196381B (zh) | 通过辅助数据对精确位置的快速确定 | |
CN106646538B (zh) | 一种基于单差滤波的变形监测gnss信号多路径改正方法 | |
CN107561568A (zh) | 基于统一模型的北斗非差非组合ppp‑rtk定位方法 | |
CN106855632B (zh) | 一种广播式vrs定位方法及系统 | |
CN102230971B (zh) | Gps多天线测姿方法 | |
CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
CN114236587A (zh) | 基于北斗地基增强的网络rtk解算方法及存储介质 | |
CN104656108B (zh) | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 | |
CN105629279B (zh) | 一种网络基准站间的宽巷模糊度固定方法 | |
CN110780323B (zh) | 一种长距离下基于北斗三频信号的实时分米级定位方法 | |
KR102180302B1 (ko) | 의사거리 측정치 및 이중차분 반송파 측정치를 이용한 전리층 보정정보 생성 시스템 및 방법 | |
CN103605145B (zh) | 基于gnss多频数据和cors实现网络实时动态定位的方法 | |
Zhang et al. | Models, methods and assessment of four-frequency carrier ambiguity resolution for BeiDou-3 observations | |
CN109597105B (zh) | 一种顾及载波系统间偏差的gps/glonass紧组合定位方法 | |
CN111352137B (zh) | 一种顾及广播星历误差的多模gnss异步rtk定位方法 | |
WO2017070732A1 (en) | A method of analysing a signal transmitted between a global satellite navigation satellite system and a receiver | |
CN104503223B (zh) | Gnss 三频高精度卫星钟差估计与服务方法 | |
Rabbou et al. | Precise point positioning using multi-constellation GNSS observations for kinematic applications | |
CN104808230A (zh) | 一种星钟与r-轨道误差混合估计的实时精密单点定位方法 | |
CN103592653B (zh) | 用于局部地区单频卫星导航用户的电离层延迟修正方法 | |
Li et al. | Regional modeling of tropospheric delay considering vertically and horizontally separation of station for regional augmented PPP | |
CN110346816B (zh) | 短到长基线场景下电离层建模约束的实时差分定位方法 | |
CN108254774A (zh) | 基于gnss多频信号的单基站长距离实时定位方法 | |
CN115144878A (zh) | 基于ppp的短距离大高差nrtk对流层延迟改正方法 |
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 |