《数字图像处理-第八章.ppt》由会员分享,可在线阅读,更多相关《数字图像处理-第八章.ppt(127页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression第第8 8章章 图图 像像 压压 缩缩各行业,各领域,每天都有大量的数据要进行存储、处理和传输。如,美国国会图书馆电子图书馆。图像压缩就是解决这样的问题:减少表示数字图像时需要的数据量。图像压缩最早可追溯到35年前,当时是建立一种模拟的方法以减少传输所需的带宽,称为带宽压缩,后来计算机技术的发展,转到数字压缩。图像压缩是一种“开放技术”,现代图像Sensor分辨率不断提高,电视广播标准不断发展,图像压缩
2、成为一种基本技术。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression第第8 8章章 图图 像像 压压 缩缩 本章主要讲述图像压缩一解压缩的理论和实践。先讲基本原理,再讲应用 。压缩技术:信息保持编码,档案保存信息有损编码:广播电视视频会议传真。最后介绍标准。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.
3、1 8.1 基础基础8.1 基础 数据压缩:减少表示给定信息量所需的数据。同一件事情可以有不同的描述版本,至少一个可能包含是不必要的数据。数据冗余是图像压缩的主要问题,例如:表示一个相同信息的两个数据集合中,携载的信息的单元数量分别为相对数据冗余 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1 基础(续)例如 第1集合有10信息单元,第2集合有1信息单元三种基本数据冗余 编码冗余像素词冗余心理视觉冗余数数 字字 图图
4、像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.1 编码冗余利用直方图的原理,减少数据量。设直方图、灰度都归一化到表示 的比特数为 ,则每像素所需平均比特为:图像,编码所需比特数MN Lavg,用m比特二进制编码,则表示灰度的可减少到m如8位。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8
5、.1 8.1 基础基础8.1.1 编码冗余例8.1说明变长编码实现数据压缩数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.1 编码冗余 一般用二进制编码时,冗余总存在。总之,利用少比待,表示出现概率大的灰度数,实现数据压缩,亦叫变长编码。有关编码我们在本章以后部分会详加介绍也可以参考一些信息编码类的书数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapte
6、r 8Image CompressionImage Compression8.1 8.1 基础基础8.1.2 像素间冗余 先看右边两幅图a,b,c,d,而e,f为相关系数45间隔联系。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.2 像素间冗余 如上页图示:说明:像素间存在相关性(尤其是相邻像素),任何给定像素可据其相邻像素,适当的预测而得,因此,单个像素所携载信息相对较少。像素间冗余:空间冗余,几何冗余,帧间冗余
7、数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.2 像素间冗余 例8.2为一个行程编码实现压缩的例子:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.3 心理视觉冗余 人眼感觉到的区域亮度不仅取决于反对光,还和其它因素有关。如马赫带(亮度不变,
8、感觉变,P32)其原因是眼睛对各种信息感受的灵敏度不同。那些不重要的信息叫心理视觉冗余,而这种冗余可以在不削减图像感知质量的情况下消除。人眼感知图像时,不是分析每个像素值,通常找特征,如边缘,纹理合并,成群大脑与书籍相联系图像理解,也叫“量化”(导致一定信息丢失)从一个范围很广的值集合有限输出值。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.3 心理视觉冗余 例8.3通过量化进行压缩。数数 字字 图图 像像 处处 理
9、理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.3 心理视觉冗余 IGS量化方法:为减少颗粒状纹路用相邻像素灰度的低位产生随机数,加到当前像素。IGS量化过程:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.4 保真度准则评估的两类准则:(1)客观;(2)主观客观:表示解压缩
10、图像主观:典型观察者+典型图像数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression8.1 8.1 基础基础8.1.4 保真度准则表8.3为绝对等级。可以并排对比,非常恶劣非常好 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.2 8.2 图像压缩模型图像压缩模型8.2 图像压缩模型。常用图像压缩系统模型。
11、数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.2 8.2 图像压缩模型图像压缩模型8.2.1 信源编码器和信源解码器信源解码器 图中信源编码目的是消除输入冗余,信道编码是增强信源编码器抗噪性。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.2 8.2 图像压缩模型图像压缩模型8.2.2 信道编
12、码器和解码器 向信源编码数据中插入冗余数据,减少信道噪声的影响。最有用的信道编码技术是,R.W.Hamming 思想:向被编码数据中加入足够位数,以确保有效的码字间变化的位数最少。如:将3位冗余码加到4位码,使得任意2个有效码字间距离为3,则1位错误可检出来。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.2 8.2 图像压缩模型图像压缩模型8.2.2 信道编码器和解码器(续)一位错误由一个非0奇偶校码字 给出。如果结果 ,解码器只要翻转码字中
13、由奇偶校验字拨出的比特位的位置(的码),然后以 解码即可。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3 信息论要素 显示一幅图像究竟需多大的数据量?8.3.1 测量信息 信息的产生可以被模拟为一个概率过程,发生概率为P(E)的随机事件E包含。的信息,也叫自信息,底数决定信息的单位,通常2,则为比特。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter
14、8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.2 信息信道离散信息系统的数学模型:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.2 信息信道(信源)信源,根据一个有限或无限可数的符号集生成一个信息的符号序列。信源符号集 叫信源字母表A 用(A,Z)描述信源。如果产生k个信源符号,据大数定理,将平均被输出 次,则根据k
15、输出得到的平均自信息:每个信源输出的平均信息:以信源的熵(或不确定度),是观察单个信源输出得到的平均信息量。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.2 信息信道(信道)为了得到具有正向信道传输矩阵Q的信道容量,先计算信源的熵(输出 前提下)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image Compression
16、Image Compression 8.3 8.3 信息论要素信息论要素8.3.2 信息信道(信道)(续)对所有 的期望值:代入整理得:叫z关于v的平均条件自信息(条件熵)其含义是:一个信源符号产生一个输出,在观察到该输出后,一个信源符号的平均信息(亦:条件熵)叫z和v的互信息。表示v为z所提供的信息量。也表示观察单个输出符号时接收到的平均信息。(v可看作编码的结果)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要
17、素8.3.2 信息信道(信道)(续)观察信道单一输出时接取到的平均信息就是信源符号概率向量z和信道矩阵Q的函数。叫信道容量:能够可靠传送信息的最大传送率,不取决于信源输入概率Q,而取决于信道的条件概率函数(Q)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.2 信息信道(信道)(续)例8.6 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Cha
18、pter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本编码定理 (扩展)的通信系统模型:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本编码定理无噪声编码定理。当信道和通信系统中不存在噪声时,香农第一定理定义了每个信源符号的最小平均码表。概念:零记忆信源;一个具有有限集合(A,Z)和在统计上独立的符号源的信息源。叫单一
19、符号或非扩充信源的n次扩充,其熵为对应单符号信源熵的n倍。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本编码定理用一整长为 的码为 进行编码 码表为超出 自信息的最小整数。可导出:(平均值)表示对应于非扩充信源的n次扩充编码的平均字长 即:(香农第一定理)说明:对无限扩充的信源进行编码可使 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter
20、8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本编码定理(续上)编码效率例8.7 噪声编码定量 当信道有噪声或易出错,我们关注的重点从尽可能紧凑的信息表示的编码尽可能稳定的通信。例8.8 噪声二值信道数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本编码定理 假设BSC有出错概率 ,我们可以将每个信息或二进
21、制符号重复几遍:0000,1111 输出端按“多数确定”的方法来解码(那个多算那个)单个错误概率 2个错 3个错 当无错误或仅一个错时,可正确解码(此时错误概率0.0003)通过加大重复传输次数,可使传输中的总误差尽可能地小。当信源包含K个符号时,可用对每个符号传r次来对信源的n阶扩展进行编码,有关键点:仅选择 个可能码,序列中的S个码字作为有效码。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.3 8.3 信息论要素信息论要素8.3.3 基本
22、编码定理 零记忆信源的信息产生率(信息单位/符号)等于信源的熵 ,其n阶扩展的信息产生率为 。如果用来编码的S个(或 )个有效码字相同概率的,可达最大信息率。尺寸为S(或 ),块长为r的码具有信息率:叫香农第二定理。说明:对信息的Rcc取RAC距离d=cc(取小的,取下边)对多数图像来说,RAC距离的概率分布是不均匀的,所以要用合适的变长码来对RAC距离进行编码。采用类似于B1码对RAC距离编码。最短的距离用最短的码字来编,而其它所有距离用如下方法编:用第1个前缀表示最短的RAC距离 第2个前缀将d赋给某个距离范围 再如上d-范围下限范围下限的二进制表示 如在上图中ec=8 cc=4 则RAC
23、码 1100 0 11 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩边界跟踪和编码边界跟踪和编码 通过跟踪二值图中的区域边界并进选编码,也可达到对常数区编码的目的,预测微分量化就是一种面向扫描线的边界跟踪方法(叫PDQ)。各参数的意义见图8.18,获得一系理差异对()()就叫双编码(DDC)加上,前一行和后一行的对应差异对连结起来,解码器可根据这些信息正确地确定边界在图中的位置。数数 字字 图图 像像 处处
24、理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩边界跟踪和编码边界跟踪和编码 PDQ和DDC编码的最后1步都是对 和 或 和 及表示新边界开始和旧边界结束的行和列的坐标用变长码编码。为了避免对每个表示新边界开始和边界结束的行和列坐标都编码,需用1个独特的码代表与区域不相交的扫描线。二值压缩技术的比较,见下页表8.8用游程编码的熵的一阶估计作为变长码所能达到的压缩性能(H)的近似。结论:(1)各种方法得到的码率都比一阶熵估计要小,说明各种方法均
25、能消除一定的像素冗余,其中游程码最好。(2)灰度编码能得到的改近约为1bit/pixel(相对于二值编码)。(3)所有5种方法压缩率都在12之间,这主要是因为它们对低位面的压缩的最差,“”代表此处数据膨胀。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩边界跟踪和编码边界跟踪和编码表8.8数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8
26、Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩预测编码无损有损数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩8.4.4 无损预测编码 其基本原理是通过对每个像素新增的信息进行提取和编码,来消除在空间上较为接近像素之间的冗余信息。一个像素的新增信息被定义为此像素实际值和预测值之间的差异。无损预测编码系统的基本组成如下:数数 字字 图图 像像 处处 理理西
27、安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩8.4.4 无损预测编码用变长代码进行 编码压缩数据流 解码时:有局部的、全局的、自适应的方法 多数根据:m:预测器的阶 round:4舍5入 :预测系数如一维线性预测编码 上式中不能对每一行的前m个像素进行求解,故这些像素必须用其他方法进行编码(如霍夫曼编码)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image Co
28、mpressionImage Compression 8.4 8.4 无误差压缩无误差压缩8.4.4 无损预测编码 例8.15数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.4 8.4 无误差压缩无误差压缩8.4.4 无损预测编码叫前像素预测器,相应的编码过程称为差分编码或前像素编码。压缩比一般为2:1左右 预测误差的概率密度函数一般用0均值不相差拉普拉斯概率密度函数表示数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院
29、电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5 有损压缩 在重构的精度上让步,而换取压缩能力的增加,有损压缩编码技术可达100:1压缩率,并且做到10:150:1之间,图像无本质区别(单色图像),与上节不同:在于是否存在量化器。8.5.1 有损预测编码有损预测编码模型 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损
30、压缩有损压缩8.5.1 有损预测编码 量化器将预测误差映射成有限范围的输出,表为 是反馈环的输入,表示为过去预测和与其对应的量化误差的函数(和),解码器的输出也用上式给出。例8.16德尔塔调制 失真对所有预测编码方法都是很普遍的,这种失真的严重性取决于所使用的量化和预测方法之间互相作用。一般,预测器和量化器在设计中独立进行的。预测器在设计中认为量化器没有误差。而量化器在设计中则需要最小化自身的误差。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8
31、.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 例8.16德尔塔调制数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 最佳预测器 将编码器的均方预测误差最小:基于这些条件的预测编码方法称为差分脉冲编码调制(DPCM)就是选择m个预测系数以使上式最小化,对上式每个子数求解,等于0,并设fn均值为0,方差为 时,得到:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学
32、院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 最佳预测器 是下列mm自相关矩阵的逆矩阵可推导出,最优时,预测误差的方差为:上述表达式形式简单,但要获得R,r所需的自相关计算需很困难。实际中逐幅图像计算预测系数的方法很少用,一般都假设1个简单的图像模型,并将其对应的自相关代入R,r,表达式中以计算全局(所有图)系数。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image Compr
33、essionImage Compression 8.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 最佳预测器 例如1个2-D马尔可夫源具有可分离自相差函数 其中 是图像的水平和垂直相关系数,用4阶线性预测器:来预测;则最优系数为:对 赋予不同的值,可得到如下4个预测器:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 最佳预测器(续上):用上述四个预测器对下副图象进行预测误差。数数 字
34、字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.1 有损预测编码 最佳预测器(续上):结论:随预测阶数增加,误差明显减少。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩最佳量化最佳量化 这是一个典型的量化函数,这个阶梯状函数 是s的奇函数
35、,该函数可用在第1象限的 个Si和ti完全描述。Si:判别电平ti:重建电平 在(Si,Si+1)内的sti 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩最佳量化最佳量化 根据以上定义,量化器的设计就是在给定优化准则和输入概率密度函数p(s)条件下,选最优的Si和ti 优化准则可以是统计的或是心理视觉准则。如果用最小均方量化误差 作准则,且p(s)为偶函数,则最小误差条件为:可看出:重建电平(ti)是所给定判别区
36、间的p(s)曲线下面积的重心,而判别值(Si)是2个重建值的中值。这种量化器叫L级Lloyd-Max(劳埃德马克斯)量化器.数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩最佳量化最佳量化 表8.10给出了单位方差的拉普拉斯概率密度函数的Lloyd-Max量化器(2、4、8级)这三个量化器分别给出1、2、3bit/pixel固定输出率。对于重建值方差 的情况,用表中给的数据乘以它们的概率密度函数的标准差即可。表中步长
37、 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩最佳量化最佳量化 尽管Lloyd-Max量化器和最优均匀量化器不是自适应的,但如果根据图像局部性质调节量化值,也能提高效率。理论上讲,可以较细地量化缓慢变化区域,而较粗地量化快速变化区域,这可同时减少颗粒噪声和斜率过载,且码率增加很少。表8.11给出了用预测器和量化器的不同组合对同一幅图像编码所得的均方根误差。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工
38、程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码 预测编码技术直接对象素在图像空间进行操作,故也称为空域方法。本节讨论基于图像变换的编码方法(也叫频域方法)。在变换偏码中我们用可逆的线性变换将图像映射成1组变换系数,然后将这些系数量化和编码。对多数自然界图像变换得到的系数值都很小,这些系数可较粗地量化,甚至完全忽略掉而只产生很少的失真。典型的变换编码系统框图如下:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系C
39、hapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码 上述任一个或全部变换编码步骤都可以根据图像局部内容调整时,叫自适应变换编码。否则(全固定)叫非自适应变换编码。变换的选择 考虑大小为NN的图像f(x,y),正向离散变换T(u,v)为:这2个也叫基函数或基图像。对h也有相同的结论。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5
40、有损压缩有损压缩8.5.2 变换编码 是z的二进制表达式的第k位(从右往左)例如:如果m=3,z=6(110)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码图8.29显示了Walsh-Hadamard的基函数(图像)。N=4时,每个块包括44=16个元素,例如,为得左上角的块,u=v=0,x,y=0,1,2,3时的g(x,y,0,0)值即为图像值。W-H所有核值均为+1或-1。数数 字字 图图
41、像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩离散余弦函数变换(离散余弦函数变换(DCT)图8.30为其基图像。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码 例8.19分别用DFT,WHT和DCT近纤变换编码。注意,压缩并不是在变换
42、步骤中取得的,而是在量化变换系数中取得的。对于一个给定的编码应用,如何选择变换取决于容许的重建误差和计算要求。变换具有将图像能量或信息集中于系数的能力,均方重建误差与所用变换的该性质直接相关。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码 一幅nn图像f(x,y)可表示成它的2-D变换T(u,v)的函数。这里用n替换N以表示子图像。上式表成矩阵形式:F是由所有f(x,y)组成的nn矩阵现在定义一
43、个变换系数模板函数:数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码(续上)为F的近似。均方误差:总的均方近似误差就是所有截除的变换系数的方差之和。一个能把最多的信息集中在最少数上的变换所产生的重建误差最小。上式最后一步是依据基函数正交性和F中的像素是由0均值和已知方差的随机过程产生的。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Ch
44、apter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码(续)一幅NN图像的 个子图像的均方误差是相同的,所以NN图像的均方误差等于子图像的均方误差。不同变换的信息集中能力不同。DCT比DFT和WHT且有更强的信息集中能力。但理论上说,KLT是所有变换中信息集中最优的变换。但KLT与输入图像数据有关。计算量极大。故不实用。非正弦类变换,如WHT计算容易,但正弦类变换(如DFT、DCT)更接近KLT的信息集中能力。用于DCT的信息集中力和计算复杂性综合得比较好而得到较多的应用,它具有使“块效应”最小的变换(相邻子图像
45、边界可见)。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 子图像尺寸选择子图像尺寸选择 子图像尺寸也是影响变换编码误差和计算复杂度的一个重要因素。多数情况下,我们将图像分成尺寸满足以下两个条件的子图像:相邻子图像之间的相关(冗余)减少到某个可接受的水平。子图像的长和宽都是2的整数次幂。第2个条件主要是为了简化对子图像变换的计算。最常用的尺寸是88和1616。数数 字字 图图 像像 处处
46、 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 子图像尺寸选择子图像尺寸选择 例8.20为变换编码中子图像尺寸的影响。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 比特分配比特分配 与 相对应的截断误差和两
47、个因素有关 截除的变换系数的数量和相对重要性 用来表示所保留系数的精度,大多数变换编码系统中,保留的系数是根据下列两个准则之一来选择确定:1)最大方差准则,叫区域编码2)最大幅度准则,叫阈值编码。对整个变换子图像的系数截断,量化和编码过程叫比特分配。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 比特分配比特分配 1区域编码:其基础是信息论中的不确定性原理,即:具有最大方差的变换系数带有
48、最多的图像信息,它们应当保留在编码过程中。方差计算:从个变换后的子图像获得基于某个图像模型(如马尔科夫自相关函数)算子。典型的区域模板如图8.36(a)所示保留系统的量化和编码:见8.36(b)数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 比特分配比特分配 区域编码:图8.36 数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8
49、Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 比特分配比特分配 两种策略:给各系数分配相同数量的比特,常将系数用它们的均方差归一化,然后均匀量化给不同系数分配几种固定数量的比特。此时,要对每个系数设计一个量化器(如Lloyd-Max量化器)。为构造量化器,将零阶或直流分量系数模型化为一个瑞利函数,其它系数模型化为拉普拉斯或高斯密度函数。因为每个系数都是子图像中象素的线性组合,所以根据中心极限定理,随子力像尺寸的增大,系数趋于高斯分配,但直流不是这样(0)率失真定理指出,一个方差为 的高斯随机
50、变量用至少 bit来表达(D:均方误差)故分配给每个量化器的级数正比于 。数数 字字 图图 像像 处处 理理西安邮电学院电子与信息工程系西安邮电学院电子与信息工程系Chapter 8Chapter 8Image CompressionImage Compression 8.5 8.5 有损压缩有损压缩8.5.2 变换编码变换编码 比特分配比特分配 2阈值编码 本质上是自适应的,为各个子图像保留的变换子数的位置随子图像的不同而不同。事实上,由于阈值编码计算简单,较常用。对任意子图像,值最大的变换系数对重建子图像的质量贡献最大。由于最大系数的位置随子图像变化,故 的元素常重新排列,构成一个1-D行