0=15*log(m+1.001);surf(m0);求出数组之后将其进行fft2得到m,然后求出其绝对值,为了避免其数值过大,求其对数,且为了避免出现0和1的情况,可以在m的基础上加上1.001,最后用surf()函数显示出3----D效果.(2)r=24;forx=1:256fory=1:256if(x-128).^2+(y-128).^2<r.^2;t(x,y)=0;endendend通过对r数值的改变,和if条件的变化来实现不同的低通和高通情况下的滤波.五.实验结果及结论1.验证二维傅里叶变换的平移性和旋转不变性;1.1产生的图像与fft2(f1)幅度谱的程序如下:clearall;closeall;a=zeros(128,128)fory=54:74forx=34:94a(x,y)=1;endendfigure(1)a1=fft2(a);subplot(1,2,1);imshow(a);subplot(1,2,2);a2=abs(a1);mesh(a2);forx=1:128fory=1:128b(x,y)=(-1).^(x+y).*a(x,y);endendfigure(2)b1=fft2(b);subplot(1,2,1);imshow(b);subplot(1,2,2);b2=abs(b1);mesh(b2);figure(3)t=imrotate(a,315,'nearest','crop')t1=fft2(t);subplot(1,2,1);imshow(t);subplot(1,2,2);t2=abs(t1);surf(t2);结果如图1.2令,则图像与的幅度谱图如下:结果分析:对比两图可以得到,1.2得的图就是将1.1得的图的频谱往中心移.即称为频谱中心化,将能量集中的4个角往中心靠拢。1.3若将f2(x,y)顺时针旋转45°得到,则及的幅度谱图如下
全文预览
