all1027otherwaygai.m

一个用matlab编写的 关于笔画提取的仿真试验程序

close all;
clear;
%A2=imread('D:\a实验用\nh1.bmp','bmp');
%figure(9),imshow(A2);
tic;
A1=imread('shuang.bmp');%%%%玉1三g处大shuang开林生徊入z点生地旧程民江曰人小日答表快主九 s shuang su1 川1  非张 xiang 强 里目1原藏1 g 广  方 jia  jh领1须学1蒙1编男全
A=rgb2gray(A1);
k=25;
hui=0;
figure(1),imshow(A),title('原字符图像');
[m,n]=size(A);
% for i=1:(m)
% for j=1:(n)
%    if A(i,j) >200
%       A(i,j)=255;
%   end
% end
% end figure(988),imshow(A),title('原字符图像');
G=[0;-2;2;0];
H=[0.125;0.375;0.375;0.125];
%[A]=filtimag1231(A1);
%[cc]=fenggfuc(A);
%imwrite(cc,'qqying.bmp');figure(3),imshow(cc);
%Ac2=edge(A,'canny');
%figure(15),imshow(Ac2);
%G1=[0;0;-2;0;2;0;0;0];
%H1=[0.125;0;0.375;0;0.375;0;0.125;0];
%G=[0.025;0.475;0.475;0.025];
%H=[0.025;0.475;0.475;0.025];huit=round(hui/(m*n));
%aaa=[huit,huit/2,huit,huit/2];
%ss=aaa*G;
%ss2=aaa*H;
%%%%边缘扩展
Ak=zeros(m+k,n+k);
for i=6:(m+6-1)
    for j=6:(n+6-1)
        Ak(i,j)=A(i-5,j-5);
    end
end
%%%对行进行1变换
[SSh,WWh]=hanb(Ak);%%%%hanb及hanb2为函数
%%%对行进行2变换
[SSh2,WWh2]=hanb2(SSh);
%%%对行进行2变换[SSh3,WWh3]=hanb3(SSh2);
%%%对列进行1变换%%%%lieb及lieb2为函数
[SSl,WWl]=lieb(Ak);
%%%对列进行2变换
[SSl2,WWl2]=lieb2(SSl);
%%%对列进行2变换[SSl3,WWl3]=lieb3(SSl2);
%%对行进行寻找变换后的最大值点c5=WWh3;   [ahmax3]=findmaxminh(c5);
% c1=WWh;
% [ahmax1]=findmaxminh(c1);
   cc3=WWh2;
   [ahmax2]=findmaxminh(cc3);
  %%寻找1列变换后的最大值点
  % c2=WWl;
 %  [almax1]=findmaxminl(c2);
   cc4=WWl2;
   [almax2]=findmaxminl(cc4);
%%%%%%%%%%%%%%%mm=1;  c6=WWl3;   [almax3]=findmaxminl(c6);
hmax2=ahmax2;
lmax2=almax2;
hmax22=ahmax2;
lmax22=almax2;
%hmax1=ahmax1;
%lmax1=almax1;
%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%
%figure(2),imshow(ccall11);
%%%%%%%%%%%%%%
%统计字符的笔画宽度
[Ahen,Ashu]=bihua22(A);
[rowabh]=widthrow(Ashu,m,n);
bb1=mean2(rowabh);
[colabh]=widthcol(Ahen,m,n);
bb2=mean2(colabh);
if bb1 >0&&bb2>0
if bb1>bb2 
    bb=round(bb2);
else 
    bb=round(bb1);
end
elseif bb1>0
    bb=bb1;
else
    bb=bb2;
end
%%% 统计行或列的不为零的边缘点,若不为零则记录其点值并存于Amh和Aml中
[mch2,nch2]=size(hmax2);
[mcl2,ncl2]=size(lmax2);
for i=6:(m+5)
    ererhm=0;
    ehm=0;
    for j=6:(n+5)
        if (abs(hmax2(i,j))>0)
           ererhm=ererhm+abs(hmax2(i,j));
           ehm=ehm+1;
       end
   end
   if ehm>0
   km=ererhm/ehm;
   Amh(1,i)=km;
else
     Amh(1,i)=0;
 end
end

for j=6:(n+5)
    ererlm=0;
    elm=0;
    for i=6:(m+5)
        if (abs(lmax2(i,j))>0)
           ererlm=ererlm+abs(lmax2(i,j));
           elm=elm+1;
       end
   end
    if elm>0
   km=ererlm/elm;
   Aml(1,j)=km;
else 
    Aml(1,j)=0; 
