%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）;??????%?取傅立葉變換的實(shí)部
II=imag（sfftI）;??????%?取傅立葉變換的虛部
A=sqrt（RR.^2+II.^2）;?%?計(jì)算頻譜幅值
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）;?%?顯示彩色圖像

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????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