CN107677267A - 基于mems‑imu的室内行人导航航向反馈修正方法 - Google Patents
基于mems‑imu的室内行人导航航向反馈修正方法 Download PDFInfo
- Publication number
- CN107677267A CN107677267A CN201710725704.9A CN201710725704A CN107677267A CN 107677267 A CN107677267 A CN 107677267A CN 201710725704 A CN201710725704 A CN 201710725704A CN 107677267 A CN107677267 A CN 107677267A
- Authority
- CN
- China
- Prior art keywords
- course
- pedestrian
- principal direction
- course angle
- angle
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/206—Instruments for performing navigational calculations specially adapted for indoor navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于MEMS‑IMU的室内行人导航航向反馈修正方法,属于行人导航技术领域。该方法结合典型建筑结构定义行人主方向上直行为航向反馈条件,构建航向差值作为状态量进行自适应无迹卡尔曼滤波(AUKF),将最优估计结果与捷联航向角之和作为反馈航向角,并运用反馈航向值参与四元数解算,从而减少航向误差的累积。与不反馈时相比,精确解算航向角的时间提高了两倍以上,对导航精度的提高有较强实用性。
Description
技术领域
本发明涉及一种基于MEMS-IMU的室内行人导航航向反馈修正方法,属于行人导航技术领域。
背景技术
行人导航技术已掀起国内外研究热潮,其中航向作为核心技术之一,对行人导航定位精度起着至关重要的作用。传统的姿态解算中利用陀螺仪获得角速率,进而通过四元数解算得到角度,因其利用积分原理,故其解算的航向角误差会随时间发散,不满足与长时间的行人定位需求。另一个较为普遍做法是利用磁力计得到航行角信息,地磁场作为地球的固有属性,可以用于计算磁航向角,但是室内地磁受到钢筋、管道等其他金属制品的影响,进而产生严重的航向误差,无法满足室内定位下高精度航向的要求。目前对航向角的修正算法中,大都需要外界信号的参与,如基于移动通信网络的辅助GPS(A-GPS)、伪卫星(Pseudolite)、无线局域网(WLAN)、射频标签(RFID)、Zigbee、蓝牙(Bluetooth,BT)等,这些定位技术都需要提前布置基础设施,成本较高,后期维护量大等缺点,不适用与如消防人员的位置跟踪等特殊的救援领域。
发明内容
针对上述问题,本发明旨在解决减少纯惯性导航中航向角误差随时间累积的问题。采用MEMS陀螺仪和加速度计等惯性器件,结合传统的捷联航姿角解算算法和PDR算法,以及室内特殊的环境构造,将行人直行的角速度和航向等特征作为航向反馈条件,并利用航向差值进行自适应无迹卡尔曼滤波,将最优结果和当前捷联航向角作为下一时刻航向角参与解算,实现的航向精准导航时间满足消防人员定位的特殊需求。
为了实现上述目的本发明采用如下技术方案:基于MEMS-IMU的室内行人导航航向反馈修正方法,包括以下步骤:
(1)通过建筑物结构特征,构建主方向和主方向航向范围。
(2)根据加速度计测得的数据进行步态识别,判断行人是否处于行走状态。
(3)根据相邻步数间航向角和相邻步数间的角速度信息判断行人是否处于直行状态。
(4)当判断行人在主方向范围内直行,用当前航向角yawi与当前行人所处主方向航向角差值作为状态量,进行自适应无迹卡尔曼滤波,获得最优估计结果
(5)将最优估计结果与当前航向角yawi之和作为此时航向输出值参与位置解算,并将该航向输出值作为下一时刻航向角参与四元数和欧拉角的捷联解算。
进一步地,步骤(1)所述构建主方向和主方向航向范围包括:将建筑物内的走廊及过道划分为四个方向,定义为建筑物的主方向;设定行人在主方向上直行的晃动角度范围,建立以主方向为基准的主方向航向角范围。
进一步地,步骤(2)所述判断行人是否处于行走状态的方法为,设置一滑动窗口求取加速度方差,然后与预设的方差阈值比较,判断行人是否处于行走状态。
进一步地,步骤(3)所述相邻步数间航向角和相邻步数间的角速度信息判断行人是否处于直行状态,分别采用三步间航向角差值和三步间隔内记录的角速度值进行判断;所述三步间航向角差值作为判断标准时,任两步的航向角绝对差值均小于预设航向角阈值时,行人直行;所述三步间隔内记录的角速度值作为判断标准时,三步间隔内的角速度绝对值小于预设阈值的个数占总个数的比例小于比例设定阈值时,为转向,大于比例设定阈值时,为弯道行走或者直行。判断所述弯道行走或者直行采用的方法是,若记录的角速度中大于零和小于零的个数分差不多时,为直行,否则为弯道行走。
本发明的有益效果如下:
采用低成本,低功耗,轻质量的MEMS惯性传感器,易于普及推广。
算法流程简单易操作,无需耗费处理器的太多资源。
与不反馈时相比,精确解算航向角的时间提高了两倍以上,在室内导航中可有效提高导航精度,实用性强。
附图说明
图1为航向反馈修正算法的流程示意图;
图2为行人在室内的主方向和主方向范围示意图;
图3为行人直行时的航向角和角速度的变化示意图;
图4为行人转向时的航向角和角速度的变化示意图;
图5为行人走弯道时的航向角和角速度的变化示意图。
具体实施方式
下面结合附图对本发明创造做进一步详细说明。
如图1所示,导航解算过程中,陀螺仪可以实时获得角速率信息,通过四元数解算,实时获得捷联航向角,解算频率一般可设置为50Hz。当行人在主方向范围内直行时,将当前捷联航向角与当前所处的主方向航向角的差值进行自适应无迹卡尔曼滤波,将最优结果与当前捷联航向角之和作为当前航向角的输出,并将其作为下一时刻的航向初始值参与四元和欧拉角解算。
基于微机电系统惯性测量单元(MEMS-IMU)的室内行人导航航向反馈修正算法主要包括以下步骤:
步骤一:四元数算法和欧拉角的捷联解算
设roll,pitch,yaw分别为横滚角,俯仰角和航向角。采用z-x-y的旋转次序,故pitch,roll,yaw为分别为绕着x,y,z轴旋转所得角度。在初始时刻利用加速度值得到俯仰和横滚角分别为
分别表示载体坐标系下x,y,z轴的加速度值。
设置初始航向角yaw=0°。由初始时刻的roll,pitch,yaw信息可以得出初始四元数为Q(t0)=[q0(t0),q1(t0),q2(t0),q3(t0)]T,公式(3)为四元数的微分方程。
其中Q(ti)=[q0(ti),q1(ti),q2(ti),q3(ti)]T为采样点i时刻的四元数,采用一阶毕卡法求解四元数微分方程,得到
其中△t为两相邻采样时刻的时间间隔。ωx,ωy,ωz为分别为载体坐标系相对于惯性坐标系的角速度信息,考虑MEMS-IMU的器件精度,故对其中地球自转角速度和导航系下载体运动引起的角速度均可忽略不计,ωx,ωy,ωz即为陀螺仪所测量的三轴角速度。
四元数更新后,为了便于对角度的直观理解,将四元数转换为欧拉角,其对应的关系为
步骤二:主方向和主方向范围设定:
从开放的卫星地图中可以发现,居民楼和商业大厦等绝大部分建筑都是矩形结构,走廊相互正交或平行,房间的形状也基本为矩形,因此可以将建筑物内的走廊朝向的四个方向定义为建筑物的主方向。如图2中实线所示。考虑到每个行人在直行过程中身体晃动情况不同,根据该行人的晃动情况设定直行晃动角度范围。设置变量range。故当前时刻的航向角与初始时刻航向角(即0°)的差值为-range(即360-range)到range,90-range到90+range,180-range到180+range,270-range到270+range的角度范围内时,即认为处于主方向范围内,如图2中虚线所示。
A1=1表示行人处于四个主方向范围,A1=0表示行人未处于四个主方向范围。
步骤三:行人直行检测
利用加速度滑动方差检测法检测步态:行人在有步态和无步态过程中加速度波动幅度差别较大,可设置一定的滑动窗口求取方差和适当的方差阈值来对行人是否处于行走状态进行区分,公式如下:
其中n为窗口中数据的个数,aj为窗口中第j个时刻的目标加速度,为第i-n+1时刻到i时刻目标加速度的均值。Var为此时窗口数据的方差。
Varmax为阈值,当Varmax=1时,能够区分出行人有无步态。A2=1,A2=0分别代表有步态和无步态。然后利用连续三步之间的角速度和航向角信息对行人直行进行判定。
行人在直行过程中,航向角变化较小,故利用三步间航向差值进行直行判定。
上式中m,n可为k,k-1,k-2,m≠n,yawk为当前步态的航向,yawk-1为前一步时的航向,yawk-2为前两步时的航向。行人在直行时身体会出现正负range的角度晃动,如果三步的任两步的航向绝对差值均小于range时,作为行人直行的判断条件之一。A3=1表示连续三步间的航向差值小于阈值,可作为直行的判定条件,而A3=0代表行人连续步数间航向变化较大,认为行人未处于直行状态。行人在室内的行走情况简单得分为直行,转向,弯道行走三种模式。图3、图4、图5分别为行人在对应情况下的z轴角速度和航向角的变化图。可看出航向角和角速度变化在三种模式下均有不同,故利用航向角和角速度对行人直行进行判定。为了与三步的航向差值判直行相对应,在角速度判直行中也采集连续运动的三步间隔内所记录的角速度值。从图3,4,5中可看出,在转向中角速度值明显大于弯道和直行。故记录在这段时间内角速度绝对值小于阈值ωth的个数占总个数的比例,记做R1,并对该比例设定阈值r1。当小于r1时,认定是转向,A4记为0,大于r1时则认定为弯道行走或者直行。A4记为1。公式如下所示:
在弯道行走中由于行人朝着某一个方向转动,故采集到的角速度值会更多得偏向正或者负的方向。如图5中,角速度分布在正值的较多。而行人在直行过程中角速度在正负的分布情况较为均匀,如图3所示。故在记录A4的同时,分别记录大于零和小于零的角速度的个数占总个数的比例,记做R2和R3。选取R2和R3中较大的比例参与直行的判定,并对此比例设定阈值r2,当大于r2时,A5记为0,当小于r2时,A5记为1。公式下所示:
步骤四:当A2,A3,A4和A5同时为1时才能说明检测到行人直行。
步骤五:当检测到行人在主方向范围内直行,将其进行逻辑“与”运算,记做:
Flag=A1&A2&A3&A4&A5 (12)
如果Flag结果不为1,则进行步骤一。如果A1,A2,A3,A4和A5同时为1时,即Flag为1时,用当前捷联航向角yawi与当前行人所处主方向航向角差值作为状态量,并对其进行自适应无迹卡尔曼滤波,首先建立6维状态向量:
式中,航向角和姿态角误差为δroll,δpitch,δyaw分别表示利用传统和修正两种算法得出的横滚角,俯仰角和航向角差值。角速度误差为δω=[δωx,δωy,δωz],δωx,δωy,δωz分别表示x,y,z轴的角速度误差。其中量测矩阵为H=[0,0,1,0,0,0],经过自适应无迹卡尔曼滤波后将最优估计结果与当前航向角yawi之和作为此刻航向输出值,公式如下:
该输出值参与行人实时的位置解算,并作为下一时刻航向初始值参与四元数和欧拉角解算。从而实现反馈修正的目的。
本发明整个算法流程如下:
1、首先利用传统的四元数和欧拉角算法对航姿角进行解算。
2、将室内环境划分为四个主方向,并根据行人直行的身体晃动划分为四个主方向范围。
3、利用步态识别、三步间航向角和三步间的角速度信息判断行人直行状态;
4、利用行人在主方向范围内直行作为航向反馈修正条件。
5、利用捷联航向角和当前主方向的航向角的差值作为状态量,与横滚角和俯仰角以及三轴角速度信息构成六维状态向量进行自适应无迹卡尔曼滤波处理,将最优结果与当前捷联航向角之和作为此刻的航向角输出并作为下一时刻的初始航向角参与四元数和欧拉角的解算。
Claims (5)
1.基于MEMS-IMU的室内行人导航航向反馈修正方法,其特征在于,包括以下步骤:
(1)通过建筑物结构特征,构建主方向和主方向航向范围;
(2)根据加速度计测得的数据进行步态识别,判断行人是否处于行走状态;
(3)根据相邻步数间航向角和相邻步数间的角速度信息判断行人是否处于直行状态:
(4)当判断行人在主方向范围内直行,用当前航向角yawi与当前行人所处主方向航向角差值作为状态量,进行自适应无迹卡尔曼滤波,获得最优估计结果
(5)将最优估计结果与当前航向角yawi之和作为此时航向输出值参与位置解算,并将该航向输出值作为下一时刻航向角参与四元数和欧拉角的捷联解算。
2.根据权利要求1所述基于MEMS-IMU的室内行人导航航向反馈修正方法,其特征在于:步骤(1)所述构建主方向和主方向航向范围包括:将建筑物内的走廊及过道划分为四个方向,定义为建筑物的主方向;设定行人在主方向上直行的晃动角度范围,建立以主方向为基准的主方向航向角范围。
3.根据权利要求1所述基于MEMS-IMU的室内行人导航航向反馈修正方法,其特征在于:步骤(2)所述判断行人是否处于行走状态的方法为,设置一滑动窗口求取加速度方差,然后与预设的方差阈值比较,判断行人是否处于行走状态。
4.根据权利要求1所述基于MEMS-IMU的室内行人导航航向反馈修正方法,其特征在于:步骤(3)所述相邻步数间航向角和相邻步数间的角速度信息判断行人是否处于直行状态,分别采用三步间航向角差值和三步间隔内记录的角速度值进行判断;所述三步间航向角差值作为判断标准时,任两步的航向角绝对差值均小于预设航向角阈值时,行人直行;所述三步间隔内记录的角速度值作为判断标准时,三步间隔内的角速度绝对值小于预设阈值的个数占总个数的比例小于比例设定阈值时,为转向,大于比例设定阈值时,为弯道行走或者直行。
5.根据权利要求4所述基于MEMS-IMU的室内行人导航航向反馈修正方法,其特征在于:判断所述弯道行走或者直行采用的方法是,若记录的角速度中大于零和小于零的个数分差不多时,为直行,否则为弯道行走。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710725704.9A CN107677267A (zh) | 2017-08-22 | 2017-08-22 | 基于mems‑imu的室内行人导航航向反馈修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710725704.9A CN107677267A (zh) | 2017-08-22 | 2017-08-22 | 基于mems‑imu的室内行人导航航向反馈修正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107677267A true CN107677267A (zh) | 2018-02-09 |
Family
ID=61134627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710725704.9A Pending CN107677267A (zh) | 2017-08-22 | 2017-08-22 | 基于mems‑imu的室内行人导航航向反馈修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107677267A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109855620A (zh) * | 2018-12-26 | 2019-06-07 | 北京壹氢科技有限公司 | 一种基于自回溯算法的室内行人导航方法 |
CN111912406A (zh) * | 2020-08-14 | 2020-11-10 | 航天科工通信技术研究院有限责任公司 | 一种基于改进式hde的室内行人导航航向修正方法及系统 |
CN112747741A (zh) * | 2020-12-07 | 2021-05-04 | 北京汽车研究总院有限公司 | 车辆的惯性导航方法、装置及车辆 |
CN113188542A (zh) * | 2021-06-16 | 2021-07-30 | 东方红卫星移动通信有限公司 | 一种行人主航向更新方法、行人航迹获取方法及系统 |
CN114322918A (zh) * | 2020-09-30 | 2022-04-12 | 广东博智林机器人有限公司 | 可移动设备状态的检测方法、装置和计算机可读存储介质 |
CN116784838A (zh) * | 2023-08-24 | 2023-09-22 | 中国科学技术大学 | 基于可穿戴惯性传感器转向识别系统、方法、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120296603A1 (en) * | 2011-05-16 | 2012-11-22 | Qualcomm Incorporated | Sensor orientation measurement with respect to pedestrian motion direction |
CN104121905A (zh) * | 2014-07-28 | 2014-10-29 | 东南大学 | 一种基于惯性传感器的航向角获取方法 |
CN105043385A (zh) * | 2015-06-05 | 2015-11-11 | 北京信息科技大学 | 一种行人自主导航定位的自适应卡尔曼滤波方法 |
CN106017461A (zh) * | 2016-05-19 | 2016-10-12 | 北京理工大学 | 基于人体/环境约束的行人导航系统三维空间定位方法 |
CN107036597A (zh) * | 2017-05-02 | 2017-08-11 | 华南理工大学 | 一种基于智能手机内置惯性传感器的室内定位导航方法 |
-
2017
- 2017-08-22 CN CN201710725704.9A patent/CN107677267A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120296603A1 (en) * | 2011-05-16 | 2012-11-22 | Qualcomm Incorporated | Sensor orientation measurement with respect to pedestrian motion direction |
CN104121905A (zh) * | 2014-07-28 | 2014-10-29 | 东南大学 | 一种基于惯性传感器的航向角获取方法 |
CN105043385A (zh) * | 2015-06-05 | 2015-11-11 | 北京信息科技大学 | 一种行人自主导航定位的自适应卡尔曼滤波方法 |
CN106017461A (zh) * | 2016-05-19 | 2016-10-12 | 北京理工大学 | 基于人体/环境约束的行人导航系统三维空间定位方法 |
CN107036597A (zh) * | 2017-05-02 | 2017-08-11 | 华南理工大学 | 一种基于智能手机内置惯性传感器的室内定位导航方法 |
Non-Patent Citations (4)
Title |
---|
JIMENEZ, A.R.: "Improved Heuristic Drift Elimination (iHDE) for pedestrian navigation in complex buildings", 《INDOOR POSITIONING AND INDOOR NAVIGATION (IPIN)》 * |
吕沧海 等: "《中远程导弹组合导航技术》", 31 March 2014, 国防工业出版社 * |
楼喜中 等: "基于角度的HDE算法在室内行人航迹定位中的研究", 《传感技术学报》 * |
赵辉 等: "基于主方向的行人自主定位航向修正算法", 《电子技术应用》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109855620A (zh) * | 2018-12-26 | 2019-06-07 | 北京壹氢科技有限公司 | 一种基于自回溯算法的室内行人导航方法 |
CN111912406A (zh) * | 2020-08-14 | 2020-11-10 | 航天科工通信技术研究院有限责任公司 | 一种基于改进式hde的室内行人导航航向修正方法及系统 |
CN114322918A (zh) * | 2020-09-30 | 2022-04-12 | 广东博智林机器人有限公司 | 可移动设备状态的检测方法、装置和计算机可读存储介质 |
CN114322918B (zh) * | 2020-09-30 | 2024-02-13 | 广东博智林机器人有限公司 | 可移动设备状态的检测方法、装置和计算机可读存储介质 |
CN112747741A (zh) * | 2020-12-07 | 2021-05-04 | 北京汽车研究总院有限公司 | 车辆的惯性导航方法、装置及车辆 |
CN113188542A (zh) * | 2021-06-16 | 2021-07-30 | 东方红卫星移动通信有限公司 | 一种行人主航向更新方法、行人航迹获取方法及系统 |
CN116784838A (zh) * | 2023-08-24 | 2023-09-22 | 中国科学技术大学 | 基于可穿戴惯性传感器转向识别系统、方法、设备及介质 |
CN116784838B (zh) * | 2023-08-24 | 2024-01-09 | 中国科学技术大学 | 基于可穿戴惯性传感器转向识别系统、方法、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107677267A (zh) | 基于mems‑imu的室内行人导航航向反馈修正方法 | |
CN109991636A (zh) | 基于gps、imu以及双目视觉的地图构建方法及系统 | |
CN105509739B (zh) | 采用固定区间crts平滑的ins/uwb紧组合导航系统及方法 | |
CN102519450B (zh) | 一种用于水下滑翔器的组合导航装置及方法 | |
CN106643715A (zh) | 一种基于bp神经网络改善的室内惯性导航方法 | |
CN106168485B (zh) | 步行航迹数据推算方法及装置 | |
CN103791902B (zh) | 适用于高机动载体的星敏感器自主导航方法 | |
CN104613965B (zh) | 一种基于双向滤波平滑技术的步进式行人导航方法 | |
CN110398245A (zh) | 基于脚戴式惯性测量单元的室内行人导航姿态估计方法 | |
CN109459028A (zh) | 一种基于梯度下降的自适应步长估计方法 | |
CN110986941B (zh) | 一种手机安装角的估计方法 | |
CN109959374B (zh) | 一种行人惯性导航全时全程逆向平滑滤波方法 | |
CN102636166A (zh) | 基于航向角的wsn/ins组合导航系统及方法 | |
Xu et al. | Predictive adaptive Kalman filter and its application to INS/UWB-integrated human localization with missing UWB-based measurements | |
CN112066980B (zh) | 一种基于人体四节点运动约束的行人导航定位方法 | |
Huang et al. | Smartphone-based indoor position and orientation tracking fusing inertial and magnetic sensing | |
CN104406592B (zh) | 一种用于水下滑翔器的导航系统及姿态角校正和回溯解耦方法 | |
CN104296741B (zh) | 采用距离平方和距离平方变化率的wsn/ahrs紧组合方法 | |
CN112362057A (zh) | 基于零速修正与姿态自观测的惯性行人导航算法 | |
CN108592907A (zh) | 一种基于双向滤波平滑技术的准实时步进式行人导航方法 | |
Qian et al. | RPNOS: Reliable pedestrian navigation on a smartphone | |
CN105115507B (zh) | 一种基于双imu的双模式室内个人导航系统及方法 | |
CN109752023A (zh) | 一种目标运动状态快速估计方法 | |
Li et al. | Real-time attitude tracking of mobile devices | |
CN104897157B (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: 20180209 |