《数字图像处理整数小波》由会员分享，可在线阅读，更多相关《数字图像处理整数小波（38页珍藏版）》请在金锄头文库上搜索。
1、张立保北京师范大学信息学院第2页第12讲第第12章章 图象编码方法图象编码方法12.1预测编码预测编码12.2变换编码变换编码12.3子带编码子带编码 12.4小波变换编码小波变换编码 章毓晋 (TH-EE-IE)第3页第12讲12.1 预测编码预测编码空域方法，消除象素间的冗余象素间的相关性使得预测成为可能仅提取每个象素中的新信息并对它们编码12.1.1无损预测编码信息保存型12.1.2有损预测编码信息损失型章毓晋 (TH-EE-IE)第4页第12讲12.1.1 无损预测编码无损预测编码系统无损预测编码系统编码器 + 解码器（有相同的预测器）章毓晋 (TH-EE-IE)第5页第12讲12.1
2、.1 无损预测编码无损预测编码过程无损预测编码过程输入序列： fn (n = 1, 2, )预测输出： （舍入成整数）预测误差：误差编码：在符号编码器中用变长码编误差解压序列：哪里取得了压缩？（消除了象素间冗余）（消除了象素间冗余） 章毓晋 (TH-EE-IE)第6页第12讲12.1.1 无损预测编码 m阶线性预测： 1-D线性预测： 一阶1-D线性预测：预测误差的概率密度函数：章毓晋 (TH-EE-IE)第7页第12讲12.1.2 有损预测编码1、有损预测编码系统、有损预测编码系统增加了1个量化器，预测器放在1个反馈环中 章毓晋 (TH-EE-IE)第8页第12讲12.1.2 有损预测编码1
3、、有损预测编码系统、有损预测编码系统输入序列： fn (n = 1, 2, )量化输出：预测输入：解压序列：编码误差：哪里又又取得了压缩？（量化，减少了（量化，减少了 心理视觉冗余）心理视觉冗余） 章毓晋 (TH-EE-IE)第9页第12讲12.1.2 有损预测编码1、有损预测编码系统、有损预测编码系统德尔塔调制（DM）预测器量化器预测系数 a 1，常数 c 0 DM方法得到的码率是1比特/象素 章毓晋 (TH-EE-IE)第10页第12讲12.1.2 有损预测编码1、有损预测编码系统、有损预测编码系统DM编码中的失真示例 章毓晋 (TH-EE-IE)第11页第12讲12.1.2 有损预测编码
4、2、最优预测、最优预测 最小化编码器的均方预测误差差值脉冲码调制法 （DPCM） 章毓晋 (TH-EE-IE)第12页第12讲12.1.2 有损预测编码2、最优预测、最优预测4阶线性预测器 章毓晋 (TH-EE-IE)第13页第12讲12.1.2 有损预测编码3、最优量化、最优量化判别重建章毓晋 (TH-EE-IE)第14页第12讲12.1.2 有损预测编码3、最优量化、最优量化最小均方量化误差 重建电平是p(s)曲线下面积的重心 判别值为2个 重建值的中值 q(s)奇函数输入概率密度函数，偶函数章毓晋 (TH-EE-IE)第15页第12讲12.2 变换编码变换编码频域方法，非信息保持型 12
5、.2.1变换编码系统12.2.2子图象尺寸选择12.2.3变换选择12.2.4比特分配章毓晋 (TH-EE-IE)第16页第12讲12.2.1 变换编码系统图象分解：减少变换的计算复杂度图象变换：解除每个子图象内部象素之间的 相关性，或者说将尽可能多的信息集中到尽可能少的变换系数上压缩不是在变换中而是在量化变换系数时取得的压缩不是在变换中而是在量化变换系数时取得的章毓晋 (TH-EE-IE)第17页第12讲12.2.2 子图象尺寸选择影响变换编码误差和计算复杂度（压缩量和计算复杂度都随子图象尺寸的增加而增加 ）两个条件： 相邻子图象之间的相关（冗余）减少到某个可接受的水平； 子图象的长和宽都是
6、2的整数次幂最常用的子图象尺寸：8 8和16 16 章毓晋 (TH-EE-IE)第18页第12讲12.2.2 子图象尺寸选择 变换编码重建误差与子图象尺寸的关系章毓晋 (TH-EE-IE)第19页第12讲12.2.3 变换选择u一个能把最多的信息集中到最少的系数上去的变换所产生的重建误差最小 u不同变换的信息集中能力不同 lKLT最优，但计算量非常大（依赖于图象）l正弦类变换（如DFT和DCT）较优l非正弦类变换（如WHT）实现简单l小波变换计算快且有局部性质（不需分解）章毓晋 (TH-EE-IE)第20页第12讲12.2.3 变换选择第5章介绍了：DFT，WHT，DCT第10章介绍了：KLT
7、信息集中能力： KLT DCT DFT WHT所需计算量：KLT DCT DFT WHTDCT是较好的（综合）选择是较好的（综合）选择章毓晋 (TH-EE-IE)第21页第12讲12.2.4 比特分配比特分配：对变换子图象的系数截断、量化和编码的全过程截断误差 截除的变换系数的数量和相对重要性 用来表示所保留系数的精度（量化）保留系数的2个准则 最大方差准则，称为分区编码 最大幅度准则，称为阈值编码章毓晋 (TH-EE-IE)第22页第12讲12.2.4 比特分配1、分区编码、分区编码具有最大方差的变换系数带有最多的图象信息事先确定模板，保留一定的系数章毓晋 (TH-EE-IE)第23页第12
8、讲12.2.4 比特分配2、阈值编码、阈值编码根据子图象特性自适应选择保留系数将系数排队，与阈值比较确定去舍63625857494836350000000061595650473734210000000060555146383322200000001054524539322319100000000053444031241811900000001434130251712830000001142292616137420000111128271514651000001011章毓晋 (TH-EE-IE)第24页第12讲12.2.4 比特分配2、阈值编码、阈值编码随子图象不同而保留不同位置的变换系数常用
9、三种对变换子图象取阈值（即产生式(12.2.4)所示模板函数）的方法：(1) 对所有子图象用一个全局阈值压缩的程度随（不同）图象而异 ？ (2) 对各个子图象分别用不同的阈值舍去同数量系数，码率是个常数 章毓晋 (TH-EE-IE)第25页第12讲12.2.4 比特分配2、阈值编码、阈值编码(3) 根据子图象中系数的位置选取阈值将取阈值和量化结合起来 章毓晋 (TH-EE-IE)第26页第12讲12.3 子带编码子带编码子带子带：由图象分解得到的一系列带限分量的集合将它们重新组合起来可以无失真地重建原始图象将图象分解为子带后进行编码的主要好处是(1)不同子带内的图象能量和统计特性不同，可以采取
10、不同的变长码甚至不同的编码方法分别进行编码，提高编码效率(2)通过频率分解，减少或消除了不同频率之间的相关性，有利于减少图象数据的冗余(3)量化等操作可在各子带内分别进行，避免了互相干扰和噪声扩散章毓晋 (TH-EE-IE)第27页第12讲12.3 子带编码子带编码1、子带分解、子带分解输出序列x*(n)是先通过将x(n)用分析滤波器h0(n)和h1(n)分解为y0(n)和y1(n)，再借助合成滤波器g0(n)和g1(n)进行重建得到的 h0(n)和h1(n)都是半带（half-band）滤波器，其中h0(n)对应一个低通滤波器，其输出是x(n)的近似部分， h1(n)对应一个高通滤波器，其输
11、出是x(n)的细节部分 章毓晋 (TH-EE-IE)第28页第12讲12.3 子带编码子带编码1、子带分解、子带分解序列x(n)，n = 0, 1, 2, 的Z-变换在时域中以2为因子的抽样 在时域中以2为因子的内插 章毓晋 (TH-EE-IE)第29页第12讲12.3 子带编码子带编码1、子带分解、子带分解序列x(n)先被抽样再被内插，得到x*(n) 根据Z-变换来考虑子带编码和解码系统系统的输出为：其中第2项（依赖于z）表达了由于抽样和内插过程而引入的混叠 章毓晋 (TH-EE-IE)第30页第12讲12.3 子带编码子带编码2、子带重建、子带重建重建无失真，x*(n) x(n)和X*(z
12、) X(z) 结合进一个矩阵表达式 章毓晋 (TH-EE-IE)第31页第12讲12.3 子带编码子带编码2、子带重建、子带重建假设Hm(z)是非奇异的（行列式不为零） 分析滤波器和合成滤波器是交叉调制的 对有限冲击响应（finite impulse response， FIR）滤波器，detHm(z) = az(2k+1) 章毓晋 (TH-EE-IE)第32页第12讲12.4 小波变换编码小波变换编码12.4.1小波变换编码系统12.4.2基于提升小波的编码章毓晋 (TH-EE-IE)第33页第12讲12.4.1 小波变换编码系统小波变换编码也是一种变换编码方式与采用正交变换（如DCT）的编
13、解码系统不同， 小波变换编解码系统中没有图象分块的模块小波变换的计算效率很高，且本质上具有局部性小波变换编码不会产生使用DCT变换在高压缩比时的块效应 章毓晋 (TH-EE-IE)第34页第12讲12.4.1 小波变换编码系统小波变换编码需考虑的几个因素小波变换编码需考虑的几个因素1.小波选择小波选择如：双正交小波2.分解层数选择分解层数选择影响小波编码计算的复杂度和重建误差3.量化设计量化设计对小波编码压缩和重建误差影响最大在不同尺度间调整量化间隔 章毓晋 (TH-EE-IE)第35页第12讲12.4.2 基于提升小波的编码可以在当前位置实现整数到整数的变换，运算速度快且节约内存。它包括三个
14、步骤：1.分裂（分裂（split）将图象数据分解成偶数部分和奇数部分 S uj(x, y) := uj1,k(x, y), vj1,k(x, y) 章毓晋 (TH-EE-IE)第36页第12讲2.预测（预测（predict） 保持偶数部分不变并用偶数部分来预测奇数部分，然后用奇数部分与预测值的差（称为细节系数）替代奇数部分vj1,k(x, y) := vj1,k(x, y) P uj 1,k(x, y) 12.4.2 基于提升小波的编码章毓晋 (TH-EE-IE)第37页第12讲3.更新（更新（update） 构造一个作用于细节函数的算子U，并叠加到偶数部分上以获得近似图象，这里要保持原始图象的一些特性 uj 1,k(x, y) := uj 1,k(x, y) +U vj 1,k(x, y) 12.4.2 基于提升小波的编码章毓晋 (TH-EE-IE)第38页第12讲12.4.2 基于提升小波的编码重建过程重建过程三个运算：三个运算：（M 合并） (1)uj1,k(x, y) := uj1,k(x, y) U vj1,k(x, y)(2)vj1,k(x, y) := vj1,k(x, y) + P uj1,k(x, y)(3)uj,k(x, y) := M uj1,k(x, y), vj1,k(x, y)章毓晋 (TH-EE-IE)
