【正文】 m(i1,j1) ~= Im(i1,j) a = a + 1。 end if a==6 %%%分叉點(diǎn)判斷 tezheng(i,j,1)=i。 tezheng(i,j,2)=j。 tezheng(i,j,3)=1。 xxx=xxx+1。 I(i,j)=0。 I(i1,j)=1。I(i1,j+1)=1。I(i,j+1)=1。I(i+1,j+1)=1。I(i+1,j)=1。I(i+1,j1)=1。I(i,j1)=1。I(i1,j1)=1。 %%%%%%%%%%%%%%1表示分叉點(diǎn) 0表示端點(diǎn) end if a==2 %%%% 端點(diǎn)判斷 tezheng(i,j,1)=i。 tezheng(i,j,2)=j。 tezheng(i,j,3)=2。 xxx=xxx+1。 I(i,j)=0。 I(i1,j)=1。I(i1,j+1)=1。I(i,j+1)=1。I(i+1,j+1)=1。I(i+1,j)=1。I(i+1,j1)=1。I(i,j1)=1。I(i1,j1)=1。 end end endendfigure,imshow(I)%%%%%%%%%%%%%%%%%%%%特征提取后處理~Izz=I。for i=1:m %%%%%%%%%%%%%%%%%%%%%%%%去除邊界點(diǎn) for j=1:n if tezheng(i,j,3)~=0 if i170||i+17m|| j170||j+17n tezheng(i,j,3)=0。 xxx=xxx1。 %%%消除端點(diǎn) Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 else a=0。b=0。c=0。d=0。 a=I(i+1,j)+I(i+2,j)+I(i+3,j)+I(i+4,j)+I(i+5,j)+I(i+6,j)+I(i+7,j)+I(i+8,j)+I(i+9,j)+I(i+10,j)+I(i+11,j)+I(i+12,j)。 b=I(i1,j)+I(i2,j)+I(i+3,j)+I(i4,j)+I(i5,j)+I(i6,j)+I(i7,j)+I(i8,j)+I(i9,j)+I(i10,j)+I(i11,j)+I(i12,j)。 c=I(i,j+1)+I(i,j+2)+I(i,j+3)+I(i,j+4)+I(i,j+5)+I(i,j+6)+I(i,j+7)+I(i,j+8)+I(i,j+9)+I(i,j+10)+I(i,j+11)+I(i,j+12)。 d=I(i,j1)+I(i,j2)+I(i,j3)+I(i,j4)+I(i,j5)+I(i,j6)+I(i,j7)+I(i,j8)+I(i,j9)+I(i,j10)+I(i,j11)+I(i,j12)。 if ( a=1 amp。amp。 b=1 amp。amp。 c=1 amp。amp。 d=1) tezheng(i,j,3)=0。 %%%消除端點(diǎn) xxx=xxx1。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 end end end endend%%%%%%%%%%%%%%%%%%%%%%%%%%消除斷點(diǎn)for i=17:m17 for j=17:n17 if tezheng(i,j,3)==2 for x=i12:i+12 for y=j12:j1 if tezheng(x,y,3)==2 tezheng(i,j,3)=0。 tezheng(x,y,3)=0。 xxx=xxx2。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 Izz(x1,y)=0。Izz(x1,y+1)=0。Izz(x,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y)=0。Izz(x+1,y1)=0。Izz(x,j1)=0。Izz(x1,y1)=0。 end end for x=i12:i+12 for y=j+1:j+12 if tezheng(x,y,3)==2 tezheng(i,j,3)=0。 tezheng(x,y,3)=0。 xxx=xxx2。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 Izz(x1,y)=0。Izz(x1,y+1)=0。Izz(x,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y)=0。Izz(x+1,y1)=0。Izz(x,j1)=0。Izz(x1,y1)=0。 end end end end end endend%figure,imshow(I)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%消除毛刺for i=10:m10 for j=10:n10 if tezheng(i,j,3)==2 for x=i5:i+5 for y=j5:j+5 if tezheng(x,y,3)==1 tezheng(i,j,3)=0。 tezheng(x,y,3)=0。 xxx=xxx2。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 Izz(x1,y)=0。Izz(x1,y+1)=0。Izz(x,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y)=0。Izz(x+1,y1)=0。Izz(x,j1)=0。Izz(x1,y1)=0。 end end end end endend%figure,imshow(I)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%消除小橋for i=10:m10 for j=10:n10 if tezheng(i,j,3)==1 for x=i6:i+6 for y=j6:j1 if tezheng(x,y,3)==1 tezheng(i,j,3)=0。 tezheng(x,y,3)=0。 xxx=xxx2。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 Izz(x1,y)=0。Izz(x1,y+1)=0。Izz(x,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y)=0。Izz(x+1,y1)=0。Izz(x,j1)=0。Izz(x1,y1)=0。 end end end for x=i6:i+6 for y=j+1:j+6 if tezheng(x,y,3)==1 tezheng(i,j,3)=0。 tezheng(x,y,3)=0。 xxx=xxx2。 Izz(i1,j)=0。Izz(i1,j+1)=0。Izz(i,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j)=0。Izz(i+1,j1)=0。Izz(i,j1)=0。Izz(i1,j1)=0。 Izz(x1,y)=0。Izz(x1,y+1)=0。Izz(x,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y)=0。Izz(x+1,y1)=0。Izz(x,j1)=0。Izz(x1,y1)=0。 end end end end endendfigure,imshow(Izz)%%%%%%%%%%%%%%%%%%%%%%%%%匹配pipei=zeros(xxx,3)。ddd=1。for i=10:m10 for j=10:n10 if ddd~=xxx+1 if tezheng(i,j,3)~=0 pipei(ddd,1)=tezheng(i,j,1)。 pipei(ddd,2)=tezheng(i,j,2)。 pipei(ddd,3)=tezheng(i,j,3)。 ddd=ddd+1。 end end endend%計(jì)算p1與格pp=zeros(xxx,xxx1)。 qqq=1。for fff=1:xxx1 for i=qqq:xxx1 d=sqrt((pipei(qqq,1)pipei(i+1,1)).^2+(pipei(qqq,2)pipei(i+1,2)).^2)。 % jiao=arctan((pipei(i+1,2)pipei(qqq,2))/(pipei(i+1,1)pipei(qqq,1))。