CN105928512A - 一种基于地磁场的室内定位方法 - Google Patents

一种基于地磁场的室内定位方法 Download PDF

Info

Publication number
CN105928512A
CN105928512A CN201610270482.1A CN201610270482A CN105928512A CN 105928512 A CN105928512 A CN 105928512A CN 201610270482 A CN201610270482 A CN 201610270482A CN 105928512 A CN105928512 A CN 105928512A
Authority
CN
China
Prior art keywords
magnetic field
point
indoor
inertial navigation
navigation system
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
Application number
CN201610270482.1A
Other languages
English (en)
Inventor
楼喜中
方俊
孙宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Xin Shengda Information Technology Co Ltd
Original Assignee
Hangzhou Xin Shengda Information Technology Co Ltd
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 Hangzhou Xin Shengda Information Technology Co Ltd filed Critical Hangzhou Xin Shengda Information Technology Co Ltd
Priority to CN201610270482.1A priority Critical patent/CN105928512A/zh
Publication of CN105928512A publication Critical patent/CN105928512A/zh
Pending legal-status Critical Current

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/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • 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

本发明公开了一种基于地磁场的室内定位方法,包括以下步骤:S1,建立地磁场数据库;S2,通过惯性导航系统输出运动载体行走过程中的磁场强度有限点集合A={a1,a2,a3…,ap}和运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)},然后根据每个相对位置与地磁数据库得出新地磁数据点集合B={b1,b2,b3,…,bp},惯性导航系统将磁场强度有限点集合A={a 1,a2,a3…,ap}与新地磁数据点集合B={b1,b2,b3,…,bp}进行计算出运动载体的整个轨迹豪斯多夫距离dH(HN,H’N);S3,惯性导航系统会根据豪斯多夫距离得出室内轨迹,最终通过惯性导航系统会通过室内轨迹遍历整个地磁场数据库得出其当前室内具体位置。本发明使得能够进行室内定位,定位的精确度更高,更加能够满足现代人们的需求。

Description

一种基于地磁场的室内定位方法
技术领域
本发明涉及地理导航技术领域,尤其涉及一种基于地磁场的室内定位方法。
背景技术
在如今的军事民生方面,导航技术俨然已经成为一种非常重要的技术,作为当前非常热门的技术之一,目前导航的手段较多,一般有无线电导航、天文导航、惯性导航和卫星导航,而现有的导航技术都存在缺陷:
无线电导航—无线电导航信号包括WLAN、GPS导航和雷达导航等,虽然克服了误差累积的问题,且具有精度高,但是布置成本高,同时在室内容易受到建筑物内部构造造成信号干扰,从而无法适合在室内进行导航。
天文导航和卫星导航—天文导航主要是利用对自然天体的测量来确定自身位置和航向的导航技术;卫星导航是指采用导航卫星对地面、海洋、空中和空间用户进行导航定位的技术,常见的GPS导航,北斗星导航等均为卫星导航,其天文导航和卫星导航主要是应用在航天器上,而在室内中其成本较高,导航容易受到建筑物的内部构造的影响,从而使得导航精度不够高。
惯性导航—惯性器件导航利用惯性元件(加速度计)来测量运载体本身的加速度,经过积分和运算得到速度和位置,从而达到对运载体导航定位的目的,惯性器件导航是目前应用最广泛的技术,它具有自主性强,抵抗外界干扰能力强等特点,但是会随着时间增长产生累积误差的问题,其中惯性导航中应用较多的是地磁匹配的惯性导航,而地磁匹配的惯性导航一般用于导弹和飞机等载体在室外且广阔的范围内进行匹配计算,但是在室内进行导航很容易受到建筑物构造和金属物等障碍物会对磁场进行干扰,使得室内导航精度得不到保障,且现有的地磁数据库很难适应室内导航的需求,其花费的成本较高,无法适用现在发展需求。
总而言之,现有的导航技术一般都是使用在室外,而将这些室外导航技术使用在室内其精度不是很高,所以不能够准确的进行定位,从而无法帮助人们进行定位导航,无法达到使用需求。
发明内容
针对现有技术中的不足,本发明的提供了一种基于地磁场的室内定位方法。
为了实现上述目的,本发明所采取的技术方案是:
一种基于地磁场的室内定位方法,包括以下步骤:
步骤1,建立地磁场数据库,利用校准完成的惯性导航系统收集所要室内的数据,所述建立地磁场数据库具体步骤为:a,惯性导航系统中的高精度磁传感器在室内划分好的网格点交点处收集磁场数据;b,惯性导航系统在室内划分好的网格点交点处停留一段时间S,其停留一段时间S优先为5秒,再将惯性导航系统中的磁场强度值截取中间若干个,若干个的磁场强度值组成一个集合M,从而计算出集合M的均值Ma=M/Ms,所述再将惯性导航系统中的磁场强度值截取中间若干个为5个,若干个的磁场强度值组成一个集合M,5个磁场强度值组成一个集合M具体表示为M={m1,m2,m3,m4,m5},其集合M的均值根据Ma=M/Ms可以得到Ma=M/5,若干个的磁场强度值组成一个集合M,Ms表示集合M的个数;c,惯性导航系统将得到集合M的均值Ma依次录入到磁场数据中的网格点交点中,从而完成自建立网格型的地磁场数据库,其地磁场数据库的单位规格是dx×dy,x和y分别是地磁数据库中网格的横坐标和纵坐标,dx表示数据库中的横轴长度,dy表示数据库中的纵轴长度,其中地磁场数据库中的网格的宽度是0.6米的间隔;
步骤2,通过惯性导航系统输出运动载体行走过程中的磁场强度有限点集合A={a1,a2,a3…,ap}和运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)},磁场强度有限点集合A={a1,a2,a3…,ap}是通过惯性导航系统中的磁传感器输出的,且输出位置{X,Y}集合中每个点的数据均是相对于初始点E(0,0)位置,然后将初始点的位置设置为W(x,y),根据初始点位置W(x,y)计算出运动载体行走轨迹位置与初始点位置W(x,y)的相对位置,然后根据每个相对位置与地磁数据库得出新地磁数据点集合B={b1,b2,b3,…,bp},惯性导航系统将磁场强度有限点集合A={a1,a2,a3…,ap}与新地磁数据点集合B={b1,b2,b3,…,bp}进行计算,然后得出运动载体的整个轨迹豪斯多夫距离dH(HN,H’N),例如其中匹配区域中新磁场有限点集合B=(b1,b2,b3,…,bp)中是由运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)}根据设置的初始点W(x,y)并按照公式{(X/0.6+x),(Y/0.6+y)}得到当前行走轨迹的网格点坐标,然后根据坐标在自建的地磁数据库中找出对应网格的磁场数据,由于每个网格都有4个数据,所以对这4个数据取平均值,最后输出得到点集B,有限点集合A={a1,a2,a3…,ap}中的数据都是载体在行走过程中,惯性导航系统自动记录下的经过的区域的磁场强度,记录频率可以根据需要制定,从而惯性导航系统计算出运动载体的整个轨迹豪斯多夫dH(HN,H’N),其中每个轨迹点所在地磁场数据库的网格坐标中;
步骤3,惯性导航系统会根据豪斯多夫距离得出室内轨迹,其公式为:dH(HN,H’N)=mindH(HN,U)(H’N∈U),也就是U表示当前匹配区域,HN表示点集合A,H’N表示室内最佳匹配位置点的集合,最终通过惯性导航系统会通过室内轨迹遍历整个地磁场数据库得出其当前室内具体位置。
本发明的有益效果:与现有技术相比,本发明具有克服现有室内定位方法中地磁场数据库中花费成本较高,使得能够进行室内定位,更加符合室内复杂的状况,也有利于避免了传统地磁匹配惯性定位中多点匹配中,点对点关系必须明确的缺点,使得精度更加准确,更加适合现代发展需求。
附图说明
图1为本发明的方法流程图;
图2为本发明的其中一处楼层地磁数据库模型。
具体实施方式
本发明结合图1和图2作出如下具体的说明。
一种基于地磁场的室内定位方法,包括以下步骤:
步骤1,建立地磁场数据库,利用校准完成的惯性导航系统收集所要室内的数据,所述建立地磁场数据库具体步骤为:a,惯性导航系统中的高精度磁传感器在室内划分好的网格点交点处收集磁场数据;b,惯性导航系统在室内划分好的网格点交点处停留一段时间S,其停留一段时间S优先为5秒,再将惯性导航系统中的磁场强度值截取中间若干个,若干个的磁场强度值组成一个集合M,从而计算出集合M的均值Ma=M/Ms,所述再将惯性导航系统中的磁场强度值截取中间若干个为5个,若干个的磁场强度值组成一个集合M,5个磁场强度值组成一个集合M具体表示为M={m1,m2,m3,m4,m5},其集合M的均值根据Ma=M/Ms可以得到Ma=M/5,若干个的磁场强度值组成一个集合M,Ms表示集合M的个数;c,惯性导航系统将得到集合M的均值Ma依次录入到磁场数据中的网格点交点中,从而完成自建立网格型的地磁场数据库,其地磁场数据库的单位规格是dx×dy,x和y分别是地磁数据库中网格的横坐标和纵坐标,dx表示数据库中的横轴长度,dy表示数据库中的纵轴长度,其中地磁场数据库中的网格的宽度是0.6米的间隔;
步骤2,通过惯性导航系统输出运动载体行走过程中的磁场强度有限点集合A={a1,a2,a3…,ap}和运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)},且输出位置{X,Y}集合中每个点的数据均是相对于初始点E(0,0)位置,然后将初始点的位置设置为W(x,y),根据初始点位置W(x,y)和相对位置有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)}计算出运动载体行走轨迹位置,然后根据每条轨迹位置与地磁数据库得出每条轨迹对应的新地磁数据点集合B={b1,b2,b3,…,bp},惯性导航系统将磁场强度有限点集合A={a1,a2,a3…,ap}与新地磁数据点集合B={b1,b2,b3,…,bp}进行计算,然后得出运动载体的整个轨迹豪斯多夫距离dH(HN,H’N),例如其中匹配区域中新磁场有限点集合B=(b1,b2,b3,…,bp)中是由运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)}根据设置的初始点W(x,y)并按照公式{(X/0.6+x),(Y/0.6+y)}得到当前行走轨迹的网格点坐标,然后根据坐标在自建的地磁数据库中找出对应网格的磁场数据,由于每个网格都有4个数据,所以对这4个数据取平均值,最后输出得到点集B,有限点集合A={a1,a2,a3…,ap}中的数据都是载体在行走过程中,惯性导航系统自动记录下的经过的区域的磁场强度,记录频率可以根据需要制定,从而惯性导航系统计算出运动载体的整个轨迹dH(HN,H’N),其中每个轨迹点所在地磁场数据库的网格坐标中;
步骤3,惯性导航系统会根据豪斯多夫距离得出室内轨迹,其公式为:dH(HN,H’N)=mindH(HN,U)(H’N∈U),其中U表示点集合B,也就是U表示当前匹配区域,HN表示点集合A,H’N表示室内最佳匹配位置点的集合,最终通过惯性导航系统会通过室内轨迹遍历整个地磁场数据库得出其当前室内具体位置。
本发明具体举例实现对室内行走轨迹定位导航具体实施如下:在第一次仿真中,测试区域即为图1所示的室内,从室内建立的坐标系(3,1)出发沿室内走廊(东西走向)向东走,根据地磁场数据库取20个数据点长度的轨迹进行最终仿真,遍历整个建筑物楼层后输出距离测度:f6=1445.9,f5=2107.25,f4=1113.25,f3=1524.55,f2=2241.75,f1=1758.4,同时也对30、40、50个数据点轨迹仿真,其中根据地磁场数据库50个数据点长的轨迹仿真后:f6=4824.05,f5=13909.3,f4=16268.8,f3=8872.4,f2=11359,f1=9792.5,50个点能更好的区分楼层,同时轨迹定位结果起始点输出为(3,2),上述仿真每组进行10次,共50次实验,可以发现采用50个点的匹配长度的序列就能够够输出误差范围内的定位位置,实验表明本发明方法能够根据自建的地磁场数据库和惯性输出模块输出信息,通过调整匹配数据序列的长度,可得到有效的匹配定位结果。
本发明的有益效果:与现有技术相比,本发明具有克服现有室内定位方法中地磁场数据库中花费成本较高,使得能够进行室内定位,更加符合室内复杂的状况,也有利于避免了传统地磁匹配惯性定位中多点匹配中,点对点关系必须明确的缺点,使得精度更加准确,更加适合现代发展需求。

Claims (4)

1.一种基于地磁场的室内定位方法,其特征在于,包括以下步骤:
步骤1,建立地磁场数据库,利用校准完成的惯性导航系统收集所要室内的数据;
步骤2,通过惯性导航系统输出运动载体行走过程中的磁场强度有限点集合A={a1,a2,a3…,ap}和运动载体行走轨迹位置的有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)},且输出位置{X,Y}集合中每个点的数据均是相对于初始点E(0,0)位置,然后将初始点的位置设置为W(x,y),根据初始点位置W(x,y)和相对位置有限点集合{X,Y}={(x1,y1),(x2,y2),(x3,y3),…,(xp,yp)}计算出运动载体行走轨迹位置,然后根据每条轨迹位置与地磁数据库得出每条轨迹对应的新地磁数据点集合B={b1,b2,b3,…,bp},惯性导航系统将磁场强度有限点集合A={a1,a2,a3…,ap}与新地磁数据点集合B={b1,b2,b3,…,bp}进行豪斯多夫距离计算,得出豪斯多夫距离;
步骤3,惯性导航系统会根据豪斯多夫距离得出室内轨迹,其公式为:dH(HN,H’N)=mindH(HN,U)(H’N∈U),其中U表示当前匹配区域,HN表示点集合A,H’N表示室内最佳匹配位置点的集合,最终通过惯性导航系统会通过室内轨迹遍历整个地磁场数据库得出其当前室内具体位置。
2.如权利要求1所述的一种基于地磁场的室内定位方法,其特征在于,所述建立地磁场数据库具体步骤为:a,惯性导航系统中的高精度磁传感器在室内划分好的网格点交点处收集磁场数据;b,惯性导航系统在室内划分好的网格点交点处停留一段时间S,再将惯性导航系统中的磁场强度值截取中间若干个,若干个的磁场强度值组成一个集合M,从而计算出集合M的均值Ma=M/Ms,若干个的磁场强度值组成一个集合M,Ms表示集合M的个数;c,惯性导航系统将得到集合M的均值Ma依次录入到磁场数据中的网格点交点中,从而完成自建立网格型的地磁场数据库,其地磁场数据库的单位规格是dx×dy,x和y分别是地磁数据库中网格的横坐标和纵坐标,dx表示数据库中的横轴长度,dy表示数据库中的纵轴长度。
3.如权利要求2所述的一种基于地磁场的室内定位方法,其特征在于,所述惯性导航系统在室内划分好的网格点交点处停留一段时间S优先为5秒。
4.如权利要求2所述的一种基于地磁场的室内定位方法,其特征在于,所述再将惯性导航系统中的磁场强度值截取中间若干个为5个,若干个的磁场强度值组成一个集合M,5个磁场强度值组成一个集合M具体表示为M={m1,m2,m3,m4,m5},其集合M的均值根据Ma=M/Ms可以得到Ma=M/5。
CN201610270482.1A 2016-04-26 2016-04-26 一种基于地磁场的室内定位方法 Pending CN105928512A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610270482.1A CN105928512A (zh) 2016-04-26 2016-04-26 一种基于地磁场的室内定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610270482.1A CN105928512A (zh) 2016-04-26 2016-04-26 一种基于地磁场的室内定位方法

Publications (1)

Publication Number Publication Date
CN105928512A true CN105928512A (zh) 2016-09-07

Family

ID=56837357

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610270482.1A Pending CN105928512A (zh) 2016-04-26 2016-04-26 一种基于地磁场的室内定位方法

Country Status (1)

Country Link
CN (1) CN105928512A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107560616A (zh) * 2017-08-17 2018-01-09 深圳市元智慧科技有限公司 磁力导航方法及其磁力导航仪
CN107621263A (zh) * 2017-08-25 2018-01-23 北京麦钉艾特科技有限公司 一种基于道路磁场特征的地磁定位方法
CN110057355A (zh) * 2019-04-18 2019-07-26 吉林大学 一种室内定位方法、装置、系统及计算设备
CN110296699A (zh) * 2019-06-06 2019-10-01 深圳数位传媒科技有限公司 一种基于地磁场的室内定位方法及装置
CN110967013A (zh) * 2018-10-01 2020-04-07 刘佳兴 一种基于室内地磁场信息和智能手机的室内区域定位系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102829782A (zh) * 2012-09-07 2012-12-19 滨州学院 一种地磁辅助惯性导航方法
CN103471589A (zh) * 2013-09-25 2013-12-25 武汉大学 一种室内行人行走模式识别和轨迹追踪的方法
CN104019813A (zh) * 2014-06-19 2014-09-03 无锡知谷网络科技有限公司 目标即时定位和构建地图的方法与系统
CN104105067A (zh) * 2013-04-09 2014-10-15 广东美晨通讯有限公司 基于地磁信息的实时定位系统及方法
CN104197934A (zh) * 2014-09-02 2014-12-10 百度在线网络技术(北京)有限公司 一种基于地磁的定位方法、装置及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102829782A (zh) * 2012-09-07 2012-12-19 滨州学院 一种地磁辅助惯性导航方法
CN104105067A (zh) * 2013-04-09 2014-10-15 广东美晨通讯有限公司 基于地磁信息的实时定位系统及方法
CN103471589A (zh) * 2013-09-25 2013-12-25 武汉大学 一种室内行人行走模式识别和轨迹追踪的方法
CN104019813A (zh) * 2014-06-19 2014-09-03 无锡知谷网络科技有限公司 目标即时定位和构建地图的方法与系统
CN104197934A (zh) * 2014-09-02 2014-12-10 百度在线网络技术(北京)有限公司 一种基于地磁的定位方法、装置及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姜浩等: "利用Hausdorff距离地磁匹配算法的室内定位方法", 《科技传播》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107560616A (zh) * 2017-08-17 2018-01-09 深圳市元智慧科技有限公司 磁力导航方法及其磁力导航仪
CN107621263A (zh) * 2017-08-25 2018-01-23 北京麦钉艾特科技有限公司 一种基于道路磁场特征的地磁定位方法
CN110967013A (zh) * 2018-10-01 2020-04-07 刘佳兴 一种基于室内地磁场信息和智能手机的室内区域定位系统
CN110057355A (zh) * 2019-04-18 2019-07-26 吉林大学 一种室内定位方法、装置、系统及计算设备
CN110057355B (zh) * 2019-04-18 2021-06-08 吉林大学 一种室内定位方法、装置、系统及计算设备
CN110296699A (zh) * 2019-06-06 2019-10-01 深圳数位传媒科技有限公司 一种基于地磁场的室内定位方法及装置

Similar Documents

Publication Publication Date Title
CN107734457B (zh) 智慧停车场导航系统及方法
Zhou et al. Activity sequence-based indoor pedestrian localization using smartphones
CN105547305B (zh) 一种基于无线定位和激光地图匹配的位姿解算方法
CN105928512A (zh) 一种基于地磁场的室内定位方法
CN105043415B (zh) 基于四元数模型的惯性系自对准方法
CN107504971A (zh) 一种基于pdr和地磁的室内定位方法及系统
CN103791902B (zh) 适用于高机动载体的星敏感器自主导航方法
CN103363988A (zh) 一种利用智能手机传感器实现地磁室内定位导航的方法
CN104897154A (zh) 基于手机的室内地磁导航方法及装置
CN105704652A (zh) 一种wlan/蓝牙定位中的指纹库采集和优化方法
CN110617795B (zh) 一种利用智能终端的传感器实现室外高程测量的方法
CN104197930A (zh) 一种基于惯性制导和射频识别的室内定位装置及方法
CN108775901A (zh) 一种实时slam场景地图构建系统、导航系统及方法
CN106950976B (zh) 基于卡尔曼和粒子滤波的室内飞艇三维定位装置及方法
Rehman et al. Accurate and direct GNSS/PDR integration using extended Kalman filter for pedestrian smartphone navigation
CN111121758B (zh) 室内磁地图的快速建模与可信定位方法
CN112923919A (zh) 基于图优化的行人定位方法及系统
CN106197418B (zh) 一种基于滑动窗口的指纹法与传感器融合的室内定位方法
CN105973224B (zh) 一种室内导航方法
CN107631727A (zh) 一种室内css/ins组合导航系统
Zhuang et al. Autonomous WLAN heading and position for smartphones
CN106971601A (zh) 一种基于WiFi的智能泊车及找回的系统与方法
Kang et al. Smartphone indoor positioning system based on geomagnetic field
Beck et al. Anchor free node tracking using ranges, odometry, and multidimensional scaling
CN115406435A (zh) 基于wlan和mems的室内电子地图构建方法、装置及存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20160907

RJ01 Rejection of invention patent application after publication