CN108419086B - 一种基于霍夫曼编码的图像压缩方法 - Google Patents

一种基于霍夫曼编码的图像压缩方法 Download PDF

Info

Publication number
CN108419086B
CN108419086B CN201810430127.5A CN201810430127A CN108419086B CN 108419086 B CN108419086 B CN 108419086B CN 201810430127 A CN201810430127 A CN 201810430127A CN 108419086 B CN108419086 B CN 108419086B
Authority
CN
China
Prior art keywords
image
offset
pixel
center reference
reference pixel
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
Application number
CN201810430127.5A
Other languages
English (en)
Other versions
CN108419086A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810430127.5A priority Critical patent/CN108419086B/zh
Publication of CN108419086A publication Critical patent/CN108419086A/zh
Application granted granted Critical
Publication of CN108419086B publication Critical patent/CN108419086B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/90Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using coding techniques not provided for in groups H04N19/10-H04N19/85, e.g. fractals
    • H04N19/91Entropy coding, e.g. variable length coding [VLC] or arithmetic coding

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)

Abstract

本发明公开的一种基于霍夫曼编码的图像压缩方法,属于数字图像处理技术领域。本发明读取图像,解析获取图像信息;进行图像分块处理;对分块后的子图像进行图像处理,生成偏移量替换后的更新图像;遍历整幅图像,概率统计得到整幅图像除中心基准像素点以外的所有像素点的偏移量;建立霍夫曼二叉树,得到霍夫曼码并存入元胞数组;根据霍夫曼码表对偏移量替换后的更新图像进行查表,得到偏移量对应的不定长霍夫曼码,并将偏移量对应的霍夫曼码储存在与原始图像相同大小的元胞数组中;反查表生成R、G、B三个通道原始图像;将生成的R、G、B三个通道原始图像合成,恢复原始图像。本发明能减少图像空间冗余度,实现图像无损压缩与恢复。

Description

