p0803.m

本文是采用matlab开发的识别数字的识别软件

I0=imread('8.jpg');% 必须为二值图像
I=im2bw(I0,0.4);
[y0 x0]=size(I);
Range=sum((~I)');
Hy=0;
for j=1:y0
    if (Range(j)>=1)
        Hy=Hy+1;
    end
end
RangeX=sum((~I));
Wx=0;
for i=1:x0
    if (RangeX(i)>=1)
        Wx=Wx+1;
    end
end
Amp=24/Hy;             % 将文字图像归一化到24像素点的高度。
I=imresize(I,Amp);
[y x]=size(I);
%I=bwmorph(~I,'skel',Inf);
%I=~I;
tic
%====== 基本结构 =======%
% 第一类：竖（V）；左斜（L）；右斜（R）；突变（P）
% 第二类：左半圆弧（C）；右半圆弧（Q）
% 的三类：结构待定（T）；
%=====================================%
Left=zeros(1,y);         % 左端轮廓检测
for j=1:y
        i=1;
    while ((i<=x)&&(I(j,i)==1))
        i=i+1;
    end
    if (i<=x)          
        Left(j)=i;
    end    
end
for j=1:y-1
    LeftD(j)=Left(j+1)-Left(j);
end
%========== 结构特征提取 =============%
j=1;
while ((Left(j)<1)&&(j<y))
       j=j+1;       
end 
Y1=j;

j=y;
while ((Left(j)<1)&&(j>1))
       j=j-1;
end 
Y2=j-1;      % 去掉急剧变化的两端

       
%============== 右边 ==================%
Right=zeros(1,y);         % 左端轮廓检测
for j=1:y
        i=x;
    while ((i>=1)&&(I(j,i)==1))
        i=i-1;
    end
    if (i>=1)          
        Right(j)=i;
    end    
end
for j=1:y-1
    RightD(j)=Right(j+1)-Right(j);
end

%=====================================%
Top=zeros(1,x);         % 顶端轮廓检测
for i=1:x
        j=1;
    while ((j<=y)&&(I(j,i)==1))
        j=j+1;
    end
    if (j<=y)          
        Top(i)=j;
    end    
end
for i=1:x-1
    TopD(i)=Top(i+1)-Top(i);
end
%==============================%
i=1;
while ((Top(i)<1)&&(i<x))
       i=i+1;       
end 
X1=i;

i=x;
while ((Top(i)<1)&&(i>1))
       i=i-1;
end 
X2=i-1;      % 去掉急剧变化的两端

%===================================%
Bottom=zeros(1,x);      % 底部轮廓检测
for i=1:x
        j=y;
    while ((j>=1)&&(I(j,i)==1))
        j=j-1;
    end
    if (j>=1)
        Bottom(i)=j;
    end    
end
for i=1:x-1
    BottomD(i)=Bottom(i+1)-Bottom(i);
end

%========== 数字 1 的宽度特征 =========%
Width=zeros(1,y);
for j=1:y
    Width(j)=Right(j)-Left(j);
end    
W=max(Width);

Po=0; % 用于检测笔划
Ne=0;
NS=0;%笔划数
for i=X1+4:X2-4
    for j=1:y-1
        if ((I(j+1,i)-I(j,i))>0) % 由黑到白
            Po=Po+1;
            if ((Po>=2)&&(j<=fix(0.7*y)))
                Po=3;
            end    
        else if  ((I(j+1,i)-I(j,i))<0) % 由白到黑 
                  Ne=Ne+1;
                  if ((Ne>=2)&&(j<=fix(0.7*y)))
                       Ne=3;
                  end       
              end
          end                   
      end
      NS=[NS max(Po,Ne)];
      Po=0; 
      Ne=0;
  end
 Comp=max(NS); 
%========== 轮廓结构特征提取 ==========%
if (min(W,Wx)>10)
    StrokeT=StrDetect01(TopD,X1,X2,3,6); % 顶部基本结构检测
    StrokeL=StrDetect01(LeftD,Y1,Y2,3,5); % 左边基本结构检测
    StrokeR=StrDetect01(RightD,Y1,Y2,3,5); % 右边基本结构检测
    StrokeB=StrDetect01(BottomD,X1,X2,3,6); % 底部基本结构检测   
%========== 识别 ==========%
    Digit=Recognition(StrokeT,StrokeL,StrokeR,StrokeB,Comp)
else
    Digit='1'
end    

t=toc
%======= 显示 ======%
px=(1:x);
py=(1:y);

figure(1);
imshow(I);

figure(2);
plot(Left);grid
title('左轮廓');

figure(3);
plot(Top);grid
title('上轮廓');

figure(4);
plot(Right);grid
title('右轮廓');

figure(5);
plot(Width);grid
title('宽度');