CN115166681B - 调频连续波信号体制穿墙雷达目标检测快速方法及系统 - Google Patents
调频连续波信号体制穿墙雷达目标检测快速方法及系统 Download PDFInfo
- Publication number
- CN115166681B CN115166681B CN202211091992.4A CN202211091992A CN115166681B CN 115166681 B CN115166681 B CN 115166681B CN 202211091992 A CN202211091992 A CN 202211091992A CN 115166681 B CN115166681 B CN 115166681B
- Authority
- CN
- China
- Prior art keywords
- sequence
- fourier transform
- preset condition
- distance
- frequency
- 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
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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种调频连续波信号体制穿墙雷达目标检测快速方法及系统,方法包括:采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息,获得对应的差拍信号;根据获得的差拍信号,获取每一差拍信号分别所对应的回波距离像;针对任意相邻两帧回波距离像,进行差分获得差分距离像的序列;筛选出满足第一预先设定条件的差分距离像;针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数;针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
Description
技术领域
本发明涉及穿墙雷达技术领域,尤其涉及一种调频连续波信号体制穿墙雷达目标检测快速方法及系统。
背景技术
穿墙雷达可穿透墙壁对墙后运动目标进行探测,是城市环境中抢险救援等应用中非常重要的探测感知手段。穿墙雷达从信号体制上分类,可以分为冲激脉冲体制和连续波体制等,其中连续波体制又可以分为线性调频连续波(FMCW)、步进频率(SFCW)、伪随机相位编码等若干种。由于连续波体制可以保证脉冲间的严格相参,有利于通过多帧积累提高目标探测性能,或通过多个接收天线提高测向能力,因此信号处理增益潜力巨大,是穿墙雷达重要的技术发展方向之一。
连续波穿墙雷达特别是FMCW穿墙雷达发射功率很低,特别是墙壁衰减很大造成目标回波微弱,对目标检测带来了很大挑战。通过前后两帧数据进行差分检测的方法,虽然可以消除场景中静止目标回波的影响,但穿透较厚墙壁后的微弱运动目标信噪比很低,特别是穿墙后室内场景中有植物、电扇等固定位置的动目标,仅仅通过运动信息难以区分,仍然会造成大量的虚警。同时,考虑到穿墙雷达往往采用便携或手持的方式工作,处理器处理速度有限,因此很多运算量需求较大的方法也无法适用。为了降低整机的体积、重量和功耗,检测算法的运算量也应该尽可能低,使其可以在单片机等运算资源受限的处理器上运行。
发明内容
(一)要解决的技术问题
鉴于现有技术的上述缺点、不足,本发明提供一种调频连续波信号体制穿墙雷达目标检测快速方法及系统,其解决了现有技术中连续波穿墙雷达特别是FMCW穿墙雷达发射功率很低,特别是墙壁衰减很大造成目标回波微弱,对目标检测不是特别精确的技术问题。
(二)技术方案
为了达到上述目的,本发明采用的主要技术方案包括:
第一方面,本发明实施例提供一种调频连续波信号体制穿墙雷达目标检测快速方法,包括:
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息,并获得与每帧发射信息分别一一对应的差拍信号;
S2、根据获得的与每帧发射信息分别一一对应的差拍信号,获取每一差拍信号分别所对应的回波距离像;
S3、针对任意相邻两帧回波距离像,进行差分获得差分距离像的序列;
所述差分距离像的序列包括:依次相邻两帧回波距离像之间的差分值;
S4、在差分距离像的序列中,筛选出满足第一预先设定条件的差分距离像;
S5、针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数;
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
优选地,S1具体包括:
将FMCW接收机按照设定频率所发射的任一帧发射信息的耦合放大一路
作为本振信号,并将该本振信号与FMCW接收机接收到的与所述发射信息对应的回波信号进行混频和低通滤波,获得与该发射信息所对应的零中频的差拍信号。
优选地,
其中,FMCW接收机的发射信息幅度归一化为1;
exp是指数;
j是虚部符号;
f 0 是线性调频起始频率;
T为扫描时间长度;
K为调频斜率,K=B/T;
B为调频信号带宽;
t为时刻;
m为发射信息的序号;
A为回波强度。
优选地,
优选地,所述S2具体包括:
f是t做傅里叶变换后的频率变量表示;
优选地,所述S3中,
所述差分距离像的序列包括:依次相邻两帧回波距离像之间的差分值;
N是预先设定值;
优选地,所述S4中的所述第一预先设定条件为:差分距离像大于Tp;
p为预先设定的正实数。
优选地,所述S5具体包括:
S51、基于满足第一预先设定条件的差分距离像,获取满足第一预先设定条件的差分距离像对应距离处的差拍信号;
M指的是满足第一预先设定条件的差分距离像的数量;
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式:
其中,傅里叶变换后的复散射值的函数为:
优选地,
所述第二预设条件为:傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点。
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
另一方面,本实施例还提供一种调频连续波信号体制穿墙雷达目标检测快速系统,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如上述任一所述的调频连续波信号体制穿墙雷达目标检测快速方法。
(三)有益效果
本发明的有益效果是:本发明的一种调频连续波信号体制穿墙雷达目标检测快速方法及系统,由于在差分距离像的序列中,筛选出满足第一预先设定条件的差分距离像;针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数;针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。在此过程中进行了两次傅里叶变换,对目标的检测更为准确。
附图说明
图1为本发明的一种调频连续波信号体制穿墙雷达目标检测快速方法流程图;
图2为本发明实施例中的H条距离维响应非相参叠加结果示意图;
图3为本发明实施例中的图2中峰值点位置为13米处结果的包络示意图。
具体实施方式
为了更好的解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
为了更好的理解上述技术方案,下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更清楚、透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
参见图1,本实施例提供一种调频连续波信号体制穿墙雷达目标检测快速方法,包括:
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息,并获得与每帧发射信息分别一一对应的差拍信号;
在本实施例的实际应用中,S1具体包括:
将FMCW接收机按照设定频率所发射的任一帧发射信息的耦合放大一路
作为本振信号,并将该本振信号与FMCW接收机接收到的与所述发射信息对应的回波信号进行混频和低通滤波,获得与该发射信息所对应的零中频的差拍信号。
其中,FMCW接收机的发射信息幅度归一化为1。
exp是指数。
j是虚部符号。
f 0 是线性调频起始频率。
T为扫描时间长度。
K为调频斜率,K=B/T。
B为调频信号带宽。
t为时刻。
m为发射信息的序号。
A为回波强度。
在本实施例的具体应用中,
在具体应用中,差拍信号是一个线性调频信号,起始频率为:
M为接收回波信号的数量,积累总时间为MT。
具体应用中,假设人行走的最大速度为7.2km/h或2m/s,雷达工作于L波段,,此时行走速度引起的多普勒频移为,远远低于,即,同时假设带宽和扫描时间,调频斜率,因此DK项不可忽略,而项可忽略,且。此时各相位
项可化简为:
S2、根据获得的与每帧发射信息分别一一对应的差拍信号,获取每一差拍信号分别所对应的回波距离像。
在本实施例的实际应用中,所述S2具体包括:
f是t做傅里叶变换后的频率变量表示。
S3、针对任意相邻两帧回波距离像,进行差分获得差分距离像的序列;
所述差分距离像的序列包括:依次相邻两帧回波距离像之间的差分值;
在本实施例的实际应用中,所述S3中,所述差分距离像的序列包括:依次相邻两帧回波距离像之间的差分值。
N是预先设定值。
S4、在差分距离像的序列中,筛选出满足第一预先设定条件的差分距离像。
所述第一预先设定条件为:差分距离像大于Tp。
p为预先设定的正实数。
参见图2,具体应用中,由于距离像目标的信杂比相对较低,因此直接通过单帧进行判断可能会造成误判,因此可以将H条(H是某个正整数,对H条回波进行非相参叠加。H可以等于M)回波做傅里叶变换后的回波距离像非相参叠加(如图2),并选择其峰值点(如图2中的圆圈处)。实际为了减少虚假峰值点,可选择某一门限(如均值的若干倍),只选择高于该门限的峰值点。图2中峰值点位置分别为13,15.5,18.5, 21.5,26.5,28米处。
S5、针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数。
所述S5具体包括:
S51、基于满足第一预先设定条件的差分距离像,获取满足第一预先设定条件的差分距离像对应距离处的差拍信号。
M指的是满足第一预先设定条件的差分距离像的数量;
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式:
其中,傅里叶变换后的复散射值的函数为:
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
所述第二预设条件为:傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点。
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
参见图3,针对上例中距离为13,15.5,18.5,21.5,26.5,28米处峰值点,沿序列方向做傅里叶变换后截取人体行走速度范围内的包络,其中距离为13米处结果的包络如图3。其速度峰值为0.5m/s,该距离上速度峰值点为所有距离处的速度峰值点的最大值,因此可判断该距离处为有效目标。
本实施例中的一种调频连续波信号体制穿墙雷达目标检测快速方法及系统,由于在差分距离像的序列中,筛选出满足第一预先设定条件的差分距离像;针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数;针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。在此过程中进行了两次傅里叶变换,对目标的检测更为准确。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例,或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。
应当注意的是,在权利要求中,不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中,这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用,仅是为了表述方便,而不表示任何顺序。可将这些词语理解为部件名称的一部分。
此外,需要说明的是,在本说明书的描述中,术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述,是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管已描述了本发明的优选实施例,但本领域的技术人员在得知了基本创造性概念后,则可对这些实施例作出另外的变更和修改。所以,权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也应该包含这些修改和变型在内。
Claims (10)
1.一种调频连续波信号体制穿墙雷达目标检测快速方法,其特征在于,包括:
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息,并获得与每帧发射信息分别一一对应的差拍信号;
S2、根据获得的与每帧发射信息分别一一对应的差拍信号,获取每一差拍信号分别所对应的回波距离像;
S3、针对任意相邻两帧回波距离像,进行差分获得差分距离像的序列;
所述差分距离像的序列包括:依次相邻两帧回波距离像之间的差分值;
S4、在差分距离像的序列中,筛选出满足第一预先设定条件的差分距离像;
S5、针对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换,得到傅里叶变换后的复散射值的函数;
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点,并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
8.根据权利要求7所述的方法,其特征在于,所述S5具体包括:
S51、基于满足第一预先设定条件的差分距离像,获取满足第一预先设定条件的差分距离像对应距离处的差拍信号;
M指的是满足第一预先设定条件的差分距离像的数量;
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式:
其中,傅里叶变换后的复散射值的函数为:
9.根据权利要求8所述的方法,其特征在于,
所述第二预设条件为:傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点;
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
10.一种调频连续波信号体制穿墙雷达目标检测快速系统,其特征在于,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1-9中任一所述的调频连续波信号体制穿墙雷达目标检测快速方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211091992.4A CN115166681B (zh) | 2022-09-07 | 2022-09-07 | 调频连续波信号体制穿墙雷达目标检测快速方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211091992.4A CN115166681B (zh) | 2022-09-07 | 2022-09-07 | 调频连续波信号体制穿墙雷达目标检测快速方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115166681A CN115166681A (zh) | 2022-10-11 |
CN115166681B true CN115166681B (zh) | 2022-12-06 |
Family
ID=83481897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211091992.4A Active CN115166681B (zh) | 2022-09-07 | 2022-09-07 | 调频连续波信号体制穿墙雷达目标检测快速方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115166681B (zh) |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10451723B2 (en) * | 2016-12-20 | 2019-10-22 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (CW) radar sensing system |
CN107153189B (zh) * | 2017-04-18 | 2021-08-03 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
CN107765240B (zh) * | 2017-10-20 | 2020-11-06 | 中国人民解放军海军七〇一工厂 | 一种运动状态的判断方法、装置和电子设备 |
CN108387894B (zh) * | 2018-04-13 | 2021-07-27 | 中南大学 | 穿墙雷达回波数据的处理方法 |
CN108776336B (zh) * | 2018-06-11 | 2022-06-03 | 电子科技大学 | 一种基于emd的自适应穿墙雷达静止人体目标定位方法 |
CN112654888A (zh) * | 2018-09-05 | 2021-04-13 | 京瓷株式会社 | 电子设备、电子设备的控制方法、以及电子设备的控制程序 |
CN110873877B (zh) * | 2019-04-25 | 2021-04-23 | 北京航空航天大学 | 目标运动轨迹的确定方法及装置 |
CN110118966B (zh) * | 2019-05-28 | 2020-10-13 | 长沙莫之比智能科技有限公司 | 基于毫米波雷达的人员检测与计数系统 |
CN110244296B (zh) * | 2019-06-17 | 2020-07-14 | 中国人民解放军国防科技大学 | 基于距离像复包络相位匹配处理的空间进动目标成像方法 |
CN110161491B (zh) * | 2019-06-28 | 2021-01-12 | 电子科技大学 | 一种针对微弱生命体的测距和呼吸频率估计方法 |
CN111352102B (zh) * | 2020-03-18 | 2023-02-24 | 重庆邮电大学 | 一种基于调频连续波雷达的多目标个数检测方法及装置 |
-
2022
- 2022-09-07 CN CN202211091992.4A patent/CN115166681B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115166681A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107966688B (zh) | 基于相位干涉技术的宽带雷达目标速度解模糊方法 | |
CN110646773B (zh) | 一种基于汽车毫米波雷达的弱目标检测方法、跟踪方法及系统 | |
Li et al. | Fast coherent integration for maneuvering target with high-order range migration via TRT-SKT-LVD | |
US8305257B2 (en) | Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents | |
US10054668B2 (en) | Probabilistic signal, detection, and track processing architecture and system | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN107656255A (zh) | 基于多径回波的超宽带雷达动目标二维定位方法 | |
US20220179062A1 (en) | Detection apparatus and method | |
CN108535730A (zh) | 一种多普勒气象雷达解速度模糊方法和系统 | |
CN110376563B (zh) | 基于gps辐射源的雷达运动目标检测方法 | |
CN109298417B (zh) | 一种基于雷达信号处理的建筑物内部结构探测方法及装置 | |
Heaney et al. | Measurements of three-dimensional propagation in a continental shelf environment | |
CN105589061A (zh) | 一种岸基雷达的信号处理算法 | |
US7942821B2 (en) | Doppler velocity detection device and ultrasonographic device using the same | |
CN115166681B (zh) | 调频连续波信号体制穿墙雷达目标检测快速方法及系统 | |
CN113419238A (zh) | 基于毫米波雷达的山体滑坡监测方法、电子设备、存储介质 | |
JP5423713B2 (ja) | アクティブソーナー装置及びその信号処理方法 | |
Kronauge | Waveform design for continuous wave radars | |
US6721678B1 (en) | Determining three dimensional velocity of a object | |
Ranney et al. | Recent MTI experiments using ARL's synchronous impulse reconstruction (SIRE) radar | |
CN110907930B (zh) | 一种基于角度估计的车载雷达目标检测估计方法及装置 | |
CN116027288A (zh) | 生成数据的方法、装置、电子设备及存储介质 | |
CN108387891B (zh) | 基于短基线一发两收的超宽带雷达干涉二维定位跟踪方法 | |
CN108983192B (zh) | 基于gps辐射源的雷达运动目标参数估计方法 | |
CN112698291B (zh) | 一种基于cnn的机载气象雷达气象目标检测方法 |
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 |