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('寬度');