基于傅立叶级数的简易有损图像压缩算法.docVIP

基于傅立叶级数的简易有损图像压缩算法 付浩 P概览 数据压缩在信息技术中有着广泛的用途。一般而言，在保证一定的图像质量前提下，使用有损压缩算法压缩数字图像可以获得比较令人满意压缩比。本文通过对数字图像信息进行恰当的处理，利用傅立叶级数近似拟合图像，从而实现了一种简易的有损图像压缩算法。 预备知识 数字图像在计算中的表示方式 我们都知道，任何可见光的色彩都是由红绿蓝(RGB)三原色按一定比例构成的。在计算机中，颜色的一种常见表示方式就是用三个字节分别表示RGB三种色彩，其中每个字节取值范围为且为整数。如此便可以表示出224色彩，已经超过了人眼所能分辨的颜色数，因此使用这种方式就已经足够了。 记任一颜色为向量，一幅宽w个像素高h个像素的图像便可用个向量表示。为了方便处理，本文只取向量的一个分量进行处理。即把一幅彩色图像分解为3幅灰度图像分别处理，且每幅灰度图像都用含项的整数数列表示。 由于人的视觉特点，一般人难以分辨压缩还原后的图像颜色与原图像的微小差别。有损图像压缩利用这一特点实现比无损压缩算法更好的压缩比。 傅立叶级数 有关傅立叶级数的详细论述请参阅有关书籍，这里仅不加证明地叙述本文所要用到相关性质。 傅立叶系数 设函数在数轴上处处连续且在有限区间上逐段光滑，且周期为2T。令： 根据迪里赫勒定理，傅立叶级数绝对一致收敛于，即 余弦级数 设为定义在有限区间上的连续函数，对其进行偶性开拓，即构造一个新函数，根据奇函数和偶函数在对称区间上积分的性质可知的傅立叶系数为： 于是，在区间上有： 贝塞尔不等式 根据贝塞尔不等式的变形，取为定义在有限区间上的连续函数，利用2的结论有，定义平方平均偏差为： 本文将用评估图像压缩后的失真程度。 理论推导 图像信息预处理 设输入为一幅宽w0像素高h0像素的灰度图像（灰度值在内），将其分解为若干矩形小块。对于每个小块，设其高为w宽为h，按照先从左到右后从上到下的顺序将其灰度值依次记为数列，其中 为简单起见，构造定义在闭区间上的一次分段函数： 易见为分段光滑的连续函数，且 计算傅立叶系数 根据的定义，下面将用有限项的傅立叶函数列近似拟合： 故有 该和式第一项等于0 由可知第二项等于 因此 失真度评估 定义 本文将使用计算恢复被压缩的图像。如果p过小，显然图像将会严重失真；如果p过大，将会占用过多的储存空间而无法达到压缩的目的。 考虑两种极端情况。若被压缩的图像小块为纯色的色块，显然只需取即可还原图像；若被压缩的图像小块相邻像素变化极为“剧烈”，显然需要取较大的p才能保证图像还原的质量。 因此，为了选择一个合适p，有必要对压缩还原后的图像失真度进行评估。 可以证明，定义的平方平均偏差等于与的平方平均偏差，即，因此可以使用来评估压缩后图像与原图像的偏差。 其中 因此 易见关于p单调递减。设q0为画质阀值，即对每一个图像小块，p的取值都应恰好使得 q可直接由用户指定，也可根据图像小块的大小以及用户对画质的要求等因素来决定，这个问题较为复杂，暂不讨论。 压缩数据储存 一个字节最多只能表示不大于255的非负整数。根据对输入图像信息的限制，对任意，都有 由得： 考虑把每一个储存在一个字节中。为了高效的利用储存空间，需要对做适当的放缩，可令 易见，从而实际储存的压缩数据为 解压还原时按照做相应计算即可。 算法描述 压缩部分 由用户指定图像小块的大小w和h，以及画质阀值q，令 输入一幅彩色或灰度图像，若为彩色图像，则将其分解为3幅灰度图像分别处理。 对某一未处理小块，按照计算，令，按照计算 若则跳转到7 令，按照计算，令 跳转到4 按照处理数据，并将以及p写入磁盘或内存缓冲区中 如果还有没处理的图像小块则跳转到3，否则算法结束。 粗略估计时间复杂度为，其中为所有p的平均值，由q决定。 解压部分 从文件读入、各个图像小块对应的p以及等数据 对某一未处理小块，按照逆向计算原来近似值 按计算小块内所有像素的灰度值，其中第i个灰度值为，将其输出到屏幕或其他文件 如果所有小块都处理完了，则算法结束，否则跳转到2 粗略估计时间复杂度为 总结 本文仅描述了算法的大致框架，且只给出了简单的理论推导，仅做启发思维之用，若要投入使用还需要许多细致的工作。事实上，此算法稍加改造亦可用于音频压缩。 参考书目 高等数学导论（下册） 中国科学技术大学出版社