一种基于小波包和矢量量化相结合的图像压缩编码方法
【专利摘要】本发明属于数字图像压缩编码和传输领域，其特征是：对图像进行多级小波变换，分解成低频系数和高频系数；对低频系数使用熵编码得到低频码流；对高频系数进一步进行小波包分解并根据分解后的各级高频子带系数在空间和频率上的对应关系，构成分类训练矢量；对训练矢量集进行能量阈值筛选并分别采用基于等失真准则的广义学习矢量量化方法来设计矢量量化器的码书，然后进一步使用熵编码进行无损压缩，从得到高频码流；最后对由低频码流和高频码流的压缩数据进行解码得到恢复图像。本发明能够在获得较大压缩比的情况下，保证图像的恢复质量和编码速度，并能较好的保留图像的高频细节信息。
【专利说明】一种基于小波包和矢量量化相结合的图像压缩编码方法
[0001]【技术领域】:本发明属于数字图像压缩编码和传输领域，涉及一种新的基于小波包变换和矢量量化相结合的图像压缩编码方法，是对传统图像有损压缩方法的创新和改进，可广泛应用于数据编码与传输、多媒体技术和地理信息系统等领域中。
[0002]【背景技术】:随着计算机技术、通讯技术和多媒体技术等领域的高速发展，数字图像因其具有信息丰富、获取便利等优势而被广泛的使用。但这也导致了数据量的极具增加，给图像数据的存储和传输带来了巨大的负担，因此研究行之有效的图像压缩编码技术有着重要的理论意义和实用价值。
[0003]矢量量化(Vector Quantization, VQ)技术是一种高效的有损数据压缩方法,它是以香农的速率失真理论作为理论基础，其基本原理是:首先基于输入矢量设计矢量量化的码书，码书由码字构成，码字的维数和输入矢量相同，而码书的规模远小于输入矢量集；然后使用矢量量化器码书中与输入矢量相匹配码字的索引来替代输入矢量进行存储和传输，而解码只需要快速的查表操作，因此具有压缩比大、编码速度快等优点。矢量量化过程如图1所示。
[0004]矢量量化最基本也是最重要的问题是设计性能良好的码书，码书的性能直接影响图像压缩编码的质量。LBG方法是码书设计的经典方法，该方法是基于局部寻优原理，但它具有计算繁琐，易产生无效码字，对初始码书敏感，适应能力差等缺点。另外一类被广泛使用的方法是基于人工神经网络的码书设计方法，利用神经网络的非线性函数逼近性能实现输入矢量集到码书的映射，其中典型的方法有竞争学习神经网络(Competitive LearningNeural Network, CLNN),它是一种硬判决的聚类算法,但存在的问题是往往导致局部最优，全局失真过大，甚至码字欠利用。
[0005]除了码书设计方法以外，影响码书质量的另外一个重要因素是训练输入矢量集的有效性，为了解决这一问题，目前很多研究工作将小波分析技术引入到矢量量化中，并取得了一定的效果。其基本思想是首先将图像进行小波变换，得到高频率子带和低频率子带；然后基于各频带系数的关系构成矢量集合；再应用矢量量化方法进行压缩编码。但现有方法还没有充分挖掘和利用小波变换后不同频带系数的在位置、方向和内容上的相关性。此外，对于高分辨率多细节图像，例如遥感图像、医学图像等，小波分解的高频系数部分其频率分辨率却较低，因此对图像细节的分辨能力不足，这往往会导致图像重建后图像的细节部分由较大的损失。
【发明内容】
:
[0006]发明目的:本发明提供一种基于小波包和矢量量化相结合的图像压缩编码方法，给出基于小波包变换系数的分类训练矢量组织方式和基于等失真广义学习矢量量化技术的矢量量化器码书设计方法，其目的在于解决现有技术的图像压缩比较低，矢量量化码书质量不高，恢复图像细节损失严重等不足，提出一种有效的基于小波包变化和矢量量化相结合的图像压缩编码方法。
[0007]技术方案:本发明是通过以下技术方案实施的:[0008]一种基于小波包和矢量量化相结合的图像压缩编码方法，其特征在于:该方法按以下步骤:
[0009]步骤(I):对于待压缩编码图像，对其进行多级小波变换，分解成低频系数和高频系数；
[0010]步骤(2):对于小波分解后的低频系数采用标量量化方法，再使用熵编码得到低频码流；
[0011]步骤(3):对于小波分解后各频带的高频系数进一步采用小波包技术分别对其进行分解；
[0012]步骤(4):根据各级高频子带系数在空间和频率上的对应关系，构成分类训练矢量;
[0013]步骤(5):根据矢量能量阈值对训练矢量进行筛选；
[0014]步骤(6):对于筛选后的分类矢量，分别采用基于等失真准则的广义学习矢量量化方法来设计矢量量化器的码书；
[0015]步骤(7):对矢量量化的码流进一步使用熵编码进行无损压缩，从得到高频码流；
[0016]步骤(8):对由低频码流和高频码流构成的压缩数据，采用上述编码的逆过程进行解码，再进行重建得到恢复图像。
[0017]在步骤(3)中，分别对小波变换后的第I级，第2级和第3级频带三个方向的高频系数进行二维小波包变换，将各频带各方向的高频系数进一步分解，分解后的系数能够更精确的描述图像各方向细节信息。
[0018]在步骤(4)中，各子带对应位置的系数代表原图像同一细节在不同分辨率下不同方向的描述，它们具有带间的相关性，根据这种空间位置和内容上的相关性来构造训练分类矢量，选择同一方向不同子带上同一位置的系数组成训练矢量，不同的级别频率子带选取不同数量的系数；再根据高频系数的方向性，将训练矢量进行分类，每类矢量用于设计各自的量化器码书，训练过程各自独立。
[0019]在步骤(6)中，将设计矢量量化器最优码书的等失真准则引入到广义学习矢量量化方法中，加入一种基于等失真广义学习矢量量化的码书设计方法，其主要步骤包括:
[0020]I)初始化码书C (O)，初始码字可从输入训练序列矢量中随机选取N个矢量；初始化码书中各码字的部分失真Di (O);设最大迭代次数为T，允许误差门限值ε >0，初始学习率为α (O)；
[0021]2)对于各学习步骤t = 1，".，Τ，计算每步的变化学习率α⑴；
[0022]①根据每个码字所属区域内的训练矢量计算各码字的部分失真Di (t)；
[0023]②基于各码字的部分失真，计算改进的失真测度为如,Ji);
[0024]③根据新的失真测度找出竞争获胜码字矢量；
[0025]④修正码书中的获胜码字和其他码字矢量；
[0026]3)调整码书中各码字分割区的所属训练矢量；
[0027]4)计算误差Et;
[0028]5)若误差小于给定阈值或达到最大学习步数，则停止，并确定为最终的码书；否则转步骤2)，进行重新计算。
[0029]优点及效果:[0030]1.本发明提出一种基于小波包变换和矢量量化相结合的图像压缩编码方法，该方法能够即保证了图像低频平滑部分低失真和稳定性，又保留了图像各方向的高频细节信息，使得压缩图像在获得高恢复质量的情况下具有较大的压缩比，同时也提高了图像的压缩编码速度。
[0031]2.本发明采用小波包对图像高频系数进行分解相对于常用的单纯小波分解压缩算法，能保留更细致的图像细节信息，这在高分辨率遥感卫星图像等需要保留图像边缘细节的领域中具有绝对的优势。同时，给出的基于小波包分解后的分类矢量组织方法，充分利用了不同频率子带的空间相关性和内容相关性，能够有效提高矢量量化码书的设计质量。此外将矢量集按照细节方向性进行分类，分别设计码书，体现了并行计算的思想，降低了算法的复杂度，极大地提高了码书生成的实时性，也解决了传统码书设计方法速度慢的问题。
[0032]3.本发明提出了一种新的矢量量化码书设计方法，该方法基于香浓的设计最优矢量量化器理论，将等失真准则引入到广义矢量量化学习算法中，并将其用于码书设计。该方法对竞争训练中的矢量评判标准进行重新定义，使得码书中各码子所属区域的部分失真趋于相等，也解决了码子的欠利用问题。该方法相对于LBG等常用的码书设计方法，能够减小码书的整体误差，极大地提高了矢量量化的有效性。同时也遵循了设计码书的最优准则，在理论上可设计最优矢量量化器。
【专利附图】
【附图说明】:
[0033]图1为矢量量化编码和解码流程框图；
[0034]图2为压缩编码流程框图；
[0035]图3为图像小波分解图；
[0036]图4为小波包分解及系数矢量构成图；
[0037]图5为基于改进的广义学习矢量量化流程图；
[0038]图6为图像解码流程框图。
【具体实施方式】:
[0039]本发明提出一种有效的基于小波包变化和矢量量化相结合的图像压缩编码方法，通过对小波变换后的低频系数进行低失真的标量量化来保证图像平滑部分的高质量恢复；通过对小波变换后的高频系数做进一步的小波包分解来有效地利用图像的高频细节部分，以详细分割各方向的高频细节，并合理利用分解各子频带系数的相关性构成了分类训练矢量集；在此基础上，根据设计最优化矢量量化码书的等失真准则，给出一种新的基于等失真准则的广义学习矢量量化(General Learning Vector Quantization, GLVQ)方法,并将其用于矢量量化码书的设计，有效地弥补了传统码书设计方法的不足。
[0040]因此，本发明公开的方法能够提高矢量量化码书的设计性能，在保证图像压缩编码质量的基础上，获得较大的图像压缩比。此外，本方法也能够更好地保留原图像的高频细节信息，并可以根据实际需求，控制图像的压缩比和重建质量。
[0041]下面结合具体实施例对本发明进行具体说明:
[0042]本发明的图像压缩编码过程如图2所示，具体步骤包括:
[0043]1.将图像进行3级小波变换，小波基函数采用DaubechieS9/7双正交滤波器。该小波基具有双正交性、紧支撑性，适合于表征平稳和非稳的各种二维数据。更为重要的是用该小波对信号进行二阶以上变换所得到的小波系数具有较好的稀疏性，适合于进行矢量量化。图像小波分解后的各频带分布如图3所示。其中LL3为低频子带，体现原图像中的平滑部分，虽然该部分数据量较小，但其能量占90%以上，对于重构图像十分重要。其他各部分均为高频子带，LH3，LH2，LHl分别代表原图像在不同分辨率下的水平方向细节；HL3，HL2，HLl分别代表原图像在不同分辨率下的垂直方向细节；HH3，HH2，HHl则分别代表原图像在不同分辨率下的垂直方向细节。
[0044]2.小波变换后低频部分与原始图像的相关系数非常接近，为了确保重构图像的质量，对其进行低失真的标量量化，利用DPCM预测去相关后再使用熵编码无损压缩方法中的算术编码进行压缩，从而得到低频码流。虽然标量量化压缩比小，但由于低频系数本身数据量低，因此对最终总体数据的压缩编码性能影响不大。
[0045]3.利用小波包技术对小波分解后的高频系数做进一步的分解，这能够实现对正交小波高频部分进行更加精细地划分，使分解序列在整个时域和频域都具有较高的时频分辨率和相同的带宽。这有利于更好地保证恢复图像中的细节信息，该方法对遥感卫星图像、医学图像等需要图像细节的应用领域更为重要。分解后各频带如图4所示，每个高频子带被进一步分解为各方向的细节部分。
[0046]4.对于分解后的多频带系数，在同一分辨率下不同方向高频子带具有相关性，体现为同一图像内容在不同方向上的细节描述；同时，在同一方向上的不同分辨率子带也具有相关性，体现为同一图像内容在不同尺度下的细节描述。因此，如何利用高频系数在空间位置和内容上的相关性，以及如何将小波包变换系数和矢量量化有效地组织在一起，是设计优化的矢量量化器的关键。本发明提出一种跨频带的分类矢量构成方法，在图4中，在同一方向不同频带HL3LL，H L2LL，HLlLL中分别取同一位置的I X I，2 X 2，4 X 4块构成21维矢量，依次取其他位置的矢量来构成该方向的训练矢量集。使用同样的方法，分别得到其他方向的训练矢量集。该方法是一种充分利用各级子带间的频率分辨率和空间分辨率带间相关性的矢量组合方法，按不同方向组织矢量量化的码书，每类矢量集具有不同的特性，可据此采用不同的量化策略。同时分类后也减少了优化码书的搜索空间，因此可以极大地提高设计码书的速度。
[0047]5.分类矢量集中的矢量具有较高的稀疏性，其中为零或接近零的元素较多，为了进一步减少训练集合、降低运算复杂度。本发明采用一种分类的阈值门限机制来筛选矢量，为每类矢量定义阈值Th。例如，对于HL3LL，HL2LL，HLlLL频带构成的矢量集Vh^设定门限阈值ThHm，该值通过对频带系数进行分析统计获得。将矢量中的第一个分量(HL3LL中的系_和阈值进行比较，若Ixl≤Th.，则直接删掉该矢量，而无需判断其他分量。这样做的原因是其他高分辨率子带为低分辨率子带进一步分解的结果，其能量将进一步减少，因此更低于门限。相反，若系数Ixl > Th11M则认为是重要细节信息。该方法可以通过对阈值的调整来控制训练矢量数，能够满足对不同压缩比和图像恢复质量的需求。
[0048]6.针对分类训练矢量集分别设计各自的码书。本发明提出一种基于等失真准则的广义矢量量化方法来训练生成码书。其原理是Ueda等学者证明的设计最优矢量量化器的等失真原则，它说明无论输入训练矢量如何分布，若要使平均失真最小，必须使码书中的每个码字在各自的VOTonoi分割区内对应的部分失真最小。方法的流程图如图5所示，设训练矢量集为
【权利要求】
1.一种基于小波包和矢量量化相结合的图像压缩编码方法，其特征在于:该方法按以下步骤: 步骤(1):对于待压缩编码图像，对其进行多级小波变换，分解成低频系数和高频系数； 步骤(2):对于小波分解后的低频系数采用标量量化方法，再使用熵编码得到低频码流； 步骤(3):对于小波分解后各频带的高频系数进一步采用小波包技术分别对其进行分解； 步骤(4):根据各级高频子带系数在空间和频率上的对应关系，构成分类训练矢量； 步骤(5):根据矢量能量阈值对训练矢量进行筛选； 步骤(6):对于筛选后的分类矢量，分别采用基于等失真准则的广义学习矢量量化方法来设计矢量量化器的码书； 步骤(7):对矢量量化的码流进一步使用熵编码进行无损压缩，从得到高频码流； 步骤(8):对由低频码流和高频码流构成的压缩数据，采用上述编码的逆过程进行解码，再进行重建得到恢复图像。
2.根据权利要求1所述的基于小波包和矢量量化相结合的图像压缩编码方法，其特征在于:在步骤(3)中，分别对小 波变换后的第I级，第2级和第3级频带三个方向的高频系数进行二维小波包变换，将各频带各方向的高频系数进一步分解，分解后的系数能够更精确的描述图像各方向细节信息。
3.根据权利要求1所述的基于小波包和矢量量化相结合的图像压缩编码方法，其特征在于:在步骤(4)中，各子带对应位置的系数代表原图像同一细节在不同分辨率下不同方向的描述，它们具有带间的相关性，根据这种空间位置和内容上的相关性来构造训练分类矢量，选择同一方向不同子带上同一位置的系数组成训练矢量，不同的级别频率子带选取不同数量的系数；再根据高频系数的方向性，将训练矢量进行分类，每类矢量用于设计各自的量化器码书，训练过程各自独立。
4.根据权利要求1所述的基于小波包和矢量量化相结合的图像压缩编码方法，其特征在于:在步骤(6)中，将设计矢量量化器最优码书的等失真准则引入到广义学习矢量量化方法中，加入一种基于等失真广义学习矢量量化的码书设计方法，其主要步骤包括: 1)初始化码书C(O)，初始码字可从输入训练序列矢量中随机选取N个矢量；初始化码书中各码字的部分失真01 (O);设最大迭代次数为T，允许误差门限值ε >0，初始学习率为α (O)； 2)对于各学习步骤t= 1，...，Τ，计算每步的变化学习率a (t)； ①根据每个码字所属区域内的训练矢量计算各码字的部分失真Di(t)； ②基于各码字的部分失真，计算改进的失真测度为; ③根据新的失真测度找出竞争获胜码字矢量； ④修正码书中的获胜码字和其他码字矢量； 3)调整码书中各码字分割区的所属训练矢量； 4)计算误差Et；5)若误差小于给定阈值或达到最大学习步数 ，则停止，并确定为最终的码书；否则转步骤(2)，进行重新计算。
【文档编号】H04N19/94GK103546759SQ201310526812
【公开日】2014年1月29日 申请日期:2013年10月29日 优先权日:2013年10月29日
【发明者】段勇, 李鹤婷, 于霞 申请人:沈阳工业大学