MATLAB绘制图像频谱程序

%filename = 'C:\Users\zhangP\Desktop\MMM\tem.jpg'; [filename, pathname]=uigetfile({'*.jpg';'*.png';'*.tif';'*.gif';'*.*'},'Select the image'); file_temp=strcat( pathname,filename); imgRgb = imread(filename); % 读入一幅彩色图像 imgray0=rgb2gray(imgRgb); imgray1=double(imgray0); mat=fft2(imgray1); sfftI=fftshift(mat); % 直流分量移到频谱中心 RR=real(sfftI); % 取傅立叶变换的实部 II=imag(sfftI); % 取傅立叶变换的虚部 A=sqrt(RR.^2+II.^2); % 计算频谱幅值 A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255; %归一化 %imshow(A); % 显示原图像的频谱 [m,n]=size(A); [x,y]=meshgrid(1:n,1:m); mesh(x,y,A) %surf(x,y,A) %或者surf %imhist(imgRgb); %imshow(imgRgb); % 显示彩色图像