CN109035309A - 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法 - Google Patents

基于立体视觉的双目摄像头与激光雷达间的位姿配准方法 Download PDF

Info

Publication number
CN109035309A
CN109035309A CN201810801144.5A CN201810801144A CN109035309A CN 109035309 A CN109035309 A CN 109035309A CN 201810801144 A CN201810801144 A CN 201810801144A CN 109035309 A CN109035309 A CN 109035309A
Authority
CN
China
Prior art keywords
camera
laser radar
coordinate
radar
coordinate 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.)
Granted
Application number
CN201810801144.5A
Other languages
English (en)
Other versions
CN109035309B (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.)
Suzhou Automotive Research Institute of Tsinghua University
Original Assignee
Suzhou Automotive Research Institute of Tsinghua 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 Suzhou Automotive Research Institute of Tsinghua University filed Critical Suzhou Automotive Research Institute of Tsinghua University
Priority to CN201810801144.5A priority Critical patent/CN109035309B/zh
Publication of CN109035309A publication Critical patent/CN109035309A/zh
Application granted granted Critical
Publication of CN109035309B publication Critical patent/CN109035309B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10044Radar image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,包括:通过将物体在激光雷达坐标系的坐标进行仿射变换,转变成物体在摄像机坐标系中的坐标,经过摄影变换转成物体在图像坐标系的位置;根据双目摄像头获取的视差图和物体在摄像机坐标系中的Z坐标的映射关系,将雷达坐标转换成深度图;将摄像机深度图以及雷达深度图进行互信息配准,对激光雷达的数据点和双目深度图进行位姿态配准。可以方便的获得激光雷达以及双目摄像头之间的位姿关系,从而可以方便快速的将视觉信息以及激光雷达的信息进行融合。

Description

基于立体视觉的双目摄像头与激光雷达间的位姿配准方法
技术领域
本发明属于激光雷达和摄像头之间的融合技术领域,具体地涉及一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法。
背景技术
自动驾驶环境感知***是一个信息冗余***,其环境感知传感器包括如激光雷达、双目摄像头、单目摄像头、超声波、毫米波雷达等。而在这个冗余感知***中,不同传感器之间的信息配准以及融合是自动驾驶的环境感知模块中的关键步骤。
激光雷达和摄像头之间的融合是目前的一个研究热点,其主要原理如图1所示。
激光雷达数据到摄像头数据的坐标变换,其实就是一个旋转平移矩阵[R,t]。为了表述方便,如图1所示,我们将摄像机坐标系定义为OcXcYcZc,雷达坐标系定义为OlXlYlZl,摄像机的图像坐标系定义为OcUcVc。对于空间任意一点M,其在摄像机坐标系下的坐标Xc,Yc,Zc,其在雷达坐标系下的坐标为(Xl,Yl,Zl),其在摄像机的图像坐标系下的坐标为(u,v),这样对于空间之间的坐标转换关系如下:
dx代表x轴方向一个像素的宽度,dy代表y轴方向上一个像素的宽度。dx、dy为摄像机的内参数。(u0,v0)称为图像平面的主点,也是摄像机的内参数。摄像机的内参数可以通过张正友的标定方法进行快速精确的标定。为了表述方便我们将内参矩阵定义为P,空间点在摄像机坐标系以及雷达坐标系的向量定义为Mc,Ml
摄像机坐标系和雷达坐标系之间满足:
R代表旋转,T代表平移。R、T与摄像机无关,所以称这两个参数为摄像机的外参数。有以上关系可以得出摄像机图像坐标系和雷达坐标系存在如下关系:
在公式中,仅有R,T两个未知矩阵,因此存在两类标定方法。
1.基于Phantom的点到点对应关系标定,此方法直接利用上述公式,通过获取空间点在雷达坐标系中的坐标以及在摄像机图像坐标系上的坐标,然后通过解线性非齐次方程组的方法进行标定。但是由于雷达数据的稀疏特性,因此需要通过曲线拟合的方法获取标定点。如图2所示,由于激光雷达点的稀疏特性,可能雷达射线并不能恰好扫描到模体的顶点A,因此我们需要先进性模型边线的直线拟合,边线拟合的交点即是模型的顶点A,其他的标定点可以同理获取。该方法操作相对复杂且该配准方法的配准精度会受到很大影响。
2.另一种方法基于雷达强度图像以及摄像头强度图像之间的互信息来进行摄像头以及激光雷达的位姿配准。激光雷达获取的强度分布统计与摄像头图像的强度分布统计之间的互信息与它们之间的位姿关系R,T存在一个一个正相关关系,当R,T靠近摄像头和激光雷达之间位姿关系的真值时,其互信息值达到最大。这种方法有一个关键的假设:对于同一个物体激光雷达的探测强度和摄像头图像的强度相同。但是这个假设对环境的要求过于苛刻,因为激光雷达是主动光探测,其探测强度受环境光照影响很小,但是摄像头图像却是被动探测,其探测强度受环境影响很大。因此此种方法的鲁棒性很差。
发明内容
为了解决上述存在的技术问题,本发明提供了一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,可以方便的获得激光雷达以及双目摄像头之间的位姿关系,从而可以方便快速的将视觉信息以及激光雷达的信息进行融合。
本发明的技术方案是:
一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,包括以下步骤:
S01:通过将物体在激光雷达坐标系的坐标进行仿射变换,转变成物体在摄像机坐标系中的坐标,经过摄影变换转成物体在图像坐标系的位置;
S02:根据双目摄像头获取的视差图和物体在摄像机坐标系中的Z坐标的映射关系,将雷达坐标转换成深度图;
S03:将摄像机深度图以及雷达深度图进行互信息配准,对激光雷达的数据点和双目深度图进行位姿态配准。
优选的技术方案中,所述Z坐标映射关系为:其中d是视差,f是摄像头焦距,b是两个摄像头之间的基线距离。
优选的技术方案中,所述互信息配准包括以下步骤:
S31:基于双目摄像头图像获取深度图I1
S32:基于初始化的旋转平移矩阵[R,t]对雷达数据进行投影变换,获取激光雷达的深度图I2
S33:通过两幅深度图进行互信息计算f_MI(I1,I2),通过梯度下降法获取更新后的旋转平移矩阵[R,t],并再次进行互信息计算s_MI(I1,I2),如果s_MI>f_MI,则接受更新后的旋转平移矩阵[R,t];否则程序终止。
与现有技术相比,本发明的优点是:
可以方便的获得激光雷达以及双目摄像头之间的位姿关系,从而可以方便快速的将视觉信息以及激光雷达的信息进行融合。可以十分方便地对位姿进行配准,避免了需求精度较高的配准仿体,同时解决了基于强度的互信息配准的受环境影响比较大的缺点。具有很强的鲁棒性。
附图说明
下面结合附图及实施例对本发明作进一步描述:
图1为激光雷达和摄像头之间的位姿变换示意图;
图2为激光雷达的稀疏点直线拟合获取标定点示意图;
图3为本发明基于立体视觉的双目摄像头与激光雷达间的位姿配准方法的流程图;
图4为双目摄像头和激光雷达的互信息校正流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面结合具体实施方式并参照附图,对本发明进一步详细说明。应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
实施例:
下面结合附图,对本发明的较佳实施例作进一步说明。
如图3所示,一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,包括以下步骤:
通过双目摄像头模块对周围环境进行视觉感知,通过相机参数对左右图像进行极线校正。获取左右视觉图像,通过深度图获取算法,如SGM,GC等算法获取计算出来的深度图,双目摄像头和PC处理终端通过USB协议或者网口协议传输。
激光雷达感知周围环境,在相同的情况下同时采集激光雷达数据点,获取周围物体在激光雷达坐标系中的坐标位置。
通过将物体在激光雷达坐标系的坐标进行仿射变换,转变成物体在摄像机坐标系中的坐标,然后经过摄影变换转成物体在图像坐标系的位置。因为双目摄像头获取的视差图和物体在摄像机坐标系中的坐标Z存在简单的映射关系因此可以十分简单的使用此公式将雷达坐标转换成深度图,从而可以对摄像机深度图以及雷达深度图进行互信息配准。
互信息配准算法的配准流程如图4所示。
基于双目摄像头图像获取深度图I1
基于初始化的旋转平移矩阵[R,t]对雷达数据进行投影变换,获取激光雷达的深度图I2
通过两幅深度图进行互信息计算f_MI(I1,I2),通过梯度下降法获取更新后的旋转平移矩阵[R,t],并再次进行互信息计算s_MI(I1,I2),如果s_MI>f_MI,就接受更新后的旋转平移矩阵[R,t]。否则程序终止。互信息关系R,T在真值邻域附近是严格的凸函数,因此梯度下降法可以收敛到真值。
可以十分方便地对位姿进行配准,避免了需求精度较高的配准仿体,同时解决了基于强度的互信息配准的受环境影响比较大的缺点。
应当理解的是,本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理,而不构成对本发明的限制。因此,在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。此外,本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (3)

1.一种基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,其特征在于,包括以下步骤:
S01:通过将物体在激光雷达坐标系的坐标进行仿射变换,转变成物体在摄像机坐标系中的坐标,经过摄影变换转成物体在图像坐标系的位置;
S02:根据双目摄像头获取的视差图和物体在摄像机坐标系中的Z坐标的映射关系,将雷达坐标转换成深度图;
S03:将摄像机深度图以及雷达深度图进行互信息配准,对激光雷达的数据点和双目深度图进行位姿态配准。
2.根据权利要求1所述的基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,其特征在于,所述Z坐标映射关系为: ,其中d是视差,f是摄像头焦距,b是两个摄像头之间的基线距离。
3.根据权利要求1所述的基于立体视觉的双目摄像头与激光雷达间的位姿配准方法,其特征在于,所述互信息配准包括以下步骤:
S31:基于双目摄像头图像获取深度图 ;
S32:基于初始化的旋转平移矩阵[R,t]对雷达数据进行投影变换,获取激光雷达的深度图;
S33:通过两幅深度图进行互信息计算f_MI(I1,I2),通过梯度下降法获取更新后的旋转平移矩阵[R,t],并再次进行互信息计算s_MI(I1,I2),如果s_MI>f_MI,则接受更新后的旋转平移矩阵[R,t];否则程序终止。
CN201810801144.5A 2018-07-20 2018-07-20 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法 Active CN109035309B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810801144.5A CN109035309B (zh) 2018-07-20 2018-07-20 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810801144.5A CN109035309B (zh) 2018-07-20 2018-07-20 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法

Publications (2)

Publication Number Publication Date
CN109035309A true CN109035309A (zh) 2018-12-18
CN109035309B CN109035309B (zh) 2022-09-27

Family

ID=64644532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810801144.5A Active CN109035309B (zh) 2018-07-20 2018-07-20 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法

Country Status (1)

Country Link
CN (1) CN109035309B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109949371A (zh) * 2019-03-18 2019-06-28 北京智行者科技有限公司 一种用于激光雷达和相机数据的标定方法
CN110322519A (zh) * 2019-07-18 2019-10-11 天津大学 一种用于激光雷达与相机联合标定的标定装置及标定方法
CN110673115A (zh) * 2019-09-25 2020-01-10 杭州飞步科技有限公司 雷达与组合导航***的联合标定方法、装置、设备及介质
CN111458721A (zh) * 2020-03-31 2020-07-28 江苏集萃华科智能装备科技有限公司 一种暴露垃圾的识别定位方法、装置及***
CN112212784A (zh) * 2020-09-01 2021-01-12 长春工程学院 一种点激光位移传感器与双目相机坐标融合的方法及***
CN112312113A (zh) * 2020-10-29 2021-02-02 贝壳技术有限公司 用于生成三维模型的方法、装置和***
CN112379390A (zh) * 2020-11-18 2021-02-19 成都通甲优博科技有限责任公司 基于异源数据的位姿测量方法、装置、***及电子设备
CN112396662A (zh) * 2019-08-13 2021-02-23 杭州海康威视数字技术股份有限公司 一种转换矩阵修正的方法及装置
CN112598729A (zh) * 2020-12-24 2021-04-02 哈尔滨工业大学芜湖机器人产业技术研究院 融合激光与相机的目标物体识别与定位方法
CN113239948A (zh) * 2021-03-12 2021-08-10 复旦大学 一种面向毫米波雷达与视频图像的数据融合方法及***
CN113379831A (zh) * 2021-06-22 2021-09-10 北京航空航天大学青岛研究院 一种基于双目摄像头和人形机器人的增强现实方法
CN114782556A (zh) * 2022-06-20 2022-07-22 季华实验室 相机与激光雷达的配准方法、***及存储介质
CN115100298A (zh) * 2022-08-25 2022-09-23 青岛杰瑞工控技术有限公司 一种用于深远海视觉养殖的光-声图像融合方法
CN116106895A (zh) * 2023-02-21 2023-05-12 郑州大学 基于视觉与无线感知融合的多点微动智能监测***与方法
CN117152274A (zh) * 2023-11-01 2023-12-01 三一重型装备有限公司 掘进机双目摄像头的位姿校正方法及***、可读存储介质
WO2024045942A1 (zh) * 2022-09-02 2024-03-07 中兴通讯股份有限公司 环境信息感知方法、装置、***、计算机设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106204731A (zh) * 2016-07-18 2016-12-07 华南理工大学 一种基于双目立体视觉***的多视角三维重建方法
CN106596856A (zh) * 2016-12-09 2017-04-26 山东理工大学 一种基于激光雷达和摄影测量的车辆污染物排放量实时监测方法
CN107886477A (zh) * 2017-09-20 2018-04-06 武汉环宇智行科技有限公司 无人驾驶中立体视觉与低线束激光雷达的融合矫正方法
CN108196535A (zh) * 2017-12-12 2018-06-22 清华大学苏州汽车研究院(吴江) 基于增强学习和多传感器融合的自动驾驶***

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106204731A (zh) * 2016-07-18 2016-12-07 华南理工大学 一种基于双目立体视觉***的多视角三维重建方法
CN106596856A (zh) * 2016-12-09 2017-04-26 山东理工大学 一种基于激光雷达和摄影测量的车辆污染物排放量实时监测方法
CN107886477A (zh) * 2017-09-20 2018-04-06 武汉环宇智行科技有限公司 无人驾驶中立体视觉与低线束激光雷达的融合矫正方法
CN108196535A (zh) * 2017-12-12 2018-06-22 清华大学苏州汽车研究院(吴江) 基于增强学习和多传感器融合的自动驾驶***

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109949371A (zh) * 2019-03-18 2019-06-28 北京智行者科技有限公司 一种用于激光雷达和相机数据的标定方法
CN110322519A (zh) * 2019-07-18 2019-10-11 天津大学 一种用于激光雷达与相机联合标定的标定装置及标定方法
CN110322519B (zh) * 2019-07-18 2023-03-31 天津大学 一种用于激光雷达与相机联合标定的标定装置及标定方法
CN112396662A (zh) * 2019-08-13 2021-02-23 杭州海康威视数字技术股份有限公司 一种转换矩阵修正的方法及装置
CN112396662B (zh) * 2019-08-13 2024-05-24 杭州海康威视数字技术股份有限公司 一种转换矩阵修正的方法及装置
CN110673115A (zh) * 2019-09-25 2020-01-10 杭州飞步科技有限公司 雷达与组合导航***的联合标定方法、装置、设备及介质
CN111458721A (zh) * 2020-03-31 2020-07-28 江苏集萃华科智能装备科技有限公司 一种暴露垃圾的识别定位方法、装置及***
CN112212784B (zh) * 2020-09-01 2022-02-08 长春工程学院 一种点激光位移传感器与双目相机坐标融合的方法及***
CN112212784A (zh) * 2020-09-01 2021-01-12 长春工程学院 一种点激光位移传感器与双目相机坐标融合的方法及***
CN112312113A (zh) * 2020-10-29 2021-02-02 贝壳技术有限公司 用于生成三维模型的方法、装置和***
CN112379390A (zh) * 2020-11-18 2021-02-19 成都通甲优博科技有限责任公司 基于异源数据的位姿测量方法、装置、***及电子设备
CN112598729B (zh) * 2020-12-24 2022-12-23 哈尔滨工业大学芜湖机器人产业技术研究院 融合激光与相机的目标物体识别与定位方法
CN112598729A (zh) * 2020-12-24 2021-04-02 哈尔滨工业大学芜湖机器人产业技术研究院 融合激光与相机的目标物体识别与定位方法
CN113239948A (zh) * 2021-03-12 2021-08-10 复旦大学 一种面向毫米波雷达与视频图像的数据融合方法及***
CN113379831B (zh) * 2021-06-22 2022-09-09 北京航空航天大学青岛研究院 一种基于双目摄像头和人形机器人的增强现实方法
CN113379831A (zh) * 2021-06-22 2021-09-10 北京航空航天大学青岛研究院 一种基于双目摄像头和人形机器人的增强现实方法
CN114782556B (zh) * 2022-06-20 2022-09-09 季华实验室 相机与激光雷达的配准方法、***及存储介质
CN114782556A (zh) * 2022-06-20 2022-07-22 季华实验室 相机与激光雷达的配准方法、***及存储介质
CN115100298B (zh) * 2022-08-25 2022-11-29 青岛杰瑞工控技术有限公司 一种用于深远海视觉养殖的光-声图像融合方法
CN115100298A (zh) * 2022-08-25 2022-09-23 青岛杰瑞工控技术有限公司 一种用于深远海视觉养殖的光-声图像融合方法
WO2024045942A1 (zh) * 2022-09-02 2024-03-07 中兴通讯股份有限公司 环境信息感知方法、装置、***、计算机设备及存储介质
CN116106895A (zh) * 2023-02-21 2023-05-12 郑州大学 基于视觉与无线感知融合的多点微动智能监测***与方法
CN116106895B (zh) * 2023-02-21 2024-01-26 郑州大学 基于视觉与无线感知融合的多点微动智能监测***与方法
CN117152274A (zh) * 2023-11-01 2023-12-01 三一重型装备有限公司 掘进机双目摄像头的位姿校正方法及***、可读存储介质
CN117152274B (zh) * 2023-11-01 2024-02-09 三一重型装备有限公司 掘进机双目摄像头的位姿校正方法及***、可读存储介质

Also Published As

Publication number Publication date
CN109035309B (zh) 2022-09-27

Similar Documents

Publication Publication Date Title
CN109035309A (zh) 基于立体视觉的双目摄像头与激光雷达间的位姿配准方法
CN108734744B (zh) 一种基于全站仪的远距离大视场双目标定方法
CN109285145B (zh) 基于智能手机的多株立木高度测量方法
KR101666959B1 (ko) 카메라로부터 획득한 영상에 대한 자동보정기능을 구비한 영상처리장치 및 그 방법
CN106960454B (zh) 景深避障方法、设备及无人飞行器
CN106408601B (zh) 一种基于gps的双目融合定位方法及装置
CN107424196B (zh) 一种基于弱标定多目相机的立体匹配方法、装置及***
CN109579695B (zh) 一种基于异构立体视觉的零件测量方法
TW201717613A (zh) 一種全景魚眼相機影像校正、合成與景深重建方法與其系統
CN106570899B (zh) 一种目标物体检测方法及装置
CN105844692B (zh) 基于双目立体视觉的三维重建装置、方法、***及无人机
García-Moreno et al. LIDAR and panoramic camera extrinsic calibration approach using a pattern plane
KR20050058085A (ko) 변이지도 및 깊이지도의 융합을 통한 3차원 장면 모델생성 장치 및 그 방법
CN109840922B (zh) 基于双目光场相机的深度获取方法和***
CN111220126A (zh) 一种基于点特征和单目相机的空间物***姿测量方法
CN109727290A (zh) 基于单目视觉三角测距法的变焦相机动态标定方法
JP4132068B2 (ja) 画像処理装置及び三次元計測装置並びに画像処理装置用プログラム
CN110779491A (zh) 一种水平面上目标测距的方法、装置、设备及存储介质
CN108279677B (zh) 基于双目视觉传感器的轨道机器人检测方法
CN103852060A (zh) 一种基于单目视觉的可见光图像测距方法
CN110827392A (zh) 具好的场景易用性的单目图像三维重建方法、***及装置
WO2018032841A1 (zh) 绘制三维图像的方法及其设备、***
CN111009030A (zh) 一种多视高分辨率纹理图像与双目三维点云映射方法
CN105513074B (zh) 一种羽毛球机器人相机以及车身到世界坐标系的标定方法
CN111798507A (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