DSP实验进度汇报组员：汪张扬、任艳波、陈雪松、谢聪、沈旭任务分配：汪张扬由于考G，上周没有任务，沈旭负责自制二值图像的处理，陈雪松和谢聪负责其他图片的处理，任艳波负责搜集图像压缩评价的相关材料以下为简要概括：读入图像进行傅里叶变换和压缩原始程序：a=imread('d:\1.jpg');b=figure;imshow(a);title('原始图像');F=fft2(a);F_mm=abs(F);figure;imshow(F);title('原始幅度谱');Fshift=fftshift(F);F_m=abs(Fshift);figure;imshow(F_m);title('幅度谱');F_p=angle(Fshift);figure;imshow(F_p);title('相位谱');T=@fft2;B1=blkproc(a,[88],T);%将图像分块为8×8矩阵进行处理figure;imshow(a);title('原始图像');mask=[1000000001000000001000000001000000000000000001000000001000000001];%与该矩阵相乘去掉中间行，即高频部分B2=blkproc(B1,[88],'P1*x',mask);fun=@ifft2;F3=blkproc(B2,[88],fun);F=mat2gray(F3);figure;imshow(F);title('压缩87.5%的图像');刚开始的原始图像：原始图像傅里叶变换后的原始幅度谱：原始幅度谱取高频点后的图像：87.5%压缩的图像取对数后的频谱图：图像进行取整后的i=imread('d:\1.jpg');figure(1);imshow(i);colorbar;j=fft2(i);k=fftshift(j);h=floor(k);n=ifft2(h)/255;figure(3);m=imresize(n,2);imshow(m,[]);colorbar-0.8-0.6-0.4-0.200.20.40.60.8取小幅值为零：i=imread('d:\1.jpg');figure(1);imshow(i);colorbar;j=fft2(i);k=fftshift(j);k(abs(k)<1)=0;g=ifft2(k)/255;figure(2);h=imresize(g,2);imshow(h);colorbar00.10.20.30.40.50.60.70.80.91结论：二维黑白图像的压缩主要有三种方式：取整，让小幅值为零，或者利用矩阵的变换只取某些变换后的值进行重建，如矩阵的左上角为低频部分，取它即可。取整重建将黑色部分变白了，小幅值变换将白色部分变黑了。二、实验前后的结构比较（一）、自然景观（1）、原图：（2）、更改DFT系数为整数，做IDFT观察图像的变化更改DFT系数为整数figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化更改小幅度频谱系数为0figure3更改后时域与原时域的误差为1471.6806242515此处取频谱虚实部中最大值的0.0005作为取零的阈值。（二）、人物照片（1）、原图（2）、更改DFT系数为整数，做IDFT观察图像的变化更改DFT系数为整数figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化更改小幅度频谱系数为0figure3更改后时域与原时域的误差为1349.14370621138此处取频谱虚实部中最大值的0.0005作为取零的阈值。（三）、卡通图片（1）、原图（2）更改DFT系数为整数，做IDFT观察图像的变化更改DFT系数为整数：figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化更改小幅度频谱系数为0figure3更改后时域与原时域的误差为2546.07530759797此处取频谱虚实部中最大值的0.0005作为取零的阈值。三、实验分析与说明（1）、实验中图片频谱中的小幅度部分需要探讨，此处取的是最大值的0.0005倍。（2）、实验中原图取自jpg格式的rgb彩色图。（3）、实验的结果图也已jpg格式存储（4）、jpg格式本身就是采用了频谱取整取零技术。四、实验结论（1）、实验中小幅度的阈值不能取得太大，否则图片会失真。本次试验去的是0.0005倍，出现较明显的失真；实验中发现取为0.0001，就几乎看不出失真。（2）、关于频谱取整，我想应该是要将4字节的double变成2字节的int，可以节省约50%的存储空间。而且从本次实验可以看出频谱取整几乎不失真。（3）、关于频谱小幅度取零，这个阈值需要根据具体图片来确定，或者寻找自适应的阈值；以为它对图片的失真度和压缩率也起主要作用。阈值太小，越不容易失真，但要压缩率也越小；阈值太大，越容易失真，但压缩率也越大。（4）、自然景观、人物照片、卡通图片这三类对同意阈值的敏感度不同；下附实验代码：下面是对对片的处理的全部代码，图片名K1.jpg为待处理的图片;%2）更改DFT系数为整数，做IDFT观察图像的变化；%3）更改小幅值的DFT系数为0，做IDFT变化观...
