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