CN113933861A - 机载激光雷达点云生成方法及*** - Google Patents
机载激光雷达点云生成方法及*** Download PDFInfo
- Publication number
- CN113933861A CN113933861A CN202111336888.2A CN202111336888A CN113933861A CN 113933861 A CN113933861 A CN 113933861A CN 202111336888 A CN202111336888 A CN 202111336888A CN 113933861 A CN113933861 A CN 113933861A
- Authority
- CN
- China
- Prior art keywords
- data
- filtering
- point cloud
- information
- current
- 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
- 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/933—Lidar systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft
-
- 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/40—Correcting position, velocity or attitude
-
- 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
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Electromagnetism (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Navigation (AREA)
Abstract
本发明机载激光雷达点云生成方法及***涉及机载激光雷达,通过当前计算周期内的GNSS数据和IMU数据,并结合前一计算周期内优化后的第一位姿信息,获得当前周期的第一位姿信息,通过将激光雷达数据和局部地图进行匹配,获得第二位姿信息,在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差,对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,解决了机载激光雷达生成的点云数据精度受卫星信号和环境因素影响过大的问题,本发明适用于机载激光雷达。
Description
技术领域
本发明涉及机载激光雷达,特别涉及机载激光雷达点云生成方法及***。
背景技术
机载激光雷达***包括三大传感器设备:GNSS单元、IMU单元和激光雷达,其中GNSS单元用于确定激光发射点的位置信息;IMU单元用于测量激光发射瞬间激光的姿态信息;激光雷达设备以脉冲激光测距作为主要技术手段,以激光束扫描的工作方式测量从传感器到地物对象的激光照射点间的距离,同时还可以获得反射率、激光脉冲回波次数等信息。目前,机载激光雷达***点云生成的具体流程为:GNSS单元使用RTK技术和IMU***进行融合,形成组合导航***,计算出位姿信息;基于组合导航***计算出的位姿信息,结合激光雷达***获取的激光点距离信息,经过一系列空间变换,得出激光脚点在当地水平坐标系下的三维坐标,即一系列离散的、空间分布不规则的三维点云数据,所述激光脚点指激光落在物体表面的光点。
整个点云生成过程是线性过程,即点云的形成是完全基于组合导航***计算出的位姿信息,如果组合导航计算出的位姿信息误差过大,那么点云精度将急剧降低,因此,现有技术存在以下问题:1、如果卫星信号出现遮挡而无法进行RTK计算,将无法生成点云信息;2、如果卫星信息变弱或者因环境因素导致IMU的精度降低,组合导航计算出的位姿信息误差偏大,那么点云精度将降低。
发明内容
本发明所解决的技术问题:提供一种机载激光雷达点云生成方法及***,解决现有的机载激光雷达生成的点云数据精度受卫星信号和环境因素影响过大的问题。
本发明解决上述技术问题采用的技术方案:机载激光雷达点云生成方法,包括以下步骤:
S01、获取当前计算周期内的GNSS数据和IMU数据,作为第一数据;
S02、根据所述第一数据和前一计算周期内优化后的第一位姿信息,获得载体在当前计算周期的第一位姿信息;
S03、获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息;
S04、在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差;
S05、根据最新的位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
S06、对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图。
进一步的,所述位姿信息包括载体姿态和位置信息。
进一步的,步骤S03中,所述第二数据和局部地图进行匹配采用NDT匹配技术。
进一步的,步骤S04中,所述融合滤波采用kalman滤波算法。
进一步的,步骤S06中,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
机载激光雷达***,包括组合导航模块、点云数据生成模块、激光雷达单元和滤波单元,所述组合导航模块包括GNSS单元和IMU单元;
所述组合导航模块在当前计算周期内,获取GNSS单元和IMU单元采集到的数据,作为第一数据,根据所述第一数据和前一计算周期内优化后的第一位姿信息,计算获得载体在当前计算周期的第一位姿信息,将所述第一位姿信息发送给滤波单元,接收滤波单元反馈的位姿信息误差,并根据滤波单元反馈的最新位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
所述激光雷达单元获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息,将所述第二数据发送给点云数据生成模块,将所述第二位姿信息发送给滤波单元;
所述滤波单元在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差,将所述第三位姿信息发送给点云数据生成模块,将所述位姿信息误差发送给组合导航模块;
所述点云数据生成模块根据对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图,将所述局部地图发送给激光雷达单元。
进一步的,所述位姿信息包括载体姿态和位置信息。
进一步的,所述激光雷达单元将所述第二数据和最新的局部地图进行匹配采用NDT匹配技术。
进一步的,所述滤波单元采用kalman滤波算法。
进一步的,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
本发明的有益效果:本发明机载激光雷达点云生成方法及***,通过当前计算周期内的GNSS数据和IMU数据,并结合前一计算周期内优化后的第一位姿信息,获得当前周期的第一位姿信息,通过将激光雷达数据和局部地图进行匹配,获得第二位姿信息,在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差,对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图,解决了机载激光雷达生成的点云数据精度受卫星信号和环境因素影响过大的问题。本发明与现有技术相比,利用位姿信息误差对第一位姿信息进行优化,利用局部地图与激光雷达数据进行匹配,减少了卫星信号和环境因素对点云数据的影响,提高了点云数据的精度。
附图说明
附图1是本发明机载激光雷达点云生成方法及***的方法流程示意图。
具体实施方式
本发明中的机载激光雷达点云生成方法如附图1所述,包括以下步骤:
S01、获取当前计算周期内的GNSS数据和IMU数据,作为第一数据;
S02、根据所述第一数据和前一计算周期内优化后的第一位姿信息,获得载体在当前计算周期的第一位姿信息;
具体的,所述位姿信息包括载体姿态和位置信息。
S03、获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息;
具体的,所述第二数据和局部地图进行匹配采用NDT匹配技术。
S04、在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差;
具体的,所述融合滤波采用kalman滤波算法。
S05、根据最新的位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
S06、对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图。
具体的,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
本发明中的机载激光雷达***包括组合导航模块、点云数据生成模块、激光雷达单元和滤波单元,所述组合导航模块包括GNSS单元和IMU单元;
所述组合导航模块在当前计算周期内,获取GNSS单元和IMU单元采集到的数据,作为第一数据,根据所述第一数据和前一计算周期内优化后的第一位姿信息,计算获得载体在当前计算周期的第一位姿信息,将所述第一位姿信息发送给滤波单元,接收滤波单元反馈的位姿信息误差,并根据滤波单元反馈的最新位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
具体的,所述位姿信息包括载体姿态和位置信息。
所述激光雷达单元获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息,将所述第二数据发送给点云数据生成模块,将所述第二位姿信息发送给滤波单元;
具体的,所述第二数据和最新的局部地图进行匹配采用NDT匹配技术。
所述滤波单元在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差,将所述第三位姿信息发送给点云数据生成模块,将所述位姿信息误差发送给组合导航模块;
具体的,所述滤波单元采用kalman滤波算法。
所述点云数据生成模块根据对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图,将所述局部地图发送给激光雷达单元。
具体的,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
对本发明中的所有周期作进一步说明,计算周期为组合导航模块计算一次位姿信息的时间,扫描周期为激光雷达单元扫描一周或多周所用的时间,滤波周期为生成点云数据的间隔时间。
Claims (10)
1.机载激光雷达点云生成方法,其特征在于,包括以下步骤:
S01、获取当前计算周期内的GNSS数据和IMU数据,作为第一数据;
S02、根据所述第一数据和前一计算周期内优化后的第一位姿信息,获得载体在当前计算周期的第一位姿信息;
S03、获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息;
S04、在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差;
S05、根据最新的位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
S06、对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图。
2.根据权利要求1所述的机载激光雷达点云生成方法,其特征在于,所述位姿信息包括载体姿态和位置信息。
3.根据权利要求1或2所述的机载激光雷达点云生成方法,其特征在于,步骤S03中,所述第二数据和局部地图进行匹配采用NDT匹配技术。
4.根据权利要求1或2所述的机载激光雷达点云生成方法,其特征在于,步骤S04中,所述融合滤波采用kalman滤波算法。
5.根据权利要求1或2所述的机载激光雷达点云生成方法,其特征在于,步骤S06中,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
6.机载激光雷达***,应用于权利要求1所述的机载激光雷达点云生成方法,其特征在于,包括组合导航模块、点云数据生成模块、激光雷达单元和滤波单元,所述组合导航模块包括GNSS单元和IMU单元;
所述组合导航模块在当前计算周期内,获取GNSS单元和IMU单元采集到的数据,作为第一数据,根据所述第一数据和前一计算周期内优化后的第一位姿信息,计算获得载体在当前计算周期的第一位姿信息,将所述第一位姿信息发送给滤波单元,接收滤波单元反馈的位姿信息误差,并根据滤波单元反馈的最新位姿信息误差对当前计算周期的第一位姿信息进行优化,得到当前计算周期内优化后的第一位姿信息;
所述激光雷达单元获取当前扫描周期内的激光雷达数据,作为第二数据,将所述第二数据和最新的局部地图进行匹配,获得载体在当前扫描周期的第二位姿信息,将所述第二数据发送给点云数据生成模块,将所述第二位姿信息发送给滤波单元;
所述滤波单元在当前滤波周期内对最新的第一位姿信息和最新的第二位姿信息进行融合滤波,获得当前滤波周期对应的第三位姿信息和位姿信息误差,将所述第三位姿信息发送给点云数据生成模块,将所述位姿信息误差发送给组合导航模块;
所述点云数据生成模块根据对所述第三位姿信息以及第二数据进行空间变换,获得当前滤波周期内的点云数据,根据多个滤波周期内的点云数据生成局部地图,将所述局部地图发送给激光雷达单元。
7.根据权利要求6所述的机载激光雷达***,其特征在于,所述位姿信息包括载体姿态和位置信息。
8.根据权利要求6或7所述的机载激光雷达***,其特征在于,所述激光雷达单元将所述第二数据和最新的局部地图进行匹配采用NDT匹配技术。
9.根据权利要求6或7所述的机载激光雷达***,其特征在于,所述滤波单元采用kalman滤波算法。
10.根据权利要求6或7所述的机载激光雷达***,其特征在于,所述多个滤波周期为当前滤波周期前连续N个滤波周期,所述N为预设正整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111336888.2A CN113933861B (zh) | 2021-11-12 | 2021-11-12 | 机载激光雷达点云生成方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111336888.2A CN113933861B (zh) | 2021-11-12 | 2021-11-12 | 机载激光雷达点云生成方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113933861A true CN113933861A (zh) | 2022-01-14 |
CN113933861B CN113933861B (zh) | 2022-06-07 |
Family
ID=79286328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111336888.2A Active CN113933861B (zh) | 2021-11-12 | 2021-11-12 | 机载激光雷达点云生成方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113933861B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108320329A (zh) * | 2018-02-02 | 2018-07-24 | 维坤智能科技(上海)有限公司 | 一种基于3d激光的3d地图创建方法 |
CN111427060A (zh) * | 2020-03-27 | 2020-07-17 | 深圳市镭神智能***有限公司 | 一种基于激光雷达的二维栅格地图构建方法和*** |
CN111522043A (zh) * | 2020-04-30 | 2020-08-11 | 北京联合大学 | 一种无人车激光雷达快速重新匹配定位方法 |
CN112051590A (zh) * | 2020-08-31 | 2020-12-08 | 广州文远知行科技有限公司 | 激光雷达与惯性测量单元的检测方法及相关装置 |
CN112230242A (zh) * | 2020-09-30 | 2021-01-15 | 深兰人工智能(深圳)有限公司 | 位姿估计***和方法 |
CN112347840A (zh) * | 2020-08-25 | 2021-02-09 | 天津大学 | 视觉传感器激光雷达融合无人机定位与建图装置和方法 |
CN113432600A (zh) * | 2021-06-09 | 2021-09-24 | 北京科技大学 | 基于多信息源的机器人即时定位与地图构建方法及*** |
-
2021
- 2021-11-12 CN CN202111336888.2A patent/CN113933861B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108320329A (zh) * | 2018-02-02 | 2018-07-24 | 维坤智能科技(上海)有限公司 | 一种基于3d激光的3d地图创建方法 |
CN111427060A (zh) * | 2020-03-27 | 2020-07-17 | 深圳市镭神智能***有限公司 | 一种基于激光雷达的二维栅格地图构建方法和*** |
CN111522043A (zh) * | 2020-04-30 | 2020-08-11 | 北京联合大学 | 一种无人车激光雷达快速重新匹配定位方法 |
CN112347840A (zh) * | 2020-08-25 | 2021-02-09 | 天津大学 | 视觉传感器激光雷达融合无人机定位与建图装置和方法 |
CN112051590A (zh) * | 2020-08-31 | 2020-12-08 | 广州文远知行科技有限公司 | 激光雷达与惯性测量单元的检测方法及相关装置 |
CN112230242A (zh) * | 2020-09-30 | 2021-01-15 | 深兰人工智能(深圳)有限公司 | 位姿估计***和方法 |
CN113432600A (zh) * | 2021-06-09 | 2021-09-24 | 北京科技大学 | 基于多信息源的机器人即时定位与地图构建方法及*** |
Non-Patent Citations (2)
Title |
---|
石鹏等: "复杂环境下微小飞行器惯性/激光雷达Robust-SLAM方法", 《导航定位与授时》 * |
石鹏等: "复杂环境下微小飞行器惯性/激光雷达Robust-SLAM方法", 《导航定位与授时》, vol. 6, no. 1, 31 January 2019 (2019-01-31), pages 14 - 21 * |
Also Published As
Publication number | Publication date |
---|---|
CN113933861B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108921947B (zh) | 生成电子地图的方法、装置、设备、存储介质以及采集实体 | |
Glennie | Rigorous 3D error analysis of kinematic scanning LIDAR systems | |
Hebert et al. | 3D measurements from imaging laser radars: how good are they? | |
Roman et al. | A self‐consistent bathymetric mapping algorithm | |
CN112684432B (zh) | 激光雷达标定方法、装置、设备及存储介质 | |
CN105264402A (zh) | 具有lidar的水下平台和相关方法 | |
US11994588B2 (en) | Three-dimensional surface scanning | |
CN106990401A (zh) | 基于全波形机载激光雷达数据二类高程误差修正方法 | |
JP7007137B2 (ja) | 情報処理装置、情報処理方法および情報処理用プログラム | |
Paffenholz et al. | Direct geo-referencing of a static terrestrial laser scanner | |
CN104849723B (zh) | 一种基于多元线阵激光雷达的模拟月表地形的识别方法 | |
JP2017072423A (ja) | 推定装置、制御方法、プログラム及び記憶媒体 | |
CN104297743A (zh) | 消除高重频机载激光雷达***的测距歧义的方法和装置 | |
CN110018491B (zh) | 激光扫描方法、装置及激光雷达 | |
JP4852006B2 (ja) | 空間情報データベース生成装置及び空間情報データベース生成プログラム | |
CN111123281B (zh) | 利用基于距离探测和3d点分析的歧义消解的空中多脉冲激光扫描*** | |
CN113933861B (zh) | 机载激光雷达点云生成方法及*** | |
KR101074277B1 (ko) | 구조물 윤곽선 추출장치 및 방법 | |
JPH10318743A (ja) | 飛行体を用いた測量方法及びその装置 | |
RU2680969C1 (ru) | Способ навигации летательных аппаратов | |
CN113238202B (zh) | 光子激光三维成像***的坐标系点云计算方法及其应用 | |
Hebel et al. | LiDAR-supported navigation of UAVs over urban areas | |
JP6470658B2 (ja) | レーザ計測システムおよびレーザ計測方法 | |
TW202137138A (zh) | 含大地座標之3d點雲地圖的產生方法和3d點雲地圖產生系統 | |
CN105652260B (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 |