傅里叶变换图像压缩

《傅里叶变换图像压缩》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文档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=ab

2、s(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=[100000000100000000100000000100000000000000000100文案大全实用标准文档0000001000000001];%与该矩阵相乘去掉中间行，

3、即高频部分B2=blkproc(B1,[88],'P1*x',mask);fun=@ifft2;F3=blkproc(B2,[88],fun);F=mat2gray(F3);figure;imshow(F);title('压缩87.5%的图像');刚开始的原始图像：傅里叶变换后的原始幅度谱：文案大全实用标准文档取高频点后的图像：取对数后的频谱图：文案大全实用标准文档图像进行取整后的i=imread('d:1.jpg');figure(1);imshow(i);colorbar;j=fft2(i);k=fftshift(

4、j);h=floor(k);n=ifft2(h)/255;figure(3);m=imresize(n,2);imshow(m,[]);colorbar文案大全实用标准文档取小幅值为零：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);colorbar文案大全实用标准文档结论：二维黑白图像的压缩主

5、要有三种方式：取整，让小幅值为零，或者利用矩阵的变换只取某些变换后的值进行重建，如矩阵的左上角为低频部分，取它即可。取整重建将黑色部分变白了，小幅值变换将白色部分变黑了。二、实验前后的结构比较（一）、自然景观（1）、原图：（2）、更改DFT系数为整数，做IDFT观察图像的变化文案大全实用标准文档更改DFT系数为整数figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化文案大全实用标准文档更改小幅度频谱系数为0figure3更改后时域与原时域的误差为1471.680624

6、2515此处取频谱虚实部中最大值的0.0005作为取零的阈值。（二）、人物照片（1）、原图（2）、更改DFT系数为整数，做IDFT观察图像的变化文案大全实用标准文档更改DFT系数为整数figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化文案大全实用标准文档更改小幅度频谱系数为0figure3更改后时域与原时域的误差为1349.14370621138此处取频谱虚实部中最大值的0.0005作为取零的阈值。（三）、卡通图片（1）、原图文案大全实用标准文档（2）更改DFT系数

7、为整数，做IDFT观察图像的变化文案大全实用标准文档更改DFT系数为整数：figure2更改后时域与原时域的误差为0（3）、更改小幅值的DFT系数为0，做IDFT变化观察图像的变化文案大全实用标准文档更改小幅度频谱系数为0figure3更改后时域与原时域的误差为2546.07530759797此处取频谱虚实部中最大值的0.0005作为取零的阈值。三、实验分析与说明（1）、实验中图片频谱中的小幅度部分需要探讨，此处取的是最大值的0.0005倍。（2）、实验中原图取自jpg格式的rgb彩色图。（3）、实验的结果图也已jpg格

8、式存储（4）、jpg格式本身就是采用了频谱取整取零技术。四、实验结论（1）、实验中小幅度的阈值不能取得太大，否则图片会失真。本次试验去的是0.0005倍，出现较明显的失真；实验中发现取为0.0001，就几乎看不出失真。文案大全实用标准文档（2）、关于频谱取整，我想应该是要将4字节的double变成2字节的int，可以