【正文】 count = count + 1。 end end end end for i=1:count 1 x(i) =txy(i, 1)。 y(i)= txy(i, 2)。end imshow(double(thin))。hold on。plot(x,y,39。.39。)。%cut函數(shù)function txy=cut(thin,txy)s(8,8)=0。delta(8,8)=0。n=size(txy,1)。for i=1:8 for j=1:8 mp{i,j}=thin(1+31*(i1):31+31*(i1),1+31*(j1):31+31*(j1))。 s(i,j)=sum(sum(mp{i,j}))/(31*31)。 mp{i,j}=(mp{i,j}s(i,j)).^2。 delta(i,j)=sum(sum(mp{i,j}))。 if delta(i,j)=70 for k=1:nif (txy(k,1)=1+31*(i1)amp。amp。txy(k,1)=31+31*(i1)amp。amp。txy(k,2)=1+31*(j1)amp。amp。txy(k,2)=31+31*(j1)amp。amp。txy(k,3)==2) txy(k,:)=[0,0,0]。 end end end endendtxy=txy(find(txy(:,1)),:)。plot(txy(:,1),txy(:,2),39。ro39。)。附錄C 圖像特征點(diǎn)代碼%single_point函數(shù)function [pxy2,error]=single_point(txy,r)error=0。x=txy(:,1)。y=txy(:,2)。n=length(x)。d(1:n,1:n)=0。for j=1:n for i=1:n if (i~=j) d(i,j)=sqrt((x(i)x(j))^2+(y(i)y(j))^2)。 else d(i,j)=2*r。 end endend[a,b]=min(d)。c=find(ar)。pxy2=txy(c,:)。pxy2=pxy2(find(pxy2(:,3)==2),:)。t=size(pxy2,1)。if t==0 error=1else plot(x,y,39。b.39。)。 hold on plot(pxy2(:,1),pxy2(:,2),39。r.39。)。 end%walk函數(shù)function [error,a,b]=walk(thin,x0,y0,num) error=0。thin(y0,x0)=0。t1=0。for n=1:num if error==1 break。 else x=x0。 y=y0。 for x=x01:x0+1 if error==1 break。 else for y=y01:y0+1 t1=sum(sum(thin(y01:y0+1,x01:x0+1)))。 if (t1==0||t1=2) error=1。 a=x0。 b=y0。 break。 else if (thin(y,x)==1amp。amp。(xx0)^2+(yy0)^2~=0) if (t1=2 ) error=1。 break 。 else thin(y,x)=0。 x0=x。 y0=y。 a=x0。 b=y0。 plot(x0,y0,39。r.39。) end end end end end end endend%last1函數(shù)function [pxy3,error2]=last1(thin,r,txy,num)error=0。[pxy2,error]=single_point(txy,r)。n=size(pxy2,1)。l=1。error2=0。for i=1:n [error,a,b]=walk(thin,pxy2(i,1),pxy2(i,2),num)。 if error~=1 pxy3(l,1)= pxy2(i,1)。 pxy3(l,2)= pxy2(i,2)。 pxy3(l,3)= pxy2(i,3)。 l=l+1。 error2=0。 plot(pxy2(i,1) ,pxy2(i,2),39。r+39。)。 endend附錄D 特征點(diǎn)匹配代碼%distance函數(shù)程序function d=distance(x0,y0,num,thin)num2=fix(num/5)。for i=1:num2 [error,a,b]=walk(thin,x0,y0,5*i)。 if error~=1 d(i)=sqrt((ax0)^2+(by0)^2)。 else break。 endend%find_point函數(shù)function pxy=find_point(x0,y0,txy,num)x=txy(:,1)。y=txy(:,2)。n=length(x)。l(1,n)=0。lnn=1。pxy(num,:)=[0,0,0]。for i=1:n l(i)=sqrt((x(i)x0)^2+(y(i)y0)^2)。endll=sort(l)。for i=1:num xiao=ll(i+lnn)。nn=find(l==xiao)。lnn=length(nn)。 pxy(i,:)=[x(nn(1)),y(nn(1)),txy(nn(1),3)]。endplot(x0,y0,39。bo39。)。x0。y0。hold onplot(pxy(:,1),pxy(:,2),39。ro39。)。%主函數(shù)mainclose all。ticclear。thin1=tuxiangyuchuli(39。39。)。thin2=tuxiangyuchuli(39。39。)。figure。txy1=point(thin1)。txy2=point(thin2)。[w1,txy1]=guanghua(thin1,txy1)。[w2,txy2]=guanghua(thin2,txy2)。thin1=w1。thin2=w2。txy1=cut(thin1,txy1)。txy2=cut(thin2,txy2)。[pxy31,error2]=last1(thin1,8,txy1,60)[pxy32,error2]=last1(thin2,8,txy2,60)error=1。 num=20。cxy1=pxy31。 cxy2=pxy32。 d1=distance(cxy1(1,1),cxy1(1,2),num,thin1)。 d2=distance(cxy2(1,1),cxy2(1,2),num,thin2)。 f=(sum(abs((d1./d2)1)))。 if ff= error=0。 else error=1。 end c11=find_point(cxy1(1,1),cxy1(1,2),txy1,1)。 c12=find_point(cxy1(1,1),cxy1(1,2),txy1,2)。 c21=find_point(cxy2(1,1),cxy2(1,2),txy2,1)。 c22=find_point(cxy2(1,1),cxy2(1,2),txy2,2)。 cxy1(2,:)=c11。 cxy1(3,:)=c12(2,:)。 cxy2(2,:)=c21。 cxy2(3,:)=c22(2,:)。 x11=cxy1(1,1)。 y11=cxy1(1,2)。 x12=cxy1(2,1)。 y12=cxy1(2,2)。 x13=cxy1(3,1)。 y13=cxy1(3,2)。 x21=cxy2(1,1)。 y21=cxy2(1,2)。 x22=cxy2(2,1)。 y22=cxy2(2,2)。 x23=cxy2(3,1)。 y23=cxy2(3,2)。 dd1(1)=juli(x11,y11,x12,y12)。 dd1(2)=juli(x12,y12,x13,y13)。 dd1(3)=juli(x13,y13,x11,y11)。 dd2(1)=juli(x21,y21,x22,y22)。 dd2(2)=juli(x22,y22,x23,y23)。 dd2(3)=juli(x23,y23,x21,y21)。 ff=(sum(abs((dd1./dd2)1))) if ff=1 error=0。 else error=1。 end cxy1(2:41,:)=find_point(pxy31(1,1),pxy31(1,2),txy1,40)。 cxy2(2:41,:)=find_point(pxy32(1,1),pxy32(1,2),txy2,40)。 f11=length(find(cxy1(:,3)==2))。 f12=length(find(cxy1(:,3)==6))。 f21=length(find(cxy2(:,3)==2))。 f22=length(find(cxy2(:,3)==6))。 fff=abs(f11f21)/(f11+f12) toc44