第一题
f=imread('E:\PICTRUE\Fig0438(a)(bld_600by600).tif');
f=tofloat(f);
F=fft2(f);
S=fftshift(log(1+abs(F)));
imshow(S,[ ])
h=fspecial('sobel')
freqz2(h)
PQ=paddedsize(size(f));
H=freqz2(h,PQ(1),PQ(2));
H1=ifftshift(H);
imshow(abs(H),[ ])
figure,imshow(abs(H1),[ ])
gs=imfilter(f,h);
gf=dftfilt(f,H1);
imshow(gs,[ ])
figure,imshow(gf,[ ])
figure,imshow(abs(gs),[ ])
figure,imshow(abs(gf),[ ])
figure,imshow(abs(gs)>0.2*abs(max(gs(:))))
figure,imshow(abs(gf)>0.*abs(max(gf(:))))
第二题
f=imread('E:\PICTRUE\Fig0438(a)(bld_600by600).tif');
f=tofloat(f);
F=fft2(f);