CN112629548A - 一种基于路侧激光雷达的三维栅格地图的创建与更新方法 - Google Patents
一种基于路侧激光雷达的三维栅格地图的创建与更新方法 Download PDFInfo
- Publication number
- CN112629548A CN112629548A CN202011585046.6A CN202011585046A CN112629548A CN 112629548 A CN112629548 A CN 112629548A CN 202011585046 A CN202011585046 A CN 202011585046A CN 112629548 A CN112629548 A CN 112629548A
- Authority
- CN
- China
- Prior art keywords
- grid
- obstacle
- updating
- laser radar
- laser
- 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
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/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
- G01C21/30—Map- or contour-matching
- G01C21/32—Structuring or formatting of map data
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/93—Lidar systems specially adapted for specific applications for anti-collision purposes
- G01S17/931—Lidar systems specially adapted for specific applications for anti-collision purposes of land vehicles
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Automation & Control Theory (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Traffic Control Systems (AREA)
Abstract
本专利公开了一种基于路侧激光雷达的三维栅格地图的创建与更新方法,所述方法包括:步骤一、基于路侧激光雷达构建三维栅格地图;步骤二、确定障碍物在空间栅格中的坐标;步骤三、基于贝叶斯公式的对空间栅格地图中的障碍物状态进行判断与更新;步骤四、在更新空间栅格地图障碍物的同时确定非障碍物的三维栅格集。通过上述方法,能够准确地快速地建立一种基于路侧激光雷达的三维栅格地图并且在场景发生变化时,可以及时准确地利用路侧数字地图更新障碍物和非障碍物的坐标位置,以支撑无人驾驶车辆的正常行驶。
Description
技术领域
本发明属于车联网和辅助驾驶技术领域,具体涉及一种基于路侧激光雷达的三维栅格地图的创建与更新方法。
背景技术
无人驾驶车辆的正常行驶依赖对道路信息的全方位感知。城市场景中,道路开放,结构复杂,交通流量大,行为差异显著,这些变化对无人驾驶车辆的环境感知能力有着更高的要求。目前,无人驾驶车载核心传感器感知距离有限,无法满足恶劣天气、遮挡区等非视距环境感知需求。
作为一种以激光束探测目标位置、速度等特征的雷达***,激光雷达具有较远的探测距离、较高的角度和速度分辨率,因此被广泛地用作车载传感器,为车辆对周边障碍物的感知提供助力。随着V2X与5G通信技术的发展,将激光雷达部署于路侧,进而构建动态更新的高精度地图成为可能。与部署在车辆上相比,将激光雷达部署在路侧一方面可以通过较高的部署位置保证视野不被复杂的城市道路场景遮挡,另一方面可以为多辆网联汽车服务,有效降低成本。
针对以上情况,本发明提出了一种基于路侧激光雷达的三维栅格地图的创建与更新方法,旨在赋能无人驾驶,辅助实现实时、精准的道路环境全方位感知。
发明内容
本发明要解决的技术问题是:提出一种基于路侧激光雷达的三维栅格地图的创建与更新方法,实现在场景发生变化时,路侧数字地图可以及时更新以支撑无人驾驶车辆的正常行驶。
为此,本发明提供的技术方案包括:
一种基于路侧激光雷达的三维栅格地图的创建与更新方法,其特征在于,所述方法包括:步骤一、基于路侧激光雷达构建三维栅格地图;其包括,通过道路两侧交替布设的激光雷达获得测量,并将所述空间设定为由若干大小相等的小立方体构成的空间栅格;步骤二、确定障碍物在空间栅格中的坐标,其包括,根据激光束的发射方向β、激光束的俯仰角γ以及激光发射点距障碍物的距离d,据此计算障碍点的实际位置,并根据障碍点实际位置与空间栅格中的栅格的对应关系确定障碍点在空间栅格中的坐标;步骤三、基于贝叶斯公式的对空间栅格地图中的障碍物状态进行判断与更新,其包括,使用A+(s)=A(s)+Mj更新单个栅格A的状态A(s),其中A+(s)为测量之后的单个栅格A的状态,Mj表示测量值m的状态,j∈0,1,M0=-1,M1=1。步骤四、在更新空间栅格地图障碍物的同时确定非障碍物的三维栅格集,当已经确定了栅格内激光束在x-y坐标系投影到障碍物之间的第I个栅格I(Xi,Yi)的位置,激光束的发射方向与x轴呈β角时,确定第I+1个栅格(Xi+1,Yi+1)的位置;当光束在Xi+1点的实际纵向位置Y=Xitanβ,l1=Y-Yi,l2=Yi+1时;若l1>l2,则Yi+1=Yi+1,第I+1个非障碍物的栅格取(Xi+1,Yi+1);若l1<l2,则Yi+1=Yi第I+1个非障碍物的栅格取(Xi+1,Yi);若l1=l2,栅格(Xi+1,Yi)和栅格(Xi+1,Yi+1)均入选下一个非障碍物的栅格。
优选地,当激光雷达中的激光器具有最大俯仰角±α、具备最小量程rmin,则,在路侧布设时的最小量程rmin=q+(h/2+H)/tanα;其中,H为激光雷达支架的高度、h为激光雷达自身高度、q为激光雷达自身半径;rmin=q+(h/2+H)/tanα。
通过上述方法,能够准确地快速地建立一种基于路侧激光雷达的三维栅格地图并且在场景发生变化时,可以及时准确地利用路侧数字地图更新障碍物和非障碍物的坐标位置,以支撑无人驾驶车辆的正常行驶。
附图说明
图1为激光雷达的激光束发射图;
图2为最短量程计算方法;
图3为路侧激光雷达的布设方法;
图4为空间栅格化示意图;
图5为坐标计算方法;
图6为坐标转化示意图;
图7为在确定了栅格I后确定栅格I+1的算法图解;
图8为可通行区域确定算法流程图。
具体实施方式
为更好地说明本发明所提出地三维栅格地图地创建与更新方法,结合专利附图对具体实施方式进行详细说明:
Step1:路侧激光雷达的布设:如图2与3所示,因激光雷达中的激光器具有最大俯仰角±α,故其具备最小量程rm,即距离r以内的障碍物将无法被检测。在路侧布设时,最小量程rmin会随着激光雷达支架的高度H而发生变化,假设激光雷达自身高度为h,半径为q即:
rmin=q+(h/2+H)/tanα
此外,路侧布设还需保证激光雷达无盲区,因此采用间隔的方式在道路两侧进行布设。假设激光雷达的最大量程为R,道路宽度为D,则如图2所示,两个相邻的激光雷达布设点位沿道路方向的距离应为:
布设位置为道路两侧交替布设,距离道路边缘距离:
r≥rmin
Step2:空间的栅格化:本发明需将激光雷达的测量空间栅格化,即将其视为由若干大小相等的小立方体构成,如图4所示。立方体的棱长由使用场景决定,棱长越小则地图精确度越高,同时占用的计算资源也越多。假设其棱长为a,则地图分辨率为
Step3:障碍物坐标计算及坐标转换:如图5所示,以激光雷达朝向道路方向的切面中点为原点建立坐标系;当激光雷达捕获障碍物时,对于其中一束激光,所获取的信息为:激光束的发射方向β,激光束的俯仰角γ以及激光发射点距障碍物的距离d,据此计算障碍点b(x1,y1,z1)的实际坐标为:
x1=dcosγcosβ;y1=dcosγsinβ;z1=dsinγ
为将障碍物完整的反映在三维栅格地图中,需将真实坐标转化为栅格坐标。如图6所示,对于真实坐标系中的一个障碍点,假设其x轴坐标为x,若na<x≤(n+1)a(即该点的x轴坐标位于第n至第n+1个栅格之间),则在栅格坐标系中,其x轴坐标X为n+1,即
故对上述障碍物点b(x1,y1,z1),在栅格坐标系内,其坐标为:
按上述逻辑,即可根据激光雷达扫描到的障碍点完成3D栅格地图的构建,并完成对障碍物轮廓的标注。
Step4:基于贝叶斯公式的栅格状态判断与更新规则:本发明基于概率分布对三维栅格是否被障碍物占据进行判断。在实现了三维栅格对道路的全覆盖后,对于足够小单个立方体A而言,其状态具有两种可能:未被障碍物占据(s=0)与被障碍物占据(s=1),用概率可以表示为:
p(s=0)+p(s=1)=1
进一步使用二者比值反映A的状态
对于该立方体,若获取了新的测量值m~{0,1},则需要更新该点的状态,仍旧使用概率表示:
根据贝叶斯公式:
带入A(s)可得:
等号两边取对数:
使用Mj表示测量值m的状态,对于一个立方体,m存在两种情况,即为0或1,即j∈0,1因此该测量值也有两种情况:
其中,M0表示未检测到障碍物,M1表示检测到障碍物。
进一步的,立方体A的状态A(s)更新规则可以写作:
A+(s)=A(s)+Mj
其中,A+(s)为测量之后的立方体A的状态,当没有测量值的初始状态下,A0(s)=0。
在实际使用中,每一台激光雷达的检测精度是固定的,因此Mj也是定值,为了便于判断立方体是否被障碍物占据,规定M0=-1,M1=1,则对于一个立方体,其状态A(s)的越大,说明该立方体处被障碍物占据的可能性越高;相反,A(s)值越小,则说明该立方体未被障碍物占据的可能性越小;当A(s)为0时,需重新对该处进行检测,并依据式x进行状态更新。
Step5:基于改进的Bresenham法的可通行三维栅格的确定:为了为道路上行驶的车辆提供参考,本发明在更新三维栅格地图障碍物的同时使用改进的Bresenham法确定可通行区域,即非障碍物的三维栅格集。图7所示为在确定了栅格I后确定栅格I+1的算法图解,假设已经确定了栅格内激光束在x-y坐标系投影到障碍物之间的第I个栅格I(Xi,Yi)的位置,要进一步确定第I+1个栅格(Xi+1,Yi+1)的位置,因已知激光束的发射方向β,在x-y坐标系内,该直线方程可以写作y=xtanβ,假设Xi+1=Xi+1,则激光束在Xi+1的实际纵向位置为Y=Xitanβ,令l1=Y-Yi,l2=Yi+1,则l1、l2分别反应的是沿着该直线,栅格(Xi+1,Yi)和栅格(Xi+1,Yi+1)距该直线的距离。比较l1、l2的大小,若:l1>l2,则Yi+1=Yi+1,下个栅格取(Xi+1,Yi+1);若:l1<l2,则Yi+1=Yi下个栅格取(Xi+1,Yi);若l1=l2,则为了保证安全,栅格(Xi+1,Yi)和栅格(Xi+1,Yi+1)均入选。此后该方法对每一束激光光束的周边的可通行栅格进行标注,直到条邻近路径之间不存在其他障碍物,则中间区域即为可通行区域,划分结果如图8所示。
Claims (4)
1.一种基于路侧激光雷达的三维栅格地图的创建与更新方法,其特征在于,所述方法包括:
步骤一、基于路侧激光雷达构建三维栅格地图;其包括,通过道路两侧交替布设的激光雷达获得测量,并将所述空间设定为由若干大小相等的小立方体构成的空间栅格;
步骤二、确定障碍物在空间栅格中的坐标,其包括,根据激光束的发射方向β、激光束的俯仰角γ以及激光发射点距障碍物的距离d,据此计算障碍点的实际位置,并根据障碍点实际位置与空间栅格中的栅格的对应关系确定障碍点在空间栅格中的坐标
步骤三、基于贝叶斯公式的对空间栅格地图中的障碍物状态进行判断与更新,其包括,使用A+(s)=A(s)+Mj更新单个栅格A的状态A(s),其中A+(s)为测量之后的单个栅格A的状态,Mj表示测量值m的状态,j∈0,1,M0=-1,M1=1。
步骤四、在更新空间栅格地图障碍物的同时确定非障碍物的三维栅格集,当已经确定了栅格内激光束在x-y坐标系投影到障碍物之间的第I个栅格I(Xi,Yi)的位置,激光束的发射方向与x轴呈β角时,确定第I+1个栅格(Xi+1,Yi+1)的位置;当光束在Xi+1点的实际纵向位置Y=Xitanβ,l1=Y-Yi,l2=Yi+1时;若l1>l2,则Yi+1=Yi+1,第I+1个非障碍物的栅格取(Xi+1,Yi+1);若l1<l2,则Yi+1=Yi第I+1个非障碍物的栅格取(Xi+1,Yi);若l1=l2,栅格(Xi+1,Yi)和栅格(Xi+1,Yi+1)均入选下一个非障碍物的栅格。
2.根据权利要求1所述的一种基于路侧激光雷达的三维栅格地图的创建与更新方法,其特征在于,当激光雷达中的激光器具有最大俯仰角±α、具备最小量程rmin,则,在路侧布设时的最小量程rmin=q+(h/2+H)/tanα;其中,H为激光雷达支架的高度、h为激光雷达自身高度、g为激光雷达自身半径;rmin=q+(h/2+H)/tanα。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011585046.6A CN112629548B (zh) | 2020-12-28 | 2020-12-28 | 一种基于路侧激光雷达的三维栅格地图的创建与更新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011585046.6A CN112629548B (zh) | 2020-12-28 | 2020-12-28 | 一种基于路侧激光雷达的三维栅格地图的创建与更新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112629548A true CN112629548A (zh) | 2021-04-09 |
CN112629548B CN112629548B (zh) | 2023-01-24 |
Family
ID=75285737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011585046.6A Active CN112629548B (zh) | 2020-12-28 | 2020-12-28 | 一种基于路侧激光雷达的三维栅格地图的创建与更新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112629548B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114253273A (zh) * | 2021-12-23 | 2022-03-29 | 南京世泽科技有限公司 | 一种基于多线激光雷达的避障方法 |
CN115493600A (zh) * | 2022-10-12 | 2022-12-20 | 江苏盛海智能科技有限公司 | 一种基于激光雷达的无人驾驶遇障路径规划方法与终端 |
WO2023071377A1 (zh) * | 2021-10-25 | 2023-05-04 | 珠海一微半导体股份有限公司 | 基于激光观测方向的栅格和体素定位方法、机器人及芯片 |
CN116381698A (zh) * | 2023-06-05 | 2023-07-04 | 蘑菇车联信息科技有限公司 | 道路遗撒物的检测方法、装置及电子设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850834A (zh) * | 2015-05-11 | 2015-08-19 | 中国科学院合肥物质科学研究院 | 基于三维激光雷达的道路边界检测方法 |
CN106529431A (zh) * | 2016-10-31 | 2017-03-22 | 武汉大学 | 基于车载激光扫描数据的路坎点自动提取及矢量化方法 |
CN109031346A (zh) * | 2018-07-09 | 2018-12-18 | 江苏大学 | 一种基于3d激光雷达的周边泊车位辅助检测方法 |
US20190353790A1 (en) * | 2018-09-12 | 2019-11-21 | Baidu Online Network Technology (Beijing) Co., Ltd. | Ranging Method Based on Laser Radar System, Device and Readable Storage Medium |
US20200116867A1 (en) * | 2018-10-12 | 2020-04-16 | Baidu Usa Llc | Automatic lidar calibration based on pre-collected static reflection map for autonomous driving |
US20200183011A1 (en) * | 2018-12-11 | 2020-06-11 | Shanghai Tusen Weilai Artificial Intelligence Technology Co., Ltd. | Method for creating occupancy grid map and processing apparatus |
CN111373337A (zh) * | 2017-08-23 | 2020-07-03 | 图森有限公司 | 使用基于相机的子地图和基于lidar的全局地图进行厘米精度定位的3d子地图重构***和方法 |
CN111912419A (zh) * | 2020-08-28 | 2020-11-10 | 清华大学苏州汽车研究院(吴江) | 基于激光雷达的高精度语义导航地图构建方法和装置 |
CN111985322A (zh) * | 2020-07-14 | 2020-11-24 | 西安理工大学 | 一种基于激光雷达的道路环境要素感知方法 |
-
2020
- 2020-12-28 CN CN202011585046.6A patent/CN112629548B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850834A (zh) * | 2015-05-11 | 2015-08-19 | 中国科学院合肥物质科学研究院 | 基于三维激光雷达的道路边界检测方法 |
CN106529431A (zh) * | 2016-10-31 | 2017-03-22 | 武汉大学 | 基于车载激光扫描数据的路坎点自动提取及矢量化方法 |
CN111373337A (zh) * | 2017-08-23 | 2020-07-03 | 图森有限公司 | 使用基于相机的子地图和基于lidar的全局地图进行厘米精度定位的3d子地图重构***和方法 |
CN109031346A (zh) * | 2018-07-09 | 2018-12-18 | 江苏大学 | 一种基于3d激光雷达的周边泊车位辅助检测方法 |
US20190353790A1 (en) * | 2018-09-12 | 2019-11-21 | Baidu Online Network Technology (Beijing) Co., Ltd. | Ranging Method Based on Laser Radar System, Device and Readable Storage Medium |
US20200116867A1 (en) * | 2018-10-12 | 2020-04-16 | Baidu Usa Llc | Automatic lidar calibration based on pre-collected static reflection map for autonomous driving |
US20200183011A1 (en) * | 2018-12-11 | 2020-06-11 | Shanghai Tusen Weilai Artificial Intelligence Technology Co., Ltd. | Method for creating occupancy grid map and processing apparatus |
CN111985322A (zh) * | 2020-07-14 | 2020-11-24 | 西安理工大学 | 一种基于激光雷达的道路环境要素感知方法 |
CN111912419A (zh) * | 2020-08-28 | 2020-11-10 | 清华大学苏州汽车研究院(吴江) | 基于激光雷达的高精度语义导航地图构建方法和装置 |
Non-Patent Citations (1)
Title |
---|
蔡自兴等: "基于激光雷达的动态障碍物实时检测", 《控制工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023071377A1 (zh) * | 2021-10-25 | 2023-05-04 | 珠海一微半导体股份有限公司 | 基于激光观测方向的栅格和体素定位方法、机器人及芯片 |
CN114253273A (zh) * | 2021-12-23 | 2022-03-29 | 南京世泽科技有限公司 | 一种基于多线激光雷达的避障方法 |
CN114253273B (zh) * | 2021-12-23 | 2024-04-12 | 南京世泽科技有限公司 | 一种基于多线激光雷达的避障方法 |
CN115493600A (zh) * | 2022-10-12 | 2022-12-20 | 江苏盛海智能科技有限公司 | 一种基于激光雷达的无人驾驶遇障路径规划方法与终端 |
CN116381698A (zh) * | 2023-06-05 | 2023-07-04 | 蘑菇车联信息科技有限公司 | 道路遗撒物的检测方法、装置及电子设备 |
CN116381698B (zh) * | 2023-06-05 | 2024-03-12 | 蘑菇车联信息科技有限公司 | 道路遗撒物的检测方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112629548B (zh) | 2023-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112629548B (zh) | 一种基于路侧激光雷达的三维栅格地图的创建与更新方法 | |
CN102222236B (zh) | 图像处理***及位置测量*** | |
CN110816548A (zh) | 传感器融合 | |
CN100371736C (zh) | 车载雷达的轴向偏差确定方法 | |
US11861784B2 (en) | Determination of an optimal spatiotemporal sensor configuration for navigation of a vehicle using simulation of virtual sensors | |
CN107850450A (zh) | 用于生成及使用定位参考数据的方法及*** | |
US20230202519A1 (en) | Route optimization for autonomous driving systems | |
JP2017194527A (ja) | 環境地図のデータ構造、環境地図の作成システム及び作成方法、並びに、環境地図の更新システム及び更新方法 | |
CN111121849B (zh) | 传感器的方位参数的自动校准方法、边缘计算单元和路侧传感*** | |
US11237269B2 (en) | Localization technique | |
US20220128995A1 (en) | Velocity estimation and object tracking for autonomous vehicle applications | |
JP2022022287A (ja) | 地図作成装置、制御方法、プログラム及び記憶媒体 | |
CN110203192A (zh) | 定位技术选择 | |
US20230065727A1 (en) | Vehicle and vehicle control method | |
US20210304602A1 (en) | Enhanced vehicle operation | |
JPWO2018180081A1 (ja) | 劣化地物特定装置、劣化地物特定方法、劣化地物特定プログラム及び劣化地物特定プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
CN113566817B (zh) | 一种车辆定位方法及装置 | |
CN113494915A (zh) | 一种车辆横向定位方法、装置及*** | |
CN117387647A (zh) | 融合车载传感器数据与道路传感器数据的道路规划方法 | |
JP2023118751A (ja) | 自己位置推定装置 | |
JP4506163B2 (ja) | 前方物体検出装置及び前方物体検出方法 | |
US11753045B2 (en) | Modeling positional uncertainty of moving objects using precomputed polygons | |
US20220194424A1 (en) | Autonomous vehicle trajectory computed using box-based probabilistic overlap model | |
CN114089376A (zh) | 一种基于单激光雷达的负障碍物检测方法 | |
CN116990832B (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 |