CN112611377A - 一种用于小车室外导航的状态预测方法、装置及存储介质 - Google Patents
一种用于小车室外导航的状态预测方法、装置及存储介质 Download PDFInfo
- Publication number
- CN112611377A CN112611377A CN202011379890.3A CN202011379890A CN112611377A CN 112611377 A CN112611377 A CN 112611377A CN 202011379890 A CN202011379890 A CN 202011379890A CN 112611377 A CN112611377 A CN 112611377A
- Authority
- CN
- China
- Prior art keywords
- trolley
- coordinate
- gps signal
- prediction
- state
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000013598 vector Substances 0.000 claims abstract description 45
- 230000008859 change Effects 0.000 claims abstract description 38
- 238000001914 filtration Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
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
- G01C21/165—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 combined with non-inertial navigation instruments
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种用于小车室外导航的状态预测方法、装置及计算机可读存储介质,所述小车室外导航的状态预测方法包括以下步骤:获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量。本发明所述的用于小车室外导航的状态预测方法,降低了室外导航中状态预测的时间复杂度。
Description
技术领域
本发明涉及小车室外导航技术领域,尤其涉及一种用于小车室外导航的状态预测方法、装置及计算机可读存储介质。
背景技术
移动小车室外定位导航一直是智能交通领域的关注热点,GPS是目前主流的实现小车室外定位的依赖工具,因为此方法有非常突出的优点,即位置信息直接在绝对帧中给出,而且不需要增加额外的设备。
GPS接收机在实际应用中存在较大定位误差,不能满足室外导航的定位精度需求,此外,在导航过程中可能存在信号接收机被树木或者建筑物等障碍物遮挡的情况,导致GPS信号强度降低甚至信号中断,使得小车无法实现定位,影响导航的精度。
现有的室外导航定位算法大都采用卡尔曼滤波算法进行传感器信息融合,主要考虑如何提高滤波算法的精准性,而缺乏对算法的时间复杂度进行分析;一般的,在树木或者建筑物等障碍物遮挡的环境下,经常发生GPS信号强度低甚至中断的情况,导致GPS定位数据无法及时更新,考虑如何保证室外导航的精准性,增加其环境适应性是必要的,但是导航的时间消耗同样值得关注,减少状态预测的计算量,降低状态预测的时间复杂度可以使室外导航任务更加高效的完成。
发明内容
有鉴于此,有必要提供一种用于小车室外导航的状态预测方法、装置及计算机可读存储介质,用以解决现有技术中室外导航中状态预测的时间复杂度较大的问题。
本发明提供了一种用于小车室外导航的状态预测方法,包括以下步骤:
获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;
确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;
根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量。
进一步地,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型,具体包括,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型
进一步地,根据所述预测模型及观测向量判断小车的位置变化是否超过设定阈值,具体包括,利用小车的位置变化判定公式确定小车的位置变化是否超过设定阈值,所述小车的位置变化判定公式为
当λk=1时,小车的位置变化超过设定阈值,当λk=0时,小车的位置变化未超过设定阈值,Zk、Zk-1为k、k-1时刻的观测向量,ε、σ分别为第一设定阈值、第二设定阈值,ξ=[1 11 0 0],为滤波估计值。
进一步地,根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标,具体包括,若信号接收机获取的GPS信号值准确,则以GPS信号值确定小车当前的经度坐标和纬度坐标,若信号接收机获取的GPS信号值不准确,则以GPS信号估计模型确定小车当前的经度坐标和纬度坐标,所述GPS信号估计模型为
其中,μ=[0 0 0 0 1 0],ω=[0 0 1 0 0 0],ET、NT分别为小车当前的经度坐标和纬度坐标。
进一步地,所述用于小车室外导航的状态预测方法还包括,根据GPS信号值判断公式判断信号接收机获取的GPS信号值是否准确,所述GPS信号值判断公式为
其中,α=[1 0 0 0 0]、β=[0 1 0 0 0]、θ=[0 0 1 0 0],当rk=1时,则信号接收机获取的GPS信号值不准确,当rk=0,则信号接收机获取的GPS信号值准确。
进一步地,根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量,具体包括,根据所述小车当前的经度坐标、纬度坐标及预测模型获取状态估计方程模型及观测方程模型,进行滤波预测及滤波矫正,获取小车的的预测状态向量。
进一步地,所述状态估计方程模型及观测方程模型分为
本发明还提供了一种用于小车室外导航的状态预测装置,包括处理器以及存储器,所述存储器上存储有计算机程序,所述计算机程序被所述处理器执行时,实现如上述任一技术方案所述的用于小车室外导航的状态预测方法。
本发明还提供了一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机该程序被处理器执行时,实现如上述任一技术方案所述的用于小车室外导航的状态预测方法。
与现有技术相比,本发明的有益效果包括:通过获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量;降低了室外导航中状态预测的时间复杂度。
附图说明
图1为本发明提供的用于小车室外导航的状态预测方法的流程示意图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
实施例1
本发明实施例提供了一种用于小车室外导航的状态预测方法,其流程示意图,如图1所示,所述方法以下步骤:
S1、获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;
S2、确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;
S3、根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量。
优选的,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型,具体包括,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型
需要说明的是,描述了小车在地球坐标系下的位置变化情况;具体实施例时,小车的状态向量为 和分别为经度坐标、纬度坐标和航向角的估计值,表示小车的前向移动速度的估计值,为偏航率的估计值,为陀螺仪偏差的估计值,为航向角的估计值,为里程计偏差的估计值;
优选的,根据所述预测模型及观测向量判断小车的位置变化是否超过设定阈值,具体包括,利用小车的位置变化判定公式确定小车的位置变化是否超过设定阈值,所述小车的位置变化判定公式为
当λk=1时,小车的位置变化超过设定阈值,当λk=0时,小车的位置变化未超过设定阈值,Zk、Zk-1为k、k-1时刻的观测向量,ε、σ分别为第一设定阈值、第二设定阈值,ξ=[1 11 0 0],为滤波估计值;
一个优选实施例中,ε、σ分别为4和8,当时,小车的位置变化超过设定阈值,再根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;
优选的,根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标,具体包括,若信号接收机获取的GPS信号值准确,则以GPS信号值确定小车当前的经度坐标和纬度坐标,若信号接收机获取的GPS信号值不准确,则以GPS信号估计模型确定小车当前的经度坐标和纬度坐标,所述GPS信号估计模型为
其中,μ=[0 0 0 0 1 0],ω=[0 0 1 0 0 0],ET、NT分别为小车当前的经度坐标和纬度坐标。
优选的,所述用于小车室外导航的状态预测方法还包括,根据GPS信号值判断公式判断信号接收机获取的GPS信号值是否准确,所述GPS信号值判断公式为
其中,α=[1 0 0 0 0]、β=[0 1 0 0 0]、θ=[0 0 1 0 0],当rk=1时,则信号接收机获取的GPS信号值不准确,当rk=0,则信号接收机获取的GPS信号值准确;
具体实施时,当λk=1、rk=1成立时,则传感器向KF(卡尔曼滤波)估计器发送信息,使KF估计器工作,KF估计器根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量。
优选的,根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量,具体包括,根据所述小车当前的经度坐标、纬度坐标及预测模型获取状态估计方程模型及观测方程模型,进行滤波预测及滤波矫正,获取小车的的预测状态向量。
优选的,所述状态估计方程模型及观测方程模型分为
一个具体实施例中,通过预测模型对状态向量以及误差协方差矩阵进行滤波预测,滤波预测过程表示为
P′k=APk-1AT+Qω
将传感器测量信息通过滤波算法与预测信息融合,计算出卡尔曼增益和最终的状态估计信息,即滤波矫正,滤波矫正过程为
Kk=P′kHT(HP′kHT+(1-γk)Rυ+(1-λk)σ2I)-1
其中,A、H、Qω、Xk、Rυ、Pk、Kk、分别为状态转移矩阵、测量矩阵、过程估计误差协方差矩、状态向量真实值、测量误差协方差矩阵、状态向量卡尔曼滤波估计值、卡尔曼误差协方差矩阵、状态向量预测值、预测误差协方差矩阵、卡尔曼增益、测量余量;状态向量预测值即为小车的的预测状态向量;
需要说明的是,Xk为k时刻的小车状态信息,Zk为由车载传感器测得的小车状态信息,为卡尔曼滤波估计值,引入λk确保估计状态信息和测量状态信息只有在小车的状态变化超过设定阈值时才被传递到KF估计器,利用γk及时的使能GPS信号估计器,使KF滤波过程不受GPS信号中断的影响,λk和γk代表算法中的触发事件,是基于事件卡尔曼滤波算法的关键;最后,将估计信息和测量信息传到KF估计器进行卡尔曼滤波迭代计算,给出小车的最新状态信息估计向量
一个具体实施例中,以导航的精准度和整体时间消耗作为评价指标,在进行室外导航之前要为小车提供目标点的经、纬度坐标以及航向角信息,导航开始之后记录整个导航过程的耗时,作为本发明实施例所述用于小车室外导航的状态预测方法的时间评价指标,构建位置精准度评价模型如下,
eψ=|ψgoal-ψ|
C=sin(NT goal)sin(NT)cos(ET goal-ET)+cos(NT goal)cos(NT)
其中,所述ep和eψ分别代表位置误差和航向角误差,其中ψgoal、NT goal和ET goal分别代表目标点的航向角、纬度坐标和经度坐标,相应的ψ、NT和ET表示当前小车的航向角、纬度坐标和经度坐标。
实施例2
本发明实施例提供了一种用于小车室外导航的状态预测装置,包括处理器以及存储器,所述存储器上存储有计算机程序,所述计算机程序被所述处理器执行时,实现如实施例1所述的用于小车室外导航的状态预测方法。
实施例3
本发明提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机该程序被处理器执行时,实现如实施例1所述的用于小车室外导航的状态预测方法。
本发明所述的用于小车室外导航的状态预测方法、装置及计算机可读存储介质,通过获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量;降低了室外导航中状态预测的时间复杂度;
本发明所述技术方案将事件触发机制加入KF算法,提供了模拟GPS信号的GPS信号估计器,将小车的状态变化作为KF的触发事件,当状态偏差超过一定阈值时,KF估计器响应并给出新的状态估计,即使GPS信号中断发生,GPS信号估计器为KF估计器提供连续可用的位置信息使滤波过程不受影响;从而在保证室外导航的精度的同时,降低算法的时间复杂度,减少完成导航任务所需的时间。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种用于小车室外导航的状态预测方法,其特征在于,包括以下步骤:
获取小车的航向角、转向角及前向移动速度和横向侧滑速度,根据航向角、转向角及前向移动速度和横向侧滑速度,建立描述小车在地球坐标系下的位置变化的预测模型;
确定小车的观测向量,根据观测向量判断小车的位置变化是否超过设定阈值,若是,则根据所述信号接收机获取的GPS信号值是否准确,确定小车当前的经度坐标和纬度坐标;
根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量。
7.根据权利要求6所述的用于小车室外导航的状态预测方法,其特征在于,根据所述小车当前的经度坐标、纬度坐标及预测模型,获取小车的的预测状态向量,具体包括,根据所述小车当前的经度坐标、纬度坐标及预测模型获取状态估计方程模型及观测方程模型,进行滤波预测及滤波矫正,获取小车的的预测状态向量。
9.一种用于小车室外导航的状态预测装置,其特征在于,包括处理器以及存储器,所述存储器上存储有计算机程序,所述计算机程序被所述处理器执行时,实现如权利要求1-8任一所述的用于小车室外导航的状态预测方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机该程序被处理器执行时,实现如权利要求1-8任一所述的用于小车室外导航的状态预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011379890.3A CN112611377B (zh) | 2020-11-30 | 2020-11-30 | 一种用于小车室外导航的状态预测方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011379890.3A CN112611377B (zh) | 2020-11-30 | 2020-11-30 | 一种用于小车室外导航的状态预测方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112611377A true CN112611377A (zh) | 2021-04-06 |
CN112611377B CN112611377B (zh) | 2024-03-12 |
Family
ID=75229756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011379890.3A Active CN112611377B (zh) | 2020-11-30 | 2020-11-30 | 一种用于小车室外导航的状态预测方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112611377B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113596104A (zh) * | 2021-07-05 | 2021-11-02 | 武汉理工大学 | 一种移动机器人和云服务器的数据传输方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108871336B (zh) * | 2018-06-20 | 2019-05-07 | 湘潭大学 | 一种车辆位置估算系统及方法 |
CN111458733B (zh) * | 2020-01-22 | 2022-05-24 | 武汉光庭科技有限公司 | 一种结合gps定位与车身信息的位姿校正方法及装置 |
CN111653122A (zh) * | 2020-05-06 | 2020-09-11 | 南京航空航天大学 | 一种车辆协同碰撞预警系统及其控制方法 |
CN111949030B (zh) * | 2020-08-17 | 2023-01-10 | 江苏常发农业装备股份有限公司 | 一种农机定位方法、农机车辆和存储介质 |
-
2020
- 2020-11-30 CN CN202011379890.3A patent/CN112611377B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113596104A (zh) * | 2021-07-05 | 2021-11-02 | 武汉理工大学 | 一种移动机器人和云服务器的数据传输方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112611377B (zh) | 2024-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106840179B (zh) | 一种基于多传感器信息融合的智能车定位方法 | |
EP3109589B1 (en) | A unit and method for improving positioning accuracy | |
US9360323B2 (en) | Systems and methods for estimating movements of a vehicle using a mobile device | |
US8452535B2 (en) | Systems and methods for precise sub-lane vehicle positioning | |
US7096116B2 (en) | Vehicle behavior detector, in-vehicle processing system, detection information calibrator, and in-vehicle processor | |
CN109752725A (zh) | 一种低速商用机器人、定位导航方法及定位导航系统 | |
CN104678415A (zh) | 车辆位置推测系统及其方法 | |
JP7272910B2 (ja) | 車載装置、状態推定方法及びプログラム | |
CN113405545B (zh) | 定位方法、装置、电子设备及计算机存储介质 | |
KR102331312B1 (ko) | 차량 내부 센서, 카메라, 및 gnss 단말기를 이용한 3차원 차량 항법 시스템 | |
CN110316197B (zh) | 倾斜估算方法、倾斜估算装置以及存储程序的非暂时性计算机可读存储介质 | |
TWI522258B (zh) | Based on electronic map, global navigation satellite system and vehicle motion detection technology Lane identification method | |
CN108387236B (zh) | 一种基于扩展卡尔曼滤波的偏振光slam方法 | |
CN110346824B (zh) | 一种车辆导航方法、系统、装置及可读存储介质 | |
CN113252048B (zh) | 一种导航定位方法、导航定位系统及计算机可读存储介质 | |
CN109631886B (zh) | 车辆定位方法、装置、电子设备、存储介质 | |
CN111751857A (zh) | 一种车辆位姿的估算方法、装置、存储介质及系统 | |
EP4345421A2 (en) | Method for calibrating sensor parameters based on autonomous driving, apparatus, storage medium, and vehicle | |
JP6248559B2 (ja) | 車両用走行軌跡算出装置 | |
CN115060257B (zh) | 一种基于民用级惯性测量单元的车辆变道检测方法 | |
CN110637209B (zh) | 估计机动车的姿势的方法、设备和具有指令的计算机可读存储介质 | |
CN115407376A (zh) | 车辆定位校准方法、装置、计算机设备、存储介质 | |
JP2016218015A (ja) | 車載センサ補正装置、自己位置推定装置、プログラム | |
WO2024149391A1 (zh) | 车载捷联组合导航方法、装置、电子设备及存储介质 | |
CN112611377B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |