CN109190434A - 一种基于亚像素级角点检测的条码识别算法 - Google Patents

一种基于亚像素级角点检测的条码识别算法 Download PDF

Info

Publication number
CN109190434A
CN109190434A CN201810663401.3A CN201810663401A CN109190434A CN 109190434 A CN109190434 A CN 109190434A CN 201810663401 A CN201810663401 A CN 201810663401A CN 109190434 A CN109190434 A CN 109190434A
Authority
CN
China
Prior art keywords
bar code
image
piecemeal
point
blocksize
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
Application number
CN201810663401.3A
Other languages
English (en)
Other versions
CN109190434B (zh
Inventor
刘怡俊
陈龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201810663401.3A priority Critical patent/CN109190434B/zh
Publication of CN109190434A publication Critical patent/CN109190434A/zh
Application granted granted Critical
Publication of CN109190434B publication Critical patent/CN109190434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06KGRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K7/00Methods or arrangements for sensing record carriers, e.g. for reading patterns
    • G06K7/10Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation
    • G06K7/14Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation using light without selection of wavelength, e.g. sensing reflected white light
    • G06K7/1404Methods for optical code recognition
    • G06K7/146Methods for optical code recognition the method including quality enhancement steps
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20164Salient point detection; Corner detection

Abstract

本发明公开了一种基于亚像素级角点检测的条码识别算法,包括下述步骤:S1,载入原始条码图像;S2,原始条码图像前期处理;S3,判断并清除非条码区域图像,具体步骤如下:首先进行分块扫描模式,所选取分块数量为6,所选取分块尺寸blockSize=(图像宽和高的最小值)/6,即blockSize分块尺寸是根据原始条码图像行和列的最小值除于分块数量来计算的,而图像的每行分块数量block_cols和每列分块数量block_rows;本发明有效提高了图像处理速度与精确度,对非条码区域的清除更加有效,图像的处理与分析结果也越接近真实,最终也达到比较满意的识别效果。

Description

一种基于亚像素级角点检测的条码识别算法
技术领域
本发明涉及条形码检测技术领域,具体涉及一种基于亚像素级角点检测的条码识别算法。
背景技术
条形码是广泛使用的一种自动识别技术,它有以下几个方面的优点:输入速度快、可靠性高、采集信息量大、灵活实用等,条码的主要功能是携带数据并实现高效的数据采集,它在数据链中发挥着至关重要的作用,在某些行业中,如果条码精度不够高,就意味着数据链被破坏,供应链也将无法进行下去;在供应链中的任何环节,若因为条码精确度不够高而导致信息读写不准确,也将会给企业带来很多烦恼,因此需要将条码精度尽可能提高。
角点是图像中亮度变化最强的地方,反映了图像的本质特征,提取图像中的角点可以有效提高图像处理速度与精准度,因此对于处理整张图像来说特别重要,角点检测与提取地越准确,图像处理与分析结果就越接近真实;同时角点检测对真实环境下的对象识别、对象匹配都起到了决定性作用。条形码传统的角点检测方法有以下几种:基于模板的角点检测、基于边缘的角点检测、基于灰度变化的角点检测等,它们的检测精度最高只能到达像素级精度。但是,随着科学技术的飞速发展,工业检测等应用对条形码精确度的要求不断提高,传统的像素级角点检测方法已经不能满足实际测量的需要。
现有技术的缺点与不足:现有条码角点检测技术仅仅只能达到像素级精度,而随着科学技术的飞速发展,工业检测等应用对条形码精确度的要求不断提高,现有的像素级角点检测方法已经不能满足实际测量的需要。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于亚像素级角点检测的条码识别算法,该方法有效提高了图像处理速度与精确度,对非条码区域的清除更加有效,图像的处理与分析结果也越接近真实,最终也达到比较满意的识别效果。
本发明的目的通过下述技术方案实现:
一种基于亚像素级角点检测的条码识别算法,包括下述步骤:
S1,载入原始条码图像;
S2,原始条码图像前期处理;
S3,判断并清除非条码区域图像,具体步骤如下:
首先进行分块扫描模式,所选取分块数量为6,所选取分块尺寸blockSize=(图像宽和高的最小值)/6,即blockSize分块尺寸是根据原始条码图像行和列的最小值除于分块数量来计算的,而图像的每行分块数量block_cols和每列分块数量block_rows,由下述公式(1)和(2)所示:
block_cols=blurmage.cols/blockSize, (1)
block_rows=blurmage.rows/blockSize, (2)
其中,blockSize为分块尺寸大小,blurmage.cols为图像列的大小,blurmage.rows为图像行的大小;
然后对所得图像进行灰度梯度扫描操作:在每个分块的两条对角线、中竖线和中横线上作扫描,即在这四条线上检测连续4个像素之间的灰度差,当灰度差的绝对值大于阈值时,标记为一次灰度变化,如下述公式(3)所示:
其中,Vi为第i个像素的梯度量化值,|fi+1+fi+fi-1+fi-2|为连续四个像素之间的灰度,t为阈值,若分块在某条扫描线上发生超过3次的灰度变化,则认为该分块可能含有条码部分,需要保留,否则清除分块;最后对所得图像进行8连通区域检查,每个分块都有8个领域分块,而分块的边缘除外,在统计每个条码区域分块的邻域分块中是条码区域分块的个数时,若达不到3个,则将该条码区域分块视为孤立分块,孤立分块是条码的可能性较小,因此将其标记为非条码区域,并清除该分块;
S4,对图像进行分块亚像素级角点检测,对每个分块的角点数进行统计,若分块的角点数小于等于预设阈值,则视为条码区域,具体步骤如下:
亚像素级角点检测的实现是基于对向量正交性的观测而实现的,即从中央点q到其领域点p的向量和p点处的图像梯度正交,如下述公式(4)所示:
其中表示在q的一个领域点pi处的图像梯度,则q的值通过最小化εi得到;设置εi为0后,公式两边同乘以即可得到关于该领域点的一个等式;当对所有领域点均进行上述变换,等式左右分别相加,即可得到下述公式(5):
设定最终可求解q=G-1·b;
通过上述公式(5)确定亚像素级角点,所设置的阈值为3,若分块的角点数小于等于预设阈值,则视为条码区域,否则即为非条码区域;
S5,条码区域分块再判断,具体步骤如下:
条码是将宽度不等的多个黑条和空白,按照一定的编码规则排列,用以表达一组信息的图形标识符,根据条码的条形排列的特点,通过霍夫变换算法进行直线检测:
运用直线的极坐标公式:ρ=хcos(θ)+уsin(θ),其中ρ,θ是一对参数,通过对图像中的像素点运用该公式处理即可得到二维的ρ,θ参数空间,在参数空间上形成峰值,从而利用该峰值即可检测出直线;最后将直线检测和角点检测的判定结果进行或运算,作为分块的最终结果,以防止出现误判;
S6,根据剩余图像轮廓定位条码位置,旋转并截取保存条码区域,具体步骤如下:
查找并绘制条码区域的所有轮廓,用矩形勾出最大的轮廓视为条码位置,计算条码黑条倾斜角,以矩形中心进行原图像旋转,截取并保存原图中所在区域;
S7,识别条码。
本发明与现有技术相比具有以下的有益效果:
本发明有效提高了图像处理速度与精确度,对非条码区域的清除更加有效,图像的处理与分析结果也越接近真实,最终也达到比较满意的识别效果。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
本发明的主要工作原理为:对摄像机采集到的图像进行前期处理、读取原图、灰度化及中值滤波,在判断并清除非条码区域中进行角点检测采取基于亚像素级角点检测方法,根据剩余图像轮廓定位条码位置,旋转并截取保存条码区域等一系列操作从而得到精确的条形码图像。
具体来说,如图1所示,一种基于亚像素级角点检测的条码识别算法,包括下述步骤:
步骤一,载入原始条码图像;
步骤二,原始条码图像前期处理;
步骤三,判断并清除非条码区域图像,具体步骤如下:
首先进行分块扫描模式,分块尺寸太大会影响定位精度,分块尺寸太小则增加计算量,选取分块尺寸不合适则易导致条码定位失败,所选取分块数量为6,所选取分块尺寸blockSize=(图像宽和高的最小值)/6,即blockSize分块尺寸是根据原始条码图像行和列的最小值除于分块数量来计算的,而图像的每行分块数量block_cols和每列分块数量block_rows,由下述公式(1)和(2)所示:
block_cols=blurmage.cols/blockSize, (1)
block_rows=blurmage.rows/blockSize, (2)
其中,blockSize为分块尺寸大小,blurmage.cols为图像列的大小,blurmage.rows为图像行的大小;
然后对所得图像进行灰度梯度扫描操作:在每个分块的两条对角线、中竖线和中横线上作扫描,即在这四条线上检测连续4个像素之间的灰度差,当灰度差的绝对值大于阈值时,标记为一次灰度变化,如下述公式(3)所示:
其中,Vi为第i个像素的梯度量化值,|fi+1+fi+fi-1+fi-2|为连续四个像素之间的灰度,t为阈值,若分块在某条扫描线上发生超过3次的灰度变化,则认为该分块可能含有条码部分,需要保留,否则清除分块;最后对所得图像进行8连通区域检查,每个分块都有8个领域分块,而分块的边缘除外,在统计每个条码区域分块的邻域分块中是条码区域分块的个数时,若达不到3个,则将该条码区域分块视为孤立分块,孤立分块是条码的可能性较小,因此将其标记为非条码区域,并清除该分块;
步骤四,对图像进行分块亚像素级角点检测,对每个分块的角点数进行统计,若分块的角点数小于等于预设阈值,则视为条码区域,具体步骤如下:
亚像素级角点检测的实现是基于对向量正交性的观测而实现的,即从中央点q到其领域点p的向量和p点处的图像梯度正交,如下述公式(4)所示:
其中表示在q的一个领域点pi处的图像梯度,则q的值通过最小化εi得到;设置εi为0后,公式两边同乘以即可得到关于该领域点的一个等式;当对所有领域点均进行上述变换,等式左右分别相加,即可得到下述公式(5):
设定G为矩阵,b为向量,最终可求解q=G-1·b;
通过上述公式(5)确定亚像素级角点,所设置的阈值为3,若分块的角点数小于等于预设阈值,则视为条码区域,否则即为非条码区域;
步骤五,条码区域分块再判断,具体步骤如下:
条码是将宽度不等的多个黑条和空白,按照一定的编码规则排列,用以表达一组信息的图形标识符,根据条码的条形排列的特点,通过霍夫变换算法进行直线检测:
运用直线的极坐标公式:ρ=хcos(θ)+уsin(θ),其中ρ,θ是一对参数,通过对图像中的像素点运用该公式处理即可得到二维的ρ,θ参数空间,在参数空间上形成峰值,从而利用该峰值即可检测出直线;最后将直线检测和角点检测的判定结果进行或运算,作为分块的最终结果,以防止出现误判;
步骤六,根据剩余图像轮廓定位条码位置,旋转并截取保存条码区域,具体步骤如下:
查找并绘制条码区域的所有轮廓,用矩形勾出最大的轮廓视为条码位置,计算条码黑条倾斜角,以矩形中心进行原图像旋转,截取并保存原图中所在区域;
步骤七,识别条码。
其中所使用的边缘检测算法为canny算子;所选滤波技术为中值滤波。
本发明有效提高了图像处理速度与精确度,对非条码区域的清除更加有效,图像的处理与分析结果也越接近真实,最终也达到比较满意的识别效果。
上述为本发明较佳的实施方式,但本发明的实施方式并不受上述内容的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (1)

1.一种基于亚像素级角点检测的条码识别算法,其特征在于,包括下述步骤:
S1,载入原始条码图像;
S2,原始条码图像前期处理;
S3,判断并清除非条码区域图像,具体步骤如下:
首先进行分块扫描模式,所选取分块数量为6,所选取分块尺寸blockSize=(图像宽和高的最小值)/6,即blockSize分块尺寸是根据原始条码图像行和列的最小值除于分块数量来计算的,而图像的每行分块数量block_cols和每列分块数量block_rows,由下述公式(1)和(2)所示:
block_cols=blurmage.cols/blockSize, (1)
block_rows=blurmage.rows/blockSize, (2)
其中,blockSize为分块尺寸大小,blurmage.cols为图像列的大小,blurmage.rows为图像行的大小;
然后对所得图像进行灰度梯度扫描操作:在每个分块的两条对角线、中竖线和中横线上作扫描,即在这四条线上检测连续4个像素之间的灰度差,当灰度差的绝对值大于阈值时,标记为一次灰度变化,如下述公式(3)所示:
其中,Vi为第i个像素的梯度量化值,|fi+1+fi+fi-1+fi-2|为连续四个像素之间的灰度,t为阈值,若分块在某条扫描线上发生超过3次的灰度变化,则认为该分块可能含有条码部分,需要保留,否则清除分块;最后对所得图像进行8连通区域检查,每个分块都有8个领域分块,而分块的边缘除外,在统计每个条码区域分块的邻域分块中是条码区域分块的个数时,若达不到3个,则将该条码区域分块视为孤立分块,孤立分块是条码的可能性较小,因此将其标记为非条码区域,并清除该分块;
S4,对图像进行分块亚像素级角点检测,对每个分块的角点数进行统计,若分块的角点数小于等于预设阈值,则视为条码区域,具体步骤如下:
亚像素级角点检测的实现是基于对向量正交性的观测而实现的,即从中央点q到其领域点p的向量和p点处的图像梯度正交,如下述公式(4)所示:
其中表示在q的一个领域点pi处的图像梯度,则q的值通过最小化εi得到;设置εi为0后,公式两边同乘以即可得到关于该领域点的一个等式;当对所有领域点均进行上述变换,等式左右分别相加,即可得到下述公式(5):
设定最终可求解q=G-1·b;
通过上述公式(5)确定亚像素级角点,所设置的阈值为3,若分块的角点数小于等于预设阈值,则视为条码区域,否则即为非条码区域;
S5,条码区域分块再判断,具体步骤如下:
条码是将宽度不等的多个黑条和空白,按照一定的编码规则排列,用以表达一组信息的图形标识符,根据条码的条形排列的特点,通过霍夫变换算法进行直线检测:
运用直线的极坐标公式:ρ=хcos(θ)+уsin(θ),其中ρ,θ是一对参数,通过对图像中的像素点运用该公式处理即可得到二维的ρ,θ参数空间,在参数空间上形成峰值,从而利用该峰值即可检测出直线;最后将直线检测和角点检测的判定结果进行或运算,作为分块的最终结果,以防止出现误判;
S6,根据剩余图像轮廓定位条码位置,旋转并截取保存条码区域,具体步骤如下:
查找并绘制条码区域的所有轮廓,用矩形勾出最大的轮廓视为条码位置,计算条码黑条倾斜角,以矩形中心进行原图像旋转,截取并保存原图中所在区域;
S7,识别条码。
CN201810663401.3A 2018-06-25 2018-06-25 一种基于亚像素级角点检测的条码识别算法 Active CN109190434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810663401.3A CN109190434B (zh) 2018-06-25 2018-06-25 一种基于亚像素级角点检测的条码识别算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810663401.3A CN109190434B (zh) 2018-06-25 2018-06-25 一种基于亚像素级角点检测的条码识别算法

Publications (2)

Publication Number Publication Date
CN109190434A true CN109190434A (zh) 2019-01-11
CN109190434B CN109190434B (zh) 2021-08-24

Family

ID=64948529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810663401.3A Active CN109190434B (zh) 2018-06-25 2018-06-25 一种基于亚像素级角点检测的条码识别算法

Country Status (1)

Country Link
CN (1) CN109190434B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109784121A (zh) * 2019-01-14 2019-05-21 新大陆数字技术股份有限公司 Dot-peen DPM码的识别方法及装置
CN110659529A (zh) * 2019-09-03 2020-01-07 珠海市杰理科技股份有限公司 条形码识别方法、装置、计算机设备和存储介质
CN113673410A (zh) * 2021-08-16 2021-11-19 北京大学 人工标记的生成方法、识别方法、装置、设备、介质和程序
CN115577728A (zh) * 2022-12-07 2023-01-06 深圳思谋信息科技有限公司 一维码定位方法、装置、计算机设备及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0887760A2 (en) * 1992-08-10 1998-12-30 United Parcel Service Of America, Inc. Method and apparatus for decoding bar code symbols using subpixel interpolation
CN101763643A (zh) * 2010-01-07 2010-06-30 浙江大学 一种结构光三维扫描仪系统自动标定方法
CN101908122A (zh) * 2010-06-01 2010-12-08 福建新大陆电脑股份有限公司 条空边界处理模块、条码识别装置及其方法
CN103279935A (zh) * 2013-06-09 2013-09-04 河海大学 基于map算法的热红外遥感图像超分辨率重建方法及系统
CN105975894A (zh) * 2016-05-11 2016-09-28 华中科技大学 一种基于自适应边缘检测和映射模型的一维码识别算法
CN107240112A (zh) * 2017-06-28 2017-10-10 北京航空航天大学 一种复杂场景下个体x角点提取方法
CN107633192A (zh) * 2017-08-22 2018-01-26 电子科技大学 一种基于机器视觉的复杂背景下条形码分割与识读方法
CN107833186A (zh) * 2017-10-26 2018-03-23 长沙全度影像科技有限公司 一种基于Encoder‑Decoder深度学习模型的单透镜空间变化图像复原方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0887760A2 (en) * 1992-08-10 1998-12-30 United Parcel Service Of America, Inc. Method and apparatus for decoding bar code symbols using subpixel interpolation
CN101763643A (zh) * 2010-01-07 2010-06-30 浙江大学 一种结构光三维扫描仪系统自动标定方法
CN101908122A (zh) * 2010-06-01 2010-12-08 福建新大陆电脑股份有限公司 条空边界处理模块、条码识别装置及其方法
CN103279935A (zh) * 2013-06-09 2013-09-04 河海大学 基于map算法的热红外遥感图像超分辨率重建方法及系统
CN105975894A (zh) * 2016-05-11 2016-09-28 华中科技大学 一种基于自适应边缘检测和映射模型的一维码识别算法
CN107240112A (zh) * 2017-06-28 2017-10-10 北京航空航天大学 一种复杂场景下个体x角点提取方法
CN107633192A (zh) * 2017-08-22 2018-01-26 电子科技大学 一种基于机器视觉的复杂背景下条形码分割与识读方法
CN107833186A (zh) * 2017-10-26 2018-03-23 长沙全度影像科技有限公司 一种基于Encoder‑Decoder深度学习模型的单透镜空间变化图像复原方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109784121A (zh) * 2019-01-14 2019-05-21 新大陆数字技术股份有限公司 Dot-peen DPM码的识别方法及装置
CN109784121B (zh) * 2019-01-14 2022-09-23 新大陆数字技术股份有限公司 Dot-peen DPM码的识别方法及装置
CN110659529A (zh) * 2019-09-03 2020-01-07 珠海市杰理科技股份有限公司 条形码识别方法、装置、计算机设备和存储介质
CN110659529B (zh) * 2019-09-03 2023-06-02 珠海市杰理科技股份有限公司 条形码识别方法、装置、计算机设备和存储介质
CN113673410A (zh) * 2021-08-16 2021-11-19 北京大学 人工标记的生成方法、识别方法、装置、设备、介质和程序
CN115577728A (zh) * 2022-12-07 2023-01-06 深圳思谋信息科技有限公司 一维码定位方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN109190434B (zh) 2021-08-24

Similar Documents

Publication Publication Date Title
US11403839B2 (en) Commodity detection terminal, commodity detection method, system, computer device, and computer readable medium
CN109190434A (zh) 一种基于亚像素级角点检测的条码识别算法
CN110286126A (zh) 一种基于视觉图像的晶圆表面缺陷分区域检测方法
Kumar Mishra et al. A review of optical imagery and airborne lidar data registration methods
US11676301B2 (en) System and method for efficiently scoring probes in an image with a vision system
CN109426814B (zh) 一种发票图片特定板块的定位、识别方法、系统、设备
US10528781B2 (en) Detection method and system for characteristic patterns of Han Xin codes
CN106446894A (zh) 一种基于轮廓识别球型目标物位置的方法
CN105426844B (zh) 一种答题卡识别方法
CN109858325B (zh) 一种表格检测方法和装置
CN107851327A (zh) 粗细搜索方法及图像处理装置
CN112233116B (zh) 基于邻域决策与灰度共生矩阵描述的凹凸痕视觉检测方法
CN106056020A (zh) 一种基于机器视觉系统的条码检测识别方法及系统
US20150302236A1 (en) Method and device for identifying a two-dimensional barcode
Xu et al. An automatic optical and SAR image registration method with iterative level set segmentation and SIFT
CN112990183B (zh) 离线手写汉字同名笔画提取方法、系统、装置
CN105046200A (zh) 基于直线检测的电子阅卷方法
Xu et al. Improved SIFT match for optical satellite images registration by size classification of blob-like structures
CN101833644A (zh) 一种基于动态模板的校正图形搜索方法
CN111563896A (zh) 一种用于接触网异常检测的图像处理方法
CN104881641B (zh) 基于移动设备的问卷和表格数字化识别方法及系统
Safdarinezhad et al. An automatic method for precise 3D registration of high resolution satellite images and Airborne LiDAR Data
Pritts et al. Minimal solvers for rectifying from radially-distorted scales and change of scales
CN113822092B (zh) 用于定位位置探测图形的方法及装置、电子设备、介质
CN113743360B (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