CN106772494A - 一种低成本gnss气压计组合rtk定位方法 - Google Patents
一种低成本gnss气压计组合rtk定位方法 Download PDFInfo
- Publication number
- CN106772494A CN106772494A CN201710027621.2A CN201710027621A CN106772494A CN 106772494 A CN106772494 A CN 106772494A CN 201710027621 A CN201710027621 A CN 201710027621A CN 106772494 A CN106772494 A CN 106772494A
- Authority
- CN
- China
- Prior art keywords
- delta
- double difference
- barometers
- dtri
- gnss
- 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
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/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
Abstract
本发明涉及一种低成本GNSS气压计组合RTK定位方法,属于全球卫星定位与导航技术领域。该方法基于低成本GNSS模块以及气压计模块,利用基站发送的单频伪距与载波相位观测值,实时解算载体的相对位置与双差实数模糊度参数,结合相对气压计提供的高程数据,采用约束LAMBDA方法固定模糊度的整数值,实时获得厘米级的相对定位结果。本发明所提供的一种低成本GNSS气压计组合RTK定位方法,可提高系统稳定性与可靠性,实现稳健的GNSS RTK定位。
Description
技术领域
本发明属于全球卫星定位与导航技术领域,涉及一种低成本GNSS气压计组合RTK(Real-time kinematic,实时动态差分)定位方法。
背景技术
全球卫星导航系统(GNSS)是我国自主建设的全球卫星导航系统,能够为地球表面和近地空间的广大用户提供全天时、全天候、高精度的定位、导航和授时服务,已广泛应用于国防、海陆空交通运输、测绘、移动通信、电力、电子金融、精细农业和减灾救灾等领域,是拓展人类活动和促进社会发展的重要空间基础设施。
GNSS精密RTK定位技术是全球卫星导航系统提供精密位置服务的关键技术之一,目前已广泛应用于驾考、精准农业、无人机、测绘等领域。它利用实时接收的载波相位和伪距数据以及参考站观测数据进行差分解算,并估计载波相位的整周模糊度,可以实时获得高精度的定位结果,极大的提高了作业效率,降低作业成本。但是基于目前的GNSS精密RTK产品一般基于采用双频方案,其硬件成本相对较高,从而限制了其应用范围以及产业化推广。
发明内容
有鉴于此,本发明的目的在于提供一种低成本GNSS气压计组合RTK定位方法,该方法基于低成本单频卫星导航模块以及气压计模块,结合基准站观测数据,实时解算载体的位置信息。
为达到上述目的,本发明提供如下技术方案:
一种低成本GNSS气压计组合RTK定位方法,该方法包括以下步骤:
S1:周跳探测;
S2:形成双差观测方程;
S3:Kalman滤波估计;
S4:气压计约束模糊度解算。
进一步,采用载波相位观测值进行解算时,由于天线周围环境影响,载波相位观测值不可避免的存在周跳现象,为了获得可靠的解算结果,需要实时的探测周跳,在步骤S1中,采用多普勒周跳探测方法,形成单差观测值探测周跳:
由于观测时间间隔较短,卫星变化不大,此方法可以有效探测到2周周跳;
δ为域值,设置为1.8,如果超过该域值,则认为有周跳发生,即重新初始化该模糊度参数。
进一步,在步骤S2中,接收的基准站(A)观测数据分别与载体模块(B)接收数据形成双差观测方程为:
式中,为A、B接收机的双差伪距观测值,表示以米为单位的A、B接收机的双差载波相位观测值,为A、B接收机到卫星的双差距离,λg为载波波长,为第一频段载波的A、B接收机双差电离层延迟,为A、B接收机双差对流层延迟,为载波的双差整周模糊度值,表示双差伪距的观测噪声,表示双差载波相位的观测噪声。
进一步,在步骤S3中,采用Kalman滤波估计对应的待估参数,其过程为:
Qk/k=(E-KkAk)HkQk (5)
式中yk为伪距与载波相位观测值,Ak为观测值的系数矩阵,为k-1历元的状态向量,为预测的状态向量,为当前k历元的状态向量,Kk为增益矩阵;
进一步,在步骤S4中,采用序贯最小二乘技术,计算模糊度浮点解及其方差协方差阵 此时采用约束LAMBDA方法解算模糊度,获得模糊度的固定解:
式中z为整数模糊度候选矢量,为最优n组模糊度解算结果;
对于每组模糊度值,更新基线向量:
式中为基线向量与模糊度的协方差,为浮点解基线向量,为固定解基线向量;
由于气压相对高度已知,此时利用已知的气压相对高度约束选取对应的模糊度整数解:
式中u为已知的气压相对高度,δu为误差范围,为固定解的气压相对高度。
本发明的有益效果在于:本发明所提供的一种低成本GNSS气压计组合RTK定位方法,可提高系统稳定性与可靠性,实现稳健的GNSS RTK定位。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明所述方法流程图。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
图1为本发明所述方法流程图,如图所示,本发明提供的一种低成本GNSS气压计组合RTK定位方法具体包括以下步骤:
步骤一:卫星导航模块与气压计实时接收数据:
使得卫星导航模块连接卫星导航天线,实时接收导航卫星的单频观测数据,并同时获取基站单频观测数据与气压计数据。
步骤二:周跳探测:
采用载波相位观测值进行解算时,由于天线周围环境影响,载波相位观测值不可避免的存在周跳现象,为了获得可靠的解算结果,需要实时的探测周跳;在本实施例中采用多普勒周跳探测方法。
形成单差观测值探测周跳:
由于观测时间间隔较短,卫星变化不大,此方法可以有效探测到2周周跳。
s为域值,一般可以设置为1.8,如果超过该域值,则认为有周跳发生,即重新初始化该模糊度参数。
步骤三:形成双差观测方程:
接收的基准站(A)观测数据分别与载体模块(B)接收数据形成双差观测方程为:
式中,为A、B接收机的双差伪距观测值,表示以米为单位的A、B接收机的双差载波相位观测值,为A、B接收机到卫星的双差距离,λg为载波波长,为第一频段载波的A、B接收机双差电离层延迟,为A、B接收机双差对流层延迟,为载波的双差整周模糊度值,表示双差伪距的观测噪声,表示双差载波相位的观测噪声。
由于两个天线距离很近,此时大气误差影响可以忽略,则原方程可以表示为:
步骤四:线性化观测方程:
公式(2)为非线性观测方程,为了解算相对基线分量,需要进行线性化处理,采用泰勒展开线性化(3)式可得:
式中:
式中为B接收机方向余玄,ΔxB,ΔyB,ΔzB为B接收机代估参数,为卫星k,j到B天线近似距离,为卫星k,j到A天线距离,A天线坐标可由单点定位获得。
步骤五:Kalman滤波估计实数模糊度参数:
Kalman滤波状态方程与过程方程如下:
采用Kalman滤波技术实时估计载体的姿态信息。
估计的参数主要包含三个位置参数、双差载波相位的模糊度:
观测值的系数矩阵为:
采用Kalman滤波估计对应的待估参数,其过程为:
Qk/k=(E-KkAk)HkQk (8)
式中yk为伪距与载波相位观测值,Ak为观测值的系数矩阵,为k-1历元的状态向量,为预测的状态向量,为当前k历元的状态向量,Kk为增益矩阵。
步骤六:高程约束模糊度解算:
采用Kalman滤波技术,计算模糊度浮点解及其方差协方差阵此时采用约束LAMBDA方法解算模糊度,获得模糊度的固定解。
式中z为整数模糊度候选矢量,为最优n组模糊度解算结果。
对于每组模糊度值,更新基线向量:
式中为基线向量与模糊度的协方差,为浮点解基线向量,为固定解基线向量。
由于气压相对高度已知,此时利用已知的气压相对高度约束选取对应的模糊度整数解。
式中u为已知的气压相对高度,δu为误差范围,为固定解的相对高程。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (5)
1.一种低成本GNSS气压计组合RTK定位方法,其特征在于:该方法包括以下步骤:
S1:周跳探测;
S2:形成双差观测方程;
S3:Kalman滤波估计;
S4:气压计约束模糊度解算。
2.根据权利要求1所述的一种低成本GNSS气压计组合RTK定位方法,其特征在于:在步骤S1中,采用多普勒周跳探测方法,形成单差观测值探测周跳:
由于观测时间间隔较短,卫星变化不大,此方法可以有效探测到2周周跳;
δ为域值,设置为1.8,如果超过该域值,则认为有周跳发生,即重新初始化该模糊度参数。
3.根据权利要求2所述的一种低成本GNSS气压计组合RTK定位方法,其特征在于:在步骤S2中,接收的基准站(A)观测数据分别与载体模块(B)接收数据形成双差观测方程为:
式中,为A、B接收机的双差伪距观测值,表示以米为单位的A、B接收机的双差载波相位观测值,为A、B接收机到卫星的双差距离,λg为载波波长,为第一频段载波的A、B接收机双差电离层延迟,为A、B接收机双差对流层延迟,为载波的双差整周模糊度值,表示双差伪距的观测噪声,表示双差载波相位的观测噪声。
4.根据权利要求3所述的一种低成本GNSS气压计组合RTK定位方法,其特征在于:在步骤S3中,采用Kalman滤波估计对应的待估参数,其过程为:
Qk/k=(E-KkAk)HkQk (5)
式中yk为伪距与载波相位观测值,Ak为观测值的系数矩阵,为k-1历元的状态向量,为预测的状态向量,为当前k历元的状态向量,Kk为增益矩阵;
5.根据权利要求4所述的一种低成本GNSS气压计组合RTK定位方法,其特征在于:在步骤S4中,采用序贯最小二乘技术,计算模糊度浮点解及其方差协方差阵此时采用约束LAMBDA方法解算模糊度,获得模糊度的固定解:
式中z为整数模糊度候选矢量,为最优n组模糊度解算结果;
对于每组模糊度值,更新基线向量:
式中为基线向量与模糊度的协方差,为浮点解基线向量,为固定解基线向量;由于气压相对高度已知,此时利用已知的气压相对高度约束选取对应的模糊度整数解:
式中u为已知的气压相对高度,δu为误差范围,为固定解的气压相对高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710027621.2A CN106772494A (zh) | 2017-01-13 | 2017-01-13 | 一种低成本gnss气压计组合rtk定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710027621.2A CN106772494A (zh) | 2017-01-13 | 2017-01-13 | 一种低成本gnss气压计组合rtk定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106772494A true CN106772494A (zh) | 2017-05-31 |
Family
ID=58945501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710027621.2A Pending CN106772494A (zh) | 2017-01-13 | 2017-01-13 | 一种低成本gnss气压计组合rtk定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772494A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085627A (zh) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星系统单频接收机定位精度的方法 |
CN109975849A (zh) * | 2017-12-28 | 2019-07-05 | 中移物联网有限公司 | 一种基线向量的确定方法、服务器和计算机存储介质 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110412638A (zh) * | 2019-08-16 | 2019-11-05 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110988948A (zh) * | 2019-11-07 | 2020-04-10 | 北京航空航天大学 | 一种基于动对动相对定位场景中完好性分析方法 |
WO2020133909A1 (zh) * | 2018-12-29 | 2020-07-02 | 北京金朋达航空科技有限公司 | 飞控与导航一体机 |
EP3746858A4 (en) * | 2018-01-30 | 2021-10-20 | Valmont Industries, Inc. | GPS ALIGNMENT SYSTEM AND METHOD USING REAL-TIME KINEMATICS (RTK) |
CN115616640A (zh) * | 2022-10-28 | 2023-01-17 | 北斗天地(北京)科技有限公司 | 一种智能手机gnss气压计组合定位方法 |
CN116106955A (zh) * | 2023-01-10 | 2023-05-12 | 北京航空航天大学 | 一种基于气压计与ppp紧组合的智能终端高精度定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008150389A1 (en) * | 2007-05-31 | 2008-12-11 | Navcom Technology, Inc. | Distance dependant error mitigation in real-time kinematic (rtk) positioning |
CN103383462A (zh) * | 2012-05-03 | 2013-11-06 | 中国科学院国家天文台 | 一种卫星导航通信终端高程定位的方法 |
CN105467415A (zh) * | 2016-01-18 | 2016-04-06 | 中国民航大学 | 一种基于差分气压高度约束的小型无人机rtk相对定位方法 |
CN106019344A (zh) * | 2016-07-19 | 2016-10-12 | 中国科学院重庆绿色智能技术研究院 | 一种通信基站天线位移监测方法 |
-
2017
- 2017-01-13 CN CN201710027621.2A patent/CN106772494A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008150389A1 (en) * | 2007-05-31 | 2008-12-11 | Navcom Technology, Inc. | Distance dependant error mitigation in real-time kinematic (rtk) positioning |
CN103383462A (zh) * | 2012-05-03 | 2013-11-06 | 中国科学院国家天文台 | 一种卫星导航通信终端高程定位的方法 |
CN105467415A (zh) * | 2016-01-18 | 2016-04-06 | 中国民航大学 | 一种基于差分气压高度约束的小型无人机rtk相对定位方法 |
CN106019344A (zh) * | 2016-07-19 | 2016-10-12 | 中国科学院重庆绿色智能技术研究院 | 一种通信基站天线位移监测方法 |
Non-Patent Citations (2)
Title |
---|
于兴旺 等: "不同 GPS接收机RTK定位精度比较分析", 《测绘信息与工程》 * |
许国昌: "《GPS理论、算法与应用》", 31 December 2011 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975849A (zh) * | 2017-12-28 | 2019-07-05 | 中移物联网有限公司 | 一种基线向量的确定方法、服务器和计算机存储介质 |
EP3746858A4 (en) * | 2018-01-30 | 2021-10-20 | Valmont Industries, Inc. | GPS ALIGNMENT SYSTEM AND METHOD USING REAL-TIME KINEMATICS (RTK) |
CN109085627B (zh) * | 2018-06-11 | 2019-07-12 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星系统单频接收机定位精度的方法 |
CN109085627A (zh) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星系统单频接收机定位精度的方法 |
WO2020133909A1 (zh) * | 2018-12-29 | 2020-07-02 | 北京金朋达航空科技有限公司 | 飞控与导航一体机 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110109163B (zh) * | 2019-04-09 | 2022-01-07 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110412638A (zh) * | 2019-08-16 | 2019-11-05 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110412638B (zh) * | 2019-08-16 | 2023-03-28 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110988948A (zh) * | 2019-11-07 | 2020-04-10 | 北京航空航天大学 | 一种基于动对动相对定位场景中完好性分析方法 |
CN110988948B (zh) * | 2019-11-07 | 2021-11-02 | 北京航空航天大学 | 一种基于动对动相对定位场景中完好性分析方法 |
CN115616640A (zh) * | 2022-10-28 | 2023-01-17 | 北斗天地(北京)科技有限公司 | 一种智能手机gnss气压计组合定位方法 |
CN115616640B (zh) * | 2022-10-28 | 2023-10-20 | 北斗天地(北京)科技有限公司 | 一种智能手机gnss气压计组合定位方法 |
CN116106955A (zh) * | 2023-01-10 | 2023-05-12 | 北京航空航天大学 | 一种基于气压计与ppp紧组合的智能终端高精度定位方法 |
CN116106955B (zh) * | 2023-01-10 | 2023-09-05 | 北京航空航天大学 | 一种基于气压计与ppp紧组合的智能终端高精度定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772494A (zh) | 一种低成本gnss气压计组合rtk定位方法 | |
CN110412638B (zh) | 一种低成本三天线gnss rtk定位及测姿方法 | |
CN106646538B (zh) | 一种基于单差滤波的变形监测gnss信号多路径改正方法 | |
CN104483691B (zh) | 一种gnss组合精密单点定位方法 | |
Han et al. | Accurate absolute GPS positioning through satellite clock error estimation | |
CN103529462B (zh) | 一种用于全球导航卫星系统的动态周跳探测与修复方法 | |
US6229479B1 (en) | Relative position measuring techniques using both GPS and GLONASS carrier phase measurements | |
CN101770033B (zh) | 连续运行参考站系统站间整周模糊度网络固定方法 | |
CN108802782B (zh) | 一种惯导辅助的北斗三频载波相位整周模糊度求解方法 | |
CN104483690A (zh) | 一种gnss三频精密单点定位模糊度固定方法 | |
CN113267796A (zh) | 一种双天线gnss、rtk定位及测向方法 | |
US11867823B2 (en) | System and method for determining GNSS positioning corrections | |
JP2007507718A (ja) | 到着時間差(tdoa)測位サービスの方法とシステム | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN111694037B (zh) | 一种终端定位方法及终端设备 | |
CN109407126A (zh) | 一种多模接收机联合定位解算的方法 | |
CN111796309A (zh) | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 | |
CN103592653B (zh) | 用于局部地区单频卫星导航用户的电离层延迟修正方法 | |
Hou et al. | Recent advances and perspectives in GNSS PPP-RTK | |
US10024974B1 (en) | Method and apparatus applicable to positioning in NLOS environment | |
CN106019344A (zh) | 一种通信基站天线位移监测方法 | |
CN105510946A (zh) | 一种bds卫星载波相位整周模糊度快速解算方法 | |
CN104991265A (zh) | 一种北斗卫星导航系统用户统一性定位方法 | |
CN116755126A (zh) | 一种基于三维模型映射匹配的北斗实时精准定位方法 | |
CN102830410B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |