资源描述：
《实验的3图像变换》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用标准文案实验三图像变换一、实验目的 1、熟悉及掌握图像的变换原理及性质，实现图像的傅里叶变换。二、实验环境MATLAB6.5以上版本、WINXP或WIN2000计算机三、相关函数l图像的变换1fft2fft2函数用于数字图像的二维傅立叶变换，如：i=imread('e:w01.tif');j=fft2(i);由于MATLAB无法显示复数图像，因此变换后的结果还需进行求模运算，即调用abs函数。之后常常还进行对数变换，即调用log函数，以减缓傅里叶谱的快速衰减，更好地显示高频信息。2ifft2ifft2函数用于数字图像的二维傅
2、立叶反变换，如：i=imread('e:w01.tif');j=fft2(i);k=ifft2(j);3fftshift用于将变换后图像频谱中心从矩阵的原点移动到矩阵的中心B=fftshift（i）4利用fft2计算二维卷积精彩文档实用标准文案利用fft2函数可以计算二维卷积,如:a=[8,1,6;3,5,7;4,9,2];b=[1,1,1;1,1,1;1,1,1];a(8,8)=0;b(8,8)=0;c=ifft2(fft2(a).*fft2(b));c=c(1:5,1:5);利用conv2(二维卷积函数)校验,如：a=[8,
3、1,6;3,5,7;4,9,2];b=[1,1,1;1,1,1;1,1,1];c=conv2(a,b);四、实验内容1、对一幅图像进行平移，显示原始图像与处理后图像，分别对其进行傅里叶变换，显示变换后结果，分析原图的傅里叶谱与平移后傅里叶频谱的对应关系。s=imread('f:1.jpg');%读入原图像i=rgb2gray(s)i=double(i)j=fft2(i);%傅里叶变换k=fftshift(j);%直流分量移到频谱中心l=log(abs(k));%对数变换m=fftshift(j);%直流分量移到频谱中心RR=re
4、al(m);%取傅里叶变换的实部II=imag(m);%取傅里叶变换的虚部精彩文档实用标准文案A=sqrt(RR.^2+II.^2);%计算频谱府幅值A=(A-min(min(A)))/(max(max(A)))*255;%归一化b=circshift(s,[800450]);%对图像矩阵im中的数据进行移位操作b=rgb2gray(b)b=double(b)c=fft2(b);%傅里叶变换e=fftshift(c);%直流分量移到频谱中心l=log(abs(e));%对数变换f=fftshift(c);%直流分量移到频谱中心WW
5、=real(f);%取傅里叶变换的实部BZZ=imag(f);%取傅里叶变换的虚部B=sqrt(WW.^2+ZZ.^2);%计算频谱府幅值B=(B-min(min(B)))/(max(max(B)))*255;%归一化subplot(2,2,1);imshow(s);title('原图像')subplot(2,2,2);imshow(uint8(b));;title('平移图像')subplot(2,2,3);imshow(A);title('离散傅里叶频谱');subplot(2,2,4);imshow(B);title('平移
6、图像离散傅里叶频谱')精彩文档实用标准文案2、对一幅图像进行旋转，显示原始图像与处理后图像，分别对其进行傅里叶变换，显示变换后结果，分析原图的傅里叶谱与旋转后傅里叶频谱的对应关系。s=imread('f:1.jpg');%读入原图像i=rgb2gray(s)i=double(i)j=fft2(i);%傅里叶变换k=fftshift(j);%直流分量移到频谱中心l=log(abs(k));%对数变换m=fftshift(j);%直流分量移到频谱中心RR=real(m);%取傅里叶变换的实部II=imag(m);%取傅里叶变换的虚部
7、A=sqrt(RR.^2+II.^2);%计算频谱府幅值A=(A-min(min(A)))/(max(max(A)))*255;%归一化精彩文档实用标准文案b=imrotate(s,-90);%对图像矩阵im中的数据进行移位操作b=rgb2gray(b)b=double(b)c=fft2(b);%傅里叶变换e=fftshift(c);%直流分量移到频谱中心l=log(abs(e));%对数变换f=fftshift(c);%直流分量移到频谱中心WW=real(f);%取傅里叶变换的实部BZZ=imag(f);%取傅里叶变换的虚部B=
8、sqrt(WW.^2+ZZ.^2);%计算频谱府幅值B=(B-min(min(B)))/(max(max(B)))*255;%归一化subplot(2,2,1);imshow(s);title('原图像')subplot(2,2,2);imshow
