证明相等的程序.txt

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

I=imread('i.bmp');
X=double(I);
[row,col]=size(X);
Y=fft2(X);
for u=1:col
I2(:,u)=fft(X(:,u));% 得到图象X的某一(u)列数据
end
for v=1:row
I3(v,:)=fft(I2(v,:));
end 

for v=1:row
I4(v,:)=fft(X(v,:));
end 
for u=1:col
I5(:,u)=fft(I4(:,u));% 得到图象X的某一(u)列数据
end

Z=Y-I3;
W=I3-I5;
for u=1:256
    for v=1:256
        if(abs(Z(v,u)<0.000001))
            Z(u,v)=0;
        end
        if(abs(W(v,u)<0.000001))
            Z(u,v)=0;
        end
        
    end
end
abs(Z(v,u))
abs(Z(v,u))
abs(Z(v,u))

>> 
ans =

     0


ans =

     0


ans =

     0

>>