一种基于霍夫曼编码的图像压缩方法
技术领域
本发明属于数字图像处理技术领域,尤其涉及一种基于霍夫曼编码的图像压缩方法。
背景技术
随着消费者对高品质图像的追求,现在图像采集设备的像素越来越高,这对图像的存储、传输、处理都带来极大的挑战。所以图像压缩技术便成为热门研究领域。图像实现压缩的意义在于减少数据存储量,节省存储空间;降低数据率以减少传输时的使用带宽,节省传输时间;压缩图像的信息量,便于特征提取。
目前业内根据人眼视觉特性和图像空间、时间特性进行压缩编码,以减少图像数据量,但是常用的压缩方式均存在局限:
JPG、PNG等图像压缩方式,压缩效率高,但是由于采用空间频域变化,对低频信息的损失较大,对图像原始数据的信息保留较少,不适合对原始图像的恢复。
JPEG2000无损压缩方式,过程繁复,实现方法复杂。
常见的视频压缩编码(H264/MPEG)在空域和时域采用多种压缩算法实现实时视频图像的压缩,但是图像信息损失较大,不适合对原始图像的恢复。
发明内容
本发明公开的一种基于霍夫曼编码的图像压缩方法,要解决的问题是提供一种图像无损压缩的方法,并解压缩完全恢复原始图像,基于霍夫曼编码实现对图像的压缩,减少图像空间冗余度,使所述的图像压缩方法具有能很好地恢复原始图像、实现方法简单等优点。
本发明目的是通过下述技术方案实现的。
本发明公开的一种基于霍夫曼编码的图像压缩方法,读取图像,解析获取图像信息;进行图像分块处理;对分块后的每一块子图像进行图像处理,生成偏移量替换后的更新图像;遍历整幅图像除中心基准像素点以外的所有像素点,概率统计得到整幅图像除中心基准像素点以外的所有像素点的偏移量;建立霍夫曼二叉树,得到霍夫曼码并存入元胞数组;根据生成霍夫曼码表对偏移量替换后的更新图像进行查表,得到偏移量对应的不定长霍夫曼码,并将偏移量对应的霍夫曼码储存在与原始图像相同大小的元胞数组中;反查表生成R、G、B三个通道原始图像;将生成的R、G、B三个通道原始图像合成,恢复原始图像。
本发明公开的一种基于霍夫曼编码的图像压缩方法,包括以下步骤:
步骤一:读取图像,解析获取图像信息,所述的图像信息包括图像宽M、高N、像素位深d和R、G、B三个通道图像。
步骤二:对步骤一获取的R、G、B通道图像分别进行m×m分块处理,每一个通道图像得到n=max(M,N)/m块子图像,每块子图像包含m2个像素点。当M、N都能够被m整除时n=max(M,N)/m;当M、N其中一个不能被m整除或者两者都不能被m整除时n=max(M,N)/m,对于剩余的行l=N%m或列c=M%m的每一个像素点与离得最近的子图像块进行处理。
为提高图像的压缩率,m×m分块处理中m优选奇数,进一步优选3×3分块处理。
步骤三:对步骤二分块后的每一块子图像进行图像处理,图像处理方法如下:保证每一块子图像中心基准像素不变,以图像中心像素点为基准,计算周围m2-1个像素点的像素偏移量,计算公式如(1)所示,其中(i,j)是中心基准像素点的像素值:
并用计算后像素点的偏移量替换原始图像对应像素点的像素值,即生成偏移量替换后的更新图像。
对于步骤二所述的当M、N其中一个不能被m整除或者两者都不能被m整除时n=max(M,N)/m,将剩余的行l=N%m或列c=M%m中的每一个像素点与离得最近的子图像块的中心基准像素计算偏移量。
步骤四:遍历整幅图像除中心基准像素点以外的所有像素点,概率统计得到整幅图像除中心基准像素点以外的所有像素点的偏移量qi,并存储在数组相应的位置si中。偏移量qi与存储位置si之间的关系如公式(2)所示:
si=qi+2d (2)
步骤五:根据步骤四概率统计整幅图像除中心基准像素点以外的所有像素点的偏移量,建立霍夫曼二叉树,得到除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码,并将对应的霍夫曼码存入元胞数组,即生成霍夫曼码表。
步骤六:根据步骤五生成霍夫曼码表对步骤三偏移量替换后的更新图像进行查表,得到除中心基准像素点以外的所有像素点的偏移量对应的不定长霍夫曼码,并将偏移量对应的霍夫曼码储存在与原始图像相同大小的元胞数组中。所述的与原始图像相同尺寸大小的元胞数组包括储存偏移量对应的霍夫曼码和中心基准像素点的像素值。所述偏移量对应的霍夫曼码指除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码。
步骤七:对R、G、B三个通道图像分别进行步骤三至步骤六所述的操作,最后得到三个与原始图像相同尺寸大小的元胞数组,计算每个与原始图像相同尺寸大小的元胞数组的平均码长di,(i=r、g、b),与原始像素位深d计算压缩比pi
压缩比pi计算公式如(3)所示:
pi=di/d (3)
步骤八:对步骤七生成的三个与原始图像相同尺寸大小的元胞数组分别进行查表处理得到除中心基准像素点以外的所有像素点的偏移量,并将查表处理得到的除中心基准像素点以外的所有像素点的偏移量分别与对应的中心基准像素点作加法处理,得到对应像素点的原始像素值,生成R、G、B三个通道原始图像。
步骤八所述的查表处理指根据步骤五生成霍夫曼码表进行查表处理。
步骤九:将步骤八生成的R、G、B三个通道原始图像合成,恢复原始图像。
有益效果:
1、本发明公开的一种基于霍夫曼编码的图像压缩方法,根据步骤四所述遍历整幅图像除中心基准像素点以外的所有像素点,概率统计整幅图像除中心基准像素点以外的所有像素点的偏移量,即仅通过概率统计遍历偏移量,而不用进行时域和频域的变换,相比现有技术中的DCT方法处理过程简单,能最大化保持原始图像信息。
2、本发明公开的一种基于霍夫曼编码的图像压缩方法,根据图像相邻像素点的相似性,对步骤一读取图像的R、G、B通道分别进行m×m分块处理,由于对图像分块操作,能够减少图像空间冗余度,满足图像局部信息的快速分析和处理要求。
3、本发明公开的一种基于霍夫曼编码的图像压缩方法,由于根据步骤五生成的霍夫曼码表进行图像压缩,又通过步骤八查表处理恢复原始图像,即通过霍夫曼编码与转换实现图像的无损压缩,能最大限度保持原始数据信息,很好地恢复原始图像。
附图说明
图1是本发明的流程图;
图2是本发明步骤四所述的偏移量概率存储数组示意图;
图3是本发明实施例的图像块分块示意图;
图4是本实施例的图像块像素偏差量计算示意图;
图5是本实施例的R通道部分偏移量霍夫曼码表;
图6是本实施例的R通道第一个图像块编码;
图7是本实施例解压缩后图像误差(扩大10倍)与JPG解压缩后图像误差(扩大10倍)对比图。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
为了验证方法的可行性,选择一幅分辨率为1024×1024未经过压缩TIFF格式图像,R、G、B每通道的位深为8bit,进行3×3分块处理。
如图1所示,本实施例公开一种基于霍夫曼编码的图像压缩方法,包括以下步骤:
步骤一:读取图像,得到图像宽M=1024、高N=1024、像素位深d=8和R、G、B三个通道图像。
步骤二:对步骤一获取的R、G、B通道图像分别进行3×3分块处理,分块方法如图3所示,每一个通道图像得到n=341块子图像,每块子图像包含m2=9个像素点。因为M=1024、N=1024都不能被m=3整除,对剩余的行l=1和列c=1中的每一个像素点与离得最近的子图像块进行处理。
为提高图像的压缩率,m×m分块处理优选3×3分块处理,并且m一般选择奇数,这里选择3×3。
步骤三:对步骤二分块后的每一块子图像进行图像处理,图像处理方法如下:如图4所示保证每一块子图像中心基准像素不变,以图像中心像素点为基准,计算周围m2-1=8个像素点的像素偏移量,每一块子图像偏移量计算公式如(1)所示:
并用计算后像素点的偏移量替换原始图像对应像素点的像素值,即生成偏移量替换后的更新图像。
对于步骤二所述的剩余的行l=1和列c=1中的每一个像素点与离得最近的子图像块的中心基准像素计算偏移量。
步骤四:遍历整幅图像除中心基准像素点以外的所有像素点,概率统计得到整幅图像除中心基准像素点以外的所有像素点的偏移量qi,并存储在数组相应的位置si中。偏移量qi与存储位置si之间的关系如公式(2)所示:
si=qi+2d (2)
其中,d=8,偏移量qi的取值范围为-255到256,si的取值范围为1到512。
步骤五:根据步骤四概率统计整幅图像除中心基准像素点以外的所有像素点的偏移量,建立霍夫曼二叉树,得到除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码,并将对应的霍夫曼码存入元胞数组,即生成霍夫曼码表,R通道部分偏移量霍夫曼码表如图5所示。
步骤六:根据步骤五生成霍夫曼码表对步骤三偏移量替换后的更新图像进行查表,得到除中心基准像素点以外的所有像素点的偏移量对应的不定长霍夫曼码,并将偏移量对应的霍夫曼码储存在与原始图像相同大小的元胞数组中,R通道中第一个图像块编码如图6所示。所述的与原始图像相同尺寸大小的元胞数组包括储存偏移量对应的霍夫曼码和中心基准像素点的像素值。所述偏移量对应的霍夫曼码指除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码。
步骤七:对R、G、B三个通道图像分别进行步骤三至步骤六所述的操作,最后得到三个与原始图像相同尺寸大小的元胞数组,计算每个与原始图像相同尺寸大小的元胞数组的平均码长di,与原始像素位深d计算压缩比pi
压缩比计算公式如(3)所示:
pi=di/d (3)
其中dr=dg=db=6,d=8,得到所以基于霍夫曼编码对这幅图像进行压缩以后的压缩率为0.75。
步骤八:对步骤七生成的三个与原始图像相同尺寸大小的元胞数组分别进行查表处理得到除中心基准像素点以外的所有像素点的偏移量,并将查表处理得到的除中心基准像素点以外的所有像素点的偏移量分别与对应的中心基准像素点作加法处理,得到对应像素点的原始像素值,生成R、G、B三个通道原始图像。
步骤八所述的查表处理指根据步骤五生成霍夫曼码表进行查表处理。
步骤九:将步骤八生成的R、G、B三个通道原始图像合成,恢复原始图像。
通过将本实施例公开一种基于霍夫曼编码的图像压缩方法解压缩后的图像与JPG格式压缩方法解压缩后的图像进行对比,说明本实施例的优点。
如图7所示为基于霍夫曼编码的图像压缩方法解压缩后的图像与原始图像的误差与JPG格式压缩方法解压缩后的图像与原始图像误差对比,为了使对比更加明显,将两者的误差都扩大了10倍,由图7看到JPG解压缩后大部分像素点都不再是原始的像素值,已经失去大部分图像信息,而基于霍夫曼编码进行图像压缩,能够保留图像的所有原始信息,只去除了图像的空间冗余度,可以完全恢复原始图像,所以与原始图像相减得到全0的图像,显示为全黑。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于霍夫曼编码的图像压缩方法,其特征在于:包括以下步骤,
步骤一:读取图像,解析获取图像信息,所述的图像信息包括图像宽M、高N、像素位深d和R、G、B三个通道图像;
步骤二:对步骤一获取的R、G、B通道图像分别进行m×m分块处理,每一个通道图像得到n=max(M,N)/m块子图像,每块子图像包含m2个像素点;当M、N都能够被m整除时n=max(M,N)/m;当M、N其中一个不能被m整除或者两者都不能被m整除时n=max(M,N)/m,对于剩余的行l=N%m或列c=M%m的所有像素点与离得最近的子图像块进行处理;
步骤三:对步骤二分块后的每一块子图像进行图像处理,图像处理方法如下:保证每一块子图像中心基准像素不变,以图像中心像素点为基准,计算周围m2-1个像素点的像素偏移量,计算公式如(1)所示,其中(i,j)是中心基准像素点的像素值:
并用计算后像素点的偏移量替换原始图像对应像素点的像素值,即生成偏移量替换后的更新图像;
对于步骤二所述的当M、N其中一个不能被m整除或者两者都不能被m整除时n=max(M,N)/m,将剩余的行l=N%m或列c=M%m中的每一个像素点与离得最近的子图像块的中心基准像素计算偏移量;
步骤四:遍历整幅图像除中心基准像素点以外的所有像素点,概率统计得到整幅图像除中心基准像素点以外的所有像素点的偏移量qi,并存储在数组相应的位置si中;偏移量qi与存储位置si之间的关系如公式(2)所示:
si=qi+2d (2)
步骤五:根据步骤四概率统计整幅图像除中心基准像素点以外的所有像素点的偏移量,建立霍夫曼二叉树,得到除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码,并将对应的霍夫曼码存入元胞数组,即生成霍夫曼码表;
步骤六:根据步骤五生成霍夫曼码表对步骤三偏移量替换后的更新图像进行查表,得到除中心基准像素点以外的所有像素点的偏移量对应的不定长霍夫曼码,并将偏移量对应的霍夫曼码储存在与原始图像相同大小的元胞数组中;所述的与原始图像相同尺寸大小的元胞数组包括储存偏移量对应的霍夫曼码和中心基准像素点的像素值;所述偏移量对应的霍夫曼码指除中心基准像素点以外的所有像素点的偏移量对应的霍夫曼码;
步骤七:重复步骤三至步骤六,分别对R、G、B三个通道图像进行处理,最后得到三个与原始图像相同尺寸大小的元胞数组,计算每个与原始图像相同尺寸大小的元胞数组的平均码长di,(i=r、g、b),与原始像素位深d计算压缩比pi
压缩比pi计算公式如(3)所示:
pi=di/d (3)
步骤八:对步骤七生成的三个与原始图像相同尺寸大小的元胞数组分别进行查表处理得到除中心基准像素点以外的所有像素点的偏移量,并将查表处理得到的除中心基准像素点以外的所有像素点的偏移量分别与对应的中心基准像素点作加法处理,得到对应像素点的原始像素值,生成R、G、B三个通道原始图像;
步骤八所述的查表处理指根据步骤五生成霍夫曼码表进行查表处理;
步骤九:将步骤八生成的R、G、B三个通道原始图像合成,恢复原始图像。
2.如权利要求1所述的一种基于霍夫曼编码的图像压缩方法,其特征在于:为提高图像的压缩率,m×m分块处理中m选3×3分块处理。
CN201810430127.5A 2018-05-08 2018-05-08 一种基于霍夫曼编码的图像压缩方法 Active CN108419086B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810430127.5A CN108419086B (zh) 2018-05-08 2018-05-08 一种基于霍夫曼编码的图像压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810430127.5A CN108419086B (zh) 2018-05-08 2018-05-08 一种基于霍夫曼编码的图像压缩方法

