资源描述：
《数字图像处理实验二 图像变换与二维数字滤波》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验二、图像变换与二维数字滤波一、实验目的1、了解图像正变换和逆变换的原理。2、理解图像变换系数的特点。3、掌握常用图像变换的实现过程。4、了解图像变换在图像数据压缩等方面的应用。5、掌握二维数字滤波器的作用以及在图像滤波中的实际应用。二、实验主要仪器设备1、计算机2、安装了图像处理工具箱的MATLAB软件。三、实验原理1、二维离散傅里叶变换、余弦变换、小波变换的正逆变换公式。2、图像的频谱分析原理。3、讨论正交变换的应用。4、二维FIR滤波器在图像滤波中的应用。四、实验内容1、在MATLAB环境中，进行图像的离散傅里叶变换和离
2、散余弦变换，观察图像的频谱并减少DCT系数，观察重建信号和误差信号，理解正交变换在压缩编码中的应用。2、在zMATLAB环境中，进行图像的近似图像和各方向的细节图像，观察重建图像，理解小波变换在图像特征检测（如边缘检测、方向检测等）中的应用。3、在MATLAB环境中，实验图像的二维FIR滤波。五、实验步骤1、选择典型图像作为研究对象（football.jpg）。2、显示原始图像。3、进行图像变换（DFT、DCT、DWT）。4、对图像进行处理（如选择不同个数的变换系数可以进行压缩）。5、对图像进行逆变换复原图像，观察重建图像和误差
3、图像进行分析。六、实验程序（1）I=imread('italy.jpg');subplot(3,3,1);imshow(I);J=rgb2gray(I);subplot(3,3,2);imshow(J);F1=fft2(J);subplot(3,3,3);imshow(log(abs(F1)+1),[010]);F2=fftshift(F1);subplot(3,3,4);imshow(log(abs(F2)+1),[010]);K=dct2(J);subplot(3,3,5),imshow(log(abs(K))+1,[010
4、]);（2）I=imread('italy.jpg');subplot(3,3,1);imshow(I);title('原图');X=rgb2gray(I);subplot(3,3,2);imshow(X);title('灰度');[cA1,cH1,cV1,cD1]=dwt2(X,'bior3.7');[cA1,cH1,cV1,cD1]=dwt2(X,'bior3.7');A1=upcoef2('a',cA1,'bior3.7',1);H1=upcoef2('h',cH1,'bior3.7',1);V1=upcoef2('v',
5、cH1,'bior3.7',1);D1=upcoef2('d',cD1,'bior3.7',1);V1=upcoef2('v',cV1,'bior3.7',1);subplot(3,3,3);image(wcodemat(A1,192));title('近似A1');subplot(3,3,5);subplot(3,3,4);image(wcodemat(H1,192));title('水平细节H1');subplot(3,3,5);image(wcodemat(V1,192));title('垂直细节V1');subplot(
6、3,3,6);image(wcodemat(D1,192));title('对角细节D1');Y=2.0*idwt2(A1,H1,V1,D1,'bior3.7');Y1=imresize(Y,0.461);subplot(3,3,7);image(Y);title('逆变换');（3）I=imread('italy.jpg');subplot(3,3,1);imshow(I);title('原图');X=rgb2gray(I);subplot(3,3,2);imshow(X);title('灰度');[cA1,cH1,cV1,c
7、D1]=dwt2(X,'bior3.7');[cA1,cH1,cV1,cD1]=dwt2(X,'bior3.7');A1=upcoef2('a',cA1,'bior3.7',1);H1=upcoef2('h',cH1,'bior3.7',1);V1=upcoef2('v',cH1,'bior3.7',1);D1=upcoef2('d',cD1,'bior3.7',1);V1=upcoef2('v',cV1,'bior3.7',1);subplot(3,3,3);image(wcodemat(A1,192));title('近似A
8、1');subplot(3,3,5);subplot(3,3,4);image(wcodemat(H1,192));title('水平细节H1');subplot(3,3,5);image(wcodemat(V1,192));title('垂直细节V1');su
