CN114675263A - 一种使用双极化探地雷达的地下管线材质识别方法 - Google Patents
一种使用双极化探地雷达的地下管线材质识别方法 Download PDFInfo
- Publication number
- CN114675263A CN114675263A CN202210372190.4A CN202210372190A CN114675263A CN 114675263 A CN114675263 A CN 114675263A CN 202210372190 A CN202210372190 A CN 202210372190A CN 114675263 A CN114675263 A CN 114675263A
- Authority
- CN
- China
- Prior art keywords
- dual
- polarization
- polarized
- underground pipeline
- metal plate
- 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
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000000463 material Substances 0.000 title claims abstract description 13
- 239000002184 metal Substances 0.000 claims abstract description 43
- 230000010287 polarization Effects 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 3
- 239000002689 soil Substances 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- 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/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- 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/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种使用双极化探地雷达的地下管线材质识别方法,包括以下步骤:步骤(1):选取适当中心频率的双极化探地雷达,采集地下管线的双极化雷达数据;步骤(2):由于双极化探地雷达的HH通道和VV通道之间存在不平衡,利用金属板反射对双极化雷达数据进行极化校准;步骤(3):计算地下管线的极化阿尔法角;步骤(4):根据拟合的地下介质相对介电常数和计算得到的极化阿尔法角,判断地下管线的类型;极化阿尔法角值小则为塑料管线,值大为金属管线,值居中则为混凝土管线。计算地下管线反射的极化阿尔法角,利用双曲线拟合计算地下土体的介电常数;根据地下管线的极化阿尔法角和地下土体介电常数,研究不同类型管线的信号特征。
Description
技术领域
本发明属于雷达识别技术领域,具体涉及一种使用双极化探地雷达的地下管线材质识别方法。
背景技术
地下管网是指城市地下管线的铺设,是将供水、供电、通信、排水等等建成一个类似网络的拉链路。它是国家为了建成较为完善的城市地下管线体系;由于其在地下深埋,所以在检测其材质时候就极为困难,但是在城市改造的时候需要对地面进行挖掘,需要一种准确的方法探测地下管道并识别其材质,为城市改造提供地下管网图,确保地下改造的顺利进行。
发明内容
本发明的目的在于提供一种使用双极化探地雷达的地下管线材质识别方法,采用双极化雷达采集地下管线的VV和HH极化雷达数据,提取管线顶点处反射信号,对该信号进行包络和叠加处理并计算对应的极化阿尔法角,结合双曲线拟合估计的地下介质介电常数,获得管线材质信息。
为实现上述目的,本发明提供如下技术方案:一种使用双极化探地雷达的地下管线材质识别方法,包括以下步骤:
步骤(1):选取适当中心频率的双极化探地雷达,采集地下管线的双极化雷达数据;
步骤(2):由于双极化探地雷达的HH通道和VV通道之间存在不平衡,利用金属板对采集的双极化雷达数据进行极化校准;
步骤(3):计算地下管线的极化阿尔法角;
步骤(4):利用双曲线拟合估计地下管线周围介质的介电常数根据计算得到的极化阿尔法角,判断地下管线的类型;极化阿尔法角值小则为塑料管线,值大为金属管线,值居中则为混凝土管线。
优选的,所述步骤(1)具体包括以下步骤:
步骤(1.1)在测距轮触发模式下,使用双极化探地雷达垂直于地下管线方向采集双极化雷达数据,确保地下管线的反曲线反射形状完好;
步骤(1.2)对获得的双极化雷达数据进行数据处理,提高信噪比,数据处理流程包括:去直流,去背景,零时校正和带通滤波。
优选的,所述步骤(2)具体包括以下步骤:
步骤(2.1):将双极化雷达正对一块金属板保持不动,在时间或测距论触发模式下采集双极化雷达数据;
步骤(2.2):金属板的理论散射矩S阵式(1)所示:
式(1)中表明金属板的HH通道和VV通道反射振幅相同,利用式(1)对后续中计算得到的SHH和SVV进行校准;
步骤(2.3):计算金属板反射信号的能量SHH和SVV。由于最大正振幅或负振幅的不稳定,不能准确反应地下管线的反射强度,采用地下管线的反射信号包络积分来代替最大振幅,对步骤(2.2)中获得的金属板反射的双极化雷达数据,取绝对值后进行希尔伯特变换:
A(t)=a(t)+b(t)i (2);
式(2)中A(t)为希尔伯特变换后的双极化雷达信号,a(t)为瞬时振幅,b(t)为相位;
步骤(2.4):通过金属板的双极化反射信号及其包络图像,确定金属板反射起始和结束时间,金属板的反射信号强度为:
公式(3)中AHH和AVV为HH通道和VV通道的反射信号包络,t1为金属板反射起始时间,t2为金属板反射结束时间;
通过式(3)求得双极化雷达的极化散射校正矩阵为:
优选的,所述步骤(3)具体包括以下步骤:
步骤(3.1):根据步骤(2.3)中的方法计算步骤(1)中采集的地下管线的双极化反射能量并进行校准:
式(5)中t1和t2分别为地下管线反射的起始时间和结束时间,AHH和AVV为希尔伯特变换后的地下管线极化雷达反射信号,δHH和δVV为步骤(2)中计算得到的校正系数,和分别为校准之后的地下管线HH通道和VV通道反射能量;
计算极化阿尔法角α为:
与现有技术相比,本发明的有益效果如下:
(1)采用双极化探地雷达,同时获得VV通道和HH通道的极化雷达数据;
(2)提取双曲线顶点位置的双极化雷达的反射信号,对该信号进行包络叠加处理;
(3)计算地下管线反射的极化阿尔法角,利用双曲线拟合计算地下土体的介电常数;
(4)根据地下管线的极化阿尔法角和地下土体介电常数,研究不同类型管线的信号特征。
附图说明
图1为本发明的双极化雷达地下管线检测示意图;
图2为本发明的金属板的反射信号示意图;
图3为本发明的金属板的双极化反射信号及其包络图像,t1和t2为确定的金属板反射起始和结束时间。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1-3所示,本发明提供一种技术方案:一种使用双极化探地雷达的地下管线材质识别方法,包括以下步骤:
步骤(1):选取适当中心频率的双极化探地雷达,采集地下管线的双极化雷达数据;
所述步骤(1)具体包括以下步骤:
步骤(1.1)在测距轮触发模式下,使用双极化雷达垂直于地下管线方向采集双极化雷达数据,确保地下管线的反曲线反射形状完好;
步骤(1.2)对获得的双极化雷达数据进行数据处理,提高信噪比,数据处理流程包括:去直流,去背景,零时校正和带通滤波。
步骤(2):由于双极化探地雷达的HH通道和VV通道之间存在不平衡,对采集的双极化雷达数据进行极化校准;
所述步骤(2)具体包括以下步骤:
步骤(2.1):将双极化雷达正对一块金属板保持不动,在时间或测距论触发模式下采集双极化雷达数据;
步骤(2.2):金属板的理论散射矩S阵式(1)所示:
式(1)中表明金属板的HH通道和VV通道反射振幅相同,利用式(1)对后续中计算得到的SHH和SVV进行校准;
步骤(2.3):计算金属板反射信号的能量SHH和SVV,由于最大正振幅或负振幅的不稳定,不能准确反应地下管线的反射强度,采用地下管线的反射信号包络积分来代替最大振幅,对步骤(2.2)中获得的金属板反射的双极化雷达数据,取绝对值后进行希尔伯特变换:
A(t)=a(t)+b(t)i (2);
式(2)中A(t)为希尔伯特变换后的双极化雷达信号,a(t)为瞬时振幅,b(t)为相位;
步骤(2.4):通过金属板的双极化反射信号及其包络图像,确定金属板反射起始和结束时间,金属板的反射信号强度为:
公式(3)中AHH和AVV为HH通道和VV通道的反射信号包络,t1为金属板反射起始时间,t2为金属板反射结束时间;
通过式(3)求得双极化雷达的极化散射校正矩阵为:
步骤(3):计算地下管线的极化阿尔法角;
所述步骤(3)具体包括以下步骤:
步骤(3.1):根据步骤(2.3)中的方法计算步骤(1)中采集的地下管线的双极化反射能量并进行校准:
式(5)中t1和t2分别为地下管线反射的起始时间和结束时间,AHH和AVV为希尔伯特变换后的地下管线极化雷达反射信号,δHH和δVV为步骤(2)中计算得到的校正系数,和分别为校准之后的地下管线HH通道和VV通道反射能量;
计算极化阿尔法角α为:
步骤(4):根据拟合的地下介质相对介电常数和计算得到的极化阿尔法角,判断地下管线的类型;极化阿尔法角值小则为塑料管线,值大为金属管线,值居中则为混凝土管线。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (4)
1.一种使用双极化探地雷达的地下管线材质识别方法,其特征在于,包括以下步骤:
步骤(1):选取适当中心频率的双极化探地雷达,采集地下管线的双极化雷达数据;
步骤(2):由于双极化探地雷达的HH通道和VV通道之间存在不平衡,利用金属板反射对双极化雷达数据进行极化校准;
步骤(3):计算地下管线的极化阿尔法角;
步骤(4):根据拟合的地下介质相对介电常数和计算得到的极化阿尔法角,判断地下管线的类型;极化阿尔法角值小则为塑料管线,值大为金属管线,值居中则为混凝土管线。
2.根据权利要求1所述的一种使用双极化探地雷达的地下管线材质识别方法,其特征在于,所述步骤(1)具体包括以下步骤:
步骤(1.1)在测距轮触发模式下,使用双极化雷达垂直于地下管线方向采集双极化雷达数据,确保地下管线的反曲线反射形状完好;
步骤(1.2)对获得的双极化雷达数据进行数据处理,提高信噪比,数据处理流程包括:去直流,去背景,零时校正和带通滤波。
3.根据权利要求1所述的一种使用双极化探地雷达的地下管线材质识别方法,其特征在于,所述步骤(2)具体包括以下步骤:
步骤(2.1):将双极化雷达正对一块金属板保持不动,在时间或测距论触发模式下采集双极化雷达数据;
步骤(2.2):金属板的理论散射矩S阵式(1)所示:
式(1)中表明金属板的HH通道和VV通道反射振幅相同,利用式(1)对后续中计算得到的SHH和SVV进行校准;
步骤(2.3):将双极化雷达正对金属板,采集金属板的反射信号用于极化校准,提取金属板的极化反射信号,计算金属板反射信号的能量SHH和SVV,由于最大正振幅或负振幅的不稳定,不能准确反应地下管线的反射强度,采用地下管线的反射信号包络积分来代替最大振幅,对步骤(2.2)中获得的金属板反射的双极化雷达数据,取绝对值后进行希尔伯特变换:
A(t)=a(t)+b(t)i (2);
式(2)中A(t)为希尔伯特变换后的双极化雷达信号,a(t)为瞬时振幅,b(t)为相位;
步骤(2.4):通过金属板的双极化反射信号及其包络图像,确定金属板反射起始和结束时间,金属板的反射信号强度为:
公式(3)中AHH和AVV为HH通道和VV通道的反射信号包络,t1为金属板反射起始时间,t2为金属板反射结束时间;
通过式(3)求得双极化雷达的极化散射校正矩阵为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210372190.4A CN114675263B (zh) | 2022-04-11 | 2022-04-11 | 一种使用双极化探地雷达的地下管线材质识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210372190.4A CN114675263B (zh) | 2022-04-11 | 2022-04-11 | 一种使用双极化探地雷达的地下管线材质识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114675263A true CN114675263A (zh) | 2022-06-28 |
CN114675263B CN114675263B (zh) | 2024-05-28 |
Family
ID=82077646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210372190.4A Active CN114675263B (zh) | 2022-04-11 | 2022-04-11 | 一种使用双极化探地雷达的地下管线材质识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114675263B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116413719A (zh) * | 2023-06-12 | 2023-07-11 | 深圳大学 | 一种月球次表层下方熔岩管道勘探方法及相关设备 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3148367A (en) * | 1958-03-24 | 1964-09-08 | Tokyo Kelki Seizojo Kk | Target identification system by radar |
JPH0552949A (ja) * | 1991-08-28 | 1993-03-02 | Toda Constr Co Ltd | 地中探知レーダ装置 |
JP2004286461A (ja) * | 2003-03-19 | 2004-10-14 | Fujitsu Ltd | 地中探知方法及び装置 |
JP2005127874A (ja) * | 2003-10-24 | 2005-05-19 | National Institute Of Information & Communication Technology | 偏波合成開口レーダ較正方法及び装置 |
CN105467222A (zh) * | 2015-12-01 | 2016-04-06 | 厦门大学 | 基于单基地测量的地表介质参数反演方法 |
CN110244295A (zh) * | 2019-07-17 | 2019-09-17 | 广州大学 | 基于极化探地雷达的地下线性目标走向测量方法及*** |
CN110857979A (zh) * | 2018-08-26 | 2020-03-03 | 吴成彬 | 一种基于相参信号的极化频选特性材料抗干扰辨识***及方法 |
CN111505628A (zh) * | 2020-04-29 | 2020-08-07 | 中国南方电网有限责任公司超高压输电公司广州局 | 基于探地雷达对地下电缆成像的探测与识别方法 |
CN111679275A (zh) * | 2020-08-06 | 2020-09-18 | 中南大学 | 基于探地雷达的地下管线识别方法 |
CN111751392A (zh) * | 2020-07-30 | 2020-10-09 | 广州大学 | 一种基于双极化探地雷达的钢筋锈蚀检测方法 |
CN112862888A (zh) * | 2021-01-29 | 2021-05-28 | 广东电网有限责任公司 | 一种电缆的定位方法、装置、计算机设备和存储介质 |
CN112882019A (zh) * | 2021-01-14 | 2021-06-01 | 长春工程学院 | 基于旋转式单极化探地雷达的全极化目标辨识分类方法 |
-
2022
- 2022-04-11 CN CN202210372190.4A patent/CN114675263B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3148367A (en) * | 1958-03-24 | 1964-09-08 | Tokyo Kelki Seizojo Kk | Target identification system by radar |
JPH0552949A (ja) * | 1991-08-28 | 1993-03-02 | Toda Constr Co Ltd | 地中探知レーダ装置 |
JP2004286461A (ja) * | 2003-03-19 | 2004-10-14 | Fujitsu Ltd | 地中探知方法及び装置 |
JP2005127874A (ja) * | 2003-10-24 | 2005-05-19 | National Institute Of Information & Communication Technology | 偏波合成開口レーダ較正方法及び装置 |
CN105467222A (zh) * | 2015-12-01 | 2016-04-06 | 厦门大学 | 基于单基地测量的地表介质参数反演方法 |
CN110857979A (zh) * | 2018-08-26 | 2020-03-03 | 吴成彬 | 一种基于相参信号的极化频选特性材料抗干扰辨识***及方法 |
CN110244295A (zh) * | 2019-07-17 | 2019-09-17 | 广州大学 | 基于极化探地雷达的地下线性目标走向测量方法及*** |
CN111505628A (zh) * | 2020-04-29 | 2020-08-07 | 中国南方电网有限责任公司超高压输电公司广州局 | 基于探地雷达对地下电缆成像的探测与识别方法 |
CN111751392A (zh) * | 2020-07-30 | 2020-10-09 | 广州大学 | 一种基于双极化探地雷达的钢筋锈蚀检测方法 |
CN111679275A (zh) * | 2020-08-06 | 2020-09-18 | 中南大学 | 基于探地雷达的地下管线识别方法 |
CN112882019A (zh) * | 2021-01-14 | 2021-06-01 | 长春工程学院 | 基于旋转式单极化探地雷达的全极化目标辨识分类方法 |
CN112862888A (zh) * | 2021-01-29 | 2021-05-28 | 广东电网有限责任公司 | 一种电缆的定位方法、装置、计算机设备和存储介质 |
Non-Patent Citations (2)
Title |
---|
张守伟;: "基于地质雷达的地下管线探测与正演模拟", 环境技术, no. 01, 25 February 2020 (2020-02-25) * |
苏兆锋;陈昌彦;肖敏;: "探地雷达探测地下管线的机理和应用研究", 岩土工程技术, no. 04, 8 August 2013 (2013-08-08), pages 192 - 194 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116413719A (zh) * | 2023-06-12 | 2023-07-11 | 深圳大学 | 一种月球次表层下方熔岩管道勘探方法及相关设备 |
CN116413719B (zh) * | 2023-06-12 | 2023-09-08 | 深圳大学 | 一种月球次表层下方熔岩管道勘探方法及相关设备 |
Also Published As
Publication number | Publication date |
---|---|
CN114675263B (zh) | 2024-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110376584B (zh) | 基于探地雷达图像特征信号识别的供水管道漏损检测方法 | |
CN109143193B (zh) | 一种基于二阶峰的高频地波雷达浪场反演方法 | |
CN114675263A (zh) | 一种使用双极化探地雷达的地下管线材质识别方法 | |
CN109343022B (zh) | 估测层间土壤含水量的方法 | |
CN109975749B (zh) | 一种校正源存在条件下的短波单站直接定位方法 | |
CN110596787A (zh) | 一种基于x波段全固态双偏振雨量雷达的降水估测方法 | |
CN110031916A (zh) | 一种基于星地链路衰减效应的降雨强度测量方法 | |
CN102401888B (zh) | 一种电磁矢量传感器阵列耦合误差的自校正方法 | |
CN113945928B (zh) | 一种基于无旋性区域的全极化sar定标方法 | |
CN109490641B (zh) | 一种中纬地区偶发e层短波场强的计算方法 | |
CN114910982B (zh) | 一种基于北斗技术的降雨预警模型构建方法 | |
CN110907904A (zh) | 一种太阳法方位波束宽度修正方法 | |
CN104820222A (zh) | 雨量雷达实时动态衰减订正与降水计算方法 | |
CN110018448A (zh) | 基于任意取向微波链路的双偏振雷达衰减订正方法 | |
CN113254850A (zh) | 一种河道淤泥量监测计算方法 | |
CN105353374A (zh) | 一种用于自旋目标的单频雷达成像方法 | |
CN106814350A (zh) | 基于压缩感知的外辐射源雷达参考信号信杂比估计方法 | |
CN115718284A (zh) | 基于圆极化反射对称性的共极化通道不平衡定标方法 | |
CN116029162A (zh) | 利用星载gnss-r数据的洪涝灾害淹没范围监测方法和*** | |
CN204315734U (zh) | 一种正交偏馈抛物面天线宽带信号极化分离装置 | |
CN115047416A (zh) | 一种基于对称区域表面扭曲成分的全极化sar定标方法 | |
JP6381856B2 (ja) | レーダ信号処理装置及びレーダ信号処理方法 | |
Li et al. | Research on methods of high coherent target extraction in urban area based on PSINSAR technology | |
CN111337551A (zh) | 双线极化天线gnss土壤湿度测量方法 | |
CN111144226A (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 |