CN107229936B - 基于球状鲁棒序列局部二值化模式的序列分类方法 - Google Patents
基于球状鲁棒序列局部二值化模式的序列分类方法 Download PDFInfo
- Publication number
- CN107229936B CN107229936B CN201710365058.XA CN201710365058A CN107229936B CN 107229936 B CN107229936 B CN 107229936B CN 201710365058 A CN201710365058 A CN 201710365058A CN 107229936 B CN107229936 B CN 107229936B
- Authority
- CN
- China
- Prior art keywords
- sequence
- binary
- frame
- value
- spherical
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/50—Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
- G06V10/507—Summing image-intensity values; Histogram projection analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了基于球状鲁棒序列局部二值化模式的序列分类方法方法,主要解决现有技术对极光序列分类准确率不高和时间复杂度不够低的问题。其技术方案是:1.选取训练数据序列和测试数据序列;2.求出训练数据序列的球状鲁棒序列局部二值化编码值,统计得到训练数据的特征向量;3.将训练数据序列的特征向量和标签输入到SVM中进行训练得到训练模型;4.求出测试数据序列的球状鲁棒序列局部二值化编码值,统计得到测试数据的特征向量;5.将测试数据序列的特征向量输入到训练好的SVM模型中得到分类结果。本发明提高了极光序列的分类准确率,降低了时间复杂度,可用于计算机对极光序列的自动分类。
Description
技术领域
本发明属于图像处理技术领域,特别涉及一种极光序列分类方法,可用于极光序列特征提取、极光序列时间分类及行星际磁场参数研究。
背景技术
极光是来自磁层的高能粒子沉降到高层大气并与中性成分碰撞激发的一种大气发光现象,是太阳风-磁层-电离层耦合造成的地球高纬活动的重要表现形式之一。人们通过极光形态及其演化的系统观测可以获得磁层和日地空间电磁活动的大量信息,有助于深入研究太阳活动对地球的影响方式与程度,对了解空间天气过程的变化规律具有重要意义。
随着不断深入地研究极光形态特征以及极光光谱学的出现,人们对不同极光的物理特征及其产生机制的研究也越来越关注。早期的极光分类只是单一地依据观察到的极光的各种视觉形式。Simmons于1998年在文章“A classification of auroral types,Journal of the British Astronomical Association,108(5),247-257.”中提出了目前不同类型的极光的分类是根据他们的地球物理的特征以及对每一个类型的大多数典型特征的详细描述是特别的参照了地磁的纬度和地磁的时间扇区和极光产生的机制。根据他提供的分类标准可以帮助研究者从地球物理的方面阐明这些极光的各种类型之间的相似和不同,对极光的不同类型都有了清楚的定义。Hu等人在文章“Synoptic distribution ofdayside aurora:Multiple-wavelength all-sky observation at Yellow RiverStation inSvalbard.Journal of Atmospheric and Solar-TerrestrialPhysics,2009,71(8):794-804.”中根据活动区域将日侧极光分为四类。Wen等人在文章“午后多重极光弧观测研究.极地研究,2005.,17(4):256-263.”中利用南极中山站的全天空摄像机观测到的极光数据对午后多重极光弧的出现率及其对应的KP指数发生规律进行了统计分析,历时175天,找出它们之间的关系。
过去的很多年中,研究者将机器学习、计算机视觉与极光数据的分类、显著性提取等相结合,并且针对极光数据的特殊性采取了不同的特征提取方式,并在此过程中发现局部二值化模式和序列局部二值化模式对极光数据具有极好的表征性。但是已有的局部二值化模式由于没有考虑时间维度和空间维度上的像素点对当前帧贡献率的不同,造成分类准确率不够高和时间复杂度不够低的不足。
发明内容
本发明的目的在于提出一种基于球状鲁棒序列局部二值化模式的序列分类方法,以提高极光序列的分类准确率,降低时间复杂度。
本发明的技术方案是在原有的序列局部二值化模式VLBP的基础上分别对权值分配和像素点选取方式进行改进,以此获得对极光序列具有更强表征性的特征,进而提高极光序列的分类效果,其实现步骤包括如下:
(1)根据极光形态特征手动从全天空极光数据中选取三种极光数据,即弧状极光数据,辐射状极光数据,帷幔状极光数据,每种数据中的序列长度为15帧到50帧之间,从每类数据中分别随机选取N段序列作为训练数据,其余的作为测试数据,N>=100;
(2)计算训练数据序列的球状鲁棒序列局部二值化编码值V:
2a)选取当前帧和前后帧的邻域像素点,并计算前后帧邻域像素点的均值g′t±1,c;
2b)用当前帧的邻域像素点计算二值化阈值TDc;
2c)用二值化阈值TDc计算当前帧的二值化编码Vt,并计算其旋转不变性二值化编码Vt rbt;
2d)用二值化阈值TDc对前后帧的像素点均值g′t±1,c计算其二值化编码Vt±1;
2e)结合当前帧的旋转不变性二值化编码Vt rbt和前后帧的二值化编码Vt±1得到球状鲁棒序列局部二值化编码V;
2f)对训练数据球状鲁棒序列局部二值化编码值V进行直方图统计,得到训练数据序列的特征向量fi j,其中i是训练数据中的序列标号,j是标签;
(3)将训练数据序列的标签j和球状鲁棒序列局部二值化编码值对应的特征向量fi j输入到SVM中进行训练得到训练模型;
(4)求出测试数据序列的球状鲁棒序列局部二值化编码值H,并进行直方图统计得到测试数据序列的特征向量fm,其中m是测试数据中的序列标号;
(5)将测试数据序列的球状鲁棒序列局部二值化编码值对应的特征向量fm输入到训练好的SVM模型中,得到测试数据序列的分类结果。
本发明由于在原有的序列局部二值化模式VLBP的基础上分别对权值分配和像素点选取方式进行改进,即分别对空间维度上相邻像素点对中心像素点的距离H和时间维度上相邻像素点对中心像素点的距离R进行量化,同时对时间维度和空间维度上选取的邻域像素点对中心像素点的表征贡献率进行调整,使其成为更合理的权重优化分配方式,以此获得了对极光序列具有更强表征性的特征,进而提高了极光序列的分类效果。
附图说明
图1是本发明的实现流程图;
图2是本发明中的序列特征提取示意图;
图4是用本发明对取不同邻域像素点个数F和维度距离比值H/R的分类准确率图。
具体实施方式
以下结合附图对本发明的内容和效果进行进一步描述。
参照图1,本实例的具体实现步骤如下:
步骤1:选取训练数据和测试数据。
根据极光形态特征手动从全天空极光数据中选取三种极光数据,即弧状极光数据,辐射状极光数据,帷幔状极光数据,每种数据中的序列长度为15帧到50帧之间;从每类数据中分别随机选取N段序列作为训练数据,其余的作为测试数据,N>=100。
步骤2:计算训练数据序列的编码值V,获得训练数据序列的特征向量fi j。
现有的计算序列的编码值方法有序列局部二值化模式VLBP,本发明采用球状鲁棒序列局部二值化模式。
参照图2,本实现的具体实现如下:
2a)选取当前帧和前后帧的邻域像素点,并计算前后帧邻域像素点的均值gt’±1,c;
2a1)以训练数据序列中第t帧Qt中的某个像素点Wc为中心像素点,设该中心像素点的像素值为gt,c,坐标为(xc,yc,t);在第t帧Qt上围绕中心像素点Wc以R为半径的圆上选取F个像素点作为邻域像素点gt,f,坐标为:
((xc+Rcos(2πf/F)),(yc+Rsin(2πf/F)),t),f=1,2,...,F;
2a2)设第t帧与其之前的第t-1帧和与其之后的第t+1帧之间在时间维度上的距离为H,根据H与空间维度上距离R的大小不同,选取不同个数的邻域像素点gt±1,c和gt±1,f’得到这两种邻域像素点的坐标,即:
gt±1,c的坐标为:(xc,yc,t±1),
2a3)分别求第t-1帧和第t+1帧上邻域像素点的均值,作为第t-1帧和第t+1帧的像素值g′t±1,c;
2b)根据第t帧中心像素点gt,c和邻域像素点gt,f求二值化阈值TDc:
2c)用二值化阈值TDc计算当前帧的二值化编码Vt,并计算其旋转不变性二值化编码Vt rbt;
2c1)用二值化阈值TDc对第t帧F个邻域像素值gt,f进行二值化编码,得到第t帧的二值化编码值Vt:
Vt=Vt(s(gt,1-TDc),...,s(gt,f-TDc),...,s(gt,F-TDc)),f=1,2,...,F
2c2)求第t帧的二值化编码值Vt对应的具有旋转不变性的二进制编码值Vt rbt:
Vt rbt=min{ROT(Vt,k)},k=0,1,...,(F-1)
其中,ROT(Vt,k)表示将第t帧的二进制编码值Vt沿顺时针方向旋转了k次后所得到的二进制编码值,min{ROT(Vt,k)}表示第t帧的二进制编码值Vt沿顺时针方向旋转了k,k=0,1,...,(F-1)次后所得到F个二进制编码值的最小值;
2d)用二值化阈值TDc对第t-1帧和第t+1帧的像素值g′t±1,c进行二值化编码,得到前后帧二值化编码值Vt±1:
Vt±1=Vt±1(s(g′t-1,c-TDc),s(g′t+1,c-TDc)),
2e)结合第t帧旋转不变性的二进制编码值Vt rbt和其前后帧二值化编码值Vt±1得到最终的球状鲁棒序列局部二值化编码值V:
V=V(Vt rbt,Vt±1)。
步骤3:训练SVM分类器。
3a)将步骤2得到的每个训练数据序列的特征向量fi j作为列向量,构成训练数据序列的特征矩阵Tr,并用fi j相应的标签j构成标签列向量r;
3b)将训练数据序列的特征矩阵Tr和标签列向量r输入到SVM分类器,对SVM分类器的参数进行拟合,得到训练好的SVM分类器。
步骤4:计算测试数据序列的编码值H,得到特征向量fm。
计算测试数据序列编码值H的过程与步骤2相同,最后得到测试数据序列的特征向量fm,其中m是测试数据中的序列标号。
步骤5:对测试数据序列分类。
5a)将步骤4得到的每个测试数据序列的特征向量fm作为列向量构成一个特征矩阵Tr;
5b)将得到的特征矩阵Tr输入到训练好的SVM分类器,得到测试数据序列的分类结果。
本发明的效果通过以下仿真实验进一步说明:
1.仿真条件与方法:
硬件平台为:Intel Core i5、2.93GHz、3.45GB RAM;
软件平台为:Windows7操作系统下的MATLAB R2012b;
2.仿真内容及结果:
实验仿真2,取不同H/R和不同F,其中H/R是时间维度上的距离H和空间维度上距离R的比值,F是邻域像素点数,用本发明对极光序列进行分类,分类准确率如图4,其中,图4(a)是取F=4,不同的H/R对应的分类准确率,图(b)是取F=4,8,16时,分类准确率最高对应H/R值,图中H的下标是邻域像素点数F的值。
从图4可以看出,邻域像素点数F越多,时间维度和空间维度的距离比值H/R越小,均值像素点的值融合邻域像素点的信息越多,特征表达性越强,分类准确率越高。
实验仿真3,分别用本发明方法与现有VLBP方法对极光序列进行分类,分类准确率和时间归一化对比结果如表1。
表1本发明与VLBP对比结果
从表1中的实验结果中可以看出,本发明的分类准确率比VLBP有所提高,且在时间复杂度上有很大改善。
实验仿真4,用本发明方法分别对训练集序列个数N为30,60,90,120,150时进行分类,其准确率如表2所示。
表2极光序列分类准确率
从表2可以看出,当训练样本数超过150时,分类准确率为88.5%,本发明分类准确率较高,能够实现计算机对极光序列的自动分类。
Claims (7)
1.基于球状鲁棒序列局部二值化模式的序列分类方法,包括:
(1)根据极光形态特征手动从全天空极光数据中选取三种极光数据,即弧状极光数据,辐射状极光数据,帷幔状极光数据,每种数据中的序列长度为15帧到50帧之间,从每类数据中分别随机选取N段序列作为训练数据,其余的作为测试数据,N>=100;
(2)计算训练数据序列的球状鲁棒序列局部二值化编码值V:
2a)选取当前帧和前后帧的邻域像素点,并分别计算前后帧邻域像素点的均值g’t-1,c和g’t+1,c;按如下步骤进行:
2a1)以训练数据序列中第t帧Qt中的某个像素点Wc为中心像素点,设该中心像素点的像素值为gt,c,坐标为(xc,yc,t);在第t帧Qt上围绕中心像素点Wc以R为半径的圆上选取F个像素点作为邻域像素点gt,f,坐标为:
((xc+R cos(2πf/F)),(yc+R sin(2πf/F)),t),f=1,2,...,F;
2a2)设第t帧与其之前的第t-1帧和与其之后的第t+1帧之间在时间维度上的距离为H,根据H与空间维度上距离R的大小不同,选取不同个数的邻域像素点gt±1,c和gt±1,f,得到这两种邻域像素点的坐标,即:
gt±1,c的坐标为:(xc,yc,t±1),
2a3)分别求第t-1帧和第t+1帧上邻域像素点的均值,作为第t-1帧和第t+1帧的像素值g’t±1,c;
2b)用当前帧的邻域像素点计算二值化阈值TDc;
2c)用二值化阈值TDc计算当前帧的二值化编码Vt,并计算其旋转不变性二值化编码Vt rbt;
2d)用二值化阈值TDc对前后帧的像素点均值g’t±1,c计算其二值化编码Vt±1;
2e)结合当前帧的旋转不变性二值化编码Vt rbt和前后帧的二值化编码Vt±1得到球状鲁棒序列局部二值化编码V;
2f)对训练数据球状鲁棒序列局部二值化编码值V进行直方图统计,得到训练数据序列的特征向量fi j,其中i是训练数据中的序列标号,j是标签;
(3)将训练数据序列的标签j和球状鲁棒序列局部二值化编码值对应的特征向量fi j输入到SVM中进行训练得到训练模型;
(4)求出测试数据序列的球状鲁棒序列局部二值化编码值E,并进行直方图统计得到测试数据序列的特征向量fm,其中m是测试数据中的序列标号;
(5)将测试数据序列的球状鲁棒序列局部二值化编码值对应的特征向量fm输入到训练好的SVM模型中,得到测试数据序列的分类结果。
3.根据权利要求1所述的方法,其中步骤2c)中用二值化阈值TDc计算当前帧的二值化编码Vt,并计算其旋转不变性二值化编码Vt rbt,按如下步骤进行:
2c1)用二值化阈值TDc对第t帧F个邻域像素值gt,f进行二值化编码,得到第t帧的二值化编码值Vt:
Vt=(s(gt,1-TDc),...,s(gt,f-TDc),...,s(gt,F-TDc)),f=1,2,...,F
2c2)求第t帧的二值化编码值Vt对应的具有旋转不变性的二进制编码值Vt rbt:
Vt rbt=min{ROT(Vt,k)},k=0,1,...,(F-1)
其中,ROT(Vt,k)表示将第t帧的二进制编码值Vt沿顺时针方向旋转了k次后所得到的二进制编码值,min{ROT(Vt,k)}表示第t帧的二进制编码值Vt沿顺时针方向旋转了k,k=0,1,...,(F-1)次后所得到F个二进制编码值的最小值。
5.根据权利要求1所述的方法,其中步骤2e)中结合当前帧的旋转不变性二值化编码Vt rbt和前后帧的二值化编码Vt±1得到球状鲁棒序列局部二值化编码V,通过下式计算:
V=(Vt rbt,Vt±1)。
7.根据权利要求1所述的方法,其中步骤(3)中进行训练得到训练模型,按如下步骤进行:
4a)将步骤(2)得到的每个训练数据序列的特征向量fi j作为列向量构成一个特征矩阵Tr,相应的标签j构成一个标签列向量r;
4b)将训练数据序列的特征矩阵Tr和标签列向量r输入到SVM分类器,对SVM分类器的中的参数进行拟合,得到训练好的SVM分类器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710365058.XA CN107229936B (zh) | 2017-05-22 | 2017-05-22 | 基于球状鲁棒序列局部二值化模式的序列分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710365058.XA CN107229936B (zh) | 2017-05-22 | 2017-05-22 | 基于球状鲁棒序列局部二值化模式的序列分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107229936A CN107229936A (zh) | 2017-10-03 |
CN107229936B true CN107229936B (zh) | 2020-06-23 |
Family
ID=59933370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710365058.XA Active CN107229936B (zh) | 2017-05-22 | 2017-05-22 | 基于球状鲁棒序列局部二值化模式的序列分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107229936B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109033944B (zh) * | 2018-06-07 | 2021-09-24 | 西安电子科技大学 | 一种全天空极光图像分类与关键局部结构定位方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971120A (zh) * | 2014-04-18 | 2014-08-06 | 西安电子科技大学 | 基于空时极向局部二值模式的极光图像序列分类方法 |
CN106649611A (zh) * | 2016-11-29 | 2017-05-10 | 哈尔滨商业大学 | 一种基于邻域旋转直角模式的图像检索方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100866792B1 (ko) * | 2007-01-10 | 2008-11-04 | 삼성전자주식회사 | 확장 국부 이진 패턴을 이용한 얼굴 기술자 생성 방법 및장치와 이를 이용한 얼굴 인식 방법 및 장치 |
US9336454B2 (en) * | 2013-05-31 | 2016-05-10 | Texas Instruments Incorporated | Vector processor calculation of local binary patterns |
-
2017
- 2017-05-22 CN CN201710365058.XA patent/CN107229936B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971120A (zh) * | 2014-04-18 | 2014-08-06 | 西安电子科技大学 | 基于空时极向局部二值模式的极光图像序列分类方法 |
CN106649611A (zh) * | 2016-11-29 | 2017-05-10 | 哈尔滨商业大学 | 一种基于邻域旋转直角模式的图像检索方法 |
Non-Patent Citations (7)
Title |
---|
Dynamic texture recognition based on completed volume local binary pattern;Tiwari, D 等;《MULTIDIMENSIONAL SYSTEMS AND SIGNAL PROCESSING》;20160430;全文 * |
Dynamic Texture Recognition Using Local Binary Patterns with an Application to Facial Expressions;Guoying Zhao 等;《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》;20070630;第29卷(第6期);论文第3节 * |
Recognizing human actions using novel space-time volume binary patterns;Baumann, F 等;《NEUROCOMPUTING》;20150115;全文 * |
人体背影序列运动特征的性别识别方法;唐浩 等;《华中科技大学学报》;20131031;全文 * |
基于空时极向 LBP 的极光序列事件检测;韩冰 等;《软件学报》;20140930;第25卷(第9期);论文第1-2节 * |
复合时空特征的双模态情感识别;王晓华 等;《中国图象图形学报》;20170131;全文 * |
局部二进制模式方法综述;刘丽 等;《中国图象图形学报》;20141231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107229936A (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110728200B (zh) | 一种基于深度学习的实时行人检测方法及系统 | |
CN110689086A (zh) | 基于生成式对抗网络的半监督高分遥感图像场景分类方法 | |
CN109325502B (zh) | 基于视频渐进区域提取的共享单车停放检测方法和系统 | |
CN104134068B (zh) | 基于稀疏编码的监控车辆特征表示及分类方法 | |
CN111967313B (zh) | 一种深度学习目标检测算法辅助的无人机图像标注方法 | |
CN108154158B (zh) | 一种面向增强现实应用的建筑物图像分割方法 | |
CN109285162A (zh) | 一种基于局部区域条件随机场模型的图像语义分割方法 | |
Liu et al. | The method of insulator recognition based on deep learning | |
CN110706235B (zh) | 一种基于双级级联分割的远红外行人检测方法 | |
CN107301382B (zh) | 基于时间依赖约束下深度非负矩阵分解的行为识别方法 | |
CN109919011A (zh) | 一种基于多时长信息的动作视频识别方法 | |
CN104239896A (zh) | 一种视频图像中人群密度等级分类方法 | |
CN111967511A (zh) | 一种基于异构特征融合网络的地基云图分类方法 | |
CN104850859A (zh) | 一种基于多尺度分析的图像特征包构建方法 | |
CN111507413A (zh) | 一种基于字典学习的城市管理案件图像识别方法 | |
CN109190666B (zh) | 基于改进的深度神经网络的花卉图像分类方法 | |
CN112507800A (zh) | 一种基于通道注意力机制和轻型卷积神经网络的行人多属性协同识别方法 | |
CN106570885A (zh) | 基于亮度和纹理融合阈值的背景建模方法 | |
CN103984963A (zh) | 一种高分辨率遥感图像场景分类方法 | |
CN112001293A (zh) | 结合多尺度信息和编解码网络的遥感影像地物分类方法 | |
CN107909073A (zh) | 多维局部二值模式和机器学习手写乐谱谱线删除方法 | |
CN113610024B (zh) | 一种多策略的深度学习遥感影像小目标检测方法 | |
CN107229936B (zh) | 基于球状鲁棒序列局部二值化模式的序列分类方法 | |
Alvarez et al. | Large-scale semantic co-labeling of image sets | |
CN110334584A (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 |