end
end
%%%%%%%%%%%%%%%%以下为边缘点值的滤波
ath=zeros(m,n);
ath1=zeros(m,n);
ath2=zeros(m,n);
atl=zeros(m,n);
atl1=zeros(m,n);atl2=zeros(m,n);
for i=1:(m)
    for j=1:(n)
        if abs(hmax2(i+5,j+3))>100 
        ath(i,j)=(hmax2(i+5,j+3));%正侧的列边缘即竖
        ath1(i,j)=abs(hmax2(i+5,j+3));%所有的列边缘即竖
    end
    end
end
    for  j=1:(n)
    for    i=1:(m)
        if abs(lmax2(i+3,j+5))>100 %&& abs(lmax2(i+3,j+5))<=255  %%%%%%%%%???????????????????
        atl(i,j)=round(lmax2(i+3,j+5));%正侧的行边缘!!!!!!!!!横
            atl1(i,j)=abs(lmax2(i+3,j+5));%所有的行边缘
        end
    end
end
% figure(22),imshow(ath),title('正侧的列边缘即竖'); 
% figure(23),imshow(ath1),title('所有的列边缘即竖'); 
% figure(24),imshow(atl),title('正侧的行边缘'); 
% figure(25),imshow(atl1),title('所有的行边缘'); 
% figure(2),imshow(alie),title('alie竖'); 
atla=zeros(m,n);
atla1=zeros(m,n);atla2=zeros(m,n);
for i=1:(m)
    for j=1:(n)
        atla(i,j)=(double(atl(i,j))+double(ath(i,j)));%所有正侧的边缘
        atla1(i,j)=double(atl1(i,j))+double(ath1(i,j)); %所有边缘
        ath2(i,j)=double(ath1(i,j))-double(ath(i,j));%负侧的列边缘!!!!!
        atl2(i,j)=double(atl1(i,j))-double(atl(i,j));%负侧的行边缘
        atla2(i,j)=double(atl2(i,j))+double(ath2(i,j));%所有的负侧边缘
        %akit(i,j)=double(atl(i,j))+(double(ath1(i,j))-double(ath(i,j)));
        apie(i,j)=double(ath(i,j))+(double(atl(i,j)));
        ana(i,j)=double(ath(i,j))+(double(atl2(i,j)));
    end
end
%%%%%%%%%%%%%列补点
% figure(25),imshow(ana),title('捺边缘'); 
% figure(26),imshow(apie),title('撇边缘'); 


dddd=0;sumt=0;
aaah=ath;%正侧的每列边缘，竖
%figure(32),imshow(ath);
for kkk=1:bb
for    j=1:(n-1)
for  i=1:(m)
       t1=aaah(i,j);
        if  t1 >0
            dddd=dddd+1;
        elseif  dddd>(6)
            sumt=dddd;
         dddd=0;
          if A(i,j) >20
             sumt=sumt+1;
             t33=aaah(i-1,j);
          aaah(i,j)=aaah(i-1,j);
          end
       else 
         dddd=0;
      end
        end
end
    end

%figure(33),imshow(aaah);


%%%%%%%%%%%滤波
dd=0;sum1=0;
 for  j=1:(n)
     i=1;
    while(i<=m)
       t1=aaah(i,j);
        if  t1 >0 
          k9=i;
          while aaah(k9,j)>0
              dd=dd+1; 
              k9=k9+1;
              if k9 > m 
                  break;
              end
          end
          if dd <bb/3
              for k8=dd:-1:1;
                  tttt=k9-k8;
                aaah(k9-k8,j)=0;
              end
             i=i+dd-1;
             dd=0;
         else
             sum1=sum1+1;
             startshui(1,sum1)=i;%起始点行坐标
             startshuj(1,sum1)=j;% 起始点列号
             i=i+dd-1;
             endshui(1,sum1)=i;
             endshuj(1,sum1)=j;
        end
        dd=0;
       
    end
     i=i+1;  
 end
end
%%%%%%%%%%%%%%%%figure(090),imshow(aaah),title('alie竖'); 

 alie=aaah;
 for i=1:(sum1-1)
     tl1=startshuj(1,i);
     tl2=endshuj(1,i+1);
     tl3=abs(tl1-tl2)+1;
     tl4=(endshui(1,i)+1);
     tl5=(startshui(1,i+1)-1);
     tttl66=abs(tl4-tl5);
        if tl3<=2
            tttt=0;
            if tl4>tl5
          for j=tl4:-1:tl5
            if A( j ,tl1) >0
                tttt=tttt+1; 
             end        
         end 
         if tttt>=tttl66
              for j=tl4:-1:tl5
             aaah(j ,tl1)=255;
             ana(j ,tl2)=0;
         end
         end
     else
         for j=tl4:tl5
            if A( j ,tl1) >0
                tttt=tttt+1; 
             end        
         end
          if tttt>=tttl66
              for j=tl4:tl5
             aaah(j ,tl1)=255;ana(j ,tl2)=0;
             end
          end
      end
     end
 end

  %%%%%%滤除小于宽度的竖
 alie=aaah;

 for i=1:sum1
