CN109544495A - 一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 - Google Patents
一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 Download PDFInfo
- Publication number
- CN109544495A CN109544495A CN201811345251.8A CN201811345251A CN109544495A CN 109544495 A CN109544495 A CN 109544495A CN 201811345251 A CN201811345251 A CN 201811345251A CN 109544495 A CN109544495 A CN 109544495A
- Authority
- CN
- China
- Prior art keywords
- image
- low
- multispectral
- resolution
- processing module
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 42
- 230000009466 transformation Effects 0.000 title claims abstract description 33
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 33
- 238000012417 linear regression Methods 0.000 claims abstract description 22
- 230000004927 fusion Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims abstract description 11
- 238000002156 mixing Methods 0.000 claims abstract description 5
- 238000007499 fusion processing Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000003786 synthesis reaction Methods 0.000 claims description 11
- 239000004615 ingredient Substances 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 6
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法,该算法首先使用SoC芯片图像加速处理模块的高斯滤波功能快速获取全色图像和多光谱的低频信息,然后对多光谱图像进行多元线性回归处理,结合全色低频信息合成低分辨率全色图像,最后利用SoC芯片多核处理器对图像进行比值变换处理获取高保真融合图像,消除了融合结果中全色和多光谱原始图像中对空间细节和光谱的干扰信息。本方法利用SoC芯片的图像加速模块的快速滤波和多核处理功能,使得算法运算效率高,鲁棒性好,解决了全色图像和多光谱图像融合中存在的光谱信息和空间细节失真问题。
Description
技术领域
本发明涉及一种基于多核SoC芯片加速模块的图像处理方法,特别是一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法。
背景技术
随着高分辨率传感器的快速开发和部署,大多数的地球观测卫星,如WorldView-2卫星IKONOS、QuickBird和geoeye-2等,可以同时提供高分辨率的全色(Panchromatic,PAN)和低分辨率的多光谱图像(Multispectral,MS)。全色图像空间分辨率高,但只有一个光谱段,多光谱空间分辨率低,但有丰富的光谱信息。为了弥补全色图像的不足,卫星上一般同时搭载一个多光谱传感器(例如WorldView-2遥感卫星多光谱图像包含红、绿、蓝、近红外四个波段)。为了满足目标判读的军事需求,将全色与多光谱图像融合成一幅图像,以便从中快速提取情报信息已成为我国遥感图像应用关注的焦点。
尽管目前已有许多有关全色和多光谱图像融合的研究,但现有方法无法有效解决融合结果中光谱失真和纹理细节模糊的问题;并且,这些方法主要是针对小尺寸图像,计算量大,难以满足大尺寸图像融合的高时效要求,部分在技术领域领先的芯片公司根据这方面需求已开发具有图像加速处理功能的芯片。“精导芯二号”是一款具有高性能图像处理能力的芯片,能够对多尺寸图像进行校正和滤波处理,满足了目前大多数领域对图像预处理的需求。为兼容常用图像处理的应用,设计了多种数据来源通道,并配有多种功能模式,能够对图像进行线性和非线性校正处理,并包含高斯、梯度和中值等多种常用滤波器,功能复杂,使用涵盖面广泛。
发明内容
本发明的目的在于提供一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法,解决宽幅图像在融合处理过程中存在的运算速率低,光谱失真和空间细节扭曲的问题。
一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法,采用多核SoC芯片上的图像加速处理***实现,包括:高斯滤波分频处理模块、多元线性回归光谱图像处理模块、合成低分辨率图像模块和比值变换算法融合处理模块。
一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法的具体步骤为:
第一步高斯滤波分频处理模块用于提取图像低频信息
高斯滤波分频处理模块对图像进行高斯滤波分频处理,高斯函数具有低通性质,模拟人眼视觉机理,通过不同的尺度参数设定,能够模拟人眼远近处所看到的图像。定义获取的新图像为,其中是原始图像的每个像素的低频信息,如式(1)所示:
其中,x表示行,y表示列,orig(x,y)是原始图像中像素值,G(x,y,σ)是高斯滤波器,low(x,y)是orig(x,y)原始图像的低频信息。
使用式(1)对原始全色图像Pan和多光谱图像MSi进行滤波处理,获取全色图像的低频信息Panlf,多光谱的低频信息高斯核窗口和尺度参数影响提取空间细节成分,进而影响全色与多光谱图像空间相关权系数的求解:尺度参数越小提取的空间细节信息越少,反之尺度越大,提取的空间细节信息越多,但存在出现空间细节提取过度的现象,从而导致光谱失真。采用标准高斯分布的值即σ,经过大量实验数据测试,选取高斯核窗口为5x5或9x9,并使用高斯滤波分频处理模块对算法进行加速处理,提高处理效率。
第二步多元线性回归光谱图像处理模块用于获取多光谱图像空间细节信息
多元线性回归光谱图像处理模块采用光谱相关的回归模型合成低分辨率全色图像,如公式(2):
式中,是多光谱图像空间相关成分,其中,MSi表示多光谱图像的各波段,其中i=1,2,3,4…;表示多光谱各波段的低频信息的均值,imghl是多光谱图像在进行多元线性回归处理和高斯滤波后得到的多光谱图像的空间细节信息。其中的系数由式(3)多元线性回归计算得出:
其中,Pan是原始全色图像。
第三步合成低分辨率图像模块用以合成低分辨率图像
合成低分辨率图像模块利用高斯滤波和多元线性回归处理结果合成低分辨率图像。合成低分辨率图像将包含多光谱图像中低分辨率空间细节成分imghl,加上全色波段中近似于多光谱波段但多余多光谱波段的光谱成分Panlf,得到低分辨率图像imgsyn,如式(4)所示:
imgsyn=imghl+Panlf (4)
第四步比值变换算法融合处理模块用以计算图像融合结果
比值变换算法融合处理模块是将高分辨率全色图像和合成的低分辨率图像进行比值运算获得图像融合结果。通过比值变换计算,将全色图像和多光谱细节差异根据计算的比值结果按比例分布到多光谱图像的各个波段,从而生成和全色图像分辨率一致的高分辨率多光谱图像。合成低分辨率图像不仅包含了全色图像和多光谱图像的细节差异,还包含了两者之间的光谱差异,在进行比值变换融合处理后,能够减少融合结果的光谱扭曲,获得高保真融合图像,公式如(6)所示:
其中,Fusi表示各波段的全色与多光谱图像融合结果,MSi是原始多光谱波段,Pan是原始全色图像,imgsyn是全色与多光谱图像融合结果。
由于计算过程中,像素点处理相互独立,利用多核SoC芯片对其进行并行优化处理,提高算法处理效率。
至此,完成了基于高斯滤波和比值变换的SoC芯片图像处理融合。
本发明针对宽幅全色与多光谱图像融合存在光谱失真和计算耗时等问题,根据遥感卫星全色与多光谱波段特性,提出了一种基于高斯滤波和比值变换的高保真融合方法,有效剔除全色与多光谱波段带来的干扰信息,避免了融合处理过程中光谱失真和空间细节扭曲问题,获得高保真融合图像结果。并分析算法模型,利用“精导芯二号”的图像加速处理模块和多核处理器特性对算法进行了优化处理,提高了处理效率,较大的减少了运算时间。
附图说明
图1一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法的流程示意图。
1.高斯滤波分频处理模块 2.多元线性回归光谱图像处理模块 3.合成低分辨率图像模块4.比值变换算法融合处理模块
具体实施方式
一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法,采用多核SoC芯片上的图像加速处理***实现,包括:高斯滤波分频处理模块1、多元线性回归光谱图像处理模块2、合成低分辨率图像模块3和比值变换算法融合处理模块4。
高斯滤波分频处理模块1用于提取图像低频信息;
多元线性回归光谱图像处理模块2用于获取多光谱图像空间细节信息;
合成低分辨率图像模块3用于合成低分辨率图像;
比值变换算法融合处理模块4用于计算图像融合结果。
一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法的具体步骤为:
第一步高斯滤波分频处理模块1用于提取图像低频信息
高斯滤波分频处理模块1对图像进行高斯滤波分频处理,高斯函数具有低通性质,模拟人眼视觉机理,通过不同的尺度参数设定,能够模拟人眼远近处所看到的图像。定义获取的新图像为,其中是原始图像的每个像素的低频信息,如式(1)所示:
其中,x表示行,y表示列,orig(x,y)是原始图像中像素值,G(x,y,σ)是高斯滤波器,low(x,y)是orig(x,y)原始图像的低频信息。
使用式(1)对原始全色图像Pan和多光谱图像MSi进行滤波处理,获取全色图像的低频信息Panlf,多光谱的低频信息高斯核窗口和尺度参数影响提取空间细节成分,进而影响全色与多光谱图像空间相关权系数的求解:尺度参数越小提取的空间细节信息越少,反之尺度越大,提取的空间细节信息越多,但存在出现空间细节提取过度的现象,从而导致光谱失真。采用标准高斯分布的值即σ,经过大量实验数据测试,选取高斯核窗口为5x5或9x9,并使用高斯滤波分频处理模块1对算法进行加速处理,提高处理效率。
第二步多元线性回归光谱图像处理模块2用于获取多光谱图像空间细节信息
多元线性回归光谱图像处理模块2采用光谱相关的回归模型合成低分辨率全色图像,如公式(2):
式中,是多光谱图像空间相关成分,其中,MSi表示多光谱图像的各波段,其中i=1,2,3,4…;表示多光谱各波段的低频信息的均值,imghl是多光谱图像在进行多元线性回归处理和高斯滤波后得到的多光谱图像的空间细节信息。其中的系数由式(3)多元线性回归计算得出:
其中,Pan是原始全色图像。
第三步合成低分辨率图像模块3用以合成低分辨率图像
合成低分辨率图像模块3利用高斯滤波和多元线性回归处理结果合成低分辨率图像。合成低分辨率图像将包含多光谱图像中低分辨率空间细节成分imghl,加上全色波段中近似于多光谱波段但多余多光谱波段的光谱成分Panlf,得到低分辨率图像imgsyn,如式(4)所示:
imgsyn=imghl+Panlf (4)
第四步比值变换算法融合处理模块4用以计算图像融合结果
比值变换算法融合处理模块4是将高分辨率全色图像和合成的低分辨率图像进行比值运算获得图像融合结果。通过比值变换计算,将全色图像和多光谱细节差异根据计算的比值结果按比例分布到多光谱图像的各个波段,从而生成和全色图像分辨率一致的高分辨率多光谱图像。合成低分辨率图像不仅包含了全色图像和多光谱图像的细节差异,还包含了两者之间的光谱差异,在进行比值变换融合处理后,能够减少融合结果的光谱扭曲,获得高保真融合图像,公式如(6)所示:
其中,Fusi表示各波段的全色与多光谱图像融合结果,MSi是原始多光谱波段,Pan是原始全色图像,imgsyn是全色与多光谱图像融合结果。
由于计算过程中,像素点处理相互独立,利用多核SoC芯片对其进行并行优化处理,提高算法处理效率。
至此,完成了基于高斯滤波和比值变换的SoC芯片图像处理融合。
Claims (2)
1.一种多核SoC芯片上的图像加速处理***,其特征在于包括:高斯滤波分频处理模块(1)、多元线性回归光谱图像处理模块(2)、合成低分辨率图像模块(3)和比值变换算法融合处理模块(4);
其中,高斯滤波分频处理模块(1)用于提取图像低频信息;
多元线性回归光谱图像处理模块(2)用于获取多光谱图像空间细节信息;
合成低分辨率图像模块(3)用于合成低分辨率图像;
比值变换算法融合处理模块(4)用于计算图像融合结果。
2.一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法,其特征在于具体步骤为:
第一步 高斯滤波分频处理模块(1)用于提取图像低频信息
高斯滤波分频处理模块(1)对图像进行高斯滤波分频处理,高斯函数具有低通性质,模拟人眼视觉机理,通过不同的尺度参数设定,能够模拟人眼远近处所看到的图像;定义获取的新图像为,其中是原始图像的每个像素的低频信息,如式(1)所示:
其中,x表示行,y表示列,orig(x,y)是原始图像中像素值,G(x,y,σ)是高斯滤波器,low(x,y)是orig(x,y)原始图像的低频信息;
使用式(1)对原始全色图像Pan和多光谱图像MSi进行滤波处理,获取全色图像的低频信息Panlf,多光谱的低频信息高斯核窗口和尺度参数影响提取空间细节成分,进而影响全色与多光谱图像空间相关权系数的求解:尺度参数越小提取的空间细节信息越少,反之尺度越大,提取的空间细节信息越多,但存在出现空间细节提取过度的现象,从而导致光谱失真;采用标准高斯分布的值即σ,经过大量实验数据测试,选取高斯核窗口为5x5或9x9,并使用高斯滤波分频处理模块(1)对算法进行加速处理,提高处理效率;
第二步 多元线性回归光谱图像处理模块(2)用于获取多光谱图像空间细节信息
多元线性回归光谱图像处理模块(2)采用光谱相关的回归模型合成低分辨率全色图像,如公式(2):
式中,是多光谱图像空间相关成分,其中,MSi表示多光谱图像的各波段,其中i=1,2,3,4…;表示多光谱各波段的低频信息的均值,imghl是多光谱图像在进行多元线性回归处理和高斯滤波后得到的多光谱图像的空间细节信息;其中的系数由式(3)多元线性回归计算得出:
其中,Pan是原始全色图像;
第三步 合成低分辨率图像模块(3)用以合成低分辨率图像
合成低分辨率图像模块(3)利用高斯滤波和多元线性回归处理结果合成低分辨率图像;合成低分辨率图像将包含多光谱图像中低分辨率空间细节成分imghl,加上全色波段中近似于多光谱波段但多余多光谱波段的光谱成分Panlf,得到低分辨率图像imgsyn,如式(4)所示:
imgsyn=imghl+Panlf (4)
第四步 比值变换算法融合处理模块(4)用以计算图像融合结果
比值变换算法融合处理模块(4)是将高分辨率全色图像和合成的低分辨率图像进行比值运算获得图像融合结果;通过比值变换计算,将全色图像和多光谱细节差异根据计算的比值结果按比例分布到多光谱图像的各个波段,从而生成和全色图像分辨率一致的高分辨率多光谱图像;合成低分辨率图像不仅包含了全色图像和多光谱图像的细节差异,还包含了两者之间的光谱差异,在进行比值变换融合处理后,能够减少融合结果的光谱扭曲,获得高保真融合图像,公式如(6)所示:
其中,Fusi表示各波段的全色与多光谱图像融合结果,MSi是原始多光谱波段,Pan是原始全色图像,imgsyn是全色与多光谱图像融合结果;
由于计算过程中,像素点处理相互独立,利用多核SoC芯片对其进行并行优化处理,提高算法处理效率;
至此,完成了基于高斯滤波和比值变换的SoC芯片图像处理融合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811345251.8A CN109544495B (zh) | 2018-11-13 | 2018-11-13 | 一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811345251.8A CN109544495B (zh) | 2018-11-13 | 2018-11-13 | 一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109544495A true CN109544495A (zh) | 2019-03-29 |
CN109544495B CN109544495B (zh) | 2023-05-23 |
Family
ID=65846961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811345251.8A Active CN109544495B (zh) | 2018-11-13 | 2018-11-13 | 一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109544495B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113570536A (zh) * | 2021-07-31 | 2021-10-29 | 中国人民解放军61646部队 | 基于cpu和gpu协同处理的全色和多光谱影像实时融合方法 |
CN113689371A (zh) * | 2021-07-27 | 2021-11-23 | 上海浦东发展银行股份有限公司 | 图像融合方法、装置、计算机设备和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104851077A (zh) * | 2015-06-03 | 2015-08-19 | 四川大学 | 一种自适应的遥感图像全色锐化方法 |
US20150288950A1 (en) * | 2013-08-16 | 2015-10-08 | University Of New Brunswick | Camera imaging systems and methods |
CN107016641A (zh) * | 2016-09-23 | 2017-08-04 | 北京航空航天大学 | 一种基于改进比值变换的全色与高光谱图像融合方法 |
-
2018
- 2018-11-13 CN CN201811345251.8A patent/CN109544495B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150288950A1 (en) * | 2013-08-16 | 2015-10-08 | University Of New Brunswick | Camera imaging systems and methods |
CN104851077A (zh) * | 2015-06-03 | 2015-08-19 | 四川大学 | 一种自适应的遥感图像全色锐化方法 |
CN107016641A (zh) * | 2016-09-23 | 2017-08-04 | 北京航空航天大学 | 一种基于改进比值变换的全色与高光谱图像融合方法 |
Non-Patent Citations (1)
Title |
---|
李旭 等: ""WorldView-2遥感图像融合新方法"", 《电子科技大学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113689371A (zh) * | 2021-07-27 | 2021-11-23 | 上海浦东发展银行股份有限公司 | 图像融合方法、装置、计算机设备和存储介质 |
CN113689371B (zh) * | 2021-07-27 | 2024-04-26 | 上海浦东发展银行股份有限公司 | 图像融合方法、装置、计算机设备和存储介质 |
CN113570536A (zh) * | 2021-07-31 | 2021-10-29 | 中国人民解放军61646部队 | 基于cpu和gpu协同处理的全色和多光谱影像实时融合方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109544495B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Pan-sharpening using an efficient bidirectional pyramid network | |
US11823363B2 (en) | Infrared and visible light fusion method | |
CN107123089B (zh) | 基于深度卷积网络的遥感图像超分辨重建方法及*** | |
US20220044375A1 (en) | Saliency Map Enhancement-Based Infrared and Visible Light Fusion Method | |
CN111080567B (zh) | 基于多尺度动态卷积神经网络的遥感图像融合方法及*** | |
US11830222B2 (en) | Bi-level optimization-based infrared and visible light fusion method | |
CN104683767B (zh) | 透雾图像生成方法及装置 | |
JP2019523509A (ja) | 暗視赤外画像における顕著性に基づく道路オブジェクト抽出方法 | |
CN112819737B (zh) | 基于3d卷积的多尺度注意力深度卷积网络的遥感图像融合方法 | |
CN103927741A (zh) | 增强目标特征的sar图像合成方法 | |
CN102509262B (zh) | 一种遥感图像薄云去除方法 | |
CN108269244B (zh) | 一种基于深度学习和先验约束的图像去雾*** | |
CN103248906B (zh) | 一种双目立体视频序列的深度图获取方法与*** | |
CN108154493B (zh) | 一种基于fpga的双波段红外图像伪彩融合算法 | |
CN105096286A (zh) | 遥感图像的融合方法及装置 | |
CN116309070A (zh) | 一种高光谱遥感图像超分辨率重建方法、装置及计算机设备 | |
CN105303542B (zh) | 基于梯度加权的自适应sfim图像融合算法 | |
CN113139902A (zh) | 一种高光谱图像超分辨率重建方法、装置及电子设备 | |
CN105160647A (zh) | 一种全色多光谱影像融合方法 | |
CN109544495A (zh) | 一种基于高斯滤波和比值变换的SoC芯片图像处理融合方法 | |
CN103955896A (zh) | 一种提升卫星图像视觉效果的真彩色增强方法 | |
CN105427268A (zh) | 一种中长波双波段红外图像特征级彩色融合方法 | |
CN107169946A (zh) | 基于非负稀疏矩阵与超球面彩色变换的图像融合方法 | |
Zhong et al. | Attention_FPNet: Two-branch remote sensing image pansharpening network based on attention feature fusion | |
CN111563866A (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 |