CN110824456B - 一种自适应分辨率三维激光扫描方法 - Google Patents
一种自适应分辨率三维激光扫描方法 Download PDFInfo
- Publication number
- CN110824456B CN110824456B CN201911068898.5A CN201911068898A CN110824456B CN 110824456 B CN110824456 B CN 110824456B CN 201911068898 A CN201911068898 A CN 201911068898A CN 110824456 B CN110824456 B CN 110824456B
- Authority
- CN
- China
- Prior art keywords
- horizontal
- frame
- scanning
- resolution
- vertical
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000010606 normalization Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 2
- 230000006978 adaptation Effects 0.000 claims 1
- 230000004927 fusion Effects 0.000 claims 1
- 230000001502 supplementing effect Effects 0.000 abstract description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/481—Constructional features, e.g. arrangements of optical elements
- G01S7/4817—Constructional features, e.g. arrangements of optical elements relating to scanning
-
- 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
-
- 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
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Mechanical Optical Scanning Systems (AREA)
Abstract
Description
技术领域
本发明涉及三维扫描技术领域,尤其是一种自适应分辨率三维激光扫描方法。
背景技术
激光视觉系统的理论研究在现实生活中被广泛应用,如无人驾驶汽车,机器人导航,建筑物重建和遥感。激光雷达系统通过估计从激光雷达发出的光和从远处物体反射的光的到达时间间隔,来测量激光发射点到物体的距离,相比于图像,其不受光照影响有其独特优势。
现有的激光雷达产品都是固定分辨率、均匀扫描的。这种方式对于结构简单的物体如墙面、平地等会浪费时间并获得大量冗余数据,对于复杂的对象则某些分辨率可能无法获得精确的细节信息。
发明内容
针对现有技术的不足,本发明提供一种自适应分辨率三维激光扫描方法,本发明利用利用人造物体本身的垂直特征,采用固定垂直分辨率、自适应水平分辨率的扫描方式,可以较快速的获得结构简单的物体信息,对于复杂物体采用高分辨率进行扫描,在扫描分辨率和扫描时间(数据量)之间取得了平衡。
本发明的技术方案为:一种自适应分辨率三维激光扫描方法,包括以下步骤:
S1)、对环境采用垂直分辨率ω、水平分辨率τ进行粗扫描,获得激光雷达数据P;
P={pi,j},i=1,2,...n,j=1,2,...,m;
其中,n是指激光雷达水平总帧数,m是指激光雷达垂直总帧数;pi,j=(αi,βj,di,j)为水平第i帧、垂直第j帧的扫描点,其中,αi为水平第i帧角度、βj为垂直第j帧角度、di,j为水平第i帧、垂直第j帧时测得的距离;
S3)、计算不同水平角度区间对应的水平扫描分辨率μx;
其中,μx为水平角度范围[αx,αx+L-1],x=1,2,...,n-L+1之间的水平扫描分辨率;
S4)、对不同水平角度区间采用垂直分辨率ω、水平扫描分辨率μx进行精细扫描得到P'={p'i,j},i=1,2,...n',j=1,2,...,m;
其中,n'是指激光雷达水平总帧数,m是指激光雷达垂直总帧数,p'i,j=(α'i,βj,d'i,j)为水平第i帧、垂直第j帧的扫描点,其中α'i为水平第i帧角度、βj为垂直第j帧角度、d'i,j为水平第i帧、垂直第j帧时测得的距离;
S5)、以精细扫描数据P'来补充粗扫描数据P,形成融合后的数据Pmix:
其中,为的离散时间傅里叶变换,是粗扫描数据水平总帧数n被长度为L的矩形窗截取的部分的距离数据组成的矩阵,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1,即水平角度范围[αx,αx+L-1],垂直帧数范围[1,m]所对应的距离,αx为水平第x帧的角度;STFT为短时傅里叶变换,DTFT为离散时间傅里叶变换,L为短时傅里叶变换窗体的宽度;
S202)、计算中高频分量所占比重高频比即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的扫描区域的物体复杂度,L为短时傅里叶变换窗体的宽度。
优选的,上述方法中,步骤S3)中,计算不同水平角度区间对应的水平扫描分辨率μx,具体如下:
S301)、构建如下函数:
其中,S为代价函数,μx为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的区域对应的水平扫描分辨率,τ是粗扫描的水平扫描分辨率,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的扫描区域的物体复杂度,c为复杂度归一化常数,L为短时傅里叶变换窗体的宽度;
S302)、令代价函数S取最小值,求出复杂度归一化常数c,然后求解出不同水平角度区间对应的水平扫描分辨率μx。
优选的,步骤S5)中,若精细扫描数据P'的水平第x帧的角度α'x与粗扫描数据P的水平第y帧的角度αy相等,则粗扫描数据P水平第y帧的距离dy,j为:
dy,j=(dy,j+d'x,j)/2,j=1,2,...,m;
若精细扫描数据P'的水平第x帧的角度α'x在粗扫描数据P中无对应角度,即α'x≠αi,i=1,2,...,n,则向粗扫描数据P中添加数据p'x,j=(α'x,βj,d'x,j),j=1,2,...,m,经过上述变化的粗扫描数据P就变为融合后的数据Pmix。
优选的,步骤S5)中,对融合后的数据Pmix进行插值滤波后处理。
本发明的有益效果为:对于未知复杂度的环境,自适应分辨率扫描可以实现对环境中简单物体的快速扫描,对于环境中的复杂物体,通过精细扫描可以获得复杂物体的细节信息,在扫描分辨率和扫描时间(数据量)之间取得了平衡。
附图说明
图1为本发明的流程示意图;
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
如图1所示,一种自适应分辨率三维激光扫描方法,包括以下步骤:
S1)、对环境采用垂直分辨率ω、水平分辨率τ进行粗扫描,获得激光雷达数据P;
P={pi,j},i=1,2,...n,j=1,2,...,m;
其中,n是指激光雷达水平总帧数,m是指激光雷达垂直总帧数;pi,j=(αi,βj,di,j)为水平第i帧、垂直第j帧的扫描点,其中,αi为水平第i帧角度、βj为垂直第j帧角度、di,j为水平第i帧、垂直第j帧时测得的距离;
其中,为的离散时间傅里叶变换,是粗扫描数据水平总帧数n被长度为L的矩形窗截取的部分的距离数据组成的矩阵,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1,即水平角度范围[αx,αx+L-1],垂直帧数范围[1,m]所对应的距离,αx为水平第x帧的角度;STFT为短时傅里叶变换,DTFT为离散时间傅里叶变换,L为短时傅里叶变换窗体的宽度;
S202)、计算中高频分量所占比重高频比即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的扫描区域的物体复杂度,L为短时傅里叶变换窗体的宽度。
S3)、计算不同水平角度区间对应的水平扫描分辨率μx;其中,μx为水平角度范围[αx,αx+L-1],x=1,2,...,n-L+1之间的水平扫描分辨率;具体如下:
S301)、构建如下函数:
其中,S为代价函数,μx为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的区域对应的水平扫描分辨率,τ是粗扫描的水平扫描分辨率,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1(即水平角度范围[αx,αx+L-1])之间的扫描区域的物体复杂度,c为复杂度归一化常数,L为短时傅里叶变换窗体的宽度;
S302)、令代价函数S取最小值,求出复杂度归一化常数c,然后求解出不同水平角度区间对应的水平扫描分辨率μx。
S4)、对不同水平角度区间采用垂直分辨率ω、水平扫描分辨率μx进行精细扫描得到P'={p'i,j},i=1,2,...n',j=1,2,...,m;
其中,n'是指激光雷达水平总帧数,m是指激光雷达垂直总帧数,p'i,j=(α'i,βj,d'i,j)为水平第i帧、垂直第j帧的扫描点,其中α'i为水平第i帧角度、βj为垂直第j帧角度、d'i,j为水平第i帧、垂直第j帧时测得的距离;
S5)、以精细扫描数据P'来补充粗扫描数据P,形成融合后的数据Pmix,对融合后的数据Pmix进行插值滤波后处理:
若精细扫描数据P'的水平第x帧的角度α'x与粗扫描数据P的水平第y帧的角度αy相等,则粗扫描数据P水平第y帧的距离dy,j为:
dy,j=(dy,j+d'x,j)/2,j=1,2,...,m;
若精细扫描数据P'的水平第x帧的角度α'x在粗扫描数据P中无对应角度,即α'x≠αi,i=1,2,...,n,则向粗扫描数据P中添加数据p'x,j=(α'x,βj,d'x,j),j=1,2,...,m,经过上述变化的粗扫描数据P就变为融合后的数据Pmix。
从上述实施例中可以看出,对于未知复杂度的环境,自适应分辨率扫描可以实现对环境中简单物体的快速扫描,对于环境中的复杂物体,通过精细扫描可以获得复杂物体的细节信息,在扫描分辨率和扫描时间(数据量)之间取得了平衡。
上述实施例和说明书中描述的只是说明本发明的原理和最佳实施例,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。
Claims (7)
1.一种自适应分辨率三维激光扫描方法,其特征在于,包括以下步骤:
S1)、对环境采用垂直分辨率ω、水平分辨率τ进行粗扫描,获得激光雷达数据P;
P={pi,j},i=1,2,...n,j=1,2,...,m;
其中,n是指激光雷达水平总帧数,m是指激光雷达垂直总帧数,pi,j为水平第i帧、垂直第j帧的扫描点;
S3)、计算不同水平角度区间对应的水平扫描分辨率μx;具体如下:
S301)、构建如下函数:
其中,S为代价函数,μx为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1,即水平角度范围[αx,αx+L-1]之间的区域对应的水平扫描分辨率,τ是粗扫描的水平扫描分辨率,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1之间的扫描区域的物体复杂度,c为复杂度归一化常数,L为短时傅里叶变换窗体的宽度;
S302)、令代价函数S取最小值,求出复杂度归一化常数c,然后求解出不同水平角度区间对应的水平扫描分辨率μx;
S4)、对不同水平角度区间采用垂直分辨率ω、水平扫描分辨率μx进行精细扫描得到P';
P'={p′i,j},i=1,2,...n',j=1,2,...,m;
其中,n'是指激光雷达水平总帧数,m是指激光雷达垂直总帧数,p′i,j为水平第i帧、垂直第j帧的扫描点;
S5)、以精细扫描数据P'来补充粗扫描数据P,形成融合后的数据Pmix:
其中,为的离散时间傅里叶变换,是粗扫描数据水平总帧数n被长度为L的矩形窗截取的部分的距离数据组成的矩阵,即为水平帧数范围[x,x+L-1],x=1,2,...,n-L+1,即水平角度范围[αx,αx+L-1],垂直帧数范围[1,m]所对应的距离,αx为水平第x帧的角度;STFT为短时傅里叶变换,DTFT为离散时间傅里叶变换,L为短时傅里叶变换窗体的宽度;
5.根据权利要求1所述的一种自适应分辨率三维激光扫描方法,其特征在于:步骤S5)中,若精细扫描数据P'的水平第x帧的角度α'x与粗扫描数据P的水平第y帧的角度αy相等,则粗扫描数据P水平第y帧的距离dy,j为:
dy,j=(dy,j+d'x,j)/2,j=1,2,...,m;
其中,d'x,j为水平第x帧、垂直第j帧时测得的距离。
6.根据权利要求5所述的一种自适应分辨率三维激光扫描方法,其特征在于:步骤S5)中,若精细扫描数据P'的水平第x帧的角度α'x在粗扫描数据P中无对应角度,即α'x≠αi,i=1,2,...,n,则向粗扫描数据P中添加数据p'x,j=(α'x,βj,d'x,j),j=1,2,...,m,经过上述变化的粗扫描数据P就变为融合后的数据Pmix;
其中,αx为水平第x帧的角度;βj为垂直第j帧角度;d'x,j为水平第x帧、垂直第j帧时测得的距离。
7.根据权利要求6所述的一种自适应分辨率三维激光扫描方法,其特征在于:步骤S5)中,对融合后的数据Pmix进行插值滤波后处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911068898.5A CN110824456B (zh) | 2019-11-05 | 2019-11-05 | 一种自适应分辨率三维激光扫描方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911068898.5A CN110824456B (zh) | 2019-11-05 | 2019-11-05 | 一种自适应分辨率三维激光扫描方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110824456A CN110824456A (zh) | 2020-02-21 |
CN110824456B true CN110824456B (zh) | 2023-05-30 |
Family
ID=69552443
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911068898.5A Active CN110824456B (zh) | 2019-11-05 | 2019-11-05 | 一种自适应分辨率三维激光扫描方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110824456B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5898483A (en) * | 1997-05-01 | 1999-04-27 | Lockheed Martin Corporation | Method for increasing LADAR resolution |
CN108169730A (zh) * | 2016-12-07 | 2018-06-15 | 岭纬公司 | 基于区域的激光雷达变密度扫描系统及方法 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09159949A (ja) * | 1995-12-08 | 1997-06-20 | Hitachi Koki Co Ltd | 多ビーム光記録装置 |
US6600168B1 (en) * | 2000-02-03 | 2003-07-29 | Genex Technologies, Inc. | High speed laser three-dimensional imager |
AU2005200937A1 (en) * | 2005-03-02 | 2006-09-21 | Maptek Pty Ltd | Imaging system |
CN101916357A (zh) * | 2010-09-03 | 2010-12-15 | 西安富立叶微电子有限责任公司 | 一种激光条码识读装置及条码识读方法 |
JP2014095758A (ja) * | 2012-11-07 | 2014-05-22 | Canon Inc | アクチュエータ及び可変形状ミラー |
CN103292699B (zh) * | 2013-05-27 | 2016-04-13 | 深圳先进技术研究院 | 一种三维扫描系统及方法 |
CN103472458B (zh) * | 2013-09-16 | 2015-04-15 | 中国科学院上海光学精密机械研究所 | 基于声光扫描的三维视频激光雷达系统 |
US10509110B2 (en) * | 2015-12-29 | 2019-12-17 | The Boeing Company | Variable resolution light radar system |
CN106289098B (zh) * | 2016-07-22 | 2018-09-14 | 武汉海达数云技术有限公司 | 三维激光扫描仪的旋转控制方法 |
KR101997095B1 (ko) * | 2016-07-22 | 2019-07-08 | 전자부품연구원 | 수평 분해능 및 영상획득 프레임이 제어되는 스캐닝 라이다 |
US10546373B2 (en) * | 2016-08-03 | 2020-01-28 | Sightline Innovation Inc. | System and method for integrated laser scanning and signal processing |
WO2018126248A1 (en) * | 2017-01-02 | 2018-07-05 | Okeeffe James | Micromirror array for feedback-based image resolution enhancement |
CN110140063B (zh) * | 2016-11-30 | 2023-10-13 | 布莱克莫尔传感器和分析有限责任公司 | 利用光学测距系统进行自适应扫描的方法和系统 |
WO2018102188A1 (en) * | 2016-11-30 | 2018-06-07 | Blackmore Sensors and Analytics Inc. | Method and system for automatic real-time adaptive scanning with optical ranging systems |
RU2679923C1 (ru) * | 2017-10-05 | 2019-02-14 | Акционерное общество "Лазерные системы" (АО "Лазерные системы") | Способ получения пространственной модели окружающей обстановки в режиме реального времени на основе данных лазерной локации и устройство для его осуществления |
CN109444855B (zh) * | 2018-12-29 | 2021-02-26 | 森思泰克河北科技有限公司 | 激光雷达及其扫描方法 |
CN110018491B (zh) * | 2019-04-16 | 2021-06-22 | 广州大学 | 激光扫描方法、装置及激光雷达 |
CN113075680B (zh) * | 2020-01-06 | 2024-07-05 | 宁波舜宇车载光学技术有限公司 | 激光雷达和激光雷达的制造方法 |
-
2019
- 2019-11-05 CN CN201911068898.5A patent/CN110824456B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5898483A (en) * | 1997-05-01 | 1999-04-27 | Lockheed Martin Corporation | Method for increasing LADAR resolution |
CN108169730A (zh) * | 2016-12-07 | 2018-06-15 | 岭纬公司 | 基于区域的激光雷达变密度扫描系统及方法 |
Non-Patent Citations (2)
Title |
---|
一种自适应分辨率的三维激光扫描测量方法;张达;余乐文;余斌;杨小聪;张元生;;矿冶(第04期) * |
地面三维激光扫描点云分辨率研究;朱凌;石若明;;遥感学报(第03期) * |
Also Published As
Publication number | Publication date |
---|---|
CN110824456A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102016107959B4 (de) | Auf strukturiertem Licht basierende Multipfadlöschung bei ToF-Bilderzeugung | |
CN109725303B (zh) | 坐标系的修正方法及装置、存储介质 | |
CN109360239B (zh) | 障碍物检测方法、装置、计算机设备和存储介质 | |
EP4109131A1 (en) | Laser radar parameter calibration method and apparatus | |
CN109178234B (zh) | 一种船舶干舷高度测量系统及其测量方法 | |
CN113111513B (zh) | 传感器配置方案确定方法、装置、计算机设备及存储介质 | |
CN112146848A (zh) | 一种确定摄像头的畸变参数的方法及装置 | |
CN116299319B (zh) | 多激光雷达的同步扫描及点云数据处理方法和雷达系统 | |
CN112014829B (zh) | 一种激光雷达扫描仪的性能指标测试方法与装置 | |
CN113777592B (zh) | 方位角标定方法和装置 | |
CN117289300A (zh) | 点云矫正方法、激光雷达及机器人 | |
Chaudhry et al. | A modeling approach for predicting the resolution capability in terrestrial laser scanning | |
CN110824456B (zh) | 一种自适应分辨率三维激光扫描方法 | |
EP3862787A1 (en) | De-jitter of point cloud data for target recognition | |
CN111679262B (zh) | 激光点云强度标定方法、装置、设备及存储介质 | |
Pan et al. | Calibration of an airborne single-photon LiDAR system with a wedge scanner | |
CN116819561A (zh) | 一种点云数据匹配方法、系统、电子设备及存储介质 | |
CN113589326B (zh) | 激光雷达的物体识别方法、装置、设备及存储介质 | |
CN113513988A (zh) | 一种激光雷达标靶检测方法、装置、车辆、存储介质 | |
CN112364741B (zh) | 无人机单目远距离障碍检测方法、装置及无人机 | |
CN114167390A (zh) | 一种车载毫米波雷达动态标定方法及系统 | |
CN111596309A (zh) | 一种基于激光雷达的车辆排队测量方法 | |
CN113706634B (zh) | 一种视觉标定方法、装置、电子设备及存储介质 | |
DE102020205470A1 (de) | Verfahren zum Verarbeiten optischer Sensordaten | |
KR101840328B1 (ko) | 3d 레이저 스캐너 |
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 |