Publications (2)

Publication Number Publication Date
CN108419086A CN108419086A (zh) 2018-08-17
CN108419086B true CN108419086B (zh) 2019-04-05

Family

ID=63138027

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810430127.5A Active CN108419086B (zh) 2018-05-08 2018-05-08 一种基于霍夫曼编码的图像压缩方法

Country Status (1)

Country Link
CN (1) CN108419086B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241898B (zh) * 2019-03-20 2020-11-10 吕衍荣 基于图像处理的数据显示方法
CN109982086B (zh) * 2019-04-10 2020-12-08 上海兆芯集成电路有限公司 图像压缩方法及其装置
CN110191345B (zh) * 2019-04-16 2020-07-03 浙江大学 基于霍夫曼编码的前景、背景分离的增量式无损压缩方法
CN110400362B (zh) * 2019-07-12 2023-02-24 重庆大学 一种基于图像的abaqus二维裂纹建模方法、***及计算机可读存储介质
CN111028135B (zh) * 2019-12-10 2023-06-02 国网重庆市电力公司电力科学研究院 一种图像文件修复方法
CN111147861A (zh) * 2020-01-02 2020-05-12 广州虎牙科技有限公司 图像压缩方法、装置、用户设备和计算机可读存储介质
CN116760952B (zh) * 2023-08-17 2023-10-20 山东欣晖电力科技有限公司 基于无人机的电力铁塔维护巡检方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882966A (zh) * 2003-11-18 2006-12-20 斯卡拉多股份公司 处理数字图像的方法以及图像表示格式
CN101715133A (zh) * 2008-09-30 2010-05-26 雅马哈株式会社 用于图像数据的无损压缩编码装置和解码装置
CN103647969A (zh) * 2011-07-06 2014-03-19 北京航空航天大学 一种基于对象的快速分形视频压缩与解压缩方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6351570B1 (en) * 1997-04-01 2002-02-26 Matsushita Electric Industrial Co., Ltd. Image coding and decoding apparatus, method of image coding and decoding, and recording medium for recording program for image coding and decoding

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882966A (zh) * 2003-11-18 2006-12-20 斯卡拉多股份公司 处理数字图像的方法以及图像表示格式
CN101715133A (zh) * 2008-09-30 2010-05-26 雅马哈株式会社 用于图像数据的无损压缩编码装置和解码装置
CN103647969A (zh) * 2011-07-06 2014-03-19 北京航空航天大学 一种基于对象的快速分形视频压缩与解压缩方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于Huffman编码的图像压缩解压研究;饶兴;《电脑知识与技术》;20110228(第4期);第887页-第889页
基于能量熵分布梯度与Huffman编码的IC图像压缩技术;梁忠伟等;《机械设计与制造》;20091231(第12期);第62页-第64页

