一种图像编码方法
【技术领域】
[0001] 本发明设及数字信号处理技术领域,尤其设及一种图像编码方法。
【背景技术】
[0002] 图像压缩是保证图像有效应用的基础,JPEGQointPhotographicExperts Group)是在国际标准化组织(ISO)领导下制定的第一个国际图像压缩标准。由于JPEG优良 的品质,使它在短短几年内获得了成功,被广泛应用于互联网和数码相机领域,网站上80% 的图像都采用了JPEG压缩标准。
[000引量化是压缩过程中减少人眼视觉冗余的有效手段,因此,根据图像内容特性,选择 合适的量化表可W提高图像的压缩效率。另一方面,人眼对于不同的图像内容中的失真具 有不同的敏感程度,但是传统的图像压缩对图像所有区域采用相同的量化方法限制了压缩 效率的提升。
【发明内容】
[0004] 在保持与JPEG压缩标准兼容的条件下,为了进一步提升图像压缩效率,本申请提 供了一图像编码方法,通过率失真优化的量化表计算方法和基于人眼视觉特性的区域自适 应系数收缩方法对图像进行编码。
[0005] 进一步地,包括如下步骤,
[0006] 步骤1)对数个图像块进行DCT变换;
[0007] 步骤2)分别统计所述数个图像块的每个频带系数在不同量化步长下的失真和码 率.
[000引步骤3)选择不超过最大失真的量化步长对应的失真相对于码率的变化率的相反 数的均值作为率失真优化中的拉格朗日因子;
[0009] 步骤4)通过计算得到量化表并进行量化;
[0010] 步骤5)在量化过程中进行基于人眼视觉特性的系数自适应收缩并得到压缩后的 图像。
[001U进一步地,步骤2)还包含,计算每个频带中量化系数在量化步长Qst巧为minQ到maxQ(其中minQ<max曲时对应的失真,其中量化步长为整数,统计所有块中每个频带量化 系数的概率分布,利用量化系数的滴估计JPEG编码的码率。
[0012] 进一步地,根据公式(1)计算每个频带在不同量化步长下的拉格朗日因子, 入i,j(Qst巧); r1 0 修邸 + 但、'邸)"、
[0013] i- , (Qstep) = -------- Cl)。 I R;.i[Qx邸)-R;.人Qs邮+ \)
[0014] 进一步地,根据用户输入的编码质量,计算图像每个频带的平均失真Di。,并选择在 保证失真不超过町。条件下,选择失真不超过D1。的最大量化步长对应的拉格朗日因子的均 值,A。,作为率失真决策中的拉格朗日因子。
[0015] 进一步地,根据公式(2)计算的量化步长作为每个频带的量化表的值
[0016] "带)+為巧:,(0'"啤)(2)。
[0017] 进一步地,据图像块像素方差将图像块分为平滑区域,边缘区域和纹理区域;不同 区域采用处理不同的系数收缩个数;
[001引进一步地,图像块方差小于给定阔值化rl时,该图像块为平滑块,图像块方差大 于给定阔值化r2时,该图像块为边缘块,图像块方差介于化rl和化r2之间时该图像块为 纹理块。
[0019] 进一步地,每个图像块按照反Zig-Zag顺序对系数进行收缩处理,其中每种块需 要处理的系数个数可W由用户指定,或者按照默认值进行。
[0020] 进一步地,进行系数收缩的默认值设置为;边缘块最小,纹理块次之,平滑块最大。
[0021] 本发明通过率失真优化的量化表计算方法和基于人眼视觉特性的区域自适应系 数收缩方法对图像进行编码,在压缩图像主观质量相同的条件下,减小图像文件大小。
【附图说明】
[0022] 为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用 的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本 领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可W根据该些附图获得其他 的附图。
[0023] 图1为本发明实施例提供的一种图像编码方法的流程图;
[0024] 图2为本发明实施例提供的一种图像编码方法的反Zig-Zag扫描顺序图;
[0025] 图3为本发明实施例采用1_JPEG压缩和自适应量化表压缩时的性能示意图;
[0026] 图4为本发明实施例采用2_JPEG压缩和自适应量化表压缩时的性能示意图。
【具体实施方式】
[0027] 下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始 至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参 考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。 [002引本技术领域技术人员可W理解,除非特意声明,该里使用的单数形式"一"、"一 个"、"所述"和"该"也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措 辞"包括"是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加 一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元 件被"连接"或"禪接"到另一元件时,它可W直接连接或禪接到其他元件,或者也可W存在 中间元件。此外,该里使用的"连接"或"禪接"可W包括无线连接或禪接。该里使用的措 辞"和/或"包括一个或更多个相关联的列出项的任一单元和全部组合。
[0029] 本技术领域技术人员可W理解,除非另外定义,该里使用的所有术语(包括技术 术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应 该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的 意义一致的意义,并且除非像该里一样定义,不会用理想化或过于正式的含义来解释。
[0030] 为便于对本发明实施例的理解,下面将结合附图W几个具体实施例为例做进一步 的解释说明,且各个实施例并不构成对本发明实施例的限定。
[0031] 一种制作图像编码方法,通过率失真优化的量化表计算方法和基于人眼视觉特性 的区域自适应系数收缩方法对图像进行编码。
[00对实施例1
[0033] 一种制作图像编码方法,包括如下步骤,
[0034] 步骤1)对数个图像块进行DCT变换;
[0035] 步骤2)分别统计所述数个图像块的每个频带系数在不同量化步长下的失真和码 率.
[0036] 步骤3)选择不超过最大失真的量化步长对应的失真相对于码率的变化率的相反 数的均值作为率失真优化中的拉格朗日因子;
[0037] 步骤4)通过计算得到量化表并进行量化;
[003引步骤5)在量化过程中进行基于人眼视觉特性的系数自适应收缩并得到压缩后的 图像。
[0039] 实施例2
[0040] 一种制作图像编码方法,其包含的步骤与实施例1相同;区别在于步骤2)还包含, 计算每个频带中量化系数在量化步长Qstep为minQ到maxQ(其中minQ<max曲时对应的失 真,其中量化步长为整数,统计所有块中每个频带量化系数的概率分布,利用量化系数的滴 估计JPEG编码的码率。
[0041] 实施例3
[0042] 一种制作图像编码方法,其包含的步骤与实施例2相同;区别在于根据公式(1)计 算每个频带在不同量化步长下的拉格朗日因子,Au(Qst巧);
[0043]
[0044] 其中,Du(Qstep)和Ru(Qstep)分别表示对应(i,j)频带在量化步长为Qst巧时 的失真和码
率,N是图像块大小,在JPEG图像中N=8。
[0045] 实施例4
[0046] 一种制作图像编码方法,其包含的步骤与实施例3相同;区别在于根据用户输入 的编码质量,计算图像每个频带的平均失真町。,并选择在保证失真不超过町。条件下,选择 失真不超过町。的最大量化步长对应的拉格朗日因子的均值,A。,作为率失真决策中的拉格 朗日因子。
[0047] 实施例5
[0048] 一种制作图像编码方法,其包含的步骤与实施例4相同;区别在于根据用户输入 的压缩质量,计算图像每个频带允许的最大失真,记作D。。对于每个频带,选择保证失真不 超过D。的最大量化步长对应的拉格朗日因子作为候选拉格朗日因子,所有频带候选拉格朗 日因子的均值作为量化步长决策中的拉格朗日因子,记作A。。根据公式(2)计算的量化步 长作为每个频带的量化表的值
[0049]
[0化0] 该样,计算得到每个频带的量化步长,从而构成JPEG编码压缩的量化矩阵输入。
[0化1] 实施例6
[0化2] -种制作图像编码方法,其包含的步骤与实施例5相同;区别在于据图像块像素 方差将图像块分为平滑区域,边缘区域和纹理区域;不同区域采用处理不同的系数收缩个 数;W图像亮度分量中8X8块的像素方差作为图像区域划分的统计量。计算亮度分量中每 个8X8块的像素方差Var(i,j)(该里i和j表示图像块坐标),并根据阔值化rl,T虹2判 断对应的8X8块所属区域。
[0化3]实施例7
[0化4] 一种制作图像编码方法,其包含的步骤与实施例6相同;区别在于图像块方差小 于给定阔值化rl时,该图像块为平滑块,图像块方差大于给定阔值化r2时,该图像块为边 缘块,图像块方差介于化rl和化r2之间时该图像块为纹理块。对于Var(i,j)小于阔值 化rl的区域记作平滑区域,Var(i,j)大于化r2为边缘区域,Var(i,j)在化rl和化r2之 间的8X8块记作纹理区域。
[0化5]实施例8
[0化6] -种制作图像编码方法,其包含的步骤与实施例7相同;区别在于每个图像块按 照反Zig-Zag顺序对系数进行收缩处理,其中每种块需要处理的系数个数由用户指定,或 者按照默认值进行。系数收缩按照公式(3)进行,
[0057]
[0化引其中,T是判断系数是否收缩的阔值,t是收缩的幅值。
[0化9] 1)计算当前编码的8X8块的像素方差,如果该方差小于化rl那么标记该块为平 滑块,如果方差大于化r2那么标记该块为边缘块,否则标记为纹理块;
[0060]。按照反Zig-Zag顺序依次对8x8块内的系数进行自适应收缩,如公式做所示, 其中不同类型的块处理的系数的个数根据用户输入的值确定,如果用户没有输入,则选择 默认值;
[0061] 为了验证提出的算法性能,表1列出了采用提出的优化应量化表和标准JPEG进行 压缩得到的编码性能对比,可W发现,提出的自适应量化表可W在客观质量一致的条件下, 节约10%左右的码率,图3和图4为采用JPEG压缩和自适应量化表压缩时的性能,其中 "x_JPEG"和"x_自适应量化表"中的数字"X"表示不同图像的编号。
[0062] 表1提出的自适应量化表和标准JPEG压缩性能比较炬DRATE)
[0063]
[0065]为了验证提出的自适应系数收缩方法,表2列出了系数收缩的JPEG压缩图像和标 准JPEG压缩的性能对比结果,可W发现,在质量接近的条件下,提出的方法可W节约1. 3% 左右的码率,其中质量评价方法是无参考图像质量评价方法[1],该分值越低代表图像质量 越好。
[0066] 表2图像区域自适应系数收缩方法和标准JPEG压缩性能对比
[0067]
[0068] 实施例9
[0069] 一种制作图像编码方法,其包含的步骤与实施例8相同;区别在于进行系数收缩 的默认值设置为;边缘块最小,纹理块次之,平滑块最大。
[0070] 本发明通过率失真优化的量化表计算方法和基于人眼视觉特性的区域自适应系 数收缩方法对图像进行编码,在压缩图像主观质量相同的条件下,减小图像文件大小。
[0071]W上,仅为本发明较佳的【具体实施方式】,但本发明的保护范围并不局限于此,任何 熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应 涵盖在本发明的保护范围之内。因此,本发明的保护范围应该W权利要求的保护范围为准。
【主权项】
1. 一种图像编码方法,其特征在于,通过率失真优化的量化表计算方法和基于人眼视 觉特性的区域自适应系数收缩方法对图像进行编码。2. 根据权利要求1所述的图像编码方法,其特征在于包括如下步骤, 步骤1)对数个图像块进行DCT变换; 步骤2)分别统计所述数个图像块的每个频带系数在不同量化步长下的失真和码率; 步骤3)选择不超过最大失真的量化步长对应的失真相对于码率的变化率的相反数的 均值作为率失真优化中的拉格朗日因子; 步骤4)通过计算得到量化表并进行量化; 步骤5)在量化过程中进行基于人眼视觉特性的系数自适应收缩并得到压缩后的图 像。3. 根据权利要求2所述的图像编码方法,其特征在于所述步骤2)还包含, 计算每个频带中量化系数在量化步长Qstep为minQ到maxQ (其中minQ〈maxQ)时对应 的失真,其中量化步长为整数,统计所有块中每个频带量化系数的概率分布,利用量化系数 的熵估计JPEG编码的码率。4. 根据权利要求3所述的图像编码方法,其特征在于, 根据公式(1)计算每个频带在不同量化步长下的拉格朗日因子,λ u (Qstep);5. 根据权利要求4所述的图像编码方法,其特征在于, 根据用户输入的编码质量,计算图像每个频带的平均失真Din,并选择在保证失真不超 过Din条件下,选择失真不超过Din的最大量化步长对应的拉格朗日因子的均值,λ ^,作为 率失真决策中的拉格朗日因子。6. 根据权利要求5所述的图像编码方法,其特征在于, 根据公式(2)计算的量化步长作为每个频带的量化表的值7. 根据权利要求1所述的图像编码方法,其特征在于, 根据图像块像素方差将图像块分为平滑区域,边缘区域和纹理区域;不同区域采用处 理不同的系数收缩个数。8. 根据权利要求7所述的图像编码方法,其特征在于, 图像块方差小于给定阈值Thrl时,该图像块为平滑块,图像块方差大于给定阈值Thr2 时,该图像块为边缘块,图像块方差介于Thrl和Thr2之间时该图像块为纹理块。9. 根据权利要求8所述的图像编码方法,其特征在于, 每个图像块按照反Zig-Zag顺序对系数进行收缩处理,其中每种块需要处理的系数个 数可以由用户指定,或者按照默认值进行。10. 根据权利要求9所述的图像编码方法,其特征在于, 进行系数收缩的默认值设置为:边缘块最小,纹理块次之,平滑块最大。
【专利摘要】本发明实施例公开了一种图像编码方法,通过率失真优化的量化表计算方法和基于人眼视觉特性的区域自适应系数收缩方法对图像进行编码;本发明通过率失真优化的量化表计算方法和基于人眼视觉特性的区域自适应系数收缩方法对图像进行编码,在压缩图像主观质量相同的条件下,减小图像文件大小。
【IPC分类】H04N19/625, H04N19/147, H04N19/124
【公开号】CN104902285
【申请号】CN201510262890
【发明人】马思伟, 张新峰, 王诗淇
【申请人】北京大学
【公开日】2015年9月9日
【申请日】2015年5月21日