基于网格特征与外围特征的汉字识别方法matlab程序

%求汉字图像外围特征 function m_peripheryft=CrossPeiphery(I) [m,n]=size(I); count=1; m_peripheryft=zeros(1,64); %搜索方向为从左向右，水平。 for k=1:8 %图像分块后，在每块的每行一次外围和二次外围特征 sum1=zeros(1,m/8); sum2=zeros(1,m/8); %每块的一次外围和二次外围特征 p1=0; p2=0; for i=(k-1)*(m/8)+1:k*(m/8) for j=1:n if I(i,j)==0 sum1(i-(k-1)*(m/8))=sum1(i-(k-1)*(m/8))+1; else while(j<n+1&I(i,j)==1) j=j+1; end while(j<n+1&I(i,j)==0) j=j+1; sum2(i-(k-1)*(m/8))=sum2(i-(k-1)*(m/8))+1; end if j>n-1 sum2(i-(k-1)*(m/8))=0; end break; end end end for i=1:m/8 p1=p1+sum1(i); end for i=1:m/8 p2=p2+sum2(i); end m_peripheryft(count)=p1/(m*n); count=count+1; m_peripheryft(count)=p2/(m*n); count=count+1; end %搜索方向是从上向下，竖直。 for k=1:8 %图像分块后，在每块的每行一次外围和二次外围特征 sum1=zeros(1,n/8); sum2=zeros(1,n/8); %每块的一次外围和二次外围特征 p1=0; p2=0; for j=(k-1)*(n/8)+1:k*(n/8) for i=1:m if I(i,j)==0 sum1(j-(k-1)*(n/8))=sum1(j-(k-1)*(n/8))+1; else while(i<m+1&I(i,j)==1) i=i+1; end while(i<m+1&I(i,j)==0) i=i+1; sum2(j-(k-1)*(n/8))=sum2(j-(k-1)*(n/8))+1; end if i>m-1 sum2(j-(k-1)*(n/8))=0; end break; end end end for i=1:n/8 p1=p1+sum1(i); end for i=1:n/8 p2=p2+sum2(i); end m_peripheryft(count)=p1/(m*n); count=count+1; m_peripheryft(count)=p2/(m*n); count=count+1; end %搜索方向为从右向左，水平。 for k=1:8 %图像分块后，在每块的每行一次外围和二次外围特征 sum1=zeros(1,m/8); sum2=zeros(1,m/8); %每块的一次外围和二次外围特征 p1=0; p2=0; for i=(k-1)*(m/8)+1:k*(m/8) for j=n:-1:1 if I(i,j)==0 sum1(i-(k-1)*(m/8))=sum1(i-(k-1)*(m/8))+1; else while(j>0&I(i,j)==1) j=j-1; end while(j>0&I(i,j)==0) j=j-1; sum2(i-(k-1)*(m/8))=sum2(i-(k-1)*(m/8))+1; end if j<2 sum2(i-(k-1)*(m/8))=0; end break; end end end for i=1:m/8 p1=p1+sum1(i); end for i=1:m/8 p2=p2+sum2(i); end m_peripheryft(count)=p1/(m*n); count=count+1; m_peripheryft(count)=p2/(m*n); count=count+1; end %搜索方向是从下向上，竖直。 for k=1:8 %图像分块后，在每块的每行一次外围和二次外围特征 sum1=zeros(1,n/8); sum2=zeros(1,n/8); %每块的一次外围和二次外围特征 p1=0; p2=0; for j=(k-1)*(n/8)+1:k*(n/8) for i=m:-1:1 if I(i,j)==0 sum1(j-(k-1)*(n/8))=sum1(j-(k-1)*(n/8))+1; else while(i>0&I(i,j)==1) i=i-1; end while(i>0&I(i,j)==0) i=i-1; sum2(j-(k-1)*(n/8))=sum2(j-(k-1)*(n/8))+1; end if i<2 sum2(j-(k-1)*(n/8))=0; end break; end end end for i=1:n/8 p1=p1+sum1(i); end for i=1:n/8 p2=p2+sum2(i); end m_peripheryft(count)=p1/(m*n); count=count+1; m_peripheryft(count)=p2/(m*n); count=count+1; end m_peripheryft=m_peripheryft./norm(m_peripheryft); end