MATLAB绘制图像频谱程序

MATLAB绘制图像频谱程序，可以绘制出选中图片的二维频谱图。本人使用MATLAB7.0b

%filename = ‘C:\Users\zhangP\Desktop\MMM\tem.jpg‘;
[filename pathname]=uigetfile（{‘*.jpg‘;‘*.png‘;‘*.tif‘;‘*.gif‘;‘*.*‘}‘Select the image‘）;  
file_temp=strcat（ pathnamefilename）;
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）;        % 显示原图像的频谱

[mn]=size（A）;
[xy]=meshgrid（1:n1:m）;
mesh（xyA）


%surf（xyA） %或者surf
%imhist（imgRgb）;
%imshow（imgRgb）; % 显示彩色图像

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      28638  2012-07-03 22:49  imageread\1.png

     文件      31715  2012-07-03 22:48  imageread\2.png

     文件      30120  2012-07-03 22:50  imageread\3.png

    I.A....   1410560  2012-07-03 10:23  imageread\matlab_图像处理工具箱.ppt

     文件        714  2012-07-03 11:40  imageread\READ.asv

     文件        714  2012-07-03 11:42  imageread\READ.m

     文件      22451  2012-06-28 15:33  imageread\tmp.jpg

     文件      20391  2012-07-03 22:32  imageread\tmp1.jpg

     文件      18761  2012-07-03 22:35  imageread\tmp2.jpg

     目录          0  2012-07-03 23:14  imageread

----------- ---------  ---------- -----  ----

              1564064                    10