请问这么求图象的频谱对吗?

shenzhi123

我想请问一下这么求图象的频谱对吗?   
I=imread('CIMG1217.JPG');   
图象是一个三维图象
I=rgb2gray(I);
I=fftshift(fft2(I));
I=double(I);
F=fft2(I);
u=1:1:800;     该图象是800*197大小的   但为什么[M,N]=size(I);u=1:1:M;v=1:1:N;替代这两行u=1:1:800;  v=1:1:197;就没有频谱了呢?
v=1:1:197;
subplot(122);surf(u,v,F);
shading interp
谢谢大家帮我看看了

yangzj

把surf(u,v,F);改成surf(v,u,F);

另外,好像你做了两次FFT

shenzhi123

如果把surf(u,v,F);改成surf(v,u,F);
再去掉I=fftshift(fft2(I));或者F=fft2(I);其中的一个，其余的不变，得出来的是图１
而且有错误：Warning: Matrix dimensions must agree, not rendering mesh.
Warning: Matrix dimensions must agree, not rendering mesh.
Warning: Matrix dimensions must agree, not rendering mesh.
Warning: Matrix dimensions must agree, not rendering mesh.
如果不改成surf(v,u,F);再去掉I=fftshift(fft2(I));或者F=fft2(I);其中的一个，得到的是图２，但是没错误

shenzhi123

请yangzj 帮忙再解决下好吗？

yangzj

I=imread('CIMG1217.jpg');
I=rgb2gray(I);
figure();
subplot(211);
imshow(I);
F=fftshift(fft2(double(I)));
F=abs(F);
[M N]=size(I);
u=1:1:M;
v=1:1:N;
subplot(212);
surf(v,u,F);
shading interp

试试

shenzhi123

谢谢啊，我试了，可是还有个问题，这个得出来的是频谱吗？我总觉得有点不像呢？

yangzj

I=zeros(800,800);
a=10;
I(400-a:400+a,400-a:400+a)=255*ones(2*a+1,2*a+1);
figure();
subplot(211);
imshow(I);
F=fftshift(fft2(double(I)));
F=abs(F);
[M N]=size(I);
u=1:1:M;
v=1:1:N;
subplot(212);
imshow(F,[]);

呵呵,那你试试这个,看看对不对

shenzhi123

真谢谢你啊，可还需要向你请教一下:@$
图好象还是有点不对,运行出来出现的应该是幅度聘谱增强图吧
如果把imshow()改成mesh()的话出现的是图6的样子,

shenzhi123

不好意思,6的第2部分截错了,是这个样子的图,

yangzj

怎么来表示这个频谱那就看你喜欢啦.

shenzhi123

哦,明白了,谢谢啊,真的很感谢你,我可以给你加分表示感谢吗?:loveliness:
还有再请教一个问题就是,最后一个是用的是自己编写的矩阵,出来的看着很明显,
但是当我用拍的图象的RGB时候时,只能先转换成灰度在做频谱,要不就画不出来,
麻烦你再帮我确认下好吗?是不是我的截图4就是频谱了,只不过没最后一个图明显?方法已经是对的了?
谢谢啊!!!!!

yangzj

原帖由 shenzhi123 于 2007-5-3 20:27 发表
哦,明白了,谢谢啊,真的很感谢你,我可以给你加分表示感谢吗?:loveliness:
还有再请教一个问题就是,最后一个是用的是自己编写的矩阵,出来的看着很明显,
但是当我用拍的图象的RGB时候时,只能先转换成灰度在做频 ...

嗯,是对的,至于要从这个频谱图看出什么问题来,就要你来分析了

[ 本帖最后由 yangzj 于 2007-5-3 20:42 编辑 ]

meteor_033

这样得到应该是圆孔衍射的频谱吧？ 不是图象的频谱