CN113538536B - 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 - Google Patents
一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 Download PDFInfo
- Publication number
- CN113538536B CN113538536B CN202110827584.XA CN202110827584A CN113538536B CN 113538536 B CN113538536 B CN 113538536B CN 202110827584 A CN202110827584 A CN 202110827584A CN 113538536 B CN113538536 B CN 113538536B
- Authority
- CN
- China
- Prior art keywords
- image block
- remote sensing
- image
- optical image
- sar
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种SAR图像信息辅助的遥感光学图像浓云检测方法及系统,方法包括S1、将遥感SAR图像和遥感光学图像均匀划分为图像块;S2、计算SAR图像块和光学图像块的HOG特征;S3、对SAR图像块和光学图像块的HOG特征进行融合;S4、将图像块的融合特征输入预先训练好的检测网络进行浓云检测;S5、若发现包含浓云的光学图像块则采用对应的SAR图像块替换。本发明借助了SAR图像信息,提升了浓云检测和去除能力。
Description
技术领域
本发明涉及遥感图像处理技术领域,更具体地说,特别涉及一种SAR图像信息辅助的遥感光学图像浓云检测方法及系统。
背景技术
受多云天气的影响,遥感光学图像通常被云层遮挡,导致无法获取清晰的遥感影像,影响对遥感影像中关键地物的判读。根据云层的薄浓,可分为薄云和浓云两种情况,薄云降低了光学图像的清晰度,浓云则使覆盖的光学图像区域完全不可见。相比于薄云而言,浓云对光学图像的影响更大,在对遥感光学图像进行判读之前必须去除。
合成孔径雷达(SAR)图像不受天气、云雾影响。基于SAR图像提供的信息能够实现光学图像浓云的检测与去除,对比文件“DarbaghshahiFN,MohammadiMR,Soryani M.Cloudremoval in remote sensing images using generative adversarial networks andSAR-to-optical image translation[J].2020.www.arXiv.org.com”提出基于生成对抗网络和SAR图像的遥感光学浓云去除方法,通过生成对抗网络联合SAR图像和光学图像获得去云后的光学图像,能够有效去除浓云。然而,上述方法在整幅图像上使用了两次生成对抗网络,计算复杂度高,运行效率低。为此,有必要对现有技术进行改进。
发明内容
本发明的目的在于提供一种SAR图像信息辅助的遥感光学图像浓云检测方法及系统,以克服现有技术所存在的缺陷。
为了达到上述目的,本发明采用的技术方案如下:
一种SAR图像信息辅助的遥感光学图像浓云检测方法,包括以下步骤:
S1、将遥感SAR图像和遥感光学图像均匀划分为图像块;
S2、计算SAR图像块和光学图像块的HOG特征;
S3、对SAR图像块和光学图像块的HOG特征进行融合;
S4、将图像块的融合特征输入预先训练好的检测网络进行浓云检测;
S5、若发现包含浓云的光学图像块则采用对应的SAR图像块替换。
进一步地,在步骤S1中遥感SAR图像Is和遥感光学图像为Io的大小为W×H,其中,W为图像的宽度,H为图像的高度,所述遥感SAR图像和遥感光学图像均有个图像块,其中为上取整运算,w为设定的图像块大小。
进一步地,所述步骤S2具体包括:
S21、计算每一个SAR图像块的HOG特征,第i个图像块的特征记为fs,i;
S22、计算每一个光学图像块的HOG特征,第i个图像块的特征记为fo,i。
进一步地,每一个SAR图像块和每一个光学图像块的HOG特征计算过程如下:
首先、根据以下公式计算每一个SAR图像块和每一个光学图像块的像素的梯度;
式中p(x,y)为图像块(x,y)点的像素值,g(x,y)为该点的梯度值,θ(x,y)为该点的梯度幅度;
然后、统计每一个SAR图像块和每一个光学图像块的梯度直方图,得到SAR图像块和光学图像块的HOG特征。
进一步地,所述步骤S3中的融合采用以下公式:
fi=[fs,i,fo,i].
进一步地,所述步骤S4中浓云检测的计算过程如下:
di=Net(fi)
其中,Net为检测网络,检测网络的输出为0和1两种结果,di为第i个图像块的浓云检测结果,di=0表示第i个光学图像块无浓云,di=1表示第i个光学图像块存在浓云。
进一步地,所述步骤S4中检测网络的建立和训练过程如下:
S41、获取多幅同个区域的遥感SAR图像和遥感光学图像的数据集;
S42、对数据集中的同个区域遥感SAR图像和遥感光学图像对进行图像块划分,并根据遥感图像的浓云区域对存在浓云的图像块添加标签;
S43、使用步骤S2-S3计算图像块的融合特征,并存储融合特征和标签数据;
S44、根据图像块大小设计浓云检测网络,采用三层的全连接网络,检测网络可以下式表示:
Net(x)=sigmoid(W2relu(W1x+b1)+b2)
其中W1、W2分别为网络的第一层和第二层权重矩阵,b1、b2分为为第一层权重矩阵和第二层权重矩阵对应的偏置项;relu为采用的修正线性单元激活函数;sigmoid()为采用的sigmoid激活函数;
S45、将交叉熵作为损失函数,使用随机梯度训练浓云检测网络。
进一步地,在步骤S42中根据遥感图像的浓云区域对存在浓云的图像块添加标签1,反之添加标签0。
本发明还提供一种用于实现上述的SAR图像信息辅助的遥感光学图像浓云检测方法的系统,包括:
划分模块,用于将遥感SAR图像和遥感光学图像均匀划分为图像块;
计算模块,用于计算SAR图像块和光学图像块的HOG特征;
融合模块,用于对SAR图像块和光学图像块的HOG特征进行融合;
浓云检测模块,用于将图像块的融合特征输入预先训练好的检测网络进行浓云检测;
替换模块,用于在发现包含浓云的光学图像块时则采用对应的SAR图像块替换;
所述划分模块、计算模块、融合模块、浓云检测模块和替换模块依次连接。
与现有技术相比,本发明的优点在于:本发明首先对同个区域的遥感SAR图像和受浓云影响的遥感光学图像进行均匀划分为图像块,之后同时提取每一图像块的SAR图像和光学图像特征,采用前馈神经网络进行浓云的检测,获得每一图像块是否包含浓云的检测结果;最后对包含浓云的光学图像区域,使用SAR图像进行替换以达到浓云去除的目的。本发明借助了SAR图像信息,提升了浓云检测和去除能力;本发明提供的浓云检测结果可作为遥感光学图像质量评估的依据以及浓云去除的引导信息。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明SAR图像信息辅助的遥感光学图像浓云检测方法流程图。
图2是本发明SAR图像信息辅助的遥感光学图像浓云检测系统框架图。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
参阅图1所示,本实施例公开了一种SAR图像信息辅助的遥感光学图像浓云检测方法,包括以下步骤:
步骤S1、将遥感SAR图像和遥感光学图像均匀划分为图像块。
具体的,记遥感SAR图像为Is,遥感光学图像为Io,图像Is和Io的大小为W×H,其中W为图像的宽度,H为图像的高度;w为设定的图像块大小;遥感SAR图像和遥感光学图像均有个图像块,其中,为上取整运算。
步骤S2、计算SAR图像块和光学图像块的HOG特征,具体为:
步骤S21、计算每一个SAR图像块的HOG特征,第i个图像块的特征记为fs,i;
步骤S22、计算每一个光学图像块的HOG特征,第i个图像块的特征记为fo,i;
其中,图像块的HOG特征计算过程如下:
首先,根据以下公式计算图像块像素的梯度:
式中,p(x,y)为图像块(x,y)点的像素值,g(x,y)为该点的梯度值,θ(x,y)为该点的梯度幅度;
然后,统计图像块的梯度直方图,得到各图像块的HOG特征。
步骤S3、对SAR图像块和光学图像块的HOG特征进行融合,公式如下:
fi=[fs,i,fo,i]
步骤S4、将图像块的融合特征输入预先训练好的检测网络进行浓云检测,检测的计算过程如下:
di=Net(fi)
其中,Net为检测网络,检测网络的输出为0和1两种结果,di为第i个图像块的浓云检测结果,di=0表示第i个光学图像块无浓云,di=1表示第i个光学图像块存在浓云。
具体的,检测网络的建立和训练过程如下:
步骤S41、获取多幅同个区域的遥感SAR图像和遥感光学图像的数据集。
步骤S42、对数据集中的同个区域遥感SAR图像和遥感光学图像对进行图像块划分,并根据遥感图像的浓云区域对存在浓云的图像块添加标签1,反之添加标签0。
步骤S43、使用步骤S2-S3计算图像块的融合特征,并存储融合特征和标签数据。
步骤S44、根据图像块大小设计浓云检测网络,采用三层的全连接网络,检测网络可以下式表示:
Net(x)=sigmoid(W2relu(W1x+b1)+b2)
其中W1、W2分别为网络的第一层和第二层权重矩阵,b1、b2分为为第一层权重矩阵和第二层权重矩阵对应的偏置项;relu为采用的修正线性单元激活函数;sigmoid()为采用的sigmoid激活函数,是深度学习里面常用的一个函数。
步骤S45、将交叉熵作为损失函数,使用随机梯度训练浓云检测网络。
步骤S5、若发现包含浓云的光学图像块则采用对应的SAR图像块替换。
参阅图2所示,本发明还提供一种用于实现上述的SAR图像信息辅助的遥感光学图像浓云检测方法的系统,包括:划分模块1,用于将遥感SAR图像和遥感光学图像均匀划分为图像块;计算模块2,用于计算SAR图像块和光学图像块的HOG特征;融合模块3,用于对SAR图像块和光学图像块的HOG特征进行融合;浓云检测模块4,用于将图像块的融合特征输入预先训练好的检测网络进行浓云检测;替换模块5,用于在发现包含浓云的光学图像块时则采用对应的SAR图像块替换;所述的划分模块1、计算模块2、融合模块3、浓云检测模块4和替换模块5依次连接。
本发明借助了SAR图像信息,增加提升了浓云检测和去除能力,本发明提供的浓云检测结果可作为遥感光学图像质量评估的依据以及浓云去除的引导信息。
虽然结合附图描述了本发明的实施方式,但是专利所有者可以在所附权利要求的范围之内做出各种变形或修改,只要不超过本发明的权利要求所描述的保护范围,都应当在本发明的保护范围之内。
Claims (4)
1.一种SAR图像信息辅助的遥感光学图像浓云检测方法,其特征在于,包括以下步骤:
S1、将遥感SAR图像和遥感光学图像均匀划分为图像块;
S2、计算SAR图像块和光学图像块的HOG特征;
S3、对SAR图像块和光学图像块的HOG特征进行融合;
S4、将图像块的融合特征输入预先训练好的检测网络进行浓云检测;
S5、若发现包含浓云的光学图像块则采用对应的SAR图像块替换;
所述步骤S2具体包括:
S21、计算每一个SAR图像块的HOG特征,第i个图像块的特征记为fs,i;
S22、计算每一个光学图像块的HOG特征,第i个图像块的特征记为fo,i;
每一个SAR图像块和每一个光学图像块的HOG特征计算过程如下:
首先、根据以下公式计算每一个SAR图像块和每一个光学图像块的像素的梯度;
式中p(x,y)为图像块(x,y)点的像素值,g(x,y)为该点的梯度值,θ(x,y)为该点的梯度幅度;
然后、统计每一个SAR图像块和每一个光学图像块的梯度直方图,得到SAR图像块和光学图像块的HOG特征;
所述步骤S3中的融合采用以下公式:
fi=[fs,i,fo,i];
所述步骤S4中浓云检测的计算过程如下:
di=Net(fi)
其中,Net为检测网络,检测网络的输出为0和1两种结果,di为第i个图像块的浓云检测结果,di=0表示第i个光学图像块无浓云,di=1表示第i个光学图像块存在浓云;
所述步骤S4中检测网络的建立和训练过程如下:
S41、获取多幅同个区域的遥感SAR图像和遥感光学图像的数据集;
S42、对数据集中的同个区域遥感SAR图像和遥感光学图像对进行图像块划分,并根据遥感图像的浓云区域对存在浓云的图像块添加标签;
S43、使用步骤S2-S3计算图像块的融合特征,并存储融合特征和标签数据;
S44、根据图像块大小设计浓云检测网络,采用三层的全连接网络,检测网络用下式表示:
Net(x)=sigmoid(W2relu(W1x+b1)+b2)
其中W1、W2分别为网络的第一层和第二层权重矩阵,b1、b2分别为第一层权重矩阵和第二层权重矩阵对应的偏置项;relu为采用的修正线性单元激活函数;sigmoid()为采用的sigmoid激活函数;
S45、将交叉熵作为损失函数,使用随机梯度训练浓云检测网络。
3.根据权利要求1所述的SAR图像信息辅助的遥感光学图像浓云检测方法,其特征在于,在步骤S42中根据遥感图像的浓云区域对存在浓云的图像块添加标签1,反之添加标签0。
4.一种用于实现权利要求1-3任意一项所述的SAR图像信息辅助的遥感光学图像浓云检测方法的系统,其特征在于,包括:
划分模块,用于将遥感SAR图像和遥感光学图像均匀划分为图像块;
计算模块,用于计算SAR图像块和光学图像块的HOG特征;
融合模块,用于对SAR图像块和光学图像块的HOG特征进行融合;
浓云检测模块,用于将图像块的融合特征输入预先训练好的检测网络进行浓云检测;
替换模块,用于在发现包含浓云的光学图像块时则采用对应的SAR图像块替换;
所述划分模块、计算模块、融合模块、浓云检测模块和替换模块依次连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110827584.XA CN113538536B (zh) | 2021-07-21 | 2021-07-21 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110827584.XA CN113538536B (zh) | 2021-07-21 | 2021-07-21 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113538536A CN113538536A (zh) | 2021-10-22 |
CN113538536B true CN113538536B (zh) | 2022-06-07 |
Family
ID=78088562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110827584.XA Active CN113538536B (zh) | 2021-07-21 | 2021-07-21 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113538536B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680542A (zh) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | 基于在线学习的遥感影像变化检测方法 |
CN107590782A (zh) * | 2017-08-21 | 2018-01-16 | 西北工业大学 | 一种基于全卷积网络的高分辨率光学图像厚云去除方法 |
CN110097101A (zh) * | 2019-04-19 | 2019-08-06 | 大连海事大学 | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 |
CN110443201A (zh) * | 2019-08-06 | 2019-11-12 | 哈尔滨工业大学 | 基于多源影像联合形状分析与多属性融合的目标识别方法 |
CN111383192A (zh) * | 2020-02-18 | 2020-07-07 | 清华大学 | 一种融合sar的可见光遥感图像去雾方法 |
CN111915592A (zh) * | 2020-08-04 | 2020-11-10 | 西安电子科技大学 | 基于深度学习的遥感图像云检测方法 |
CN112381815A (zh) * | 2020-11-26 | 2021-02-19 | 浙江易智信息技术有限公司 | 多维融合遥感大数据去云方法 |
CN112819688A (zh) * | 2021-02-01 | 2021-05-18 | 西安研硕信息技术有限公司 | 一种sar图像转类光学图像的转换方法及系统 |
CN112861672A (zh) * | 2021-01-27 | 2021-05-28 | 电子科技大学 | 一种基于光学-sar异源遥感图像匹配方法 |
CN113012209A (zh) * | 2021-02-23 | 2021-06-22 | 中国科学院空天信息创新研究院 | 融合sar图像信息的光学图像去云方法 |
-
2021
- 2021-07-21 CN CN202110827584.XA patent/CN113538536B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680542A (zh) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | 基于在线学习的遥感影像变化检测方法 |
CN107590782A (zh) * | 2017-08-21 | 2018-01-16 | 西北工业大学 | 一种基于全卷积网络的高分辨率光学图像厚云去除方法 |
CN110097101A (zh) * | 2019-04-19 | 2019-08-06 | 大连海事大学 | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 |
CN110443201A (zh) * | 2019-08-06 | 2019-11-12 | 哈尔滨工业大学 | 基于多源影像联合形状分析与多属性融合的目标识别方法 |
CN111383192A (zh) * | 2020-02-18 | 2020-07-07 | 清华大学 | 一种融合sar的可见光遥感图像去雾方法 |
CN111915592A (zh) * | 2020-08-04 | 2020-11-10 | 西安电子科技大学 | 基于深度学习的遥感图像云检测方法 |
CN112381815A (zh) * | 2020-11-26 | 2021-02-19 | 浙江易智信息技术有限公司 | 多维融合遥感大数据去云方法 |
CN112861672A (zh) * | 2021-01-27 | 2021-05-28 | 电子科技大学 | 一种基于光学-sar异源遥感图像匹配方法 |
CN112819688A (zh) * | 2021-02-01 | 2021-05-18 | 西安研硕信息技术有限公司 | 一种sar图像转类光学图像的转换方法及系统 |
CN113012209A (zh) * | 2021-02-23 | 2021-06-22 | 中国科学院空天信息创新研究院 | 融合sar图像信息的光学图像去云方法 |
Non-Patent Citations (2)
Title |
---|
光学遥感图像中复杂海背景下的舰船检测;王慧利等;《光学精密工程》;20180315;第26卷(第03期);第723-732页 * |
改进CGAN网络的光学遥感图像云去除方法;裴傲等;《农业工程学报》;20200723;第36卷(第14期);第194-202页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113538536A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Multi-level fusion and attention-guided CNN for image dehazing | |
CN110119703B (zh) | 一种安防场景下融合注意力机制和时空图卷积神经网络的人体动作识别方法 | |
CN108805002B (zh) | 基于深度学习和动态聚类的监控视频异常事件检测方法 | |
CN113592736B (zh) | 一种基于融合注意力机制的半监督图像去模糊方法 | |
WO2019136591A1 (zh) | 基于弱监督时空级联神经网络的显著目标检测方法及系统 | |
CN111723693B (zh) | 一种基于小样本学习的人群计数方法 | |
CN108133188A (zh) | 一种基于运动历史图像与卷积神经网络的行为识别方法 | |
CN108564025A (zh) | 一种基于可变形卷积神经网络的红外图像物体识别方法 | |
CN109558811B (zh) | 一种基于运动前景关注及非监督的关键帧提取的动作识别方法 | |
CN113065645B (zh) | 孪生注意力网络、图像处理方法和装置 | |
CN110097028B (zh) | 基于三维金字塔图像生成网络的人群异常事件检测方法 | |
CN112990077B (zh) | 基于联合学习与光流估计的面部动作单元识别方法及装置 | |
CN114972976B (zh) | 基于频域自注意力机制的夜间目标检测、训练方法及装置 | |
CN113128360A (zh) | 基于深度学习的司机驾驶行为检测与识别方法 | |
WO2021243947A1 (zh) | 对象再识别方法及装置、终端和存储介质 | |
CN110458115B (zh) | 一种基于时序的多帧集成目标检测算法 | |
CN117727069A (zh) | 基于多尺度信息交互网络的文本-图像行人重识别方法 | |
CN115359407A (zh) | 一种视频中的多车辆跟踪方法 | |
CN111192206A (zh) | 一种提高图像清晰度的方法 | |
CN116363535A (zh) | 基于卷积神经网络的无人机航拍影像中的船舶检测方法 | |
CN112396036A (zh) | 一种结合空间变换网络和多尺度特征提取的遮挡行人重识别方法 | |
CN111507353A (zh) | 一种基于文字识别的中文字段检测方法及系统 | |
CN113538536B (zh) | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 | |
CN112613564A (zh) | 一种剔除重叠框的目标检测后处理方法 | |
CN114219734A (zh) | 一种基于深度学习的sar辅助下遥感光学影像去云方法 |
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 |