95-04-15.m

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

clear all;
f=imread('d:\image\coins.tif');
PQ=paddedsize(size(f));
D0=0.05*PQ(1);
liw=lpfilter('ideal',PQ(1),PQ(2),D0);
gli=dftfilt(f,liw);
lbw=lpfilter('btw',PQ(1),PQ(2),D0,2);
glb=dftfilt(f,lbw);
lgw=lpfilter('gaussian',PQ(1),PQ(2),D0);
glg=dftfilt(f,lgw);
hiw=hpfilter('ideal',PQ(1),PQ(2),D0);
ghi=dftfilt(f,hiw);
hbw=hpfilter('btw',PQ(1),PQ(2),D0,2);
ghb=dftfilt(f,hbw);
hgw=hpfilter('gaussian',PQ(1),PQ(2),D0);
ghg=dftfilt(f,hgw);
hhiw=hpfilter('ideal',PQ(1),PQ(2),D0);
hhiw2=0.5+2*hhiw;
gei=dftfilt(f,hhiw);
hhbw=hpfilter('btw',PQ(1),PQ(2),D0,2);
hhbw2=0.5+2*hhbw;
geb=dftfilt(f,hhbw2);
hhgw=hpfilter('gaussian',PQ(1),PQ(2),D0);
hhgw2=0.5+2*hhgw;
geg=dftfilt(f,hhgw2);
figure
subplot(3,4,1);
imshow(f)
subplot(3,4,2);
gli=gscale(gli);
imshow(gli,[])
subplot(3,4,3);
glb=gscale(glb);
imshow(glb,[])
subplot(3,4,4);
glg=gscale(glg);
imshow(glg,[])
subplot(3,4,6);
ghi=gscale(ghi);
imshow(ghi,[])
subplot(3,4,7)
ghb=gscale(ghb);
imshow(ghb,[])
subplot(3,4,8);
ghg=gscale(ghg);
imshow(ghg,[])
subplot(3,4,10);
gei=gscale(gei);
imshow(gei,[])
subplot(3,4,11);
geb=gscale(geb);
imshow(geb,[])
subplot(3,4,12);
geg=gscale(geg);
imshow(geg,[])