071041219彭金金实验五数字图像的变换技术及其MATLAB实现.docVIP

数字图象处理实验 实验报告 学号：071041219 姓名：彭金金 班级：071041B 指导老师：石庆研 实验五 数字图像的变换技术及其MATLAB实现 一实验目的： 1、了解图像变换的意义和手段； 2、加深对傅里叶变换、离散余弦变换、沃尔什变换和小波变换基本原理的理解； 3、熟悉常用变换在图像增强、图像恢复、图像压缩以及特征提取等方面的应用； 4、掌握常用变换的MATLAB实现。 二实验内容： 利用MATLAB工具箱实现二维傅立叶正反变换，并利用傅里叶变换实现快速卷积： M文件： clc clear I=imread(lena.gif); J=fft2(I); K = fftshift(J); A=conv2(J,J); subplot(1,3,1);imshow(log(abs(J)),[],InitialMagnification,fit) title(原图傅立叶变换频谱) subplot(1,3,2);imshow(log(abs(K)),[],InitialMagnification,fit) title(原图傅立叶变换频谱中心变换) subplot(1,3,3);imshow(log(abs(A)),[],InitialMagnification,fit) title(卷积频谱) B=ifft2(A); figure subplot(1,2,1);imshow(I),title(原图) subplot(1,2,2);imshow(B),title(卷积后的图像) 2、利用MATLAB工具箱实现二维离散余弦正反变换，并利用DCT实现图像压缩： M文件： clc clear I=imread(lena.gif); I=im2double(I); K=dctmtx(8);%余弦变换后生成一个8*8矩阵 J1=blkproc(I,[8 8],P1*x*P2,K,K); M=[1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; J2=blkproc(J1,[8 8],P1.*x,M); I2=blkproc(J2,[8 8],P1*x*P2,K,K); subplot(1,2,1);imshow(J1);title(原始图像余弦变换后矩阵) subplot(1,2,2);imshow(J2);title(压缩后的图像余弦变换后矩阵) figure; subplot(1,2,1);imshow(I);title(原始图像) subplot(1,2,2);imshow(I2);title(压缩后的图像) 3、利用MATLAB提供的哈达玛变换函数实现图像的压缩： clc close clear I=imread(lena.gif); I1=im2double(I); [M,N]=size(I1); s=64; K=hadamard(8); J1=blkproc(I1,[8 8],P1*x*P2,K,K); L=im2col(J1,[8 8],distinct); L2=L; [Y,P]=sort(L); [m,n]=size(L); s=m-s; for i=1:n L2(P(1:s),i)=0; end J2=col2im(L2,[8 8],[M N],distinct); I2=blkproc(J2,[8 8],P1*x*P2,K,K); I2=I2/255; subplot(1,2,1);imshow(J1);title(原始图像哈达玛变换后矩阵) subplot(1,2,2);imshow(J2);title(压缩后的图像哈达玛变换后矩阵) figure; subplot(1,2,1);imshow(I);title(原始图像) subplot(1,2,2);imshow(I2);title(变换后的图像) 熟悉MALAB提供的一维和二维离散小波变换函数利用小波变换实现图像的去噪和压缩： (1)一维小波变换： Clc clear close [x,fs,n]=wavread(1.wav); [c1,d1]=dwt(x,db1); subplot(1,3,1);plot(x);title(原始信号) subplot(1,3,2);plot(c1);title(小波分解近似分量) subplot(1,3,3);plot(d1);title(小波分解细节分量) （2）二维小波变换实现图像去噪 clc clear close I=imread(lena.gif); In=imnoise(I,gaussian,0,0.