DIP 例4.15 空间域滤波得到频率域滤波

做了整整一天才完成这个例子。完成的代码也存在一个BUG，在最后提一下好了。

J=imread('Fig1016(a)(building_original).tif');imshow(J)
p=double(J);
P=834+2;Q=1114+2;%size(J)
transform=ones(P,Q);
for i=1:P
    for j=1:Q
        transform(i,j)=(-1)^(i+j);
    end