车牌识别的matlab程序--(详细注释-并有使用注意点)

ysite=1； k1=1； for l=1:n2

for k=1:markrow3(ysite)%从顶边至第一个峰下降点扫描 if sbw1(k,l)==1 xdata(k1)=l； ydata(k1)=k； k1=k1+1； break； end end end

else %检测下边 ysite=n1；

if markrow4(n1) ==0

if markrow4(n1-1) ==maxw ysite= 0； %无下边 else

ysite= n1-1； end end if ysite ~=0 k1=1； for l=1:n2 k=m2；

while k>=markrow(ysite) %从底边至最后一个峰的上升点扫描

if sbw1(k,l)==1 xdata(k1)=l； ydata(k1)=k； k1=k1+1； break； end k=k-1；

end end end end

%(2)线性拟合，计算与x夹角

fresult = fit(xdata',ydata','poly1')； %poly1 Y = p1*x+p2 p1=fresult.p1；

angle=atan(fresult.p1)*180/pi； %弧度换为度，360/2pi, pi=3.14 %(3)旋转车牌图象

subcol = imrotate(subcol1,angle,'bilinear','crop')； %旋转车牌图象 sbw = imrotate(sbw1,angle,'bilinear','crop')；%旋转图像

figure,subplot(2,1,1),imshow(subcol)；title('车牌灰度子图')；%输出车牌旋转后的灰度图像标题显示车牌灰度子图

subplot(2,1,2),imshow(sbw)；title('')；%输出车牌旋转后的灰度图像 title(['车牌旋转角: ',num2str(angle),'度'] ,'Color','r')；%显示车牌的旋转角度

%Step8 旋转车牌后重新计算车牌水平投影，去掉车牌水平边框，获取字符高度

histcol1=sum(sbw)； %计算垂直投影 histrow=sum(sbw')； %计算水平投影

figure,subplot(2,1,1),bar(histcol1)；title('垂直投影(旋转后)')； subplot(2,1,2),bar(histrow)； title('水平投影(旋转后)')； figure,subplot(2,1,1),bar(histrow)； title('水平投影(旋转后)')； subplot(2,1,2),imshow(sbw)；title('车牌二值子图(旋转后)')； %去水平(上下)边框,获取字符高度 maxhight=max(markrow2)； findc=find(markrow2==maxhight)； rowtop=markrow(findc)；

rowbot=markrow(findc+1)-markrow1(findc+1)；

sbw2=sbw(rowtop:rowbot,:)； %子图为(rowbot-rowtop+1)行 maxhight=rowbot-rowtop+1； %字符高度(rowbot-rowtop+1)

%Step9 计算车牌垂直投影，去掉车牌垂直边框，获取车牌及字符平均宽度

histcol=sum(sbw2)； %计算垂直投影

figure,subplot(2,1,1),bar(histcol)；title('垂直投影(去水平边框后)')；%输出车牌的垂直投影图像

subplot(2,1,2),imshow(sbw2)； %输出垂直投影图像

title(['车牌字符高度： ',int2str(maxhight)],'Color','r')；%输出车牌字符高度

%对垂直投影进行峰谷分析

meancol=mean(histcol)；%求垂直投影的平均值 mincol=min(histcol)；%求垂直投影的平均值 levelcol=(meancol+mincol)/4；%求垂直投影的1/4 count1=0； l=1； for k=1:width

if histcol(k)<=levelcol count1=count1+1； else

if count1>=1

markcol(l)=k； %字符上升点

markcol1(l)=count1； %谷宽度(下降点至下一个上升点) l=l+1； end count1=0； end end

markcol2=diff(markcol)；%字符距离(上升点至下一个上升点) [m1,n1]=size(markcol2)； n1=n1+1； markcol(l)=width； markcol1(l)=count1；

markcol2(n1)=markcol(l)-markcol(l-1)；

%Step10 计算车牌上每个字符中心位置，计算最大字符宽度maxwidth l=0； for k=1:n1

markcol3(k)=markcol(k+1)-markcol1(k+1)；%字符下降点 markcol4(k)=markcol3(k)-markcol(k)； %字符宽度(上升点至下降点)

markcol5(k)=markcol3(k)-double(uint16(markcol4(k)/2))；%字符中心位置

end

markcol6=diff(markcol5)； %字符中心距离(字符中心点至下一个字符中心点)

maxs=max(markcol6)； %查找最大值，即为第二字符与第三字符中心距离

findmax=find(markcol6==maxs)； markcol6(findmax)=0；

maxwidth=max(markcol6)；%查找最大值，即为最大字符宽度

%Step11 提取分割字符,并变换为22行?14列标准子图 l=1；

[m2,n2]=size(subcol)； figure；

for k=findmax-1:findmax+5

cleft=markcol5(k)-maxwidth/2； cright=markcol5(k)+maxwidth/2-2； if cleft<1 cleft=1；

cright=maxwidth； end if cright>n2 cright=n2； cleft=n2-maxwidth； end

SegGray=sbw(rowtop:rowbot,cleft:cright)； SegBw1=sbw(rowtop:rowbot,cleft:cright)；

SegBw2 = imresize(SegBw1,[22 14])；%变换为22行?14列标准子图

subplot(2,n1,l),imshow(SegGray)； if l==7

title(['车牌字符宽度： ',int2str(maxwidth)],'Color','r')； end

subplot(2,n1,n1+l),imshow(SegBw2)； fname=strcat('c:\\work\\sam\\image',int2str(k),'.jpg')；%这块儿是自建模板，网上也有下载的，注意文件路径

imwrite(SegBw2,fname,'jpg') l=l+1； end

%Step12 将计算计算获取的字符图像与样本库进行匹配，自动识别出字符代码。

liccode=char(['0':'9' 'A':'Z' '粤桂海云贵川京津沪'])； %建立自动识别字符代码表

SubBw2=zeros(22,14)； l=1；

[m2,n2]=size(sbw)； for k=findmax-1:findmax+5

cleft=markcol5(k)-maxwidth/2； cright=markcol5(k)+maxwidth/2-2； if cleft<1 cleft=1；

cright=maxwidth； end if cright>n2 cright=n2；

cleft=n2-maxwidth； end

SegBw1=sbw(rowtop:rowbot,cleft:cright)；