资源描述:
《实验二图像几何变换与傅里叶变换》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验二图像几何变换与傅里叶变换课程名称数字图像处理实验名称—图像儿何变换与傅里叶变换姓名学号专业班级实验日期2013年―月24日成绩指导教师一、实验目的与要求1.掌握图像平移、缩放、旋转与镜像变换;。2.掌握傅里叶变换和反变换的方法,观察图像的频域表示的特点二、实验例程1、几何变换图像'卜移,H写平移函数immove.infunctionoutImagc=immovc(inImagc,Tx,Ty,i)Im,n,i]=size(inlmage);%计算图象大小outlmage=uint8(zeros(m+abs(Ty),n+abs(Tx),i
2、));if(Tx>0&&Ty>0)%往右下角平移的情况ou(Image(l+Ty:m+Ty,l+Tx:n+Tx,l:i)=inImage;elseif(Tx<0&&Ty<0)%往左I:角平移的情况outlmage(1:m,1:n,1:i)=inImage;elseif(Tx<0&&Ty>0)%往左下角平移的情况outImage(l+Ty:m+Ty,1:n,l:i)=inImage;else%往右上角平移的情况outlmage(l:m,l+Tx:n+Tx,l:i)=inImage;enclendend2、傅里叶变换1、matIab计算和形象化
3、2DDFT傅立叶变换函数和逆向傅立叶变换函数F=fft2(f)仁real(ifft2(F))例子:I=imread('dragfl.bmp');仁real(ifft2(F));F=fft2(I);imshow(F,[])fl=real(ifft2(F));imshow(f,[J)傅立叶频谱计算函数s二abs(F)例子:频谱可视化I=imrcad('dragfl.bmp');F=fft2(I);S=abs(F);imshow(S,[])中心频谱函数:把变换的原点移到频域区域中心Fc=fftshift(F)把变换的频域区域中心移到原点:F=if
4、ftshift(F)例子:中心频谱可视化I=imread('dragfl.bmp1);F=fft2(I);Fc=fftshift(F);S=abs(Fc);imshow(S,[])f2=ifftshift(Fc);imshow(abs(f),[])频谱动态范围压缩(频谱增强):例子:^imreadCdragfl.bmp');F=fft2(I);S=abs(Fc);S2=log(l+abs(Fc));imshow(S2,[])三、实验内容及步骤1.读入彩色图像,实现图像的平移、缩放、旋转和平移I=imread('d:autumn.tif);
5、H=immove(I,80,-80);figure;imshow(H);title('平移后图像);BIJ.mmove(Ibo0「50);figure;imshow(B)>le(■卡陵迅函謹);-H抿5RH藏GN.mresize(Io・5);figura.mshow(G)s-A^」EEls^)_C=imrcsizc(I,2);figure;imshow(C);title(缩放后图像);缩放后图像D=imrotate(l,30);figure;imshow(D);title('旋转后图像');旋转后图像E=imrotatc(I,54);fig
6、ure;imshow(E);title(旋转后图像');旋转后图像F=flipdim(IJ);figure;imshow(F);title(唾直镜像图像);垂直镜像图像G=flipdim(L2);figure;imshow(G);title(,水平镜像图像')小十坦I豕QI豕imcrop(LL80,80,80,80J)imcrop(I,[45,75,65,61])2、读入一幅灰度图像,实现图像的平移、缩放、旋转和平移J=rgb2gray(I);imshow(J);title('M图像');原图像H=immove(J,80,-80);figu
7、re;imshow(H);titleC平移后图像');平移后图像B=immove(J,50,-50);figure;imshow(B);title(,A
8、z移后图像);平移后图像G=imrcsizc(J,0.3);figure;imshow(G);title(缩放后图像');缩放后图像c".mrcsizc(J2)八Hgul.e;imshow(c)三tA^^Gs^)-DJ.mroul
9、像F=flipdim(J,l);figure;imshow(F);title(,垂直镜像图像J;垂直镜像图像G=flipdim(J,2);figure;imshow(G);tit