多媒体与数字图像处理第三次实验(2)

%相位

JP = angle(J); iJP = ifft2(JP);

figure, imshow(abs(iJP)*100); title('相位傅立叶反变换');

B) 利用MATLAB软件实现数字图像傅立叶变换的程序 I=imread(‘at3_1m4_04.tif'); %读入原图像文件 imshow(I); %显示原图像

fftI=fft2(I); %二维离散傅立叶变换 sfftI=fftshift(fftI); %直流分量移到频谱中心 RR= real(sfftI); %取傅立叶变换的实部 II= imag(sfftI); %取傅立叶变换的虚部 A=sqrt(RR.^2+II.^2);%计算频谱幅值

A=（A-min(min(A))）/(max(max(A))-min(min(A)))*225; %归一化 figure; %设定窗口

imshow(A); %显示原图像的频谱

C) 绘制一个二值图像矩阵,并将其傅立叶函数可视化。 f=zeros(30,30); f(5:24,13:17)=1;

imshow(f,'notruesize')

F=fft2(f);

F2=log(abs(F));

figure,imshow(F2,[-1 5],'notruesize');colormap(jet);