CN107290755A - 基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 - Google Patents
基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 Download PDFInfo
- Publication number
- CN107290755A CN107290755A CN201710487373.XA CN201710487373A CN107290755A CN 107290755 A CN107290755 A CN 107290755A CN 201710487373 A CN201710487373 A CN 201710487373A CN 107290755 A CN107290755 A CN 107290755A
- Authority
- CN
- China
- Prior art keywords
- signal
- target
- photon counting
- echo
- local oscillation
- 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
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/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/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,涉及激光雷达技术领域。解决了现有目标距离和目标强度的获取方法因距离解算计算量大,导致的耗时长及解算强度准确率低的问题。本发明首先,通过对本振信号改造,使其与回波信号具有相同的形式,再将改造后的本振信号和Gm‑APD探测结果混频处理,获得相关峰谱。一方面通过搜索相关峰谱中波峰极大值的位置,解算出距离信息;一方面通过信号发生器输出的本振信号中光子计数脉冲的个数M和Gm‑APD单光子探测器探测到的回波信号中光子计数脉冲的个数M′进行处理,解算出强度信息。本发明主要用于对目标进行距离和强度进行探测。
Description
技术领域
本发明涉及激光雷达技术领域。
背景技术
由于Gm-APD的极高响应灵敏度和极快的光电响应速度,因此基于Gm-APD探测器的激光雷达系统在远距离和微弱信号的应用领域中发挥了重要的作用。但是目前军事和民用的各个领域应用中,仅靠距离像已经很难满足很多复杂的目标的探测要求了,然而由于Gm-APD雪崩效应的0和1的逻辑输出,只能响应信号到达的时间,无法响应信号的强度信息。这极大的限制了基于Gm-APD探测器的激光雷达系统的应用前景。人们为了解决这一问题,探索了很多方法。
其中,比较主要的一种方法是复合调制脉冲编码的4D成像光子计数激光雷达系统,该方法采用了独特的复合调制脉冲编码发射,在脉冲位置调制的基础上加入了脉冲幅度的调制,然后通过回波脉冲序列和发射脉冲序列的相关处理可以得到目标的距离信息,通过不同幅度的脉冲探测结果的统计解算可以获得目标的强度信息。
这种方法虽然实现了距离像和强度像的同时获取,但这个数据处理方法需要将回波信号中不同强度的脉冲进行分别统计,这一方面计算量大,比较耗费时间,另一方法由于分辨的错误也会造成解算强度的误差。因此,亟需提供一种耗时短和解算强度准确率高的距离像和强度像的获取方法。
发明内容
本发明是为了解决现有目标距离和目标强度的获取方法因距离解算计算量大,导致的耗时长及解算强度准确率低的问题,本发明提供了一种基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法。
基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,所述的4D成像光子计数激光雷达系统包括信号发生器、激光器、发射光学系统、单向反射器、扫描器、全反射镜、接收光学系统、窄带滤波片、Gm-APD单光子探测器和信号处理模块;
所述的信号发生器的复合调制随机脉冲信号输出端同时与激光器的复合信号输入端和信号处理模块的复合信号输入端连接,激光器的脉冲信号输出端与发射光学系统的脉冲信号输入端连接,发射光学系统输出的光信号经单向反射器透射后入射至扫描器,扫描器输出光探测信号至目标,经目标反射后的回波信号入射至该扫描器,扫描器输出的回波信号依次经单向反射器和全反射镜的反射后,入射至接收光学系统进行汇聚回波信号,被汇聚的回波信号经过窄带滤波片滤波后,最后由Gm-APD单光子探测器进行光子计数探测,探测结果输入到信号处理模块,信号处理模块对接收的信号进行处理,从而得到目标的距离和强度信息;
其中,信号发生器的复合调制随机脉冲信号输出端输出的信号为本振信号;
目标距离和目标强度的获取方法包括如下步骤:
步骤一,信号处理模块对本振信号采用回波高斯脉冲波形函数进行处理,获得改造后的本振信号;
步骤二,对改造后的本振信号和Gm-APD单光子探测器输出的光子计数探测结果进行时间上的自相关处理,使改造后的本振信号与光子计数探测结果进行混频,从而获得相关峰谱;
步骤三,采用峰值估算方法对相关峰谱中的波峰极大值进行搜索,获得相关峰谱中波峰极大值所对应的时间,该时间为回波信号的往返时间τ,将回波信号的往返时间τ代入雷达距离方程中,从而获得目标的距离;
步骤四,信号处理模块对信号发生器输出的本振信号中光子计数脉冲的个数M和Gm-APD单光子探测器探测到的回波信号中光子计数脉冲的个数M′进行处理,获得本振信号的响应概率P=M′/M;
对本振信号的响应概率P进行强度解算,从而获得目标的强度信息Ns。
所述的强度解算采用如下公式实现:
Ns=-ln(1-P)。
所述的雷达距离方程的表达式为:R=cτ/2;其中,R为目标的距离值,c为光速。
原理分析:首先,通过对本振信号改造,使其与回波信号具有相同的形式,再将改造后的本振信号和Gm-APD探测结果混频处理,获得相关峰谱。一方面通过搜索相关峰谱中波峰的位置,解算出距离信息;一方面通过信号发生器输出的本振信号中光子计数脉冲的个数M和Gm-APD单光子探测器探测到的回波信号中光子计数脉冲的个数M′进行处理,解算出强度信息。
本发明带来的有益效果是,通过对本振信号的改造,使其与回波信号具有相同的形式,从而提高相关处理的效率,提高距离解算算法的精度,目标距离的测量结果的精度提高了20%以上。通过对强度解算算法进行改进,极大的简化了解算的过程,从而避免了多余环节内引入的误差,有助于强度解算精度的提高,强度解算精度提高了30%以上。
附图说明
图1为本发明所述的4D成像光子计数激光雷达系统的原理示意图;
图2为本发明所述的基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法的流程图;
图3为距离信息解算的相关峰谱示意图;其中,Tpluse表示发射脉冲的宽度,I表示信号脉冲的强度,t表示时间,τ表示回波信号的往返延迟时间,图3a为本振信号为矩形脉冲的信号波形,图3b为改造后的本振信号的波形,图3c为Gm-APD单光子探测器输出的探测结果的波形,图3d为相关峰谱的波形。
具体实施方式
具体实施方式一:参见图1说明本实施方式,本实施方式所述的基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,所述的4D成像光子计数激光雷达系统包括信号发生器1、激光器2、发射光学系统3、单向反射器4、扫描器5、全反射镜6、接收光学系统7、窄带滤波片8、Gm-APD单光子探测器9、信号处理模块10;
所述的信号发生器1的复合调制随机脉冲信号输出端同时与激光器2的复合信号输入端和信号处理模块10的复合信号输入端连接,激光器2的脉冲信号输出端与发射光学系统3的脉冲信号输入端连接,发射光学系统3输出的光信号经单向反射器4透射后入射至扫描器5,扫描器5输出光探测信号至目标,经目标反射后的回波信号入射至该扫描器5,扫描器5输出的回波信号依次经单向反射器4和全反射镜6的反射后,入射至接收光学系统7进行汇聚回波信号,被汇聚的回波信号经过窄带滤波片8滤波后,最后由Gm-APD单光子探测器9进行光子计数探测,探测结果输入到信号处理模块10,信号处理模块10对接收的信号进行处理,从而得到目标的距离和强度信息;
其中,信号发生器1的复合调制随机脉冲信号输出端输出的信号为本振信号;
目标距离和目标强度的获取方法包括如下步骤:
步骤一,信号处理模块10对本振信号采用回波高斯脉冲波形函数进行处理,获得改造后的本振信号;
步骤二,对改造后的本振信号和Gm-APD单光子探测器9输出的光子计数探测结果进行时间上的自相关处理,使改造后的本振信号与光子计数探测结果进行混频,从而获得相关峰谱;
步骤三,采用峰值估算方法对相关峰谱中的波峰极大值进行搜索,获得相关峰谱中波峰极大值所对应的时间,该时间为回波信号的往返时间τ,将回波信号的往返时间τ代入雷达距离方程中,从而获得目标的距离;
步骤四,信号处理模块10对信号发生器1输出的本振信号中光子计数脉冲的个数M和Gm-APD单光子探测器9探测到的回波信号中光子计数脉冲的个数M′进行处理,获得本振信号的响应概率P=M′/M;
对本振信号的响应概率P进行强度解算,从而获得目标的强度信息Ns。
具体通过图3说明本实施方式,本振信号的矩形脉冲(如图3a)改造成和回波信号具有相同形式的本振信号(如图3b);然后使用改造后的本振信号和Gm-APD单光子探测器9的探测结果进行相关,从而可以有效的提高相关的效果,有效的缩窄相关峰,从而有效的提高测距精度。
具体实施方式二:参见图1说明本实施方式,本实施方式与具体实施方式一所述的一种基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法的区别在于,所述的强度解算采用如下公式实现:
Ns=-ln(1-P)。
本实施方式,如图3(a)本振信号由多种等差强度的脉冲随机分布组成,我们使用1,2…i…n来表示不同强度的单个脉冲,1,2…i…n叫做编码,表示不同强度的脉冲。一次探测使用的是脉冲序列,包括M个脉冲,这个脉冲序列里包含着n种不同强度调制的脉冲随机分布组成。n种不同强度的脉冲是等概率出现的,因此,对于一次探测的脉冲序列,M个各种强度脉冲,编码为1,2…i…n的脉冲个数是相等的分别为M/n。我们设最大强度的回波脉冲的光电子形式信号强度为Ns,也就是说编码为n的脉冲强度为Ns,这样编码为i(i∈[1,n])的回波脉冲,信号的强度为(iNs)/n,那么现在来计算编码为i的脉冲触发Gm-APD产生雪崩事件概率
由于本振信号中包括n种不同强度的脉冲,n种不同强度的脉冲中,每种强度的脉冲出现的概率是均等的,均为因此,总的响应概率应该为本振信号中n种不同强度脉冲响应概率的加权平均:
将公式一代入公式二中,并进行整理,从而获得Ns=-ln(1-P)。
具体实施方式三:参见图1说明本实施方式,本实施方式与具体实施方式一所述的一种基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法的区别在于,所述的雷达距离方程的表达式为:R=cτ/2;其中,R为目标的距离值,c为光速。
Claims (3)
1.基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,所述的4D成像光子计数激光雷达系统包括信号发生器(1)、激光器(2)、发射光学系统(3)、单向反射器(4)、扫描器(5)、全反射镜(6)、接收光学系统(7)、窄带滤波片(8)、Gm-APD单光子探测器(9)和信号处理模块(10);
所述的信号发生器(1)的复合调制随机脉冲信号输出端同时与激光器(2)的复合信号输入端和信号处理模块(10)的复合信号输入端连接,激光器(2)的脉冲信号输出端与发射光学系统(3)的脉冲信号输入端连接,发射光学系统(3)输出的光信号经单向反射器(4)透射后入射至扫描器(5),扫描器(5)输出光探测信号至目标,经目标反射后的回波信号入射至该扫描器(5),扫描器(5)输出的回波信号依次经单向反射器(4)和全反射镜(6)的反射后,入射至接收光学系统(7)进行汇聚回波信号,被汇聚的回波信号经过窄带滤波片(8)滤波后,最后由Gm-APD单光子探测器(9)进行光子计数探测,探测结果输入到信号处理模块(10),信号处理模块(10)对接收的信号进行处理,从而得到目标的距离和强度信息;
其中,信号发生器(1)的复合调制随机脉冲信号输出端输出的信号为本振信号;
其特征在于,目标距离和目标强度的获取方法包括如下步骤:
步骤一,信号处理模块(10)对本振信号采用回波高斯脉冲波形函数进行处理,获得改造后的本振信号;
步骤二,对改造后的本振信号和Gm-APD单光子探测器(9)输出的光子计数探测结果进行时间上的自相关处理,使改造后的本振信号与光子计数探测结果进行混频,从而获得相关峰谱;
步骤三,采用峰值估算方法对相关峰谱中的波峰极大值进行搜索,获得相关峰谱中波峰极大值所对应的时间,该时间为回波信号的往返时间τ,将回波信号的往返时间τ代入雷达距离方程中,从而获得目标的距离;
步骤四,信号处理模块(10)对信号发生器(1)输出的本振信号中光子计数脉冲的个数M和Gm-APD单光子探测器(9)探测到的回波信号中光子计数脉冲的个数M′进行处理,获得本振信号的响应概率P=M′/M;
对本振信号的响应概率P进行强度解算,从而获得目标的强度信息Ns。
2.根据权利要求1所述的基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,其特征在于,所述的强度解算采用如下公式实现:
Ns=-ln(1-P)。
3.根据权利要求1所述的基于4D成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法,其特征在于,所述的雷达距离方程的表达式为:R=cτ/2;其中,R为目标的距离值,c为光速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710487373.XA CN107290755B (zh) | 2017-06-23 | 2017-06-23 | 基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710487373.XA CN107290755B (zh) | 2017-06-23 | 2017-06-23 | 基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107290755A true CN107290755A (zh) | 2017-10-24 |
CN107290755B CN107290755B (zh) | 2020-06-09 |
Family
ID=60098251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710487373.XA Active CN107290755B (zh) | 2017-06-23 | 2017-06-23 | 基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107290755B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445471A (zh) * | 2018-03-26 | 2018-08-24 | 武汉大学 | 一种单光子激光雷达多探测器条件下的测距精度评估方法 |
CN110058255A (zh) * | 2019-05-10 | 2019-07-26 | 中国工程物理研究院电子工程研究所 | 方位探测装置及方法 |
CN110161521A (zh) * | 2019-06-11 | 2019-08-23 | 中国科学院光电技术研究所 | 一种基于真随机编码的光子计数激光雷达 |
CN110554404A (zh) * | 2019-08-08 | 2019-12-10 | 华中科技大学 | 一种强背景噪声下Gm-APD阵列激光雷达成像方法及系统 |
CN111060887A (zh) * | 2019-11-05 | 2020-04-24 | 哈尔滨工业大学 | 一种基于凹凸搜索的Gm-APD激光雷达低信噪比回波数据信号提取方法 |
CN112965078A (zh) * | 2021-02-04 | 2021-06-15 | 郜键 | 一种基于msm探测的amcw远距离激光成像装置 |
CN115840215A (zh) * | 2023-02-13 | 2023-03-24 | 宁波飞芯电子科技有限公司 | 探测方法和探测设备 |
CN117590353A (zh) * | 2024-01-19 | 2024-02-23 | 山东省科学院海洋仪器仪表研究所 | 一种光子计数激光雷达的弱回波信号快速提取及成像方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137450A (zh) * | 2015-08-10 | 2015-12-09 | 哈尔滨工业大学 | 低虚警双Gm-APD探测器光子计数激光雷达 |
CN105182361A (zh) * | 2015-08-06 | 2015-12-23 | 哈尔滨工业大学 | 一种基于复合调制脉冲编码的4d成像光子计数激光雷达 |
-
2017
- 2017-06-23 CN CN201710487373.XA patent/CN107290755B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182361A (zh) * | 2015-08-06 | 2015-12-23 | 哈尔滨工业大学 | 一种基于复合调制脉冲编码的4d成像光子计数激光雷达 |
CN105137450A (zh) * | 2015-08-10 | 2015-12-09 | 哈尔滨工业大学 | 低虚警双Gm-APD探测器光子计数激光雷达 |
Non-Patent Citations (3)
Title |
---|
ZIJING ZHANG ET AL.: "Premixing photon-counting chirped amplitude modulation lidar for range and velocity measurement in photon starved scenes", 《OPTIK》 * |
张子静: "基于概率统计的光子激光雷达性能提高的理论与实验研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
赵远: "高灵敏度的光子偏振激光雷达系统", 《红外与激光工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445471B (zh) * | 2018-03-26 | 2019-08-02 | 武汉大学 | 一种单光子激光雷达多探测器条件下的测距精度评估方法 |
CN108445471A (zh) * | 2018-03-26 | 2018-08-24 | 武汉大学 | 一种单光子激光雷达多探测器条件下的测距精度评估方法 |
CN110058255A (zh) * | 2019-05-10 | 2019-07-26 | 中国工程物理研究院电子工程研究所 | 方位探测装置及方法 |
CN110161521B (zh) * | 2019-06-11 | 2022-07-19 | 中国科学院光电技术研究所 | 一种基于真随机编码的光子计数激光雷达 |
CN110161521A (zh) * | 2019-06-11 | 2019-08-23 | 中国科学院光电技术研究所 | 一种基于真随机编码的光子计数激光雷达 |
CN110554404A (zh) * | 2019-08-08 | 2019-12-10 | 华中科技大学 | 一种强背景噪声下Gm-APD阵列激光雷达成像方法及系统 |
CN111060887A (zh) * | 2019-11-05 | 2020-04-24 | 哈尔滨工业大学 | 一种基于凹凸搜索的Gm-APD激光雷达低信噪比回波数据信号提取方法 |
CN111060887B (zh) * | 2019-11-05 | 2022-02-22 | 哈尔滨工业大学 | 一种基于凹凸搜索的Gm-APD激光雷达低信噪比回波数据信号提取方法 |
CN112965078A (zh) * | 2021-02-04 | 2021-06-15 | 郜键 | 一种基于msm探测的amcw远距离激光成像装置 |
CN112965078B (zh) * | 2021-02-04 | 2024-02-09 | 苏州奥瑞图光电科技有限公司 | 一种基于msm探测的amcw远距离激光成像装置 |
CN115840215A (zh) * | 2023-02-13 | 2023-03-24 | 宁波飞芯电子科技有限公司 | 探测方法和探测设备 |
CN117590353A (zh) * | 2024-01-19 | 2024-02-23 | 山东省科学院海洋仪器仪表研究所 | 一种光子计数激光雷达的弱回波信号快速提取及成像方法 |
CN117590353B (zh) * | 2024-01-19 | 2024-03-29 | 山东省科学院海洋仪器仪表研究所 | 一种光子计数激光雷达的弱回波信号快速提取及成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107290755B (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107290755A (zh) | 基于4d成像光子计数激光雷达系统实现的目标距离和目标强度的获取方法 | |
CN105182361A (zh) | 一种基于复合调制脉冲编码的4d成像光子计数激光雷达 | |
CN109343069B (zh) | 可实现组合脉冲测距的光子计数激光雷达及其测距方法 | |
CN105425244B (zh) | 前置混频的啁啾调制光子计数激光雷达 | |
CN108089194B (zh) | 一种基于复合伪随机编码的光子计数激光雷达 | |
CN108415010B (zh) | 一种基于梯形lfmcw调制的雷达多目标检测方法 | |
JP5138854B2 (ja) | 光学距離測定 | |
US7312856B2 (en) | Programmable pulse capture device with automatic gain control | |
CN103616696A (zh) | 一种激光成像雷达装置及其测距的方法 | |
CN108008402B (zh) | 一种用于激光测距的单光子雪崩二极管探测器阵列 | |
CN104457452A (zh) | 一种伪随机码体制激光引信系统及其目标识别方法 | |
CN105487067A (zh) | 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 | |
CN110221308A (zh) | 一种相干脉冲激光测距的方法、相关装置及存储介质 | |
CN108519604B (zh) | 一种基于伪随机码调制解调的固态面阵激光雷达测距方法 | |
CN110646779A (zh) | 相位编码不饱和调制方法、装置、激光雷达测距测速方法及激光雷达系统 | |
CN110161521A (zh) | 一种基于真随机编码的光子计数激光雷达 | |
CN107807353A (zh) | 采用n重脉冲编码的三维成像光子计数系统及计数方法 | |
CN106371108A (zh) | 一种全光纤脉冲相干多普勒激光测风方法及雷达系统 | |
Feng et al. | FPGA-based digital chaotic anti-interference lidar system | |
CN113640752B (zh) | 一种基于脉间相位频谱双捷变的波形设计方法 | |
US20230050937A1 (en) | Detection method and detection apparatus | |
CN112034483B (zh) | 基于相干探测关联成像技术的多目标距离--速度解耦方法 | |
CN117590418A (zh) | 一种基于扩频体制的激光回波探距方法和系统 | |
CN106646422A (zh) | 增强相干测风雷达多普勒频移信号信噪比的预处理系统 | |
CN116148815A (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 |