机器视觉及matlab仿真研究.docVIP

机器视觉技术及应用 学院：电气与信息工程学院 班级：自动化072 姓名：景斌 学号：200706070217 2011-01-04 题目：一 产生下图所示图像 f1(m,n)，其中图像大小为 256×256，中间亮条为128×32，暗处=0，亮处=100。对其进行FFT： ① 同屏显示原图f1(m,n)和FFT(f1)的幅度谱图； ② 若令，重复以上过程，比较二者幅度谱的异同，简述理由； ③ 若将 f2(m,n)顺时针旋转 90 度得到 f3(m,n)，试显示 FFT(f3)的幅度谱，并与 FFT(f2) 的幅度谱进行比较； ④ 若将f1(m,n) 顺时针旋转90度得到f4(m,n)，令f5(m,n)=f1(m,n)+f4(m,n)，试显示FFT(f5) 的幅度谱，并指出其与 FFT(f1)和FFT(f4)的关系； ⑤ 若令f6(m,n)=f2(m,n)+f3(m,n)，试显示FFT(f6)的幅度谱，并指出其与 FFT(f2)和FFT(f3) 的关系，比较FFT(f6)和FFT(f5)的幅度谱。 1、显示图像 clear,close all *\清除所有工作平台变量，关闭打开的图形窗口\* a=zeros(256,256); *\设置256*256大小的图形\* a(64:192,112:144)=100; *\设置128*32大小的亮条\* imshow(a) *\在窗口中显示要求图像\* 2.对其进行FFT: clear,close all a=zeros(256,256); a(64:192,112:144)=100; *\设置128*32大小的亮条\* imshow(a) F=fft2(a); *\对a进行傅立叶变换\* b=log(abs(F)); *\对F取对数\* imshow(b,[.1 5]);colormap(jet);colorbar *\显示b图像\* 3.同屏显示原图f1（m,n）和FFT（f1）： clear,close all *\清除所有工作平台变量，关闭打开的图形窗口\* a=zeros(256,256); *\设置256*256大小的图形\* a(64:192,112:144)=100; subplot(1,2,1),imshow(a) *\ 在窗口中显示原图像(第一行第一列) \* F=fft2(a); *\对a进行傅立叶变换\* b=log(abs(F)); *\对F取对数\* subplot(1,2,2),imshow(b) *\在窗口中显示b图像(第一行第二列) \* 4.实现 clear,close all a=zeros(256,256); a(64:2:192,112:2:144)=100; *\隔行隔列显示\* imshow(a) a(65:2:191,113:2:143)=100; *\隔行隔列显示\* imshow(a) 5.对f2进行二维快速傅里叶变换并图原图同屏显示： clear,close all *\清除所有工作平台变量，关闭打开的图形窗口\* a=zeros(256,256); *\设置256*256大小的图形\* a(64:2:192,112:2:144)=100; *\ 隔行隔列显示\* a(65:2:191,113:2:143)=100; *\隔行隔列显示\* subplot(1,2,1),imshow(a); F=fft2(a); b=log(abs(F)); *\傅立叶变换\* subplot(1,2,2),imshow(b) 6.将 f2(m,n)顺时针旋转 90 度得到 f3(m,n)，试显示 FFT(f3)的幅度谱 clear,close all a=zeros(256,256); a(64:2:192,112:2:144)=100; *\隔行隔列显示\* a(65:2:191,113:2:143)=100; *\隔行隔列显示\* c=imrotate(a,90); *\旋转90度\* subplot(1,2,1),imshow(a); *\一行两列显示\* subplot(1,2,2),imshow(c) 7