【正文】 論文終于得以定稿。 四年寒窗，所收獲的不僅僅是愈加豐厚的知識，更重要的是在閱讀、實踐中所培養(yǎng)的思維方式、表達能力和廣闊視野。在此，衷心的感謝他們。 感謝對我論文進行評審的各位老師，感謝對論文的指導和提出的寶貴意見。最后，再次感謝所有關心我、幫助過我的老師、同學、朋友和親人。附錄主掃描程序clcclear allclose alloriginI=imread(39。)。I=originI。endI=double(I)。var=0。 endendM1=M/(m*n)。 endendvar1=var/(m*n)。for x=1:m for y=1:n if I(x,y)=M1 I(x,y)=150+sqrt(100*(I(x,y)M1)*(I(x,y)M1)/var1)。 end endendfigure, imshow(uint8(I))%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 分割M = 10。 L= n/M。var1=zeros(H,L)。 for y=1:L。var=0。 for j=1:M。 end end aveg1(x,y)=aveg/(M*M)。 for j=1:M。 end end var1(x,y)=var/(M*M)。Vmean=0。 Vmean=Vmean+var1(x,y)。Vmean1=Vmean/(H*L)。gtotle=0。vtemp=0。 gtotle=gtotle+aveg1(x,y)。 vtotle=vtotle+var1(x,y)。V1=vtotle/vtemp。gtotle1=0。vtemp1=0。 gtotle1=gtotle1+aveg1(x,y)。 vtotle1=vtotle1+var1(x,y)。V2=vtotle1/vtemp1。T1=G2。T3=G1100。for x=1:H for y=1:L if aveg1(x,y)T1 amp。 var1(x,y)T2 moban(x,y)=1。amp。 end endendfor x=2:H1 for y=2:L1 if moban(x,y)==1 if moban(x1,y) + moban(x1,y+1) +moban(x,y+1) + moban(x+1,y+1) + moban(x+1,y) + moban(x+1,y1) + moban(x,y1) + moban(x1,y1) =4 moban(x,y)=0。for x=1:H for y=1:L if moban(x,y)==1 for i=1:M for j=1:M I(i+(x1)*M,j+(y1)*M)=G1。 end end end endendfigure, imshow(uint8(I)) %%%%%%%%%%%%%%%%%%%%%%二值化temp=(1/9)*[1 1 1。1 1 1]。 In=zeros(m,n)。 for b=2:n1。 endendI=In。for x=5:m5。 sum1=I(x,y4)+I(x,y2)+I(x,y+2)+I(x,y+4)。 sum3=I(x2,y+2)+I(x4,y+4)+I(x+2,y2)+I(x+4,y4)。 sum5=I(x2,y)+I(x4,y)+I(x+2,y)+I(x+4,y)。 sum7=I(x4,y4)+I(x2,y2)+I(x+2,y+2)+I(x+4,y+4)。 sumi=[sum1,sum2,sum3,sum4,sum5,sum6,sum7,sum8]。 summin=min(sumi)。 b=summ/8。 else sumf = summax。 else Im(x,y)=255。 endend for i=1:m for j =1:n if (Icc(i,j)==128) Icc(i,j)=0。 end。 In=Im。 end end if Im(i,j)==0 if Im(i1,j) + Im(i1,j+1) +Im(i,j+1) + Im(i+1,j+1) + Im(i+1,j) + Im(i+1,j1) + Im(i,j1) + Im(i1,j1) =7 In(i,j)=1。 end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%細化figure,imshow(In)。figure,imshow(I)。Im=I。for i=2:m1 for j=2:n1 if Im(i,j)==1 a = 0。 end if Im(i1,j+1) ~= Im(i,j+1) a = a + 1。 end if Im(i+1,j+1) ~= Im(i+1,j) a = a + 1。 end if Im(i+1,j1) ~= Im(i,j1) a = a + 1。 end if Im(i1,j1) ~= Im(i1,j) a = a + 1。 tezheng(i,j,2)=j。 xxx=xxx+1。 I(i1,j)=1。I(i,j+1)=1。I(i+1,j)=1。I(i,j1)=1。 %%%%%%%%%%%%%%1表示分叉點 0表示端點 end if a==2 %%%% 端點判斷 tezheng(i,j,1)=i。 tezheng(i,j,3)=2。 I(i,j)=0。I(i1,j+1)=1。I(i+1,j+1)=1。I(i+1,j1)=1。I(i1,j1)=1。for i=1:m %%%%%%%%%%%%%%%%%%%%%%%%去除邊界點 for j=1:n if tezheng(i,j,3)~=0 if i170||i+17m|| j170||j+17n tezheng(i,j,3)=0。 %%%消除端點 Izz(i1,j)=0。Izz(i,j+1)=0。Izz(i+1,j)=0。Izz(i,j1)=0。 else a=0。c=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)。 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)。 if ( a=1 amp。 b=1 amp。 c=1 amp。 d=1) tezheng(i,j,3)=0。 Izz(i1,j)=0。Izz(i,j+1)=0。Izz(i+1,j)=0。Izz(i,j1)=0。 end end end endend%%%%%%%%%%%%%%%%%%%%%%%%%%消除斷點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。 xxx=xxx2。Izz(i1,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j1)=0。Izz(i1,j1)=0。Izz(x1,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y1)=0。Izz(x1,y1)=0。 tezheng(x,y,3)=0。 Izz(i1,j)=0。Izz(i,j+1)=0。Izz(i+1,j)=0。Izz(i,j1)=0。 Izz(x1,y)=0。Izz(x,y+1)=0。Izz(x+1,y)=0。Izz(x,j1)=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。 xxx=xxx2。Izz(i1,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j1)=0。Izz(i1,j1)=0。Izz(x1,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y1)=0。Izz(x1,y1)=0。 tezheng(x,y,3)=0。 Izz(i1,j)=0。Izz(i,j+1)=0。Izz(i+1,j)=0。Izz(i,j1)=0。 Izz(x1,y)=0。Izz(x,y+1)=0。Izz(x+1,y)=0。Izz(x,j1)=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。 xxx=xxx2。Izz(i1,j+1)=0。Izz(i+1,j+1)=0。Izz(i+1,j1)=0。Izz(i1,j1)=0。Izz(x1,y+1)=0。Izz(x+1,y+1)=0。Izz(x+1,y1)=0。Izz(x1,y1)=0。ddd=1。 pipei(ddd,2)=tezheng(i,j,2)。 ddd=ddd+1。 qqq=1。 % jiao=arctan((pipei(i+1,2)pipei(qqq,2))/(pipei(i+1,1)pipei(qqq,1)