数字图像处理-实验二..docVIP

图像压缩实验 实验目的及要求： 学会利用一种算法进行数据压缩，增进对图像压缩理论的感性认识。 二、实验内容： 用MATLAB编写二个子函数，一个函数用于对图像进行压缩，并将压缩后的图像以文件的方式保存进来，同时还要计算压缩比：压缩前的文件尺寸与压缩后的文件尺寸之比；另一个函数用于进行解压缩，它读取经压缩程序压缩后的图像文件，把图像恢复出来。显然恢复的图像与原图像之间会存在偏差，请计算峰值信噪比。实验采用512x512的LENA图像，该图像每象素由8比特表示。算法如下： 压缩算法：先把图像按16x16进行分块，对每一小块进行二维沃尔什变换，对每一个变换后的块按表1进行比特分配，对变换后的系数进行量化。对表1中每个系数，其取值的范围可以这样确定：对某个位置上的系数，将变换后的所有16x16的小块中同一位置的系数值取出来，再进行统计，得出该位置上的系数值的上下限。将系数值的上下限也存储到最后的压缩文件中，用于解压缩时从编码值恢复系数值。由于量化后表达每个象素所需要的比特数变少了，在存储数据时，把不足8比特的像素与相邻的像素组合起来形成一个字节。压缩后数据的存放格式自行确定。 表1（空着的部分全为0，表示把这些地方的数据丢掉不用） 8 8 8 8 8 8 8 8 4 4 4 4 4 4 4 4 8 8 4 4 4 4 4 4 4 4 2 2 2 2 2 2 8 4 4 4 4 4 4 4 4 4 2 2 2 2 2 2 8 4 4 4 2 2 2 2 2 2 2 2 8 4 4 2 8 4 4 2 8 4 4 2 8 4 4 2 4 4 4 2 4 4 4 2 4 2 2 2 4 2 2 2 4 2 2 4 2 2 4 2 2 4 2 2 解压缩算法：实现与压缩相反的一个过程。 三、实验步骤： 1 ．创建一个GUI用户界面窗口（如图1所示） 2 ．编辑步骤1创建的GUI用户界面窗口所对应的m文件（如图2所示） 图1 图2 3 ．输入如下程序，实现图像压缩功能 新建walsh.m文件，输入以下代码并保存： function [k]=walsh(x) b=[0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1, 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1, 0 0 0 0 1 1 1 1 0 0 0 0 1 1 1 1, 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1]; n=0; g=zeros(16,16); for u=1:16, for v=1:16, for i=1:4, m=b(i,v)*b(5-i,u); n=n+m; end g(u,v)=(-1)^n; n=0; end end k=g*x*g./256; 新建 dewalsh.m文件，输入以下代码并保存： function [k]=dewalsh(x) b=[0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1, 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1, 0 0 0 0 1 1 1 1 0 0 0 0 1 1 1 1, 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1]; n=0; g=zeros(16,16); for u=1:16, for v=1:16, for i=1:4, m=b(i,v)*b(5-i,u); n=n+m; end g(u,v)=(-1)^n; n=0; end end k=g*x*g; 新建 一个M文件，输入以下代码保存并运行： function code img=imread(lena.bmp); % lena.bmp为512x512像素，256色图 subplot(1,2,1),imshow(img);title(原始图像);