CN111352137B - 一种顾及广播星历误差的多模gnss异步rtk定位方法 - Google Patents
一种顾及广播星历误差的多模gnss异步rtk定位方法 Download PDFInfo
- Publication number
- CN111352137B CN111352137B CN202010336832.6A CN202010336832A CN111352137B CN 111352137 B CN111352137 B CN 111352137B CN 202010336832 A CN202010336832 A CN 202010336832A CN 111352137 B CN111352137 B CN 111352137B
- Authority
- CN
- China
- Prior art keywords
- error
- broadcast ephemeris
- satellite
- rtk
- asynchronous
- 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
Images
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/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
-
- 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
Abstract
本发明公开了一种顾及广播星历误差的多模GNSS异步RTK定位方法,首先参考精密星历和精密钟差评定异步RTK模型中每颗卫星的广播星历轨道误差和钟差误差在通信延迟时间内变化的统计值;然后,将这些统计值作为附加产品添加到RTCM差分信号格式文件中,播发给RTK终端设备;其次,用户在接收到广播星历轨道误差和钟差在通信延迟时间变化的统计值后,结合终端RTK定位时的差分数据龄期(通信延迟)优化计算每颗卫星的随机模型误差,然后再利用双差观测值模型进行多模GNSS异步RTK定位解算,从而提高终端RTK的定位效果。
Description
技术领域
本发明属于GNSS导航定位技术领域,特别涉及一种顾及广播星历误差的多模GNSS异步RTK定位方法。
背景技术
实时动态差分(RTK)技术是实时高精度全球卫星导航定位系统(GNSS)应用中使用最为广泛的厘米级定位技术,RTK终端一般是通过无线网络接入高精度GNSS服务产品/基准站数据。在实际应用中不可避免会发生通信延迟或中断的情况,当基准站数据发生延迟时,为了获得流动站当前时刻的位置,终端必须使用基准站的历史观测值,此时,基准站和流动站的数据是异步的,称为异步RTK。RTK终端一般是基于广播星历,异步双差观测值中的广播星历轨道和钟差误差难以消除,RTK定位效果会随着通信延迟的增长而显著下降。GNSS系统已经从80年代的单GPS发展为如今的GPS/GLONASS/BDS/Galileo四大卫星导航定位系统。在进行异步RTK定位时,由于不同系统轨道类型和卫星钟特性的差异,不同广播星历误差在同样的通信延迟时间内会呈现明显的差异,因此采用多模GNSS异步RTK定位时需考虑不同卫星广播星历误差差异的影响。
发明内容
针对上述问题,本文提出一种顾及广播星历误差的多模GNSS异步RTK定位方法,通过评定GPS/GLONASS/BDS/Galileo卫星的广播星历误差在通信延迟时间的变化优化多模GNSS异步RTK定位的随机模型,提高多模GNSS异步RTK的定位精度及可靠性。
本发明首先参考精密星历和精密钟差评定异步RTK模型中每颗卫星的广播星历轨道误差和钟差误差在通信延迟时间内变化的统计值;然后,将这些统计值作为附加产品添加到RTCM差分信号格式文件中,播发给RTK终端设备;其次,用户在接收到广播星历轨道误差和钟差在通信延迟时间变化的统计值后,结合终端RTK定位时的差分数据龄期(通信延迟)优化计算每颗卫星的随机模型误差,然后再利用双差观测值模型进行多模GNSS异步RTK定位解算,从而提高终端RTK的定位效果。
本发明的具体技术方案如下:
一种顾及广播星历误差的多模GNSS异步RTK定位方法,其特征在于,包括以下步骤:
步骤1:参考精密星历和精密钟差,评定不同卫星的广播星历误差在通信延迟时间内变化的统计值,并将这些值发送给用户;其中时刻ti某卫星广播星历轨道和钟差误差在通信延迟dt时间内的变化值计算公式如下:
式中Δr(ti)、Δa(ti)、Δc(ti)分别表示广播星历轨道在时刻ti的径向、切向和法向误差,Δclk(ti)表示卫星钟差在时刻ti的误差,Δr(ti+dt)、Δa(ti+dt)、Δc(ti+dt)分别表示广播星历轨道在时刻ti+dt的径向、切向和法向误差,Δclk(ti+dt)表示卫星钟差在时刻ti+dt的误差,其中广播星历的轨道和钟差误差是以精密星历为参考计算的;
Δr(ti,ti+dt)、Δa(ti,ti+dt)、Δc(ti,ti+dt)和Δclk(ti,ti+dt)分别表示广播星历轨道径向、切向、法向以及卫星钟差误差从时刻ti到时刻ti+dt的变化值;卫星的广播星历轨道误差和钟差误差在通信延迟时间内变化的RMS统计值可用以下公式计算:
其中N表示采样点的总个数;σr(dt)、σa(dt)、σc(dt)、σclk(dt)分别表示广播星历轨道径向、切向、法向以及卫星钟差误差在dt通信延迟时间内变化的RMS统计结果;
步骤2:将步骤1评定的GPS/GLONASS/BDS/Galileo每颗卫星在特定通信延迟时间的广播星历轨道误差以及钟差误差变化的统计值进行编码,并作为附加产品添加到RTCM格式文件中,播发给用户端;
步骤4:异步RTK参数估计,在步骤3确定了异步RTK定位的随机模型后,结合函数模型
y=Aa+Bb+ε (5)
用卡尔曼滤波或者序贯最小二乘估计模糊度、坐标参数;其中a是双差模糊度参数向量,b是坐标等实参数向量;A、B是与之对应的系数矩阵;y是多模GNSS卫星的双差观测值向量,包括载波相位和伪距和观测值,ε是随机噪声向量;
步骤5:RTK固定解更新,在求解出模糊度浮点向量后固定双差模糊度,然后利用条件最小二乘/卡尔曼滤波更新即可得到RTK的固定解坐标。
上述顾及广播星历误差的多模GNSS异步RTK定位方法中,步骤2中将和通信延迟时间相关的广播星历轨道和钟差误差作为附加产品添加到RTCM格式文件中,将GPS,GLONASS,BDS,Galileo等系统所有卫星的广播星历轨道径向、切向和法向误差和卫星钟差误差在指定通信延迟时间内变化的统计值编码后发送给用户。
本发明的有益技术效果如下:
本发明给出了一种GNSS异步RTK模型中与通信延迟相关的广播星历轨道和钟差误差的评估方法,定量计算GPS/GLONASS/BDS/Galileo卫星的广播星历误差在特定通信延迟时间变化的统计值并播发给用户,通过充分考虑不同GNSS卫星广播星历误差变化的差异性,为异步RTK定位模型多模GNSS卫星观测值分配更为合理的权重,进而降低广播星历误差变化较大的卫星对异步RTK的定位不良影响,克服了GNSS RTK定位在通信短时中断/延迟等情况下定位精度低、可靠性差的问题。
附图说明
图1为本发明实施例的流程图。
图2为本发明GNSS卫星广播星历轨道径向误差随通信延迟时间变化的统计结果。
图3为本发明GNSS卫星广播星历钟差误差随通延迟信时间变化的统计结果。
图4为本发明不同时间段GNSS卫星广播星历误差在15秒延迟时间内变化的统计结果。
图5为本发明多模GNSS异步RTK随机模型优化前和优化后的比对效果说明图。
具体实施方式
以下结合附图和实施例详细说明本发明技术方案。参见图1,实施例的各步骤具体实现如下:
一、计算广播星历误差在通信延迟时间内变化的统计值
对于同步观测值,基准站和流动站对应的卫星发射时刻非常接近,由于广播星历轨道误差和钟差误差具有较强的时间相关性,因而站间做差后广播星历轨道误差和钟差误差基本可以消除;然而对于异步观测值,基准站和流动站对应的卫星发射时刻相差较大,站间做差后广播星历误差难以消除,此外不同系统不同卫星的广播星历误差在同样的通信延迟时间内的变化可能存在明显差异。通过参考精密星历和精密钟差,可以评定不同卫星的广播星历误差在通信延迟时间内变化的统计值,并将这些值发送给用户。需要说明的,不同时段的广播星历误差在给定通信延迟时间内的变化统计值具有一致性,因此,我们无需实时计算异步RTK模型的广播星历误差,当卫星更新或者广播星历出现异常时不定期更新与通信延迟时间相关的广播星历误差即可。
某时刻ti卫星广播星历轨道和钟差误差在通信延迟dt时间内的变化值计算公式如下:
式中Δr(ti)、Δa(ti)、Δc(ti)分别表示广播星历轨道在时刻ti的径向、切向和法向误差,Δclk(ti)表示卫星钟差在时刻ti的误差,Δr(ti+dt)、Δa(ti+dt)、Δc(ti+dt)和广播星历轨道在时刻ti+dt的径向、切向和法向误差,Δclk(ti+dt)表示卫星钟差在时刻ti+dt的误差,广播星历的轨道和钟差误差是以精密星历为参考计算的;
Δr(ti,ti+dt)、Δa(ti,ti+dt)、Δc(ti,ti+dt)和Δclk(ti,ti+dt)分别表示广播星历轨道径向、切向、法向以及卫星钟差误差从时刻ti到时刻ti+dt的变化值;卫星的广播星历轨道误差或者钟差误差在通信延迟dt时间内变化的统计值σ可用以下公式计算:
其中N表示采样点的总个数。
二、生成RTCM格式文件
附加产品包括和通信延迟时间相关的广播星历轨道或者钟差误差产品。主要包括广播星历轨道径向、切向、法向误差以及钟差误差在指定通信延迟时间内(如1s、2s.....30s)的变化统计值。
具体实施时,本领域技术人员可自行指定广播星历轨道、钟差误差在通信延迟时间内变化的统计值σ的传播方式。本发明提出一种改进的RTCM格式文件,将广播星历轨道和钟差误差在通信延迟时间内变化的统计值作为附加产品添加到RTCM数据流中,如表1所示。把该RTCM格式文件播发给RTK终端定位用户。用户可以根据当前定位时刻的差分数据龄期计算每颗卫星的广播星历轨道误差和钟差误差。
表1 广播星历误差在通信延迟时间内变化的统计值
电文号 | Message NO. | Bits(10) |
卫星系统 | System type | Bits(3) |
卫星号 | PRN | Bits(6) |
通信延迟时间 | t<sub>d</sub> | Bits(6) |
轨道径向 | σ<sub>r</sub>(t<sub>d</sub>) | Bits(8) |
轨道切向 | σ<sub>a</sub>(t<sub>d</sub>) | Bits(8) |
轨道法向 | σ<sub>c</sub>(t<sub>d</sub>) | Bits(8) |
卫星钟差 | σ<sub>clk</sub>(t<sub>d</sub>) | Bits(8) |
格式说明:第1行是RTCM的电文号,在实际使用中可自定义,应区别于现有的电文号;第2行是系统类型,包括GPS,GLONASS,BDS,Galileo等;第3行是卫星的PRN号;第4行是通信延迟时间,单位是秒,其范围可以设定为1-30秒;第5到7行分别是广播星历卫星轨道径向、切向和法向误差在指定通信延迟时间内变化的统计值,单位是毫米(mm);第5到7行分别是广播星历钟差误差在指定通信延迟时间内变化的统计值,单位是毫米(mm);第3列是不同字段所占用的比特数,在实际使用中可适当调整。
三、用户进行RTK定位
传统RTK随机模型中的先验误差仅包含是接收机端观测值的噪声,一般根据卫星高度角确定;是指空间信号测距误差,主要包括广播星历轨道和钟差误差,异步RTK模型中,这项误差主要和差分数据龄期/通信延迟时间dt相关;公式(4)给出了异步RTK模型中空间信号误差的计算方法,其中ωr,ωa,ωc分别表示轨道径向、切向、法向误差在卫星测距信号方向的投影系数。根据公式(3)和(4)即可计算出多模GNSS异步RTK模型中每颗卫星观测值的等效随机误差。
在确定了异步RTK定位的随机模型后,结合函数模型
y=Aa+Bb+ε (5)
即可用卡尔曼滤波或者序贯最小二乘估计模糊度、坐标参数等。其中a是双差双差模糊度参数向量,b是坐标等实参数向量等;A、B是与之对应的系数矩阵;y是双差观测值向量,包括载波相位和伪距和观测值,ε是随机噪声向量。
在求解出模糊度浮点向量及其方差协方差阵后,可以用LAMBDA方法等固定双差模糊度,然后利用条件最小二乘更新即可得到RTK的固定解坐标。
图2和图3给出了几颗GPS/BDS/GLONASS/Galileo卫星在特定通信延迟时间(1-15s)内广播星历径向误差以及钟差误差变化的RMS统计值结果,从图中可以看出随着通信延迟时间的增长广播星历误差变化的统计值逐渐变大,但是不同GNSS卫星在同样的通信延迟时间内广播星历径向以及钟差误差变化的RMS统计值差异较大,因此根据不同GNSS卫星广播星历的误差变化优化多模GNSS异步RTK的随机模型显得尤为必要。
图4是根据不同天的星历数据评估的几颗GNSS卫星在15秒通信延迟内广播星误差变化的RMS统计结果,从图中可以看出不同天的广播星历误差在给定通信延迟时间内变化的统计值较为一致,因此,在实际操作中无需实时评估异步RTK模型的广播星历误差,当卫星更新或者广播星历出现异常时不定期更新与通信延迟时间相关的广播星历误差即可。
图5是某超短基线多模GNSS异步RTK随机模型优化前和优化后的定位误差结果,异步RTK差分数据龄期为0-15s,从图中可以看出随机模型优化后多模GNSS异步RTK的定位精度有了显著提升。随机模型优化前东(E)、北(N)、高(U)方向误差的标准差(STD)分别为1.2,1.2和2.4cm,随机模型优后三个方向的STD分别为0.6、1.0和1.5cm,在东方向和高程方向定位误差显著降低。
以上实施例仅供说明本发明之用,而非对本发明的限制,有关技术领域的技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变换或变型,因此所有等同的技术方案,都落入本发明的保护范围。
Claims (1)
1.一种顾及广播星历误差的多模GNSS异步RTK定位方法,其特征在于,包括以下步骤:
步骤1:参考精密星历和精密钟差,评定不同卫星的广播星历误差在通信延迟时间内变化的统计值,并将这些值发送给用户;其中时刻ti某卫星广播星历轨道和钟差误差在通信延迟dt时间内的变化值计算公式如下:
式中Δr(ti)、Δa(ti)、Δc(ti)分别表示广播星历轨道在时刻ti的径向、切向和法向误差,Δclk(ti)表示卫星钟差在时刻ti的误差,Δr(ti+dt)、Δa(ti+dt)、Δc(ti+dt)分别表示广播星历轨道在时刻ti+dt的径向、切向和法向误差,Δclk(ti+dt)表示卫星钟差在时刻ti+dt的误差,其中广播星历的轨道和钟差误差是以精密星历为参考计算的;
Δr(ti,ti+dt)、Δa(ti,ti+dt)、Δc(ti,ti+dt)和Δclk(ti,ti+dt)分别表示广播星历轨道径向、切向、法向以及卫星钟差误差从时刻ti到时刻ti+dt的变化值;卫星的广播星历轨道误差和钟差误差在通信延迟时间内变化的RMS统计值采用以下公式计算:
其中N表示采样点的总个数;σr(dt)、σa(dt)、σc(dt)、σclk(dt)分别表示广播星历轨道径向、切向、法向以及卫星钟差误差在dt通信延迟时间内变化的RMS统计结果;
步骤2:将步骤1评定的GPS/GLONASS/BDS/Galileo每颗卫星在特定通信延迟时间的广播星历轨道误差以及钟差误差变化的统计值进行编码,并作为附加产品添加到RTCM格式文件中,播发给用户端;
步骤4:异步RTK参数估计,在步骤3确定了异步RTK定位的随机模型后,结合函数模型
y=Aa+Bb+ε (5)
用卡尔曼滤波或者序贯最小二乘估计模糊度、坐标参数;其中a是双差模糊度参数向量,b是坐标等实参数向量;A、B是与之对应的系数矩阵;y是多模GNSS卫星的双差观测值向量,包括载波相位和伪距的观测值,ε是随机噪声向量;
步骤5:RTK固定解更新,在求解出模糊度浮点向量后固定双差模糊度,然后利用条件最小二乘/卡尔曼滤波更新即可得到RTK的固定解坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010336832.6A CN111352137B (zh) | 2020-04-26 | 2020-04-26 | 一种顾及广播星历误差的多模gnss异步rtk定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010336832.6A CN111352137B (zh) | 2020-04-26 | 2020-04-26 | 一种顾及广播星历误差的多模gnss异步rtk定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111352137A CN111352137A (zh) | 2020-06-30 |
CN111352137B true CN111352137B (zh) | 2021-08-31 |
Family
ID=71193548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010336832.6A Active CN111352137B (zh) | 2020-04-26 | 2020-04-26 | 一种顾及广播星历误差的多模gnss异步rtk定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111352137B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113156477B (zh) * | 2021-04-15 | 2022-08-26 | 中国科学院精密测量科学与技术创新研究院 | 一种基于Android智能手机的高精度RTK定位方法 |
CN116133114A (zh) * | 2021-11-12 | 2023-05-16 | 中国移动通信有限公司研究院 | 定位方法、装置、设备及可读存储介质 |
CN114779301B (zh) * | 2022-03-30 | 2023-05-05 | 江苏城乡建设职业学院 | 一种基于广播星历的卫星导航实时精密单点定位方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104714244B (zh) * | 2015-03-31 | 2017-11-17 | 东南大学 | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 |
CN104808230A (zh) * | 2015-05-04 | 2015-07-29 | 东南大学 | 一种星钟与r-轨道误差混合估计的实时精密单点定位方法 |
CN106371115B (zh) * | 2016-06-27 | 2019-01-25 | 上海达华测绘有限公司 | 基于北斗短报文技术多模多频海上精密定位方法 |
CN107229061B (zh) * | 2017-07-18 | 2019-09-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
-
2020
- 2020-04-26 CN CN202010336832.6A patent/CN111352137B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111352137A (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108415050B (zh) | 一种基于低轨星座导航增强系统的ppp-rtk定位方法 | |
CN111352137B (zh) | 一种顾及广播星历误差的多模gnss异步rtk定位方法 | |
US10802160B2 (en) | Rapid determination of precise position by aiding data | |
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
EP1654559B1 (en) | Method for generating clock corrections for a wide-area or global differential gps system | |
EP2202532B1 (en) | Supporting a relative positioning | |
CN108076662B (zh) | 具有使用未组合公式来解算模糊度的能力的gnss接收机 | |
WO2008034728A1 (en) | Integrated mobile-terminal navigation | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
WO2017070732A1 (en) | A method of analysing a signal transmitted between a global satellite navigation satellite system and a receiver | |
CN111694030A (zh) | 一种基于格网虚拟观测值的bds局域差分方法及系统 | |
CN107966722B (zh) | 一种gnss钟差解算方法 | |
CN110780323B (zh) | 一种长距离下基于北斗三频信号的实时分米级定位方法 | |
CN109597105B (zh) | 一种顾及载波系统间偏差的gps/glonass紧组合定位方法 | |
CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及系统 | |
CN115933356B (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
CN111983641B (zh) | 一种用于实时生成北斗星基增强系统完好性参数的方法 | |
CN114879239B (zh) | 一种增强瞬时ppp固定解的区域三频整数钟差估计方法 | |
CN115373005A (zh) | 卫星导航信号间高精度产品转化方法 | |
CN114924295A (zh) | 一种载波相位平滑伪距定位方法、装置以及存储介质 | |
CN109884679B (zh) | 一种单模gnss系统跨频点混合双差rtk解算方法 | |
CN111948682B (zh) | 一种基于bds三频组合的伪距和载波相位随机模型构建方法 | |
CN114563806A (zh) | 一种Android移动设备的PPP-RTK实时定位方法和系统 | |
CN110618435B (zh) | 观测数据的生成方法与提高vrs稳定性的电文数据生成方法 | |
CN110045634B (zh) | 一种gnss基准站的非差误差建模方法 |
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 |