matlab车牌识别程序代码.docx
《matlab车牌识别程序代码.docx》由会员分享，可在线阅读，更多相关《matlab车牌识别程序代码.docx（58页珍藏版）》请在冰点文库上搜索。
matlab车牌识别程序代码
%车牌识别程序主体
clc;
closeall;
clearall;
%==========================================================
%说明：
%
%
%===========================================================
%==============测定算法执行的时间，开始计时=================
tic%%%%%记录程序运行时间
%=====================读入图片================================
[fn,pn,fi]=uigetfile（'*.jpg','选择图片'）;
I=imread（[pnfn]）;
figure;
imshow（I）;
title（'原始图像'）;%显示原始图像
chepailujing=[pnfn]
I_bai=I;
[PY2,PY1,PX2,PX1]=caitu_fenge（I）;
%I=rgb2hsv（I）;
%[PY2,PY1,PX2,PX1]=caitu_tiqu（I,I_bai）;%用HSI模型识别蓝色，用rgb模型识别白色
%================分割车牌区域=================================
%===============车牌区域根据面积二次修正======================
[PY2,PY1,PX2,PX1,threshold]=SEC_xiuzheng（PY2,PY1,PX2,PX1）;
%==============更新图片=============================
Plate=I_bai（PY1:
PY2,PX1:
PX2,:
）;%使用caitu_tiqu
%==============考虑用腐蚀解决蓝色车问题=============
bw=Plate;figure,imshow（bw）;title（'车牌图像'）;%hsv彩图提取图像
%==============这里要根据图像的倾斜度进行选择这里选择的图片20090425686.jpg
bw=rgb2gray（bw）;figure,imshow（bw）;title（'灰度图像'）;
%================倾斜校正======================
qingxiejiao=rando_bianhuan（bw）
bw=imrotate（bw,qingxiejiao,'bilinear','crop'）;figure,imshow（bw）;title（'倾斜校正'）;%取值为负值向右旋转
%==============================================
bw=im2bw（bw,graythresh（bw））;%figure,imshow（bw）;
bw=bwmorph（bw,'hbreak',inf）;%figure,imshow（bw）;
bw=bwmorph（bw,'spur',inf）;%figure,imshow（bw）;title（'擦除之前'）;
'）;闭合运算bw=bwmorph（bw,'open',5）;%figure,imshow（bw）;title（'
bw=bwareaopen（bw,threshold）;figure,imshow（bw）;title（'擦除'）;
%==================加入进度条================================
%h=waitbar（0,'程序运行中，请稍等......'）
%fori=1:
10000
%waitbar（i/5000,h）
%end
%close（h）;
%wavplay（wavread（'程序运行中.wav'）,22000）;
%==========================================================
bw=~bw;figure,imshow（bw）;title（'擦除反色'）;
%=============对图像进一步裁剪，保证边框贴近字体===========
bw=touying（bw）;figure;imshow（bw）;title（'Y方向处理'）;
bw=~bw;
bw=bwareaopen（bw,threshold）;
bw=~bw;%figure,imshow（bw）;title（'二次擦除'）;
[y,x]=size（bw）;%对长宽重新赋值
%=================文字分割=================================
fenge=shuzifenge（bw,qingxiejiao）
[m,k]=size（fenge）;
%=================显示分割图像结果=========================
figure;
fors=1:
2:
k-1
subplot（1,k/2,（s+1）/2）;imshow（bw（1:
y,fenge（s）:
fenge（s+1）））;
end
%================给七张图片定位===============桂AV6388
han_zi=bw（1:
y,fenge
（1）:
fenge
（2））;
zi_mu=bw（1:
y,fenge（3）:
fenge（4））;
zm_sz_1=bw（1:
y,fenge（5）:
fenge（6））;
zm_sz_2=bw（1:
y,fenge（7）:
fenge（8））;
shuzi_1=bw（1:
y,fenge（9）:
fenge（10））;
shuzi_2=bw（1:
y,fenge（11）:
fenge（12））;
shuzi_3=bw（1:
y,fenge（13）:
fenge（14））;
%==========================识别====================================
%======================把修正数据读入==============================
xiuzhenghanzi=imresize（han_zi,[11055],'bilinear'）;
xiuzhengzimu=imresize（zi_mu,[11055],'bilinear'）;
xiuzhengzm_sz_1=imresize（zm_sz_1,[11055],'bilinear'）;
xiuzhengzm_sz_2=imresize（zm_sz_2,[11055],'bilinear'）;
xiuzhengshuzi_1=imresize（shuzi_1,[11055],'bilinear'）;
xiuzhengshuzi_2=imresize（shuzi_2,[11055],'bilinear'）;
xiuzhengshuzi_3=imresize（shuzi_3,[11055],'bilinear'）;
%============把0-9,A-Z以及省份简称的数据存储方便访问====================
hanzishengfen=duquhanzi（imread（'cpgui.bmp'）,imread（'cpguizhou.bmp'）,imread（'cpjing.bmp'）,imread（'cpsu.bmp'）,imread（'cpyue.bmp'））;
%因数字和字母比例不同。
这里要修改
shuzizimu=duquszzm（imread（'0.bmp'）,imread（'1.bmp'）,imread（'2.bmp'）,imread（'3.bmp'）,imread（'4.bmp'）,...
imread（'5.bmp'）,imread（'6.bmp'）,imread（'7.bmp'）,imread（'8.bmp'）,imread（'9.bmp'）,...
imread（'10.bmp'）,imread（'11.bmp'）,imread（'12.bmp'）,imread（'13.bmp'）,imread（'14.bmp'）,...
imread（'15.bmp'）,imread（'16.bmp'）,imread（'17.bmp'）,imread（'18.bmp'）,imread（'19.bmp'）,...
imread（'20.bmp'）,imread（'21.bmp'）,imread（'22.bmp'）,imread（'23.bmp'）,imread（'24.bmp'）,...
imread（'25.bmp'）,imread（'26.bmp'）,imread（'27.bmp'）,imread（'28.bmp'）,imread（'29.bmp'）,...
imread（'30.bmp'）,imread（'31.bmp'）,imread（'32.bmp'）,imread（'33.bmp'））;
zimu=
duquzimu（imread（'10.bmp'）,imread（'11.bmp'）,imread（'12.bmp'）,imread（'13.bmp'）,imread（'14.bmp'）,...
imread（'15.bmp'）,imread（'16.bmp'）,imread（'17.bmp'）,imread（'18.bmp'）,imread（'19.bmp'）,...
imread（'20.bmp'）,imread（'21.bmp'）,imread（'22.bmp'）,imread（'23.bmp'）,imread（'24.bmp'）,...
imread（'25.bmp'）,imread（'26.bmp'）,imread（'27.bmp'）,imread（'28.bmp'）,imread（'29.bmp'）,...
imread（'30.bmp'）,imread（'31.bmp'）,imread（'32.bmp'）,imread（'33.bmp'））;
shuzi=
duqushuzi（imread（'0.bmp'）,imread（'1.bmp'）,imread（'2.bmp'）,imread（'3.bmp'）,imread（'4.bmp'）,...
imread（'5.bmp'）,imread（'6.bmp'）,imread（'7.bmp'）,imread（'8.bmp'）,imread（'9.bmp'））;
%============================识别结果================================
i=1;%shibiezm_sz该函数识别数字有问题
jieguohanzi=shibiehanzi（hanzishengfen,xiuzhenghanzi）;shibiejieguo（1,i）=jieguohanzi;i=i+1;
jieguozimu=shibiezimu（zimu,xiuzhengzimu）;shibiejieguo（1,i）=jieguozimu;
i=i+1;
jieguozm_sz_1=shibiezm_sz（shuzizimu,xiuzhengzm_sz_1）;shibiejieguo（1,i）
=jieguozm_sz_1;i=i+1;
jieguozm_sz_2=shibiezm_sz（shuzizimu,xiuzhengzm_sz_2）;shibiejieguo（1,i）
=jieguozm_sz_2;i=i+1;
jieguoshuzi_1=shibieshuzi（shuzi,xiuzhengshuzi_1）;shibiejieguo（1,i）=jieguoshuzi_1;i=i+1;
jieguoshuzi_2=shibieshuzi（shuzi,xiuzhengshuzi_2）;shibiejieguo（1,i）=jieguoshuzi_2;i=i+1;
jieguoshuzi_3=shibieshuzi（shuzi,xiuzhengshuzi_3）;shibiejieguo（1,i）=jieguoshuzi_3;i=i+1;
%==========================对话框显示显示=============================================
shibiejieguo
===============后面的注释掉了（语音读出这一部分）%%%%========
msgbox（shibiejieguo,'识别结果'）;
%%=====================导出文本==================
%fid=fopen（'Data.xls','a+'）;
%fprintf（fid,'%s\r\n',shibiejieguo,datestr（now））;
%fclose（fid）;
%%===================读出声音===================
%duchushengyin（shibiejieguo）;
%================读取计时==========================
t=toc%%%记录程序运行时间
%============================================
自己编写的子函数
functionbw_fir=touying（imane_bw）
X_yuzhi=1;
[y,x]=size（imane_bw）;
Y_touying=（sum（（~imane_bw）'））';%往左边投影统计黑点
X_touying=sum（（~imane_bw））;%往下面投影
%找黑体边缘
Y_up=fix（y/2）;
Y_yuzhi=mean（Y_touying（（fix（y/2）-10）:
（fix（y/2）+10）,1））/1.6;
while（（Y_touying（Y_up,1）>=Y_yuzhi）&&（Y_up>1））%找到图片上边界
Y_up=Y_up-1;
end
Y_down=fix（y/2）;
while（（Y_touying（Y_down,1）>=Y_yuzhi）&&（Y_downY_down=Y_down+1;
end
%去除左边边框干扰
X_right=1;
if（X_touying（1,fix（x/14）））<=X_yuzhi
X_right=fix（x/14）
end
%找黑体边缘
bw_fir=imane_bw（Y_up:
Y_down,X_right:
x）;
自己编写的子函数
functionfenge=shuzifenge（imfenge,qingxiejiao）
[y,x]=size（imfenge）;
%===============用函数设定分割阈值========================
SS=x*y
ifSS<=20000
shedingyuzhi=4;
elseifSS>20000&&SS<=30000
shedingyuzhi=4;
elseifSS>30000&&SS<=50000
shedingyuzhi=4;
elseifSS>50000&&SS<=80000
shedingyuzhi=4;
else
shedingyuzhi=4;
end
ganrao=SS/100;
%=========================================================
%定义数组histogram存储垂直方向的黑点数
histogram=sum（~imfenge）;
%=================文字分割=============================
k=1;
forh=1:
x-1
if
（（histogram（1,h）<=shedingyuzhi）&&（histogram（1,h+1）>shedingyuzhi））||（（h==1）&&histogram（1,h）>shedingyuzhi）
fenge（1,k）=h;
k=k+1;
elseif
（（histogram（1,h）>shedingyuzhi）&&（histogram（1,h+1）<=shedingyuzhi））||（（h==x-1）&&histogram（1,h）>shedingyuzhi）
fenge（1,k）=h+1;
k=k+1;
end
end
k=k-1;%去掉多产生的一个K值
ifk<10
msgbox（'提取出错','警告'）;
wavplay（wavread（'提取出错.wav'）,22000）;
pause;
end
%==============================================
if
（sum（histogram（1,fenge（1,1）:
fenge（1,2）））fori=3:
k
fenge（1,i-2）=fenge（1,i）;
end
end
[mn]=size（fenge）;
ifn<14
msgbox（'提取出错','警告'）;
wavplay（wavread（'提取出错.wav'）,22000）;
pause;
end
fenge=fenge（1,1:
14）;
自己编写的子函数
functionshibiezm_sz=shibiezm_sz（shuzizimu,xiuzhengzm_sz）;
[y,x,z]=size（xiuzhengzm_sz）;
fork=1:
34
sum=0;
fori=1:
y
forj=1:
x
ifshuzizimu（i,j,k）==xiuzhengzm_sz（i,j）%统计黑白
sum=sum+1;
end
end
end
baifenbi（1,k）=double（sum/（x*y））;
end
chepai=find（baifenbi>=max（baifenbi））;
chepai=chepai（1,1）;
%===================数字直接返回=========================
if（chepai>=1）&&（chepai<=10）
zm_sz=chepai-1;
zm_sz=num2str（zm_sz）;
elseif（chepai>=11）&&（chepai<=34）
%==================字母对应序号转字母====================
%=========11-15ABCDE16-20FGHIJ21-25KLMNP26-30QRSTU31-35VWXYZ
ifchepai==11
zm_sz='A';
elseifchepai==12
zm_sz='B';
elseifchepai==13
zm_sz='C';
elseifchepai==14
zm_sz='D';
elseifchepai==15
zm_sz='E';
elseifchepai==16
zm_sz='F';
elseifchepai==17
zm_sz='G';
elseifchepai==18
zm_sz='H';
elseifchepai==19
zm_sz='J';
elseifchepai==20
zm_sz='K';
elseifchepai==21
zm_sz='L';
elseifchepai==22
zm_sz='M';
elseifchepai==23
zm_sz='N';
elseifchepai==24
zm_sz='P';
elseifchepai==25
zm_sz='Q';
elseifchepai==26
zm_sz='R';
elseifchepai==27
zm_sz='S';
elseifchepai==28
zm_sz='T';
elseifchepai==29
zm_sz='U';
elseifchepai==30
zm_sz='V';
elseifchepai==31
zm_sz='W';
elseifchepai==32
zm_sz='X';
elseifchepai==33
zm_sz='Y';
elseifchepai==34
zm_sz='Z';
end
end
shibiezm_sz=zm_sz;
自己编写的子函数
functionshibiezm_sz=shibiezm_sz（shuzizimu,xiuzhengzm_sz）;
[y,x,z]=size（xiuzhengzm_sz）;
fork=1:
34
sum=0;
fori=1:
y
forj=1:
x
ifshuzizimu（i,j,k）==xiuzhengzm_sz（i,j）%统计黑白
sum=sum+1;
end
end
end
baifenbi（1,k）=double（sum/（x*y））;
end
chepai=fi
升级会员 