2023年計(jì)算機(jī)視覺實(shí)驗(yàn)報(bào)告Eperiment

Experiment3:EdgeDetection

Class：電子1203班StudentID:Name:王影

I.Aim

Theaimofthis1aboratorysessionistolearntodealwithimagedata

byMat1ab.Bytheendofthissession,youshouldbeableto

performimagepreprocessingofedgedetectioninspatialdomain

andfrequencydomain.

II.KnowledgerequiredintheExperiment

i.Youaresupposedtohave1earnedthebasicski1IsofusingMatlab;

ii.YouneedtoreviewMat1abprogramminglanguageandM-filefor

mat.

iii.Youshou1dhavestudiededgedetectionmethods.

IH.ExperimentContents

Demand:Pleaseshowthefigureonthe1eftandlistthecodeson由erightres

pectivelybel1oweachquestion.(請(qǐng)將運(yùn)營結(jié)果(圖片)和程序代碼貼在每題下方)

i.Read“car?jpg”fi1e(todothisbyimreadfunction),convertthecolor

imageintograysealeimage,andthenperfbrmedgedetectionusingRo

terts,Prewitt,Sobe1operatorseparate1yinspatia1domainanddisp1

aytheresultsinaMatlabwindow.

程序：

c1ear;

im=imread('car.jpg');

I=rgb2gray(im);

subplot(3Z2Z1)；imshow(I)；

title('Grayimage*);

[Y,X]=size(I)；

im_edge=zeros(Y,X)；

T=30;

fork=2:Y-l

forkk=2:X-l

im_edge(k,kk)=abs(I(k+lzkk+1)-I(kzkk))+abs(I(k,kk+

1)-1(k+lzkk))；

if(im_edge(k,kk)>T)

im_edge(k,kk)=1;

else

im_edge(k,kk)=0；

end

end

end

subplot(3,2,2)；

imshow(im_edge,[])；%[]EAEy^YxQSl-Ed0~255pA-SII§AU\￡

tit1e('Robertimage');

[Yx]=size(I)；

imedge=zeros(Y,X);

fork=2:Y-1

forkk=2:X-l

imedge(k,kk)=abs(I(k-1,kk+l)-I(k-lzkk-1))+abs(I(k,kk

+1)-im(kzkk-1))+abs(I(k+1zkk+1)-I(k+1,kk-1))+??.

abs(I(k+1,kk-1)-1(k-lzkk-1))+abs(I(k+lzkk)-I(k

-1,kk))+abs(I(k+l,kk+l)-工(k-l,kk+1));

end

end

subp1ot(3,2,3);

imshow(imedge,[])；

title(1Prewitimage*)；

[YX]=size(I);

im_edge=zeros(YzX);

fork=2:Y-l

forkk=2:X—1

im_edge(k,kk)=abs(I(k-1,kk+1)-I(k-1,kk-1))+2*abs(工(k,k

k+1)-I(k,kk-1))+abs(I(k+1,kk+1)-1(k+1,kk-1))+...

abs(工(k+l,kk-1)-I(k-lzkk-1))+2*abs(I(k+lzkk)

-I(k-1,kk))+abs(I(k+1zkk+1)-I(k-1,kk+l));

end

end

subplot(3,2,4)；

imshow(im_edgez[]);

title('Sobelimage1);

圖像如下：

ii.Read“car.jpg”fi1e(todothisbyimreadfunetion),convertthecolor

imageintograyscaleimage,thenperformedgedetectioninfrequencydoma

inusingGaussianHighpassfiIteranddisplaytheresuItinaMatlab

window.

第二題程序：

%頻域邊沿檢測(cè)，運(yùn)用Gassian高通濾波器進(jìn)行濾波，進(jìn)行邊沿檢測(cè)

%頻域邊沿檢測(cè)，運(yùn)用Gassian高通濾波器進(jìn)行濾波，進(jìn)行邊沿檢測(cè)

clear;

im=imread('car.jpg*);

I=rgb2gray(im)；

subplot(1,2,1);

imshow(I);

title(/grayimage);

%shiftingimage(mu1tiplytheimageby(-l)x+y)

[row,col]=size(I);

[Y,X]=meshgrid(1:co1,1：row);

II=double(I).*(/).A(X+Y);

F=fft2(11);

%subplot(2,2,2);

%title('Fourierspectrum');

%creathighpassfiIter

D=zeros(row,col);

u0=floor(row/2);%傅立葉變換中心

v0=floor(col/2);

D0=40;%截止頻率

n=2;

fori=l：row

forj=l:col

d=((i-uO)-2+(j-v0)A2)A0.5；

%D(i,j)=1/(1+(D0/d)A(2*n));

D(i,j)=l-exp((-d八2)/(2*(D0)人2));

end

end

%fi1tering

G=F.*D;

%Inverttheresu1tandshifting

g=rea1(ifft2(G))；

im=g.*(-l).A(X+Y);

im=im>40;%閾值擬定edge

subplot(l,2,2);

imshow(im);

t