CN109753996B - 基于三维轻量化深度网络的高光谱图像分类方法 - Google Patents
基于三维轻量化深度网络的高光谱图像分类方法 Download PDFInfo
- Publication number
- CN109753996B CN109753996B CN201811542534.1A CN201811542534A CN109753996B CN 109753996 B CN109753996 B CN 109753996B CN 201811542534 A CN201811542534 A CN 201811542534A CN 109753996 B CN109753996 B CN 109753996B
- Authority
- CN
- China
- Prior art keywords
- data
- dimensional
- network
- hyperspectral image
- layer
- 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
Landscapes
- Image Analysis (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明提出的一种基于三维轻量化深度网络的高光谱图像分类方法,通过构建三位轻量化深度网络,实现了有限样本条件下,高光谱图像深度特征的自主提取,高精度的分类。本发明与现有的基于深度学习的高光谱图像分类方法相比,网络模型更深,精度更高,且参数量更少。
Description
技术领域
本发明涉及一种基于三维轻量化深度网络的高光谱图像分类方法,属遥感图像处理领域。
背景技术
高光谱图像同时包含光谱信息和空间信息,在军事及民用领域都有着重要的应用。然而,高光谱图像的高维特性、波段间高度相关性、光谱混合等使得高光谱图像分类面临巨大挑战。近年来,随着深度学习新技术的出现,基于深度学习的高光谱图像分类在方法和性能上得到了突破性的进展。但是,深度学习的模型参数较多,需要大量的训练样本。在深度学习相关技术中,一般来说,有效的增加网络深度对于提升网络的性能非常重要。然而高光谱图像有标注的样本相对较少,难以完全满足深层次深度学习模型的训练,容易出现过拟合问题。因此,研究适于高光谱图像高精度分类的深度模型一直是一项具有挑战性的任务,迫切需要一种算法能够在有限样本条件下,实现高精度的高光谱图像分类。
高光谱图像分类问题旨在给定一副具有部分标注像素的图像,通过相关算法,预测出图像中所有像素对应的具体地物类别。传统的高光谱图像分类方法一般利用人工预设特征,例如SIFT,HOG,PHOG等,从高光谱图像中提取特征,然后借由多层感知器,支持向量机等模型来分类。但是这些人工预设特征的设计及选取依赖专业性知识,且很难选取一种有通用性的特征。
近年来,随着深度学习的兴起,完全数据驱动和不需要先验知识的深度神经网络在图像处理和计算机视觉等领域展现出了突出的优势,其应用范围涵盖了高级图像识别、中低级图像处理等诸多方面,例如目标识别、检测、分类和图像去噪、动态去模糊、重建等等。高光谱图像分类领域也引入了深度学习的相关技术,且取得了明显优于传统方法的分类效果。但是,受高光谱图像训练样本数量的限制,应用在高光谱图像分类中的深度学习模型相对较浅,尽管计算机视觉方面大量实验已经表明有效的增加深度对于提升分类表现非常有益。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,针对高光谱图像分类问题,结合深度学习相关技术,设计一种在有限样本条件下能够实现高精度分类的算法。
技术方案
一种基于三维轻量化深度网络的高光谱图像分类方法,其特征在于步骤如下:
步骤1:数据预处理:将待处理的高光谱图像数据进行最大最小归一化,归一化公式如下:
步骤2:数据划分:统计待处理的高光谱图像中,各个类别标注样本的个数,然后按照比例,从每个类别中,抽取5%-10%的标注样本作为训练数据,剩余的部分作为测试数据;具体做法如下,对于一个尺寸为M×N×L的三维高光谱图像数据,M,N分别表示高光谱图像的高和宽,L表示数据的波段数;抽取样本时,以待处理像素为中心,抽取S×S×L的数据块作为中心像素的样本数据,S表示邻域大小;
步骤3:构建三维轻量化深度网络模型:设计网络依次包含了两部分结构:1)特征提取部分,该部分包含了基于不对称三维卷积网络的数据处理模块及多组轻量化模块构成的特征提取模块,所述的数据处理模块包含一个不对称的三维卷积层,归一化层,激励函数及池化层,其中不对称三维卷积层采用结构不对称的三维卷积核,归一化采用了batchnormalization,激励函数采用了ReLU,池化层采用了三维最大池化操作;数据经不对称卷积层处理之后,再利用轻量化模块构成的结构进一步进行特征提取;
2)分类部分,该部分由一个三维的自适应平均池化层和宽度为256的全连接层构成,三维自适应评价池化层会根据输入数据调整池化操作的核的尺寸及步长,能够将任何维度的输入数据处理成固定尺寸的输出数据;因此,在处理光谱维度不同的高光谱图像数据时,全连接层的宽度无需针对数据进行调整变化;
步骤4:训练网络模型:将训练数据批量地输入到构建好的三维轻量化深度网络中,以标注的类别为指导信息,利用梯度下降算法对网络参数进行训练,直至网络收敛;训练过程中,每次随机不重复的从训练集里抽取10-20个样本为一批训练数据,将该数据输入到网络,抽取特征并计算预测结果,以预测结果与实际结果之间的交叉熵为损失函数,计算网络权值的偏导数,并利用梯度下降算法,更新网络参数;训练过程遍历整个训练集一次为一轮训练;
步骤5:生成分类结果:基于训练过的模型,对高光谱图像中所有像素进行类别预测,然后将预测的类别放到对应的位置,得到分类结果图。
步骤2中的S取27。
步骤3中的多组轻量化模块为四组轻量化模块,依次包含1个宽度为32的轻量化模块2,2个宽度均为64的轻量化模块1和轻量化模块2,两个宽度均为128的轻量化模块1和轻量化模块2,一个宽度为256的轻量化模块1;所述的轻量化模块1的结构:从输入端到输出端,轻量化模块1右侧主干部分依次包含一个point-wise卷积层,三维的depth-wise卷积层,还有另一个point-wise卷积层;第一个point-wise卷积层及三维depth-wise卷积层后均接着一个batch normaliaztion层和ReLU激励层;第二个point-wise卷积层后边仅有一个batch normalization层,左侧分枝结构依次包含一个窗口宽度为3,步长为2的平均池化层和一个point-wise卷积层,在右侧主干和左侧分枝结构通过按元素逐元素相加合并后会经过一个ReLU激励层得到模块的输出;所述的轻量化模块2的结构:右侧主干部分与轻量化模块1结构类似,左侧分枝部分不包含任何参数,是一个恒等映射。
步骤4中整个训练过程进行60轮,前50轮学习率设为0.01,最后10轮,学习率衰减到0.001,整个训练过程中,动量项设为0.9。
有益效果
本发明提出的一种基于三维轻量化深度网络的高光谱图像分类方法,通过构建三位轻量化深度网络,实现了有限样本条件下,高光谱图像深度特征的自主提取,高精度的分类。本发明与现有的基于深度学习的高光谱图像分类方法相比,网络模型更深,精度更高,且参数量更少。
附图说明
图1基于三维轻量化深度网络的高光谱图像分类方法流程图
图2轻量化模块1,有特征尺度压缩功能的轻量化模块
图3轻量化模块2,特征提取轻量化模块
图4三维轻量化深度网络结构示意图
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明的技术方案是基于三维轻量化深度网络的高光谱图像分类方法。该方法从待处理的高光谱图像中,提取少量有标注的样本对本技术方案中所提出的三维轻量化深度网络进行训练,然后利用训练好的网络模型对整副图像进行分类。
本技术方案的具体措施如下:
步骤1:数据预处理。待处理的高光谱图像数据进行最大最小归一化。
步骤2:数据划分。统计待处理的高光谱图像中,各个类别标注样本的个数,然后按照比例,从每个类别中,抽取5%-10%的标注样本作为训练数据,剩余的部分作为测试数据。
步骤3:构建网络模型。本发明所构建的网络结构依次包含两部分,1)特征提取部分,由不对称的三维卷积层,可学习的归一化层(batch normalization),激励函数(ReLU)及三维的最大池化层(max pooling)及四组不同宽度的轻量化模块构成。2)三维自适应平均池化层及全连接层构成的分类部分。
步骤4:训练网络模型。将训练数据批量的输入到构建好的三维轻量化深度网络中,以标注的类别为指导信号,利用梯度下降算法对网络参数进行训练,直至网络收敛。在测试数据上评估网络精度。
步骤5:生成分类结果。基于训练过的模型,对高光谱图像中所有像素进行类别预测,然后将预测的类别放到对应的位置,得到分类结果图。
下面结合附图对本发明进一步说明:
基于三维轻量化网络的高光谱图像分类方法
步骤1:数据预处理。待处理的高光谱图像数据进行最大最小归一化,归一化公式如下:
步骤2:数据划分。统计待处理的高光谱图像中,各个类别标注样本的个数,然后按照比例,从每个类别中,抽取5%-10%的标注样本作为训练数据,剩余的部分作为测试数据。具体做法如下,对于一个尺寸为M×N×L的三维高光谱图像数据,M,N分别表示高光谱图像的高和宽,L表示数据的波段数。抽取样本时,以待处理像素为中心,抽取S×S×L的数据块作为中心像素的样本数据,S表示邻域大小,一般取27。
步骤3:构建网络模型。本发明设计的网络依次包含了两部分结构:1)特征提取部分,该部分包含了基于不对称三维卷积网络的数据处理模块及多组轻量化模块构成的特征提取模块。数据处理模块主要包含了一个不对称的三维卷积层,归一化层,激励函数及池化层。其中不对称三维卷积层采用结构不对称的三维卷积核,该卷积核在光谱维尺寸大于空间维尺寸,从而起到在该处理模块处理数据过程中更注重光谱维信息的作用,例如,卷积层可采用光谱维尺度为8,空间维尺度为3×3的卷积核,并将卷积层宽度设为32。该模块中,归一化采用了batch normalization,激励函数采用了ReLU,池化层采用了三维最大池化操作。数据经不对称卷积层处理之后,再利用轻量化模块构成的结构进一步进行特征提取。四组轻量化模块依次包含1个宽度为32的轻量化模块,2个宽度为64的轻量化模块,两个宽度为128的轻量化模块,一个宽度为256的轻量化模块。根据具体数据的不同,该结构也可以适当的进行增减。轻量化模块包含两种不同的结构,有特征尺度压缩功能的轻量化模块1。从输入端到输出端,轻量化模块1右侧主干部分依次包含一个point-wise卷积层,三维的depth-wise卷积层,还有另一个point-wise卷积层。第一个point-wise卷积层及三维depth-wise卷积层后均接着一个batch normaliaztion层和ReLU激励层。第二个point-wise卷积层后边仅有一个batch normalization层,在主干和左侧的分枝结构通过按元素逐元素相加合并后会经过一个ReLU激励层得到模块的输出。左侧分枝结构依次包含一个窗口宽度为3,步长为2的平均池化层和一个point-wise卷积层,具体结构如图2所示。轻量化模块2,右侧主干部分与轻量化模块1结构类似,左侧分枝部分不包含任何参数,是一个恒等映射,具体结构如图3所示。模型整体结构参照图4。
2)分类部分。该部分由一个三维的自适应平均池化层和宽度为256的全连接层构成。三维自适应评价池化层会根据输入数据调整池化操作的核的尺寸及步长,能够将任何纬度的输入数据处理成固定尺寸的输出数据。因此,在处理光谱纬度不同的高光谱图像数据时,全连接层的宽度无需针对数据进行调整变化。网络整体结构请参照图4。
步骤4:训练网络模型。将训练数据批量地输入到构建好的三维轻量化深度网络中,以标注的类别为指导信息,利用梯度下降算法对网络参数进行训练,直至网络收敛。训练过程中,每次随机不重复的从训练集里抽取10-20样本为一批训练数据,将该数据输入到网络,抽取特征并计算预测结果,以预测结果与实际结果之间的交叉熵为损失函数,计算网络权值的偏导数,并利用梯度下降算法,更新网络参数。训练过程便利整个训练集一次为一轮训练。整个训练过程进行60轮,前50轮学习率设为0.01,最后10轮,学习率衰减到0.001。整个训练过程中,动量项设为0.9。
步骤5:生成分类结果。基于训练过的模型,对高光谱图像中所有像素进行类别预测,然后将预测的类别放到对应的位置,得到分类结果图。
Claims (3)
1.一种基于三维轻量化深度网络的高光谱图像分类方法,其特征在于步骤如下:
步骤1:数据预处理:将待处理的高光谱图像数据进行最大最小归一化,归一化公式如下:
步骤2:数据划分:统计待处理的高光谱图像中,各个类别标注样本的个数,然后按照比例,从每个类别中,抽取5%-10%的标注样本作为训练数据,剩余的部分作为测试数据;具体做法如下,对于一个尺寸为M×N×L的三维高光谱图像数据,M,N分别表示高光谱图像的高和宽,L表示数据的波段数;抽取样本时,以待处理像素为中心,抽取S×S×L的数据块作为中心像素的样本数据,S表示邻域大小;
步骤3:构建三维轻量化深度网络模型:设计网络依次包含了两部分结构:1)特征提取部分,该部分包含了基于不对称三维卷积网络的数据处理模块及多组轻量化模块构成的特征提取模块,所述的数据处理模块包含一个不对称的三维卷积层,归一化层,激励函数及池化层,其中不对称三维卷积层采用结构不对称的三维卷积核,归一化采用了batchnormalization,激励函数采用了ReLU,池化层采用了三维最大池化操作;数据经不对称卷积层处理之后,再利用轻量化模块构成的结构进一步进行特征提取;
2)分类部分,该部分由一个三维的自适应平均池化层和宽度为256的全连接层构成,三维自适应评价池化层会根据输入数据调整池化操作的核的尺寸及步长,能够将任何维度的输入数据处理成固定尺寸的输出数据;因此,在处理光谱维度不同的高光谱图像数据时,全连接层的宽度无需针对数据进行调整变化;
步骤3中的多组轻量化模块为四组轻量化模块,依次包含1个宽度为32的轻量化模块2,2个宽度均为64的轻量化模块1和轻量化模块2,两个宽度均为128的轻量化模块1和轻量化模块2,一个宽度为256的轻量化模块1;所述的轻量化模块1的结构:从输入端到输出端,轻量化模块1右侧主干部分依次包含一个point-wise卷积层,三维的depth-wise卷积层,还有另一个point-wise卷积层;第一个point-wise卷积层及三维depth-wise卷积层后均接着一个batch normaliaztion层和ReLU激励层;第二个point-wise卷积层后边仅有一个batchnormalization层,左侧分枝结构依次包含一个窗口宽度为3,步长为2的平均池化层和一个point-wise卷积层,在右侧主干和左侧分枝结构通过按元素逐元素相加合并后会经过一个ReLU激励层得到模块的输出;所述的轻量化模块2的结构:右侧主干部分与轻量化模块1结构类似,左侧分枝部分不包含任何参数,是一个恒等映射;
步骤4:训练网络模型:将训练数据批量地输入到构建好的三维轻量化深度网络中,以标注的类别为指导信息,利用梯度下降算法对网络参数进行训练,直至网络收敛;训练过程中,每次随机不重复的从训练集里抽取10-20个样本为一批训练数据,将该数据输入到网络,抽取特征并计算预测结果,以预测结果与实际结果之间的交叉熵为损失函数,计算网络权值的偏导数,并利用梯度下降算法,更新网络参数;训练过程遍历整个训练集一次为一轮训练;
步骤5:生成分类结果:基于训练过的模型,对高光谱图像中所有像素进行类别预测,然后将预测的类别放到对应的位置,得到分类结果图。
2.据权利要求1所述的一种基于三维轻量化深度网络的高光谱图像分类方法,其特征在于步骤2中的S取27。
3.根据权利要求1所述的一种基于三维轻量化深度网络的高光谱图像分类方法,其特征在于步骤4中整个训练过程进行60轮,前50轮学习率设为0.01,最后10轮,学习率衰减到0.001,整个训练过程中,动量项设为0.9。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811542534.1A CN109753996B (zh) | 2018-12-17 | 2018-12-17 | 基于三维轻量化深度网络的高光谱图像分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811542534.1A CN109753996B (zh) | 2018-12-17 | 2018-12-17 | 基于三维轻量化深度网络的高光谱图像分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109753996A CN109753996A (zh) | 2019-05-14 |
CN109753996B true CN109753996B (zh) | 2022-05-10 |
Family
ID=66403790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811542534.1A Active CN109753996B (zh) | 2018-12-17 | 2018-12-17 | 基于三维轻量化深度网络的高光谱图像分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109753996B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110210439A (zh) * | 2019-06-10 | 2019-09-06 | 西北工业大学 | 基于轻量化三维卷积网络的行为识别方法 |
CN110222773B (zh) * | 2019-06-10 | 2023-03-24 | 西北工业大学 | 基于不对称分解卷积网络的高光谱图像小样本分类方法 |
CN110322414B (zh) * | 2019-07-05 | 2021-08-10 | 北京探境科技有限公司 | 一种基于ai处理器的图像数据在线量化矫正方法及系统 |
CN111091550A (zh) * | 2019-12-12 | 2020-05-01 | 创新奇智(北京)科技有限公司 | 多尺寸自适应的pcb锡膏区域检测系统及检测方法 |
CN115607170B (zh) * | 2022-11-18 | 2023-04-25 | 中国科学技术大学 | 一种基于单通道脑电信号的轻量级睡眠分期方法及应用 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105320965A (zh) * | 2015-10-23 | 2016-02-10 | 西北工业大学 | 基于深度卷积神经网络的空谱联合的高光谱图像分类方法 |
CN106022355A (zh) * | 2016-05-09 | 2016-10-12 | 西北工业大学 | 基于3dcnn的高光谱图像空谱联合分类方法 |
CN107844828A (zh) * | 2017-12-18 | 2018-03-27 | 北京地平线信息技术有限公司 | 神经网络中的卷积计算方法和电子设备 |
CN108345831A (zh) * | 2017-12-28 | 2018-07-31 | 新智数字科技有限公司 | 基于点云数据的道路图像分割的方法、装置以及电子设备 |
CN108388917A (zh) * | 2018-02-26 | 2018-08-10 | 东北大学 | 一种基于改进深度学习模型的高光谱图像分类方法 |
CN108596213A (zh) * | 2018-04-03 | 2018-09-28 | 中国地质大学(武汉) | 一种基于卷积神经网络的高光谱遥感图像分类方法及系统 |
CN108717732A (zh) * | 2018-05-21 | 2018-10-30 | 电子科技大学 | 一种基于MobileNets模型的表情追踪方法 |
CN108830796A (zh) * | 2018-06-20 | 2018-11-16 | 重庆大学 | 基于谱空结合和梯度域损失的高光谱图像超分辨重构方法 |
-
2018
- 2018-12-17 CN CN201811542534.1A patent/CN109753996B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105320965A (zh) * | 2015-10-23 | 2016-02-10 | 西北工业大学 | 基于深度卷积神经网络的空谱联合的高光谱图像分类方法 |
CN106022355A (zh) * | 2016-05-09 | 2016-10-12 | 西北工业大学 | 基于3dcnn的高光谱图像空谱联合分类方法 |
CN107844828A (zh) * | 2017-12-18 | 2018-03-27 | 北京地平线信息技术有限公司 | 神经网络中的卷积计算方法和电子设备 |
CN108345831A (zh) * | 2017-12-28 | 2018-07-31 | 新智数字科技有限公司 | 基于点云数据的道路图像分割的方法、装置以及电子设备 |
CN108388917A (zh) * | 2018-02-26 | 2018-08-10 | 东北大学 | 一种基于改进深度学习模型的高光谱图像分类方法 |
CN108596213A (zh) * | 2018-04-03 | 2018-09-28 | 中国地质大学(武汉) | 一种基于卷积神经网络的高光谱遥感图像分类方法及系统 |
CN108717732A (zh) * | 2018-05-21 | 2018-10-30 | 电子科技大学 | 一种基于MobileNets模型的表情追踪方法 |
CN108830796A (zh) * | 2018-06-20 | 2018-11-16 | 重庆大学 | 基于谱空结合和梯度域损失的高光谱图像超分辨重构方法 |
Non-Patent Citations (4)
Title |
---|
A Survey of Deep Learning: Platforms, Applications and Emerging Research Trends;William Grant Hatcher等;《IEEE Access 》;20180427;第24411-24432页 * |
S-NET: A LIGHTWEIGHT CONVOLUTIONAL NEURAL NETWORK FOR N-DIMENSIONAL SIGNALS;Yingxuan Cui等;《2018 IEEE International Conference on Multimedia & Expo Workshops (ICMEW)》;20181129;第1-4页 * |
三维卷积神经网络模型联合条件随机场优化的高光谱遥感影像分类;李竺强等;《光学学报》;20180831;第38卷(第8期);第0828001-1-0828001-10页 * |
基于OpenCV的精量喷雾图像处理技术;陈钊等;《农机化研究》;20180630(第6期);第194-223页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109753996A (zh) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109754017B (zh) | 基于可分离的三维残差网络和迁移学习高光谱图像分类方法 | |
CN109753996B (zh) | 基于三维轻量化深度网络的高光谱图像分类方法 | |
CN109584248B (zh) | 基于特征融合和稠密连接网络的红外面目标实例分割方法 | |
CN110378381B (zh) | 物体检测方法、装置和计算机存储介质 | |
CN108615010B (zh) | 基于平行卷积神经网络特征图融合的人脸表情识别方法 | |
CN106683048B (zh) | 一种图像超分辨率方法及设备 | |
CN111985543B (zh) | 一种高光谱图像分类模型的构建方法、分类方法及系统 | |
CN105320965B (zh) | 基于深度卷积神经网络的空谱联合的高光谱图像分类方法 | |
CN110163813B (zh) | 一种图像去雨方法、装置、可读存储介质及终端设备 | |
CN110222773B (zh) | 基于不对称分解卷积网络的高光谱图像小样本分类方法 | |
CN107563433B (zh) | 一种基于卷积神经网络的红外小目标检测方法 | |
CN112446270A (zh) | 行人再识别网络的训练方法、行人再识别方法和装置 | |
CN112308152B (zh) | 基于光谱分割与同质区域检测的高光谱图像地物分类方法 | |
CN107832797B (zh) | 基于深度融合残差网的多光谱图像分类方法 | |
CN107808138B (zh) | 一种基于FasterR-CNN的通信信号识别方法 | |
CN107944483B (zh) | 基于双通道dcgan和特征融合的多光谱图像分类方法 | |
CN110619352A (zh) | 基于深度卷积神经网络的典型红外目标分类方法 | |
CN110136162B (zh) | 无人机视角遥感目标跟踪方法及装置 | |
CN110852369B (zh) | 联合3d/2d卷积网络和自适应光谱解混的高光谱图像分类方法 | |
CN104866868A (zh) | 基于深度神经网络的金属币识别方法和装置 | |
CN111695407A (zh) | 基于多光谱融合的性别识别方法、系统、存储介质及终端 | |
Peng et al. | Litchi detection in the field using an improved YOLOv3 model | |
CN116188943A (zh) | 太阳射电频谱爆发信息检测方法及装置 | |
CN116563683A (zh) | 一种基于卷积神经网络和多层感知机的遥感影像场景分类方法 | |
CN110866552B (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 |