CN114814778B - 一种基于毫米波雷达的载体速度解算方法 - Google Patents
一种基于毫米波雷达的载体速度解算方法 Download PDFInfo
- Publication number
- CN114814778B CN114814778B CN202210744334.4A CN202210744334A CN114814778B CN 114814778 B CN114814778 B CN 114814778B CN 202210744334 A CN202210744334 A CN 202210744334A CN 114814778 B CN114814778 B CN 114814778B
- Authority
- CN
- China
- Prior art keywords
- doppler
- carrier
- speed
- millimeter wave
- wave radar
- 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
- 238000004364 calculation method Methods 0.000 title claims description 10
- 239000011159 matrix material Substances 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000003068 static effect Effects 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000012163 sequencing technique Methods 0.000 claims abstract description 3
- 101100083446 Danio rerio plekhh1 gene Proteins 0.000 claims 5
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 claims 3
- 238000003491 array Methods 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000004075 alteration Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000000779 smoke Substances 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了基于毫米波雷达的载体速度解算方法,包括:毫米波雷达对回波信号进行二维FFT处理,生成距离多普勒热图矩阵;毫米波雷达视野内,计算所有目标在距离多普勒热图矩阵中的回波能量值;通过计算地面静止目标相对雷达的速度;在距离多普勒热图矩阵上按多普勒单元数生成数组,并对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组;对数组进行排序计算,最大的多普勒索引即为当前帧载体速度单元索引;根据当前帧载体速度单元索引与已知的速度分辨率,计算出当前帧的载体速度并滤波输出。本发明能有效区分地面静止目标与运动目标;载体速度估算精度高,能有效滤除雷达视野内运动目标的干扰。
Description
技术领域
本发明属于雷达技术领域,尤其涉及一种基于毫米波雷达的载体速度解算方法。
背景技术
毫米波雷达体积小、重量轻、测量精度高,且穿透烟雾、灰尘能力强,具有全天时全天候工作能力等优点,已广泛使用于汽车辅助驾驶、交通流量检测、无人机定高避障、智慧交通等领域,主要应用于雷达视场内典型目标的检测与跟踪。在毫米波雷达安装在机器人、汽车、无人机等运动载体上时,往往需要通过接入到运动载体系统来获取平台速度信息,以此来辅助雷达区分相对地面静止与运动目标;而应用在摩托车、电动车上时,往往难以获取到载体的速度信息。
现有毫米波雷达解算载体速度的方法,都是基于雷达检测的点云结果,来进行载体速度的解算。
专利公开号CN112946623A的中国专利,通过判断点云数据集对应X轴正负两侧的点云数据的数量同时大于一个点且速度相等,则确定该点云数据为有效点云数据,利用该有效点云数据进行载体速度的计算。
专利公开号CN113075661A的中国专利,通过探测区域内的所在物体与车辆的相对速度 ,通过众数运算将该探测区域内的所在物体中相同速度最多的判断为静止物。
专利公开号CN113911173A的中国专利,通过依靠毫米波雷达和加速度计两种传感器实现列车速度的可靠测量。
专利公开号CN113911174A的中国专利,通过相邻两帧点云图像时间差内的运行距离,来计算列车的运行速度;融合点云匹配测速与毫米波雷达测速实现对列车速度的确定。
以上专利中,利用毫米波雷达进行载体速度的估算可以分为两大类,一是利用毫米波雷达与其他传感器融合进行载体速度估算;一是利用毫米波雷达单传感器,通过毫米波雷达检测到的目标点云信息,进行载体速度的估算。这些方法,对静止目标的识别率较低,容易受运动目标干扰。在特定场景下,在雷达视野内存在其他运动目标,且地面较为空旷,雷达探测到地面静止目标点云较少,甚至没有地面静止目标点云时,会存在将雷达视野内其他运动目标相对于雷达的速度,计算成载体运动的速度,导致出现明显的速度偏差。
发明内容
有鉴于此,本发明提供了一种毫米波雷达的载体速度解算方法,通过距离多普勒热图进行能量积累,完成对地面静止目标的识别;通过识别的地面静止目标,与其对应的载体速度单元索引,直接计算当前帧的载体速度;采用滑动平均滤波法对载体速度进行滤波;可实现准确识别地面静止目标,降低运动目标的干扰,较为准确的估算毫米波雷达的载体速度,解决特定场景下出现明显的速度偏差的问题。
本发明公开的一种基于毫米波雷达的载体速度解算方法,包括以下步骤:
毫米波雷达对回波信号进行二维FFT处理,生成距离多普勒热图矩阵;
毫米波雷达视野内,计算所有目标在距离多普勒热图矩阵中的回波能量值;
通过计算地面静止目标相对雷达的速度,此速度即为毫米波雷达载体相对于地面的速度;
地面静止目标的判别,在距离多普勒热图矩阵上按多普勒单元数生成数组,并对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组;对数组进行排序计算,最大的多普勒索引即为当前帧载体速度单元索引;根据当前帧载体速度单元索引,与已知的速度分辨率,计算出当前帧的载体速度;
采用滑动平均滤波法对载体速度进行滤波输出。
进一步的,S1步骤中所述距离多普勒热图矩阵如下:
进一步的,S2步骤中对于不同RCS的目标,在距离多普勒热图矩阵中计算其回波能量值,对于高RCS目标,提取点云目标。
进一步的,S4步骤具体包括如下:
对距离多普勒热图矩阵,按多普勒单元数生成数组doppler[n],数组长度为多普勒单元数n;
在距离多普勒热图矩阵中,对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组doppler[n]:
j为第j个多普勒单元;
根据当前帧载体速度单元索引,计算当前帧的载体速度:
进一步的,S5步骤采用滑动平均滤波法对载体速度进行滤波,具体包括如下:
连续取N帧估算的载体速度,组成一个循环队列,采用先进先出的原则,当前帧估算的载体速度放入队尾,并舍弃掉原来队首的数据,对当前队列中的N个数据取算术平均值
本发明的有益效果如下:
1)毫米波雷达能够有效区分地面静止目标与运动目标;
2)毫米波雷达载体速度估算精度高,场景适用性强,鲁棒性好。
3)在雷达视野内存在其他运动目标,且地面较为空旷,雷达探测到地面点云目标较少,甚至没有地面目标点云时,也能够有效估算毫米波雷达载体速度,能有效滤除雷达视野内运动目标的干扰。
4)无需根据目标相对与雷达的位置进行极坐标转换,计算简单,运算速度快,准确性高,实时性好,能够很好得应用于产品中。
附图说明
图1本发明毫米波雷达的载体速度解算方法流程图;
图2只有地面静止目标时的距离多普勒热图;
图3运动目标与地面目标距离的多普勒热图;
图4载体速度计算结果示例图。
具体实施方式
下面结合附图对本发明作进一步的说明,但不以任何方式对本发明加以限制,基于本发明教导所作的任何变换或替换,均属于本发明的保护范围。
本发明提供的基于毫米波雷达的载体速度解算方法技术方案如下:
S1毫米波雷达对回波信号进行二维FFT处理,生成距离多普勒热图矩阵;
S2毫米波雷达视野内,所有目标都在距离多普勒热图矩阵中有回波能量值;
S3通过计算地面静止目标相对雷达的速度,此速度即为毫米波雷达载体相对于地面的速度;
S4地面静止目标的判别,在距离多普勒热图矩阵上进行;
首先,对距离多普勒热图矩阵,按多普勒单元数生成数组;
其次,在距离多普勒热图矩阵中,对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组;
然后,对数组进行排序计算,得到数组值从大到小排序中,最大的多普勒索引值;最大的多普勒索引即为当前帧载体速度单元索引;
根据当前帧载体速度单元索引,与已知的速度分辨率,即可计算出当前帧的载体速度;
S5采用滑动平均滤波法对载体速度进行滤波输出。
在一个实施例中,毫米波雷达的载体速度解算方法,包括以下步骤:
S1毫米波雷达对回波信号进行二维FFT处理,生成距离多普勒热图矩阵:
S2毫米波雷达视野内,不同RCS的目标,都在距离多普勒热图矩阵中有不同的回波能量值。低RCS目标回波能量较低,不能通过cfar算法形成点云目标,高RCS目标回波能量较高,可以形成点云目标。
S3计算地面静止目标相对雷达的速度,此速度即为毫米波雷达载体相对于地面的速度。
S4在距离多普勒热图矩阵上进行地面静止目标的判别。
首先,对距离多普勒热图矩阵,按多普勒单元数生成数组doppler[n],数组长度为多普勒单元数n。
其次,在距离多普勒热图矩阵中,对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组doppler[n]:
j为第j个多普勒单元。
根据当前帧载体速度单元索引,计算当前帧的载体速度:
S5采用滑动平均滤波法对载体速度进行滤波。
具体实现方式为连续取N帧估算的载体速度,组成一个循环队列,采用先进先出的原则,当前帧估算的载体速度放入队尾,并舍弃掉原来队首的数据,对当前队列中的N个数据取算术平均值。
参考图2和图3,图2为只有地面静止目标时的距离多普勒热图,静止目标的分布,随着距离的增加(即图中距离索引越来越大),静止目标分布在同一个多普勒索引附近;图3为运动目标与地面目标的距离多普勒热图,静止目标的分布,随着距离的增加(即图中距离索引越来越大),静止目标分布在同一个多普勒索引附近,而运动目标分布在其他远离静止目标的多普勒索引附近;图4载体速度的计算结果,虚线为当前帧计算的载体速度,实线为经过滑动平均滤波法对载体速度进行滤波后的载体速度。
本发明的有益效果如下:
1)毫米波雷达能够有效区分地面静止目标与运动目标;
2)毫米波雷达载体速度估算精度高,场景适用性强,鲁棒性好。
3)在雷达视野内存在其他运动目标,且地面较为空旷,雷达探测到地面点云目标较少,甚至没有地面目标点云时,也能够有效估算毫米波雷达载体速度,能有效滤除雷达视野内运动目标的干扰。
4)无需根据目标相对与雷达的位置进行极坐标转换,计算简单,运算速度快,准确性高,实时性好,能够很好得应用于产品中。
本文所使用的词语“优选的”意指用作实例、示例或例证。本文描述为“优选的”任意方面或设计不必被解释为比其他方面或设计更有利。相反,词语“优选的”的使用旨在以具体方式提出概念。如本申请中所使用的术语“或”旨在意指包含的“或”而非排除的“或”。即,除非另外指定或从上下文中清楚,“X使用A或B”意指自然包括排列的任意一个。即,如果X使用A;X使用B;或X使用A和B二者,则“X使用A或B”在前述任一示例中得到满足。
而且,尽管已经相对于一个或实现方式示出并描述了本公开,但是本领域技术人员基于对本说明书和附图的阅读和理解将会想到等价变型和修改。本公开包括所有这样的修改和变型,并且仅由所附权利要求的范围限制。特别地关于由上述组件(例如元件等)执行的各种功能,用于描述这样的组件的术语旨在对应于执行所述组件的指定功能(例如其在功能上是等价的)的任意组件(除非另外指示),即使在结构上与执行本文所示的本公开的示范性实现方式中的功能的公开结构不等同。此外,尽管本公开的特定特征已经相对于若干实现方式中的仅一个被公开,但是这种特征可以与如可以对给定或特定应用而言是期望和有利的其他实现方式的一个或其他特征组合。而且,就术语“包括”、“具有”、“含有”或其变形被用在具体实施方式或权利要求中而言,这样的术语旨在以与术语“包含”相似的方式包括。
本发明实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以多个或多个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。上述提到的存储介质可以是只读存储器,磁盘或光盘等。上述的各装置或系统,可以执行相应方法实施例中的存储方法。
综上所述,上述实施例为本发明的一种实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何背离本发明的精神实质与原理下所做的改变、修饰、代替、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (5)
1.一种基于毫米波雷达的载体速度解算方法,其特征在于,包括以下步骤:
S1:毫米波雷达对回波信号进行二维FFT处理,生成距离多普勒热图矩阵;
S2:毫米波雷达视野内,计算所有目标在距离多普勒热图矩阵中的回波能量值;
S3:通过计算地面静止目标相对雷达的速度,此速度即为毫米波雷达载体相对于地面的速度;
S4:地面静止目标的判别,在距离多普勒热图矩阵上按多普勒单元数生成数组,并对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组;对数组进行排序计算,计算得到当前帧载体速度单元索引;根据当前帧载体速度单元索引,与已知的速度分辨率,计算出当前帧的载体速度;
S5:对载体速度进行滤波输出;
其中S4步骤中计算得到当前帧载体速度单元索引,包括如下:
对距离多普勒热图矩阵,按多普勒单元数生成数组doppler[n],数组长度为多普勒单元数n;
在距离多普勒热图矩阵中,对所有多普勒单元对应的所有距离单元回波能量求和,并赋值给数组doppler[n]:
j为第j个多普勒单元;
对数组doppler[n]的值进行排序计算,得到数组值从大到小排序中,最大的三个值的多普勒索引值Idxmax1,Idxmax2,Idxmax3,且
doppler[Idxmax1]≥doppler[Idxmax2]≥doppler[Idxmax3];
当三个多普勒索引值相邻时,判断Idxmax1为地面静止目标在距离多普勒热图上的多普勒索引,取Idxmax=Idxmax1,Idxmax为当前帧载体速度单元索引;当三个多普勒索引值不相邻时,Idxmax取Idxmax1,Idxmax2,Idxmax3中最接近上一帧载体速度单元索引的值。
2.根据权利要求1所述的基于毫米波雷达的载体速度解算方法,其特征在于,S1步骤中所述距离多普勒热图矩阵如下:
RD=(aij)m×n
其中aij对应距离多普勒热图矩阵第i个距离单元第j个多普勒单元;m为距离单元数,n为多普勒单元数。
3.根据权利要求1所述的基于毫米波雷达的载体速度解算方法,其特征在于,S2步骤中对于不同RCS的目标,在距离多普勒热图矩阵中计算其回波能量值,对于高RCS目标,提取点云目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210744334.4A CN114814778B (zh) | 2022-06-29 | 2022-06-29 | 一种基于毫米波雷达的载体速度解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210744334.4A CN114814778B (zh) | 2022-06-29 | 2022-06-29 | 一种基于毫米波雷达的载体速度解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114814778A CN114814778A (zh) | 2022-07-29 |
CN114814778B true CN114814778B (zh) | 2022-09-20 |
Family
ID=82522400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210744334.4A Active CN114814778B (zh) | 2022-06-29 | 2022-06-29 | 一种基于毫米波雷达的载体速度解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114814778B (zh) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190011547A1 (en) * | 2017-07-07 | 2019-01-10 | Veoneer Us, Inc. | Systems and methods for resolving velocity ambiguity in an automotive radar system |
CN111045008B (zh) * | 2020-01-15 | 2023-06-09 | 深圳市华讯方舟微电子科技有限公司 | 基于展宽计算的车载毫米波雷达目标识别方法 |
CN111308437B (zh) * | 2020-02-27 | 2022-03-22 | 南京慧尔视智能科技有限公司 | 一种毫米波mimo交通雷达求熵解速度模糊的方法 |
CN111398943B (zh) * | 2020-04-02 | 2022-07-29 | 森思泰克河北科技有限公司 | 目标姿态的确定方法及终端设备 |
US20210364616A1 (en) * | 2020-05-21 | 2021-11-25 | Hertzwell Pte. Ltd. | Radar system and computer-implemented method for radar target detection |
CN112098990B (zh) * | 2020-11-19 | 2021-02-02 | 长沙莫之比智能科技有限公司 | 车载高分辨毫米波雷达对于中高速车辆的检测与跟踪方法 |
CN112505648B (zh) * | 2020-11-19 | 2023-06-30 | 西安电子科技大学 | 基于毫米波雷达回波的目标特征提取方法 |
CN112690270A (zh) * | 2020-12-21 | 2021-04-23 | 西安京维智网科技有限公司 | 基于毫米波雷达的高压输电线路杆塔探鸟驱鸟系统及方法 |
CN114002654A (zh) * | 2021-06-30 | 2022-02-01 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种毫米波雷达解速度模糊方法 |
CN114384489A (zh) * | 2022-01-17 | 2022-04-22 | 长沙莫之比智能科技有限公司 | 一种浴室跌倒检测雷达的水雾杂波滤除方法 |
-
2022
- 2022-06-29 CN CN202210744334.4A patent/CN114814778B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114814778A (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144839B (zh) | 通过传感器融合检测长对象 | |
CN110058239B (zh) | 一种车载毫米波雷达装置及目标探测方法 | |
CN111144432B (zh) | 在传感器融合系统中消除模糊检测的方法 | |
Ogawa et al. | Pedestrian detection and tracking using in-vehicle lidar for automotive application | |
CN112526521B (zh) | 一种汽车毫米波防撞雷达的多目标跟踪方法 | |
CN110379178B (zh) | 基于毫米波雷达成像的无人驾驶汽车智能泊车方法 | |
CN107103275B (zh) | 使用雷达和视觉基于车轮进行的车辆检测和追踪 | |
CN111798700B (zh) | 盲区监测报警方法和装置 | |
CN111123212A (zh) | 一种基于复杂杂波背景下的场面监视雷达的信号处理方法 | |
CN112810619A (zh) | 基于雷达的辅助驾驶系统前方目标车辆识别方法 | |
CN111157995A (zh) | 一种汽车雷达侧视角度超分辨成像方法 | |
CN114137526A (zh) | 基于标签的车载毫米波雷达多目标检测方法和系统 | |
CN115128571B (zh) | 一种基于毫米波雷达的多人与非机动车识别方法 | |
CN114814778B (zh) | 一种基于毫米波雷达的载体速度解算方法 | |
CN111123269B (zh) | 用于无人机避障雷达的地面杂波抑制方法、模块及装置 | |
CN115033026B (zh) | 搭载云台的斜侧装无人机毫米波雷达避障和定高方法 | |
CN107256382A (zh) | 基于图像识别的虚拟保险杠控制方法和系统 | |
CN115965847A (zh) | 交叉视角下多模态特征融合的三维目标检测方法和系统 | |
CN114443636B (zh) | 一种基于毫米波雷达的隧道轨迹数据镜像去除方法 | |
CN114137523A (zh) | 车载毫米波雷达原始点云的噪点排除方法和系统 | |
CN115629385A (zh) | 基于毫米波雷达和相机关联的车辆排队长度实时检测方法 | |
CN114839628A (zh) | 利用多个距离和分辨率的物体检测 | |
Kruse et al. | Target classification based on near-distance radar sensors | |
CN117250595B (zh) | 一种车载毫米波雷达金属井盖目标误报抑制方法 | |
CN116577767B (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 |