【正文】 FP_OrientationMap = DrawOrientationMap(FP_Orientation)。 figure。imshow(FP_OrientationMap,39。InitialMagnification39。,39。fit39。)。FP_Enhancement = Enhancement(FP_Normalized,FP_Orientation,W)。imshow(FP_Enhancement)。figure。 FP_Binarized = Binarized(FP_Enhancement,W)。imshow(FP_Binarized)。 figure。 FP_Thined = Thined(FP_Binarized)。imshow(FP_Thined)。 FP_FeaturePointsMatrix = GetFeaturePoints(FP_Thined)。FP_DelFeaturePointsMatrix = DelFeaturePoints(FP_FeaturePointsMatrix,FP_Orientation,W)。附錄B 提取特征點仿真程序function FP_FeaturePointsMatrix = GetFeaturePoints(FP_Thining)Z = double(FP_Thining)/255。[len,wid] = size(FP_Thining)。FP_FeaturePointsMatrix = zeros(len,wid)。Cn = zeros(len,wid)。Sn = zeros(len,wid)。 P = zeros(1,9)。for i = 2:len1 for j = 2:wid1 P(1) = Z(i1,j1)。 P(2) = Z(i1,j)。 P(3) = Z(i1,j+1)。 P(4) = Z(i,j+1)。 P(5) = Z(i+1,j+1)。 P(6) = Z(i+1,j)。 P(7) = Z(i+1,j1)。 P(8) = Z(i,j1)。 P(9) = P(1)。 for m = 1:8 Cn(i,j) = Cn(i,j)+abs(P(m+1)P(m))。 Sn(i,j) = Sn(i,j)+P(m)。 end endend for i = 1:len for j = 1:wid if Z(i,j) == 1 if Cn(i,j) == 2 amp。amp。 Sn(i,j) == 1 FP_FeaturePointsMatrix(i,j) = 1。 % endpoint elseif Cn(i,j) == 6 amp。amp。 Sn(i,j) == 3 FP_FeaturePointsMatrix(i,j) = 3。 % points on line else FP_FeaturePointsMatrix(i,j) = 2。 % forkpoint end end endend figure。hold on。for i = 1:len for j = 1:wid if FP_FeaturePointsMatrix(i,j) == 1 plot(j,len+1i,39。ro39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif FP_FeaturePointsMatrix(i,j) == 2 plot(j,len+1i,39。k39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif FP_FeaturePointsMatrix(i,j) == 3 plot(j,len+1i,39。bx39。,39。MarkerSize39。,15,39。LineWidth39。,2)。 end endend附錄C 剔除偽特征點仿真程序function FP_DelFeaturePointsMatrix = DelFeaturePoints(FP_FeaturePointsMatrix,FP_Orientation,W)[len,wid] = size(FP_FeaturePointsMatrix)。M = FP_FeaturePointsMatrix。 B = 20。for i = 1:B for j = 1:wid if M(i,j) == 1 || M(i,j) == 3 M(i,j) = 2。 end endendfor i = lenB+1:len for j = 1:wid if M(i,j) == 1 || M(i,j) == 3 M(i,j) = 2。 end endendfor i = 1:len for j = 1:B if M(i,j) == 1 || M(i,j) == 3 M(i,j) = 2。 end endendfor i = 1:len for j = widB+1:wid if M(i,j) == 1 || M(i,j) == 3 M(i,j) = 2。 end endendfigure。hold on。for i = 1:len for j = 1:wid if M(i,j) == 1 plot(j,len+1i,39。ro39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif M(i,j) == 2 plot(j,len+1i,39。k39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif M(i,j) == 3 plot(j,len+1i,39。bx39。,39。MarkerSize39。,15,39。LineWidth39。,2)。 end endend D1 = 5。D2 = 10。D3 = 5。D4 = 10。for i = 1:len for j = 1:wid if M(i,j) == 1 a = i。 % delete short line b = j。 N = FP_FeaturePointsMatrix。 flag = 1。 while(flag) if (N(a1,b1) == 2) || (N(a1,b1) == 1) N(a,b) = 0。 a = a1。 b = b1。 flag = 1。 elseif (N(a1,b) == 2) || (N(a1,b) == 1) N(a,b) = 0。 a = a1。 b = b。 flag = 1。 elseif (N(a1,b+1) == 2) || (N(a1,b+1) == 1) N(a,b) = 0。 a = a1。 b = b+1。 flag = 1。 elseif (N(a,b+1) == 2) || (N(a,b+1) == 1) N(a,b) = 0。 a = a。 b = b+1。 flag = 1。 elseif (N(a+1,b+1) == 2) || (N(a+1,b+1) == 1) N(a,b) = 0。 a = a+1。 b = b+1。 flag = 1。 elseif (N(a+1,b) == 2) || (N(a+1,b) == 1) N(a,b) = 0。 a = a+1。 b = b。 flag = 1。 elseif (N(a+1,b1) == 2) || (N(a+1,b1) == 1) N(a,b) = 0。 a = a+1。 b = b1。 flag = 1。 elseif (N(a,b1) == 2) || (N(a,b1) == 1) N(a,b) = 0。 a = a。 b = b1。 flag = 1。 else flag = 0。 end end T = (ai)/(bj)。 if T = 0 theta1 = atan(T)。 else theta1 = atan(T)+pi。 end I = floor((i1)/W)+1。 J = floor((j1)/W)+1。 theta2 = (FP_Orientation(I,J)1)*pi*。 delta = abs(theta1theta2)。 if (sqrt((ia)^2+(jb)^2) D1) amp。amp。 (delta *pi) M(i,j) = 2。 M(a,b) = 2。 end end endendfigure。hold on。for i = 1:len for j = 1:wid if M(i,j) == 1 plot(j,len+1i,39。ro39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif M(i,j) == 2 plot(j,len+1i,39。k39。,39。MarkerSize39。,10,39。LineWidth39。,2)。 elseif M(i,j) == 3 plot(j,len+1i,39。bx39。,39。MarkerSize39。,15,39。LineWidth39。,2)。 end endend N = M。D = D3。for i = 1:len for j = 1:wid if M(i,j) == 3 Block = N(iD:i+D,jD:j+D)。