Also Published As

Publication number Publication date
CN108419086A (zh) 2018-08-17

Similar Documents

Publication Publication Date Title
CN108419086B (zh) 一种基于霍夫曼编码的图像压缩方法
US7162091B2 (en) Intra compression of pixel blocks using predicted mean
US7181072B2 (en) Intra compression of pixel blocks using predicted mean
CN105578190B (zh) 应用于视频硬解码的无损压缩方法及***
US8086052B2 (en) Hybrid video compression method
US8073269B2 (en) Image decoding apparatus and image decoding method
EP1917813A2 (en) Image data processing
EP1886503A1 (en) Compression and decompression of media data
US8923389B1 (en) Unified scaling with differential coding for internal bit depth extension and reference frame compression
CN111402380A (zh) 一种gpu压缩纹理处理方法
CN112422989A (zh) 一种视频编码方法
US6584226B1 (en) Method and apparatus for implementing motion estimation in video compression
CN101511025B (zh) 图像压缩/解压缩的方法、装置
US8594194B2 (en) Compression method using adaptive field data selection
US8379715B2 (en) System and method for video compression using non-linear quantization and modular arithmetic computation
CN105635731B (zh) 高效视频编码的帧内预测参考点预处理方法
CN110611509B (zh) 一种红外温度数据压缩方法及***
CN107707917A (zh) 一种基于信息熵的视频自适应采样率设定方法
Tang et al. Feature Fusion Enhanced Super Resolution for Low Bitrate Screen Content Compression
CN107124614B (zh) 一种具有超高压缩比的图像数据压缩方法
JP5533302B2 (ja) 動画像符号化装置
CN110620922A (zh) 一种视频编码带宽的优化电路及其优化方法
WO2023168712A1 (zh) 零游程值编解码方法及视频编解码方法、装置和***
CN102905129B (zh) 静止图像的分布式编码方法
Kumar et al. Comparative analysis of image coding methods A State-of-the-Art Survey.

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Li Dong

Inventor after: Yang Kai

Inventor after: Gao Jianyin

Inventor after: Li Jianguo

Inventor after: Wu Yu

Inventor before: Yang Kai

Inventor before: Gao Jianyin

Inventor before: Li Jianguo

Inventor before: Wu Yu

GR01 Patent grant
GR01 Patent grant