tss=startshui(1,i);
tes=endshui(1,i);
tssj=startshuj(1,i);
tlon=abs(tes-tss)+1;
     if  (tlon ==round(bb))||tlon <=round(bb*3/4)
        for s1=tss:tes
         alie( s1,tssj)=0;
     end
 end
end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure(44),imshow(alie),title('alie竖'); 
 for i=1:sum1-1
     for j=(i+1):(sum1)
        tss1=startshuj(1,i);
         tss2=startshuj(1,j); 
         tbin=abs(tss1-tss2);
         if tbin<=2
          temp1=0;
         temp2=0;
         temp3=0;
         temp4=0;  
              tssi1=startshui(1,i);
         tssi2=startshui(1,j);
          tesi1=endshui(1,i);
         tesi2=endshui(1,j);
             if tssi1 <tssi2
                temp1=tssi1;
                temp3=tssi2;
            else
                temp1=tssi2;
                temp3=tssi1;
            end
             if tesi1 >tesi2
                temp2=tesi1;
                temp4=tesi2;
            else
                temp2=tesi2;
                temp4=tesi1;
            end
       if temp1>0 && temp2 > 0&&(temp4>=temp1&&temp4>=temp3)%  &&(tesi1>=tesi2&&tssi1>=tssi2&&tssi1<=tesi2)
            for kc1=temp1:temp2
             alie(kc1,tss1)=255;
             alie(kc1,tss2)=0;
         end
       end
     end
     end
 end
%aaah=alie;
figure(2),imshow(alie),title('alie竖'); 
    %%%%%%%%%%%行补点
aaa=atl2;%负侧的每行边缘
dddd=0;sum2=0;
for kkkk=1:bb
 for  i=1:(m)
    for    j=n:-1:(2)
        
        t1=aaa(i,j);
        if  t1>0
            dddd=dddd+1;
        else
         if  dddd>=bb/2
             dddd=0;
             if A(i,j) >100
               aaa(i,j)=aaa(i,j+1);
              end
          else 
             dddd=0;  
         
     end  
    
       end

end
    end
end

end
%%%%%%%%%%%滤波

dd=0;
 for  i=1:(m)
     j=1;
    while(j<=n)
       t1=aaa(i,j);
        if  t1 >0 
          k9=j;
          while aaa(i,k9)>0
              dd=dd+1; 
              k9=k9+1;
              if k9 > n 
                  k9=k9-1;
                  break;
              end
          end
          if dd < bb 
              for k8=dd:-1:1;
                  tttt=k9-k8;
                aaa(i,k9-k8)=0;
            end
             j=j+dd-1;
             dd=0;
         else
             sum2=sum2+1;
             startheni(1,sum2)=i;%起始点行坐标
             starthenj(1,sum2)=j;% 起始点列号
            j=j+dd-1; 
             endheni(1,sum2)=i;
             endhenj(1,sum2)=j;
        end
        dd=0;
       
    end
     j=j+1;  
 end
  
end

 for i=1:(sum2-1)
     tl1=endheni(1,i);
     tl2=endheni(1,i+1);
     tl3=abs(tl1-tl2);
        if tl3<=1  
            tl4=(endhenj(1,i)+1);
           tl5=(starthenj(1,i+1)-1);
           tl6=abs(tl4-tl5);
           if tl6<=bb
          for j=tl4:tl5
            if A(tl1,j) >0
                 aaa(tl1,j )=aaa(tl1,j-1 );
             end        
         end
     end
     end
 end
 
 ahen=aaa;
for i=1:sum2
tsh=starthenj(1,i);
teh=endhenj(1,i);
tshi=startheni(1,i);
tlon=abs(teh-tsh)+1;
     if   tlon <bb
        for s1=tsh:teh
         ahen( tshi,s1)=0;
     end
 end
             
 end

figure(3),imshow(ahen),title('ahen横'); 
 %%%%%%%%%%%撇的提取

for    j=1:(n)
for  i=1:(m)
    na1(i,j)=ana(i,j)-aaah(i,j)-ahen(i,j);
end
end
%figure(4),imshow(na1),title('初捺'); 
 cunna=zeros(m,n);
 na=na1;
 tnc=0;longn(1,1)=0;
 for  i=1:(m)
     j=1;
    while(j<=n)
        t1=na(i,j);
        kk8=i;kk10=i;
         kk9=j;kk11=j;
        if  t1 >0 
            tb=1; cunna(kk8,kk9)=t1;
          [t2,kk8,kk9]=napos2(na,kk8,kk9);
        while t2>0 && kk9<n && kk8<m
          cunna(kk8,kk9)=t2; 
          na(kk8,kk9)=0;