先行后列移动后的频谱图2程序.txt

图像处理 证明二维夫里叶变换等价于两步一位傅里叶变换的程序

I=imread('i.bmp');
X=double(I);
[row,col]=size(X);
for v=1:row
I2(v,:)=fft(X(v,:));
end
for u=1:col
I3(:,u)=fft(I2(:,u));% 得到图象X的某一(u)列数据
end
I4=fftshift(I3);
RR=real(I4);
II=imag(I4);
C=sqrt(RR.^2+II.^2);
C=(C-min(min(C)))/(max(max(C))-min(min(C)))*255;
figure;
imshow(C)