CN109816731A - 一种rgb和深度信息精确配准的方法 - Google Patents

一种rgb和深度信息精确配准的方法 Download PDF

Info

Publication number
CN109816731A
CN109816731A CN201711170866.7A CN201711170866A CN109816731A CN 109816731 A CN109816731 A CN 109816731A CN 201711170866 A CN201711170866 A CN 201711170866A CN 109816731 A CN109816731 A CN 109816731A
Authority
CN
China
Prior art keywords
rgb
depth
coordinate
space coordinate
follows
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
Application number
CN201711170866.7A
Other languages
English (en)
Other versions
CN109816731B (zh
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201711170866.7A priority Critical patent/CN109816731B/zh
Publication of CN109816731A publication Critical patent/CN109816731A/zh
Application granted granted Critical
Publication of CN109816731B publication Critical patent/CN109816731B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明实施例公开了一种RGB和深度信息精确配准的方法,包括准备深度传感器和RGB传感器设备,对深度传感器做预处理,设置同步拍摄信号,搭建半自动化的棋盘格数据采集平台,同步拍摄棋盘格的深度图和RGB图进行标定,计算RGB图相应的像素坐标(iy,jy),通过计算出的IR空间坐标、RGB空间坐标、RGB像素值通过选择器选择输出想用坐标数据。通过本申请的方案,能够实现RGB和深度信息精确配准。

Description

一种RGB和深度信息精确配准的方法
技术领域
本发明属于图像处理、自然交互和集成电路技术领域,具体涉及一种RGB 和深度信息(Depth)精确配准的方法及装置。
背景技术
深度感知是计算机视觉领域重要的研究方向之一,能够让智能设备获得类似人类视觉感知的能力。现有的深度感知技术可分为主动式方法和被动式方法。基于双目立体视觉的深度感知技术是一种被动式方法,受环境光的影响很大且算法和硬件实现复杂。主动式方法主要有结构光法(coding light)和飞行时间法(time offilght)。其中基于结构光法的深度感知技术可以获取较为准确的深度信息,且不受环境光的影响、算法和硬件实现相对简单便捷。在三维重建、增强现实、虚拟现实、人机交互等应用领域,结合了depth深度信息和RGB图像信息的RGBD数据是不可或缺的,RGBD数据可以使深度学习算法更智能。
而RGBD数据的获取首先要解决深度传感器和RGB传感器之间的精确标定问题。而现有的摄像机标定方法或多或少存在一些问题,如传统的摄像机标定方法,标定精度较低;非线性的标定方法算法复杂,难以实现实时标定;张正友提出的基于单平面棋盘格模板的摄像机标定法方便快捷,但是深度传感器直接拍摄的深度图像不能生成棋盘格标定板的特征,无法直接用于标定。
发明内容
鉴于此,本发明提供了一种RGB和Depth精确配准的方法及装置,通过融合深度传感器和RGB传感器的标定与配准方法,实现RGBD快速精确配准,具有精度高、稳定性好的特点。
为实现以上发明目的,本发明采用以下技术方案:
第一方面,本发明实施例提供了一种RGB和深度信息精确配准的方法,包括以下步骤:
步骤1、准备深度传感器和RGB传感器设备,对深度传感器做预处理,设置同步拍摄信号;
步骤2、搭建半自动化的棋盘格数据采集平台,将棋盘格标定板置于深度传感器和RGB传感器设备之前,标定板和传感器设备分别以一定的方式固定在水平地面的滑轨上,其中传感器设备和标定板都可以一定的自由度前后移动或旋转,以便两者之间有一定的角度;
步骤3、同步拍摄棋盘格的深度图和RGB图进行标定,计算两个传感器的内外参数;
步骤4、同步拍摄深度图和RGB图,取深度图的像素坐标(id,jd),计算出点(id,jd)的距离值dis,利用标定出的深度传感器的内参结合dis,将深度图的像素坐标转换到空间坐标Pd(Xd,Yd,Zd),通过两个传感器的外参,将深度图的空间坐标Pd(Xd,Yd,Zd)转换到RGB图的空间坐标Pr(Xr,Yr,Zr),再利用 RGB传感器的内参结合Pr(Xr,Yr,Zr)计算出RGB图相应的像素坐标(ir,jr);
步骤5、通过步骤4计算出的IR空间坐标、RGB空间坐标、RGB像素值通过选择器选择输出想用坐标数据。
第二方面,本发明实施例提供了一种RGB和深度信息精确配准的装置,包括:
1、预处理模块。准备深度传感器和RGB传感器设备,对深度传感器做预处理。具体步骤:
(1)关闭深度传感器的红外投射器,避免其发出的散斑在红外图像中产生很多亮点,而不利于棋盘角点的检测。
(2)添加额外的红外辅助光源。
2、标定模块。主要是使用公开的彩色摄像机标定算法(如张正友标定算法) 计算出深度传感器的内参Ind{fx,d,fy,d,cx,d,cy,d}、RGB传感器的内参 Inr{fx,r,fy,r,cx,r,cy,r}、两个传感器的外参Ex{R,T}。其中fx,d,fy,d为深度传感器的焦距,cx,d,cy,d为深度传感器的中心点,fx,r,fy,r为RGB传感器的焦距,cx,r,cy,r为RGB传感器的中心点,R为3*3的旋转矩阵,T为3*1的平移矩阵。
3、点对点映射模块。用于实现深度图和RGB图像素点的一一映射。
将深度图的像素坐标转换到空间坐标。结合Pd=dis·Ind -1·pd其中 pd=[id,jd,1]T表示深度图像素点的齐次坐标,Ind表示深度传感器的内参矩阵,dis表示每个点到摄像头的距离值。
(3)将深度图的空间坐标转换到RGB图的空间坐标。因为两个传感器的坐标系不同,通过两个传感器的外参Ex{R,T},即其旋转平移关系将两个传感器的坐标系联系起来。Pr=dis·R·Ind -1·pd+T。其中,R、T分别表示旋转矩阵和平移矩阵。
(4)将RGB图的空间坐标转换到RGB图的像素坐标。假设RGB图上一个像素点用齐次 坐标表示为:pr=[ir,jr,1]T,RGB传感器的内参矩阵表示为 Inr,根据小孔成像模型有像素 坐标到空间坐标的映射关系为:z′·pr=Inr·Pr,其中zt是缩放因子,令z′=Zr,从而得到 RGB图的像素坐标为:
4、输出模块。通过选择器选择输出点云数据以及r、g、b数据。
5、寄存器控制模块。用于保存参数。f焦距、s基线、d参考距离、R、T 参数、μ像素大小。
6、深度计算模块。针对红外摄像头输入的红外散斑图像通过深度计算模块生成红外散斑图像对应的运动矢量图,每个像素点值用偏移量offset表示。
采用本发明技术方案的有益效果将通过以下实施例的阐述而得到具体的体现。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例提供方法的整体流程图;
图2是本发明实施例对深度传感器预处理后采集的棋盘格标定板效果图;
图3是本发明实施例配准算法具体实现流程图;
图4是本发明实施例中坐标转换算法的示意图。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
下面结合附图对本发明进行进一步的详细说明。
图1是本发明的整体流程图,本发明通过张正友标定算法计算深度传感器和RGB传感器的内外参数,再通过坐标转换将深度图和RGB图像素点一一映射,同时获取RGBD数据。下面对本发明的各个流程进行详细说明:
步骤一、准备深度传感器和RGB传感器设备,对深度传感器做预处理:
①关闭深度传感器的红外投射器,避免其发出的散斑在红外图像中产生很多亮点,而不利于棋盘角点的检测。
本发明在这里根据深度传感器只对红外光成像的特性,添加外置红外光源对物体照射,并关闭投射器,可以呈现清晰的棋盘格特征,在此基础上对深度传感器与RGB传感器进行标定,计算出两者的内外参数,使深度传感器的标定方便且容易实现。这也是本发明的优点。
②添加额外的红外辅助光源。
步骤2、搭建半自动化的棋盘格数据采集平台:
将棋盘格标定板置于深度传感器和RGB传感器设备之前,标定板和传感器设备分别固定在水平地面的滑轨上,其中传感器设备和标定板都可以前后移动或旋转,以便两者之间形成角度。
经上述步骤1、步骤2采集到的深度传感器的棋盘格标定板效果如图2 所示,图中可以清晰的呈现棋盘格的特征。
步骤3、同步拍摄棋盘格的深度图和RGB图进行标定,计算两个传感器的内外参数:
①调整标定板和传感器的相对角度,发送同步拍摄信号,同步采集标定板的深度图像和RGB图像。重复该过程直到采集到一定数量的标定板图像。
②使用现有的彩色摄像机标定算法(如张正友标定算法)计算深度传感器的内参Ind{fx,d,fy,d,cx,d,cy,d}、RGB传感器的内参Inr{fx,r,fy,r,cx,r,cy,r}、两个传感器的外参Ex{R,T}。其中fx,d,fy,d为深度传感器的焦距,cx,d,cy,d为深度传感器的中心点,fx,r,fy,r为RGB传感器的焦距,cx,r,cy,r为RGB传感器的中心点,R为3*3的旋转矩阵,T为3*1的平移矩阵。
步骤四、深度图和RGB图像素点的配准。
图3和图4实现同步拍摄深度图和RGB图,将深度图的像素坐标(id,jd) 经坐标转换计算出RGB图相应的像素坐标(ir,jr)。
①将深度图像的像素坐标转换到空间坐标。假设深度图上一个像素点用齐次坐标表示为:pd=[id,jd,1]T,其空间坐标表示为:pr=[ir,jr,1]T,深度传感器的内参矩阵表示为:
根据小孔成像模型有像素坐标到空间坐标的映射关系为:
其中z是缩放因子,因为深度传感器能够提供每个点到摄像头的距离值dis,所以令缩放因子z=dis。通过上式可计算出深度图的空间坐标:
简写为:Pd=dis·Ind -1·pd
②将深度图的空间坐标转换到RGB图的空间坐标。因为两个传感器的坐标系不同,通过两个传感器的外参Ex{R,T},即其旋转平移关系将两个传感器的坐标系联系起来。假设RGB图相应的空间坐标表示为:Pr=[Xr,Yr,Zr]T,两者空间坐标的映射关系为:
即:
简写为:Pr=dis·R·Ind -1·pd+T。
③将RGB图的空间坐标转换到RGB图的像素坐标。假设RGB图上一个像素点用齐次坐标表示为:pr=[ir,jr,1]T,RGB传感器的内参矩阵表示为:
根据小孔成像模型有像素坐标到空间坐标的映射关系为:
简写为:z′·pr=Inr·Pr。其中z′是缩放因子,令z′=Zr,从而得到RGB 图的像素坐标为:
步骤五、通过步骤4计算出的IR空间坐标、RGB空间坐标、RGB像素值通过选择器选择输出相关数据。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (4)

1.一种RGB和深度信息精确配准的方法,包括以下步骤:
步骤1、准备深度传感器和RGB传感器设备,对深度传感器做预处理,设置同步拍摄信号;
步骤2、搭建半自动化的棋盘格数据采集平台,将棋盘格标定板置于深度传感器和RGB传感器设备之前,标定板和传感器设备分别以一定的方式固定在水平地面的滑轨上,其中传感器设备和标定板都可以一定的自由度前后移动或旋转,以便两者之间有一定的角度;
步骤3、同步拍摄棋盘格的深度图和RGB图进行标定,计算两个传感器的内外参数;
步骤4、同步拍摄深度图和RGB图,取深度图的像素坐标(id,jd),计算出点(id,jd)的距离值dis,利用标定出的深度传感器的内参结合dis,将深度图的像素坐标转换到空间坐标Pd(Xd,Yd,Zd),通过两个传感器的外参,将深度图的空间坐标Pd(Xd,Yd,Zd)转换到RGB图的空间坐标Pr(Xr,Yr,Zr),再利用RGB传感器的内参结合Pr(Xr,Yr,Zr)计算出RGB图相应的像素坐标(ir,jr);
步骤5、通过步骤4计算出的IR空间坐标、RGB空间坐标、RGB像素值通过选择器选择输出想用坐标数据。
2.根据权利要求1所述的方法,其特征在于,优选的,所述准备深度传感器和RGB传感器设备,对深度传感器做预处理,包括:
关闭深度传感器的红外投射器,避免其发出的散斑在红外图像中产生很多亮点,而不利于棋盘角点的检测;
添加额外的红外辅助光源。
3.根据权利要求1中所述的方法,其特征在于,所述同步拍摄棋盘格的深度图和RGB图进行标定,计算两个传感器的内外参数,包括:
调整标定板和传感器的相对角度,发送同步拍摄信号,同步采集标定板的深度图像和RGB图像,重复该过程直到采集到一定数量的标定板图像;
使用现有的彩色摄像机标定算法(如张正友标定算法)计算深度传感器的内参Ind{fx,d,fy,d,cx,d,cy,d}、RGB传感器的内参Inr{fx,r,fy,r,cx,r,cy,r}、两个传感器的外参Ex{R,T},其中fx,d,fy,d为深度传感器的焦距,cx,d,cy,d为深度传感器的中心点,fx,r,fy,r为RGB传感器的焦距,cx,r,cy,r为RGB传感器的中心点,R为3*3的旋转矩阵,T为3*1的平移矩阵。
4.根据权利要求1中所述的方法,其特征在于,所述步骤4还包括:
将深度图像的像素坐标转换到空间坐标,将深度图上一个像素点用齐次坐标表示为:pd=[id,jd,1]T,其空间坐标表示为:Pd=[Xd,Yd,Zd]T,深度传感器的内参矩阵表示为:
根据小孔成像模型有像素坐标到空间坐标的映射关系为:
其中z是缩放因子,深度传感器能够提供每个点到摄像头的距离值dis,令缩放因子z=dis,通过上式可计算出深度图的空间坐标:
标记为:Pd=dis·Ind -1·pd
将深度图的空间坐标转换到RGB图的空间坐标,通过两个传感器的外参Ex{R,T},即其旋转平移关系将两个传感器的坐标系联系起来,RGB图相应的空间坐标表示为:Pr=[Xr,Yr,Zr]T,两者空间坐标的映射关系为:
即:
标记为:Pr=dis·R·Ind -1·pd+T;
(3)将RGB图的空间坐标转换到RGB图的像素坐标,RGB图上一个像素点用齐次坐标表示为:pr=[ir,jr,1]T,RGB传感器的内参矩阵表示为:
根据小孔成像模型有像素坐标到空间坐标的映射关系为:
标记为:z′·pr=Inr·Pr,其中z′是缩放因子,令z′=Zr,得到RGB图的像素坐标为:
CN201711170866.7A 2017-11-21 2017-11-21 一种rgb和深度信息精确配准的方法 Active CN109816731B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711170866.7A CN109816731B (zh) 2017-11-21 2017-11-21 一种rgb和深度信息精确配准的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711170866.7A CN109816731B (zh) 2017-11-21 2017-11-21 一种rgb和深度信息精确配准的方法

Publications (2)

Publication Number Publication Date
CN109816731A true CN109816731A (zh) 2019-05-28
CN109816731B CN109816731B (zh) 2021-08-27

Family

ID=66600169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711170866.7A Active CN109816731B (zh) 2017-11-21 2017-11-21 一种rgb和深度信息精确配准的方法

Country Status (1)

Country Link
CN (1) CN109816731B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110415303A (zh) * 2019-07-29 2019-11-05 华天慧创科技(西安)有限公司 一种相机3d成像标定装置及3d成像标定方法
CN110599550A (zh) * 2019-09-09 2019-12-20 香港光云科技有限公司 Rgb-d模组的校准系统及其设备和方法
CN110827230A (zh) * 2019-09-09 2020-02-21 香港光云科技有限公司 一种利用tof提升rgb图像质量的方法和装置
CN111654626A (zh) * 2020-06-05 2020-09-11 合肥泰禾光电科技股份有限公司 一种包含深度信息的高分辨率相机
CN112150372A (zh) * 2019-06-28 2020-12-29 深圳创想未来机器人有限公司 一种深度图矫正方法、深度图矫正装置和机器人
CN112509023A (zh) * 2020-12-11 2021-03-16 国网浙江省电力有限公司衢州供电公司 一种多源相机系统及rgbd配准方法
CN112911091A (zh) * 2021-03-23 2021-06-04 维沃移动通信(杭州)有限公司 多点激光器的参数调整方法、装置和电子设备
CN113099121A (zh) * 2021-04-15 2021-07-09 西安交通大学 一种基于弱监督学习的isp实现方法
CN114862923A (zh) * 2022-07-06 2022-08-05 武汉市聚芯微电子有限责任公司 一种图像配准方法及装置、存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103607584A (zh) * 2013-11-27 2014-02-26 浙江大学 一种kinect拍摄的深度图与彩色摄像机拍摄视频的实时配准方法
JP2016157188A (ja) * 2015-02-23 2016-09-01 株式会社デンソーアイティーラボラトリ 姿勢推定モデル生成装置及び姿勢推定装置
CN107016704A (zh) * 2017-03-09 2017-08-04 杭州电子科技大学 一种基于增强现实的虚拟现实实现方法
CN107067470A (zh) * 2017-04-05 2017-08-18 东北大学 基于红外热像仪与深度相机的便携式三维温度场重建系统
CN107103626A (zh) * 2017-02-17 2017-08-29 杭州电子科技大学 一种基于智能手机的场景重建方法
CN107301665A (zh) * 2017-05-03 2017-10-27 中国科学院计算技术研究所 具有可变焦光学摄像头的深度摄像头及其控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103607584A (zh) * 2013-11-27 2014-02-26 浙江大学 一种kinect拍摄的深度图与彩色摄像机拍摄视频的实时配准方法
JP2016157188A (ja) * 2015-02-23 2016-09-01 株式会社デンソーアイティーラボラトリ 姿勢推定モデル生成装置及び姿勢推定装置
CN107103626A (zh) * 2017-02-17 2017-08-29 杭州电子科技大学 一种基于智能手机的场景重建方法
CN107016704A (zh) * 2017-03-09 2017-08-04 杭州电子科技大学 一种基于增强现实的虚拟现实实现方法
CN107067470A (zh) * 2017-04-05 2017-08-18 东北大学 基于红外热像仪与深度相机的便携式三维温度场重建系统
CN107301665A (zh) * 2017-05-03 2017-10-27 中国科学院计算技术研究所 具有可变焦光学摄像头的深度摄像头及其控制方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112150372A (zh) * 2019-06-28 2020-12-29 深圳创想未来机器人有限公司 一种深度图矫正方法、深度图矫正装置和机器人
CN112150372B (zh) * 2019-06-28 2024-02-02 深圳创想未来机器人有限公司 一种深度图矫正方法、深度图矫正装置和机器人
CN110415303A (zh) * 2019-07-29 2019-11-05 华天慧创科技(西安)有限公司 一种相机3d成像标定装置及3d成像标定方法
CN110599550A (zh) * 2019-09-09 2019-12-20 香港光云科技有限公司 Rgb-d模组的校准系统及其设备和方法
CN110827230A (zh) * 2019-09-09 2020-02-21 香港光云科技有限公司 一种利用tof提升rgb图像质量的方法和装置
CN111654626A (zh) * 2020-06-05 2020-09-11 合肥泰禾光电科技股份有限公司 一种包含深度信息的高分辨率相机
CN112509023B (zh) * 2020-12-11 2022-11-22 国网浙江省电力有限公司衢州供电公司 一种多源相机系统及rgbd配准方法
CN112509023A (zh) * 2020-12-11 2021-03-16 国网浙江省电力有限公司衢州供电公司 一种多源相机系统及rgbd配准方法
CN112911091A (zh) * 2021-03-23 2021-06-04 维沃移动通信(杭州)有限公司 多点激光器的参数调整方法、装置和电子设备
CN112911091B (zh) * 2021-03-23 2023-02-24 维沃移动通信(杭州)有限公司 多点激光器的参数调整方法、装置和电子设备
CN113099121A (zh) * 2021-04-15 2021-07-09 西安交通大学 一种基于弱监督学习的isp实现方法
CN114862923B (zh) * 2022-07-06 2022-09-09 武汉市聚芯微电子有限责任公司 一种图像配准方法及装置、存储介质
CN114862923A (zh) * 2022-07-06 2022-08-05 武汉市聚芯微电子有限责任公司 一种图像配准方法及装置、存储介质

Also Published As

Publication number Publication date
CN109816731B (zh) 2021-08-27

Similar Documents

Publication Publication Date Title
CN109816731A (zh) 一种rgb和深度信息精确配准的方法
CN106097348B (zh) 一种三维激光点云与二维图像的融合方法
CN105931240B (zh) 三维深度感知装置及方法
CN105534694B (zh) 一种人体特征可视化装置及方法
CN106651794B (zh) 一种基于虚拟相机的投影散斑校正方法
CN103903222B (zh) 三维感测方法与三维感测装置
CN104034269B (zh) 一种单目视觉测量方法与装置
CN108256504A (zh) 一种基于深度学习的三维动态手势识别方法
CN100449459C (zh) 一种单摄像机虚拟鼠标系统的校准方法
CN104539928B (zh) 一种光栅立体印刷图像合成方法
CN104504671A (zh) 一种用于立体显示的虚实融合图像生成方法
CN105395252A (zh) 具有人机交互的可穿戴式血管介入手术三维立体图像导航装置
WO2023071834A1 (zh) 用于显示设备的对齐方法及对齐装置、车载显示系统
CN110074788B (zh) 一种基于机器学习的身体数据获取方法和装置
CN110264504A (zh) 一种用于增强现实的三维配准方法和系统
CN109615601A (zh) 一种融合彩色与灰阶深度图像的方法
CN105496556B (zh) 一种用于手术导航的高精度光学定位系统
CN109920000B (zh) 一种基于多相机协同的无死角的增强现实方法
CN103399629B (zh) 获取手势屏幕显示坐标的方法和装置
CN108388341A (zh) 一种基于红外摄像机-可见光投影仪的人机交互系统及装置
CN107221000A (zh) 基于增强现实的针灸穴位可视化平台及其图像处理方法
CN103445863A (zh) 基于平板电脑的手术导航和增强现实系统
US20130057656A1 (en) System and method for managing face data
CN107808412A (zh) 一种基于低成本的三维热源环境模型确定环境信息方法
CN107462174A (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