【求助】关于MATLAB图像的傅立叶变换

[zukaldl]

物体图像的生成可以直接由矩阵运算生成，也可利用Windows下的画图工具，生成一幅黑白图像，并调用命令函数imread（）输入图像，输入的图像是一个巨大的二维矩阵，利用MATLAB函数库中的fft2（）命令对该矩阵进行二位离散傅立叶变换，得到图像的频谱，该频谱是一个复数矩阵，然后用取模函数abs（）对该复数矩阵取模，得到振幅谱矩阵，利用函数fftshift（）对取模后的矩阵进行频谱位移，这是因为变换后的二维矩阵的直流分量位于图像的周边角，该函数交换矩阵的1、3象限和2、4象限，使直流分量移到频谱中心，从而使FFT频谱可视效果与实际图像相吻合。最后利用imshow（）函数将图像显示出来。
编写程序如下：

clear
a=imread('E:\1.bmp');
grid on
figure(1)
imshow(a,[])
afft=fft2(a);
aabs=abs(afft);
aabss=fftshift(aabs);
figure(2)
imshow(aabss,[])
colormap(gray);
colorbar
figure(3)
plot(aabss)
colormap(gray);
figure(4)
meshc(aabss)
maxx1=max(max(aabss));
但是程序实现不了，请高手改正指点，谢谢～～！

[zukaldl]

指点一下~~

[zukaldl]

怎样搞的~~

[zukaldl]

~~~~~~~~~~~~~~~~~~~~~

[zukaldl]

```````````````````````````````````

[ninja]

aabs=abs(afft); ------> aabs=abs(afft(:,:,1));
bmp图像为M*N*3的矩阵，只需取颜色表的一维来求解就行了。

[zukaldl]

我试试
先谢谢了

[zukaldl]

试过了
真谢谢了
~~
没想到等了1个多星期真来高手了
谢谢~~
这是我论文上的一个点
我想用MATLAB来证明的
但是我又没学过
所以才来求救
没想到真来救星了
再次谢谢~~
:lovely: :cry:

[ninja]

晕，今天忽然发现这个论坛的。刚才是新注册ID发的第一帖。呵呵

[zukaldl]

``````````````````````