96-0408.m

Matlab數位影音處理 影像方面的處理

%96-04-08
clear all;
f=imread('d:\image\color\House.bmp');
fb=rgb2gray(f);
F=fft2(double(fb));
fd=fftshift(log(1+abs(F)));
fa=gscale(fd);
h=fspecial('sobel');
PQ=paddedsize(size(F));
H=freqz2(h,PQ(1),PQ(2));
H=ifftshift(H);
fs=imfilter(double(fb),h);
fs1=dftfilt(fb,H);
D0=0.05*PQ(1);
HL=lpfilter('ideal',PQ(1),PQ(2),D0);
%HL=fftshift(HL);
f1=dftfilt(fb,HL);
HB=lpfilter('btw',PQ(1),PQ(2),D0,2);
%HB=fftshift(HB);
f2=dftfilt(fb,HB);
HG=lpfilter('gaussian',PQ(1),PQ(2),D0);
%HG=fftshift(HG);
f3= real(dftfilt(fb,HG));
figure,imshow(fb)
figure
subplot(3,3,1)
imshow(abs(fa),[])
subplot(3,3,2)
imshow(abs(fs),[])
subplot(3,3,3)
imshow(abs(fs1),[])
subplot(3,3,4)
imshow(f1,[])
subplot(3,3,5)
imshow(f2,[])
subplot(3,3,6)
imshow(f3,[])