CN106767789B - 一种基于自适应卡尔曼滤波的行人航向最优融合方法 - Google Patents

一种基于自适应卡尔曼滤波的行人航向最优融合方法 Download PDF

Info

Publication number
CN106767789B
CN106767789B CN201710022551.1A CN201710022551A CN106767789B CN 106767789 B CN106767789 B CN 106767789B CN 201710022551 A CN201710022551 A CN 201710022551A CN 106767789 B CN106767789 B CN 106767789B
Authority
CN
China
Prior art keywords
magnetic
axis
angle
epsilon
adaptive kalman
Prior art date
Application number
CN201710022551.1A
Other languages
English (en)
Other versions
CN106767789A (zh
Inventor
黄欣
熊智
许建新
徐丽敏
孔雪博
赵宣懿
万众
李一博
Original Assignee
南京航空航天大学
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 南京航空航天大学 filed Critical 南京航空航天大学
Priority to CN201710022551.1A priority Critical patent/CN106767789B/zh
Publication of CN106767789A publication Critical patent/CN106767789A/zh
Application granted granted Critical
Publication of CN106767789B publication Critical patent/CN106767789B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments

Abstract

本发明公开了一种基于自适应卡尔曼滤波器的行人航向最优融合方法,属于行人导航技术领域。本方法构建了9阶自适应卡尔曼滤波算法模型,每到量测时刻,融合捷联航向角与磁航向角,利用磁传感器实时统计信息,建立磁异常辨识模型,进而在不同的辨识状态下实时修正自适应卡尔曼滤波量测噪声阵,实现捷联航向角与磁航向角的最优融合,保证组合航向角的高精度。本发明采用低精度的消费级传感器芯片,无论室内外,均能有效的保证行人航向角高精度,实用性较强。

Description

一种基于自适应卡尔曼滤波的行人航向最优融合方法

技术领域

本发明涉及一种基于自适应卡尔曼滤波的行人航向最优融合方法,属于行人导航技术领域。

背景技术

行人导航技术已掀起国内外研究热潮,其中航向作为核心技术之一,对行人导航定位精度起着至关重要的作用。目前较为普遍的做法是利用陀螺仪获得角速率,进而四元素解算得到角度,但是低成本的惯性传感器,其误差较大且发散严重,长时间内将会导致行人定位失效。地磁场作为地球的固有属性,可以用于计算磁航向角,但是室内地磁受到钢筋、管道等其他金属制品的影响,会产生较大的误差,无法满足所有路段高精度航向的要求。目前,国内外主要考察磁航向角与陀螺解算出的角度的差值大小,当这一差值大于某一预定阈值时,即为判定有磁场异常情况。但是这种判别方法受限于磁航向角的解算精度,本身磁传感器尚未标定完全会带有误差,投影解算航向时又会引入二次误差,判别精度较低。同时针对磁航向角与捷联解算出来的航向角融合方面,其仅仅只有完全融合与完全不融合两种情况,其灵活性与适应性较差,当磁异常情况微弱时,可以尝试提取有效信息进行航向融合。

发明内容

本发明采用低精度的陀螺传感器与磁传感器,提出了一种基于自适应卡尔曼滤波的行人航向最优融合方法,解决不同磁环境下的航向融合问题。利用磁异常辨识算法有效区分磁异常环境,进而在不同环境下恰当的选取量测噪声阵,构建自适应卡尔曼滤波器灵活实现捷联航向角与磁航向角的最优融合,实现室内外不同环境下行人航向的高精度。

本发明为解决其技术问题采用如下技术方案:

一种基于自适应卡尔曼滤波的行人航向最优融合方法,包括如下步骤:

步骤1,首先建立自适应卡尔曼滤波状态方程,选用“东北天”地理坐标系,构建9阶状态模型,如下式所示

其中为东向平台角误差;为北向平台角误差;为天向平台角误差;εbx为x轴陀螺随机常数;εby为y轴陀螺随机常数;εbz为z轴陀螺随机常数;εrx为x轴陀螺一阶马尔科夫过程;εry为y轴陀螺一阶马尔科夫过程;εrz为z轴陀螺一阶马尔科夫过程;W为系统随机过程噪声序列;A为系统矩阵;G为系统噪声矩阵;W为系统噪声序列;X为状态量;为状态量导数;wgx为x轴随机白噪声驱动;wgy为y轴随机白噪声驱动;wgz为z轴随机白噪声驱动;wrx为x轴马尔科夫白噪声驱动;wry为y轴马尔科夫白噪声驱动;wrz为z轴马尔科夫白噪声驱动;

步骤2,在步骤1自适应卡尔曼滤波状态方程建立好的基础之上,开始导航;利用陀螺仪每0.005秒采集一次数据,经误差修正后,通过四元素解算得出当前捷联航向角,利用误差修正后的磁传感器信息解算出当前磁航向角;

步骤3,在步骤2的基础上,判断当前解算时间是否达到1秒,无则返回步骤2,有则进行步骤4;

步骤4,在步骤3的基础上,研究利用磁传感器实时统计信息进行磁异常辨识,再利用二维椭圆标定算法修正磁传感器信息之后,按照下式构造磁异常辨识模型

上式代表第k时刻的磁环境辨识情况,其中,σ3为第一种磁异常判别参数;σ4为第二种磁异常判别参数;var()代表方差函数;min()代表求取最小值函数;mag1为第一个时刻总磁场强度;代表窗口大小为N的滑动数组;N为数组大小;

数组存储总磁场强度并实时更新,总磁场强度如下式所示,由三轴磁传感器测量值的平方和开根号组成。

其中:magi为第i时刻总磁场强度;magix为x轴磁场强度;magiy为y轴磁场强度;magiz为z轴磁场强度;

步骤5:在步骤4的基础上,进行自适应卡尔曼滤波量测建模,构建一维量测模型,量测周期为1秒,量测方程如下式所示:

其中,为捷联解算出的航向角;为磁航向角;为姿态误差角;为航向噪声;Hv(t)为量测矩阵;X(t)为状态量;Vv(t)为量测噪声;磁航向角解算公式如下式所示,其中X、Y分别为前向和横侧向所测量得到的磁信息。

姿态误差角与平台误差角存在如下关系:

其中,θ为俯仰角,φx、φy、φy为三个平台误差角。因此:

同时利用步骤4的结果,依据下式对卡尔曼滤波器量测噪声阵Vv(t)进行修正。

其中λ为设定的经验值;η为事先人为设定的增益常量;ε3与ε4为通过统计学方法获得的参数值,即在磁正常环境下由式(2)获得的均值。

步骤6:在步骤5的基础上,利用上述自适应卡尔曼滤波器对组合航向角进行修正,同时将陀螺估计误差反馈给陀螺误差修正模型,并返回步骤2。

本发明的有益效果如下:

1.采用低精度的消费级传感器芯片,无论室内外,均能有效的保证行人航向角高精度,实用性较强。

2.解决了磁异常环境下磁航向角误用的问题。

3.传感器成本低且普及较广,算法的实用性与推广性较强。

附图说明

图1为自适应卡尔曼滤波的行人航向最优融合流程示意图。

具体实施方案

下面结合附图对本发明创造做进一步详细说明。

如图1所示,导航解算过程中,陀螺仪可以实时获得角速率信息,通过四元素解算,实时获得捷联航向角,解算频率一般可设置为50Hz。磁传感器实时采集磁信息,主要分为两部分,一部分用于解算出磁航向角,一部分用于磁异常辨识,实时灵活调整量测噪声阵。以捷联航向角为主,每过1秒,利用自适应卡尔曼滤波器,利用磁航向角建立观测方程,实时修正捷联航向角误差,提高航向角精度,并返回估计误差对惯性传感器进行误差修正。

基于自适应卡尔曼滤波的行人航向最优融合算法主要包括以下三个步骤:

步骤一:自适应卡尔曼滤波器状态建模

选用“东北天”地理坐标系,构建9阶状态模型,如式(1)所示。

其中为东向平台角误差;为北向平台角误差;为天向平台角误差;εbx为x轴陀螺随机常数;εby为y轴陀螺随机常数;εbz为z轴陀螺随机常数;εrx为x轴陀螺一阶马尔科夫过程;εry为y轴陀螺一阶马尔科夫过程;εrz为z轴陀螺一阶马尔科夫过程;W为系统随机过程噪声序列;A为系统矩阵;G为系统噪声矩阵;W为系统噪声序列;X为状态量;为状态量导数;wgx为x轴随机白噪声驱动;wgy为y轴随机白噪声驱动;wgz为z轴随机白噪声驱动;wrx为x轴马尔科夫白噪声驱动;wry为y轴马尔科夫白噪声驱动;wrz为z轴马尔科夫白噪声驱动。

步骤二:磁异常辨识模型

研究利用磁传感器实时统计信息进行磁异常辨识,磁异常环境下,无论是磁场强度模值还是磁场强度变化趋势,均与磁正常环境有着较大差异。再利用二维椭圆标定算法修正磁传感器信息之后,按照(2)式构造磁异常辨识模型。

上式代表第k时刻的磁环境辨识情况。其中,σ3为第一种磁异常判别参数;σ4为第二种磁异常判别参数;var()代表方差函数;min()代表求取最小值函数;mag1为第一时刻总磁场强度;代表窗口大小为N的滑动数组;N为数组大小;数组存储总磁场强度并实时更新,总磁场强度如式(3)所示,由三轴磁传感器测量值的平方和开根号组成。

其中:magi为第i时刻总磁场强度;magix为x轴磁场强度大小;magiy为y轴磁场强度大小;magiz为z轴磁场强度大小。

步骤三:自适应卡尔曼滤波器量测建模

构建一维量测模型,量测周期为1秒,量测方程如下式所示:

其中,为捷联解算出的航向角;为磁航向角;为姿态误差角;为航向噪声;Hv(t)为量测矩阵;X(t)为状态量;Vv(t)为量测噪声。磁航向角解算公式如下式(5)所示,其中X、Y分别为前向和横侧向所测量得到的磁信息。

姿态误差角与平台误差角存在如下关系:

其中,θ为俯仰角,φx、φy、φy为三个平台误差角。因此:

卡尔曼滤波器量测噪声阵Vv(t)为

其中λ为设定的经验值,η为事先人为设定的增益常量,ε3与ε4为通过统计学方法获得的参数值(即在磁正常环境下由式(2)获得的均值)。

至此构建了自适应卡尔曼滤波器,实时的针对当前环境进行磁异常辨识,在不同的辨识结果下,适应性的调整卡尔曼滤波器量测噪声阵,实现捷联航向角与磁航向角的最优融合。

整个解算流程如下:

1、首先依据状态方程和量测方程构建自适应卡尔曼滤波器模型;

2、开始导航:利用陀螺仪采集数据,经误差修正后,通过四元素解算得出当前捷联航向角,利用误差修正后的磁传感器信息解算出当前磁航向角;

3、是否达到1秒,无则返回步骤2,是则进行步骤4、5、6;

4、利用磁异常辨识算法辨识出磁正常状态和异常状态,根据辨识结果实时调整自适应卡尔曼滤波器中的量测噪声阵;

5、滤波,获得修正过后的组合航向角,同时将陀螺估计误差反馈给陀螺误差修正模型;

6、返回步骤2。

Claims (1)

1.一种基于自适应卡尔曼滤波的行人航向最优融合方法,其特征在于,包括如下步骤:
步骤1,首先建立自适应卡尔曼滤波状态方程,选用“东北天”地理坐标系,构建9阶状态模型,如下式所示
其中为东向平台角误差;为北向平台角误差;为天向平台角误差;εbx为x轴陀螺随机常数;εby为y轴陀螺随机常数;εbz为z轴陀螺随机常数;εrx为x轴陀螺一阶马尔科夫过程;εry为y轴陀螺一阶马尔科夫过程;εrz为z轴陀螺一阶马尔科夫过程;W为系统随机过程噪声序列;A为系统矩阵;G为系统噪声矩阵;W为系统噪声序列;X为状态量;为状态量导数;wgx为x轴随机白噪声驱动;wgy为y轴随机白噪声驱动;wgz为z轴随机白噪声驱动;wrx为x轴马尔科夫白噪声驱动;wry为y轴马尔科夫白噪声驱动;wrz为z轴马尔科夫白噪声驱动;
步骤2,在步骤1自适应卡尔曼滤波状态方程建立好的基础之上,开始导航;利用陀螺仪每0.005秒采集一次数据,经误差修正后,通过四元素解算得出当前捷联航向角,利用误差修正后的磁传感器信息解算出当前磁航向角;
步骤3,在步骤2的基础上,判断当前解算时间是否达到1秒,无则返回步骤2,有则进行步骤4;
步骤4,在步骤3的基础上,研究利用磁传感器实时统计信息进行磁异常辨识,再利用二维椭圆标定算法修正磁传感器信息之后,按照下式构造磁异常辨识模型
上式代表第k时刻的磁环境辨识情况,其中,σ3为第一种磁异常判别参数;σ4为第二种磁异常判别参数;var()代表方差函数;min()代表求取最小值函数,mag1为第一个时刻总磁场强度;代表窗口大小为N的滑动数组;N为数组大小;
数组存储总磁场强度并实时更新,总磁场强度如下式所示,由三轴磁传感器测量值的平方和开根号组成
其中:magi为第i时刻总磁场强度;magix为x轴磁场强度;magiy为y轴磁场强度;magiz为z轴磁场强度;
步骤5:在步骤4的基础上,进行自适应卡尔曼滤波量测建模,构建一维量测模型,量测周期为1秒,量测方程如下式所示:
其中,为捷联解算出的航向角;为磁航向角;为姿态误差角;为航向噪声;Hv(t)为量测矩阵;X(t)为第t时刻状态量;Vv(t)为量测噪声;磁航向角解算公式如下式所示,其中X、Y分别为前向和横侧向所测量得到的磁信息
姿态误差角与平台误差角存在如下关系:
其中,θ为俯仰角,φx、φy、φy为三个平台误差角;因此:
同时利用步骤4的结果,依据下式对卡尔曼滤波器量测噪声阵Vv(t)进行修正
其中λ为设定的经验值;η为事先人为设定的增益常量;ε3与ε4为通过统计学方法获得的参数值,即在磁正常环境下由式(2)获得的均值;
步骤6:在步骤5的基础上,利用上述自适应卡尔曼滤波器对组合航向角进行修正,同时将陀螺估计误差反馈给陀螺误差修正模型,并返回步骤2。
CN201710022551.1A 2017-01-12 2017-01-12 一种基于自适应卡尔曼滤波的行人航向最优融合方法 CN106767789B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710022551.1A CN106767789B (zh) 2017-01-12 2017-01-12 一种基于自适应卡尔曼滤波的行人航向最优融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710022551.1A CN106767789B (zh) 2017-01-12 2017-01-12 一种基于自适应卡尔曼滤波的行人航向最优融合方法

Publications (2)

Publication Number Publication Date
CN106767789A CN106767789A (zh) 2017-05-31
CN106767789B true CN106767789B (zh) 2019-12-24

Family

ID=58947965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710022551.1A CN106767789B (zh) 2017-01-12 2017-01-12 一种基于自适应卡尔曼滤波的行人航向最优融合方法

Country Status (1)

Country Link
CN (1) CN106767789B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107300697A (zh) * 2017-06-07 2017-10-27 南京航空航天大学 基于无人机的运动目标ukf滤波方法
CN109147389A (zh) * 2018-08-16 2019-01-04 大连民族大学 自主汽车或者辅助驾驶系统规划路线的方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103115624A (zh) * 2013-01-24 2013-05-22 南京航空航天大学 一种基于地磁匹配的地磁日变修正方法
CN103759730A (zh) * 2014-01-16 2014-04-30 南京师范大学 一种基于导航信息双向融合的行人与智能移动载体的协同导航系统及其导航方法
CN103776446A (zh) * 2013-10-29 2014-05-07 哈尔滨工程大学 一种基于双mems-imu的行人自主导航解算算法
JP2015055543A (ja) * 2013-09-11 2015-03-23 株式会社フジクラ 磁気素子制御装置及び磁気素子制御方法
US9037224B1 (en) * 2010-08-02 2015-05-19 Chi Yung Fu Apparatus for treating a patient
CN104864874A (zh) * 2015-06-19 2015-08-26 北京理工大学 一种低成本单陀螺航位推算导航方法及系统
CN104897172A (zh) * 2015-06-18 2015-09-09 南京航空航天大学 基于运动捕捉系统的旋转mems惯导磁航向角误差补偿方法
CN104931028A (zh) * 2015-06-30 2015-09-23 北京联合大学 一种基于深度学习的三轴磁电子罗盘误差补偿方法
CN105043385A (zh) * 2015-06-05 2015-11-11 北京信息科技大学 一种行人自主导航定位的自适应卡尔曼滤波方法
CN105353428A (zh) * 2015-12-11 2016-02-24 吉林大学 一种地面参考区磁场延拓的地空协同电磁数据校正方法
CN105378431A (zh) * 2013-07-12 2016-03-02 微软技术许可有限责任公司 利用磁场异常进行户内定位-寻找
WO2016048566A1 (en) * 2014-09-26 2016-03-31 Intel Corporation Virtual gyroscope using dual magnetometers for electronic devices
CN106123900A (zh) * 2016-06-20 2016-11-16 南京航空航天大学 基于改进型互补滤波的室内行人导航磁航向解算方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9037224B1 (en) * 2010-08-02 2015-05-19 Chi Yung Fu Apparatus for treating a patient
CN103115624A (zh) * 2013-01-24 2013-05-22 南京航空航天大学 一种基于地磁匹配的地磁日变修正方法
CN105378431A (zh) * 2013-07-12 2016-03-02 微软技术许可有限责任公司 利用磁场异常进行户内定位-寻找
JP2015055543A (ja) * 2013-09-11 2015-03-23 株式会社フジクラ 磁気素子制御装置及び磁気素子制御方法
CN103776446A (zh) * 2013-10-29 2014-05-07 哈尔滨工程大学 一种基于双mems-imu的行人自主导航解算算法
CN103759730A (zh) * 2014-01-16 2014-04-30 南京师范大学 一种基于导航信息双向融合的行人与智能移动载体的协同导航系统及其导航方法
WO2016048566A1 (en) * 2014-09-26 2016-03-31 Intel Corporation Virtual gyroscope using dual magnetometers for electronic devices
CN105043385A (zh) * 2015-06-05 2015-11-11 北京信息科技大学 一种行人自主导航定位的自适应卡尔曼滤波方法
CN104897172A (zh) * 2015-06-18 2015-09-09 南京航空航天大学 基于运动捕捉系统的旋转mems惯导磁航向角误差补偿方法
CN104864874A (zh) * 2015-06-19 2015-08-26 北京理工大学 一种低成本单陀螺航位推算导航方法及系统
CN104931028A (zh) * 2015-06-30 2015-09-23 北京联合大学 一种基于深度学习的三轴磁电子罗盘误差补偿方法
CN105353428A (zh) * 2015-12-11 2016-02-24 吉林大学 一种地面参考区磁场延拓的地空协同电磁数据校正方法
CN106123900A (zh) * 2016-06-20 2016-11-16 南京航空航天大学 基于改进型互补滤波的室内行人导航磁航向解算方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
"Clinical Predictors of Abnormal Magnetic Resonance Imaging Findings in Patients With Asymmetric Sensorineural Hearing Loss";Ahsan SF等;《JAMA Otolaryngol Head Neck Surg》;20151231;第141卷(第5期);451-456 *
"Efficient computer aided diagnosis of abnormal parts detection in magnetic resonance images using hybrid abnormality detection algorithm";C. Lakshmi Devasena等;《Central European Journal of Computer Science》;20131231;第3卷(第3期);117-128 *
"High Accuracy Acquisition of the Magnetic Heading Signal";L Zhang;《International Conference on Mechatronics and Automation》;20071231;2304-2308 *
"低阶卡尔曼滤波器在低成本SIAH RS 中的实现研究";赖际舟等;《中国空间科学技术》;20040229(第1期);61-65 *
"基于惯性元件与磁罗盘信息融合的步行者定位系统";刘瑜;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160415(第4期);I136-765 *
"基于惯性系下陀螺误差在线估计修正的惯性与星光组合导航方法";赵慧等;《兵工学报》;20161231;第37卷(第12期);2259-2267 *
"磁罗盘与光纤陀螺的组合航向测量系统设计";李俊等;《计算机仿真》;20080229;第25卷(第2期);13-15,19 *
"航空磁异常探测关键技术研究";韩磊;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150215(第02期);A011-410 *
"非高斯背景噪声下的微弱磁异常信号检测算法";张坚等;《海军工程大学学报》;20110831;第23卷(第4期);22-26 *

Also Published As

Publication number Publication date
CN106767789A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
US9551561B2 (en) Determining location using magnetic fields from AC power lines
CN104736963B (zh) 测绘系统和方法
CN104101908B (zh) 一种用球形偏置线圈测量地磁场矢量的装置
CN106031263A (zh) 用于位置估计的方法及系统
US6711517B2 (en) Hybrid inertial navigation method and device
CN103379619B (zh) 一种定位方法和系统
AU703786B2 (en) Method to determine the direction of the Earth's magnetic field
US7451549B1 (en) Automatic calibration of a three-axis magnetic compass
EP1886517B1 (fr) Procédé et dispositif de localisation d'un terminal dans un réseau local sans fil
AU2005201584B2 (en) Method for estimating the accuracy of azimuthal orientations and portable sighting device
CN105043380A (zh) 基于微机电传感器、WiFi定位、磁场匹配的室内导航方法
CN104703130A (zh) 基于室内定位的定位方法及其装置
US20170307403A1 (en) Methods for improved heading estimation
TWI262317B (en) Method for controlling a magnetic flux sensor, a control device, and a mobile terminal device
US20110066392A1 (en) Systems and methods for calibration of gyroscopes and a magnetic compass
CN103630137B (zh) 一种用于导航系统的姿态及航向角的校正方法
CN103941309B (zh) 地磁传感器校准设备及其方法
CN104121905A (zh) 一种基于惯性传感器的航向角获取方法
CN105607093B (zh) 一种组合导航系统及获取导航坐标的方法
EP3364153A1 (en) Method for updating all attitude angles of agricultural machine on the basis of nine-axis mems sensor
CN105652306A (zh) 基于航迹推算的低成本北斗与mems紧耦合定位系统及方法
CN105222788B (zh) 基于特征匹配的飞行器航路偏移误差的自校正方法
CN103679711A (zh) 一种遥感卫星线阵推扫光学相机在轨外方位参数标定方法
CN102034238B (zh) 基于光学成像测头和视觉图结构的多摄像机系统标定方法
AU2014309335B2 (en) Drilling methods and systems with automated waypoint or borehole path updates based on survey data corrections

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