车牌识别附源代码

1、车牌识别电子 1301 孙洪江 2013234020113一、目的与要求车牌定位系统的目的在于正确获取整个图像中车牌的区域，并识别出车牌号。通过设计实现车牌识别系统，能够提高学生分析问题和解决问题的能力， 还能培养一定的科研能力二、设计原理：牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、 摄像设备、照明设 备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、 车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图 像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包 括

2、车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时 触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位 出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。三、详细设计步骤：为了进行牌照识别，需要以下几个基本的步骤：a. 牌照定位，定位图片中的牌照位置；b. 牌照字符分割，把牌照中的字符分割出来；c. 牌照字符识别，把分割好的字符进行识别，最终组成牌照号码。牌照识别过程中，牌照颜色的识别依据算法不同，可能在上述不同步骤实现， 通常与牌照识别互相配合、互相验证。(1 ) 牌照定位：自然环境下，汽车图像背景复杂、光照不均匀，如何在自然背景中准

3、确地确定牌照区域是整个识别过程的关键。首先对采集到的视频图像进行 大范围相关搜索，找到符合汽车牌照特征的若十区域作为候选区， 然后对 这些侯选区域做进一步分析、评判，最后选定一个最佳的区域作为牌照区 域，并将其从图象中分割出来。(2)牌照字符分割完成牌照区域的定位后，再将牌照区域分割成单个字符，然后进行识别。字 符分割J股采用垂直投影法。由丁字符在垂直方向上的投影必然在字符问或字符 内的间隙处取得局部最小值的附近，并且这个位置应满足牌照的字符书写格式、 字符、尺寸限制和一些其他条件。利用垂直投影法对复杂环境下的汽车图像中的 字符分割有较好的效果。(3)牌照字符识别：字符识别方法目前主要有基丁模

4、板匹配算法和基丁人工神经网络算法。基丁模板匹配算法首先将分割后的字符二值化，并将其尺寸大小缩放为字符数据库中 模板的大小，然后与所有的模板进行匹配，最后选最佳匹配作为结果。基丁人工 神经元网络的算法有两种：一种是先对待识别字符进行特征提取，然后用所获得特征来训练神经网络分配器；另一种方法是直接把待处理图像输入网络，由网络自动实现特征提取直至识别出结果。 实际应用中，牌照识别系统的识别率与牌照 质量和拍摄质量密切相关。牌照质量会受到各种因素的影响，如生锈、污损、油 漆剥落、字体褪色、牌照被遮挡、牌照倾斜、高亮反光、多牌照、假牌照等等； 实际拍摄过程也会受到环境亮度、拍摄亮度、车辆速度等等因素的影

5、响。这些影 响因素不同程度上降低了牌照识别的识别率，也正是牌照识别系统的困难和挑战 所在。为了提高识别率，除了不断的完善识别算法，还应该想办法克服各种光照 条件，使采集到的图像最利丁识别。四、各个功能模块的主要实现程序1.载入车牌图像：I=imread('car1.jpg');figure(1),imshow(I);title('original image'); %务车牌的原图显示出来，结果如下：2.将彩图转换为灰度图并绘制直方图:I1=rgb2gray(I);%将彩图转换为灰度图figure(2),subplot(1,2,1),imshow(I1);titl

6、e('gray image');figure(2),subplot(1,2,2),imhist(I1) ;title('结果如下所示：灰度图直方图');%绘制灰度图的直方图3.用roberts算子进行边缘检测：I2=edge(I1,'roberts',0.18,'both');% 选择阈值 0.18,用 roberts 算子进行边缘检测figure(3),imshow(I2);title('roberts operator edge detection image');结果如下：4. 图像实施腐蚀操作：se=1;1

7、;1;I3=imerode(I2,se);%对图像实施腐蚀操作，即膨胀的反操作 figure(4),imshow(I3);title('corrosion image');5. 平滑图像se=strel('rectangle',25,25);%构造结构元素以正方形构造一个seI4=imclose(I3,se);figure(5),imshow(I4);title('smothing image');结果如下所示：D BW、giMk 口 Wig 日理*'! 4 -r trH以 n 6. 删除二值图像的小对象I5=bwareaopen(I4,

8、2000);%去除聚团灰度值小丁 2000的部分figure(6),imshow(I5);title('remove the small objects'); %用 imshow 函数显小滤波后图像结果如下所示：7. 车牌定位y,x,z=size(I5);% 返回I5各维的尺寸，存储在 x,y,z中myI=double(I5);% 将I5转换成双精度tic %tic表示计时的开始，toc表示计时的结束Blue_y=zeros(y,1);%产生一个 y*1 的零阵for i=1:yfor j=1:xif(myI(i,j,1)=1)Blue_y(i,1)= Blue_y(i,1)+

9、1;%蓝色像素点统计endendendtemp MaxY=max(Blue_y);%Y 方向车牌区域确定PY1=MaxY;while (Blue_y(PY1,1)>=5)&&(PY1>1)PY1=PY1-1;endPY2=MaxY;while (Blue_y(PY2,1)>=5)&&(PY2<y)PY2=PY2+1;endIY=I(PY1:PY2,:,:);Blue_x=zeros(1,x);%进一步确定 x方向的车牌区域for j=1:xfor i=PY1:PY2if(myI(i,j,1)=1)Blue_x(1,j)= Blue_x(1

10、,j)+1;endendendPX1=1;while (Blue_x(1,PX1)<3)&&(PX1<x)PX1=PX1+1;endPX2=x;while (Blue_x(1,PX2)<3)&&(PX2>PX1)PX2=PX2-1;endPX1=PX1-1;%对车牌区域的校正PX2=PX2+1;dw=I(PY1:PY2-8,PX1:PX2,:);t=toc;figure(7),subplot(1,2,1),imshow(IY),title('Linedirectionareas');figure(7),subplot(1,

11、2,2),imshow(dw),title('positioning color images');8. 字符分割与识别对分割出的彩色车牌图像进行灰度转换、二值化、均值滤波、腐蚀膨胀以及字符分割以从车牌图像中分离出组成车牌号码的单个字符图像，对分割出来的字符进行预处理(二值化、归一化)，然后分析提取,对分割出的字符图像进行识别给出文本形式的车牌号码。代码如下:imwrite(dw,'dw.jpg');%将彩色车牌写入dw文件中a=imread('dw.jpg');%读取车牌文件中的数据b=rgb2gray(a);% 将车牌图像转换为灰度图imwr

12、ite(b,'gray licence plate.jpg');%将灰度图像写入文件中figure(8);subplot(3,2,1),imshow(b),title('车牌灰度图像')g_max=double(max(max(b);g_min=double(min(min(b);T=round(g_max-(g_max-g_min)/3); % T为二值化的阈值m,n=size(b);d=(double(b)>=T); % d:.值图像imwrite(d,'binary licence plate.jpg');subplot(3,2,2)

13、,imshow(d),title('before filtering binary licence plate')h=fspecial('average',3);d=im2bw(round(filter2(h,d);%使用指定的滤波器h对h进行d即均值滤波imwrite(d,'after average licence plate.jpg');单位矩阵subplot(3,2,3),imshow(d),title('after average licence plate') se=eye(2); % eye(n) returns t

14、he n-by-n identity matrixm,n=size(d);%返回矩阵b的尺寸信息，并存储在 m,n中if bwarea(d)/m/n>=0.365 %计算二值图像中对象的总面积与整个面积的比是否大于0.365d=imerode(d,se);%如果大于0.365则图像进行腐蚀elseif bwarea(d)/m/n<=0.235 %计算二值图像中对象的总面积与整个面积的比是否小于0.235d=imdilate(d,se);%如果小于则实现膨胀操作endimwrite(d,'expansion or corrosion the licence plate.jpg

15、');subplot(3,2,4),imshow(d),title('expansion or corrosion the licence plate');运行结果如下所示：9. 字符分割在汽车牌照自动识别过程中，字符分割有承前启后的作用。它在前期牌照定位的基础 上进行字符的分割，然后再利用分割的结果进行字符识别。字符识别的算法很多，因为车 牌字符间间隔较大，不会出现字符粘连情况，所以此处采用的方法为寻找连续有文字的块,若长度大于某阈值，则认为该块有两个字符组成，需要分割。一般分割出来的字符要进行 进一步的处理，以满足下一步字符识别的需要。但是对于车牌的识别，并不需要太

16、多的处 理就已经可以达到正确识别的目的。在此只进行了归一化处理，然后进行后期处理。%寻找连续有文字的块，若长度大于某阈值，则认为该块有两个字符组成，需要分割9首先创建子函数 qiege与getword ,而后调用子程序，将车牌的字符分割开并且进行归一 化处理d=qiege(d);m,n=size(d);subplot(3,2,5),imshow(d),title(n)k1=1;k2=1;s=sum(d)；j=1;while j=nwhile s(j)=0j=j+1;endk1=j;while s(j)=0 && j<=n-1j=j+1;endk2=j-1;if k2-k1

17、>=round(n/6.5)val,num=min(sum(d(:,k1+5:k2-5);d(:,k1+num+5)=0; % 分割endendd=qiege(d);y1=10;y2=0.25;flag=0;word1=;while flag=0m,n=size(d);left=1;wide=0;while sum(d(:,wide+1)=0wide=wide+1;endif wide<y1d(:,1:wide)=0;d=qiege(d);elsetemp=qiege(imcrop(d,1 1 wide m);m,n=size(temp);all=sum(sum(temp);two

18、_thirds=sum(sum(temp(round(m/3):2*round(m/3),:);if two_thirds/all>y2flag=1;word1=temp; % WORD 1endd(:,1:wide)=0;d=qiege(d);endendword2,d=getword(d);word3,d=getword(d);word4,d=getword(d);word5,d=getword(d);word6,d=getword(d);word7,d=getword(d);figure(9);subplot(2,7,1),imshow(word1),title('1

19、9;);subplot(2,7,2),imshow(word2),title('2');subplot(2,7,3),imshow(word3),title('3');subplot(2,7,4),imshow(word4),title('4');subplot(2,7,5),imshow(word5),title('5');subplot(2,7,6),imshow(word6),title('6');subplot(2,7,7),imshow(word7),title('7');m,n=size

20、(word1);word1=imresize(word1,40 20);word2=imresize(word2,40 20);word3=imresize(word3,40 20);word4=imresize(word4,40 20);word5=imresize(word5,40 20);word6=imresize(word6,40 20);word7=imresize(word7,40 20);subplot(2,7,8),imshow(word1),title('1');subplot(2,7,9),imshow(word2),title('2');

21、subplot(2,7,10),imshow(word3),title('3');subplot(2,7,11),imshow(word4),title('4');subplot(2,7,12),imshow(word5),title('5');subplot(2,7,13),imshow(word6),title('6');subplot(2,7,14),imshow(word7),title('7');imwrite(word1,'1.jpg');imwrite(word2,'2.jpg

22、');imwrite(word3,'3.jpg');imwrite(word4,'4.jpg');imwrite(word5,'5.jpg');imwrite(word6,'6.jpg');imwrite(word7,'7.jpg');运行结果如下：s W茶®D 3I10. 车牌识别：模板匹配是图象识别方法中最具代表性的基本方法之一，它是将从待识别的图象或图象区域f(i,j)中提取的若干特征量与模板T(i,j)相应的特征量逐个进行比较，计算它们之间规格化的互相关量，其中互相关量最大的一个就表示期间

23、相似程度最高，可将图象归于相应的类。也可以计算图象与模板特征量之间的距离，用最小距离法判定所属类。此处采用相减的方法来求得字符与模板中哪一个字符最相似，然后找到相似度最大的输出。汽车牌照的字符一般有七个，大部分车牌第一位是汉字，通常代表车辆所属省份，紧接其后的为字母与数字。车牌字符识别与一般文字识别在于它的字符数有限，汉字共约50多个，大写英文字母 26个，数字10个。为了实验方便，结合本次设计所选汽车牌照的特点，只建立了7个数字26个字母与10个数字的模板。其他模板设计的方法与此相同。首先取字符模板，接着依次取待识别字符与模板进行匹配，将其与模板字符相减，得到的0越多那么就越匹配。把每一幅相

24、减后的图的0值个数保存，即为识别出来的结果源代码如下：liccode=char('0':'9' 'A':'Z''苏豫陕鲁京辽浙');% 建立自动识别字符代码表SubBw2=zeros(40,20);l=1;for I=1:7ii=int2str(I);t=imread(ii,'.jpg');SegBw2=imresize(t,40 20,'nearest');SegBw2=double(SegBw2)>20;if l=1%kmin=37;kmax=43;elseif l=2%k

25、min=11;kmax=36;else l>=3%kmin=1;第一位汉字识别第二位AZ字母识别第三位以后是字母或数字识别kmax=36;endfor k2=kmin:kmaxfname=strcat('字符模板 ',liccode(k2),'.jpg');SamBw2 = imread(fname);SamBw2=double(SamBw2)>1;for i=1:40for j=1:20SubBw2(i,j)=SegBw2(i,j)-SamBw2(i,j);endend%以上相当于两幅图相减得到第三幅图Dmax=0;for k1=1:40for l

26、1=1:20if ( SubBw2(k1,l1) > 0 | SubBw2(k1,l1) <0 )Dmax=Dmax+1;endendendError(k2)=Dmax;endError1=Error(kmin:kmax);MinError=min(Error1);findc=find(Error1=MinError);Code(l*2-1)=liccode(findc(1)+kmin-1);Code(l*2)=''l=l+1;endfigure(10),imshow(dw),title ('车牌号码:',Code,'Color',&

27、#39;b');File Edit View Insert Too Is Desktop Window Help车牌号码一湾A 0 4 S 8 90W五、设计体会以前老课程设计不仅是对前面所学知识的一种检验，而且也是对自己能力的一种提高。在以后的工作、 基础比较差，一但是在设计过程是觉得自己什么东西都不会，什么东西都不懂，而且又急于求成，结果造成什么都没学好, 还是什么都不会。通过这次课程设计，我才明白学习是一个长期积累的过程， 生活中都应该不断的学习，努力提高自己知识和综合素质，特别是对于我， 定不能太过于心急，要静下心来慢慢的研究。虽然这个设计做的也不太好， 中所学到的东西是这次课

28、程设计的最大收获和财富，使我终身受益六、参考文献：1 杨万山等，基于BP神经网络的工程图纸图形符号的识别，微型电脑应用，Vol.16 ,No.2,20002 王年、李婕、任彬、汪炳权.多层次汽车车牌照定位分割方法.安徽大学学报，1999(6)Vol.23.No.23 崔 江、王友仁.车牌自动识别方法中的关键技术研究.计算机测量与控制，2003.11(4)4 许志影、李晋平.MATLA敬其在图像处理中的应用.计算机与现代化，2004(4)源代码如下:新建文件夹0.jpg');figure(1),imshow(I);title('原始车牌');%将车牌的原图显示出来I1=r

29、gb2gray(I);%将彩图转换为灰度图原始车牌');灰度图直方图');选择阈值0.18 ,用roberts算子进行边缘检测 边缘检测图像');figure(2),subplot(1,2,1),imshow(I1);title(' figure(2),subplot(1,2,2),imhist(I1);title(' I2=edge(I1,'roberts',0.18,'both');% figure(3),imshow(I2);title(' se=1;1;1; I3=imerode(I2,se);% 对图像实

30、施腐蚀操作，即膨胀的反操作 figure(4),imshow(I3);title('腐蚀后图像');se=strel('rectangle',25,25);%构造结构元素以正方形构造一个seI4=imclose(I3,se);%图像聚类、填充图像figure(5),imshow(I4);title('闭运算后图像');I5=bwareaopen(I4,2000);% 去除聚团灰度值小于 2000的部分figure(6),imshow(I5);title('去除后图像');y,x,z=size(I5);% 返回I5各维的尺寸，存储在

31、x,y,z中myI=double(I5);% 将I5转换成双精度tic %tic表示计时的开始，toc表示计时的结束Blue_y=zeros(y,1);% 产生一个 y*1 的零阵for i=1:yfor j=1:xif(myI(i,j,1)=1)%如果myI(i,j,1) 即myI的图像中坐标为(i,j) 的点值为1,即该点为车牌背景颜色蓝色% 则 Blue_y(i,1) 的值加 1Blue_y(i,1)= Blue_y(i,1)+1;%蓝色像素点统计endendendtemp MaxY=max(Blue_y);%Y 方向车牌区域确定%temp为向量yellow_y的元素中的最大值，Max洌

32、该值的索引PY1=MaxY;while (Blue_y(PY1,1)>=5)&&(PY1>1)PY1=PY1-1;endPY2=MaxY;while (Blue_y(PY2,1)>=5)&&(PY2<y)PY2=PY2+1;endIY=I(PY1:PY2,:,:);%行方向车牌区域确定Blue_x=zeros(1,x);% 进一步确定x方向的车牌区域for j=1:xfor i=PY1:PY2if(myI(i,j,1)=1)Blue_x(1,j)= Blue_x(1,j)+1;endendendPX1=1;while (Blue_x(1,

33、PX1)<3)&&(PX1<x)PX1=PX1+1;endPX2=x;while (Blue_x(1,PX2)<3)&&(PX2>PX1)PX2=PX2-1;end PX1=PX1-1;%对车牌区域的校正 PX2=PX2+1; dw=I(PY1:PY2-8,PX1:PX2,:); t=toc; figure(7),subplot(1,2,1),imshow(IY),title('车牌行方向');% 行方向车牌区域确定figure(7),subplot(1,2,2),imshow(dw),title('车牌区域

34、9;);%定位后的车牌区域如下所示：imwrite(dw,'dw.jpg');%将彩色车牌写入dw文件中a=imread('dw.jpg');%读取车牌文件中的数据b=rgb2gray(a);% 将车牌图像转换为灰度图 figure(8);subplot(3,2,1),imshow(b),title('车牌灰度图像')g_max=double(max(max(b); g_min=double(min(min(b); T=round(g_max-(g_max-g_min)/3); % T为二值化的阈值m,n=size(b); d=(double(b

35、)>=T); % d:二值图像subplot(3,2,2),imshow(d),title('二值化图像')%匀值滤波前 %滤波 h=fspecial('average',3); 9建立预定义的滤波算子，average为均值滤波，模板的尺寸为 3*3 d=im2bw(round(filter2(h,d);%使用指定的滤波器h对h进行d即均值滤波subplot(3,2,3),imshow(d),title('滤波后图像')se=eye(2); % eye(n) returns the n-by-n identity matrix单位矩阵m,n

36、=size(d);% 返回矩阵b的尺寸信息，并存储在m,n中if bwarea(d)/m/n>=0.365 %计算二值图像中对象的总面积与整个面积的比是否大于0.365d=imerode(d,se);% 如果大于0.365则图像进行腐蚀elseif bwarea(d)/m/n<=0.235 %计算二值图像中对象的总面积与整个面积的比是否小于0.235 d=imdilate(d,se);%如果小于则实现膨胀操作end subplot(3,2,4),imshow(d),title('膨胀或腐蚀图像');d=qiege(d); m,n=size(d); subplot(3

37、,2,5),imshow(d),title(n) k1=1;k2=1;s=sum(d);j=1; while j=n while s(j)=0 j=j+1; end k1=j; while s(j)=0 && j<=n-1 j=j+1;endk2=j-1;if k2-k1>=round(n/6.5)val,num=min(sum(d(:,k1+5:k2-5);d(:,k1+num+5)=0; % 分割endend%再切割d=qiege(d);%切割出7个字符y1=10;y2=0.25;flag=0;word1=;while flag=0m,n=size(d);lef

38、t=1;wide=0;while sum(d(:,wide+1)=0wide=wide+1;endif wide<y1 %认为是左侧干扰d(:,1:wide)=0;d=qiege(d);elsetemp=qiege(imcrop(d,1 1 wide m);m,n=size(temp);all=sum(sum(temp);two_thirds=sum(sum(temp(round(m/3):2*round(m/3),:);if two_thirds/all>y2flag=1;word1=temp; % WORD 1endd(:,1:wide)=0;d=qiege(d);endend

39、%分割出第二个字符word2,d=getword(d);%分割出第三个字符word3,d=getword(d);%分割出第四个字符word4,d=getword(d);%分割出第五个字符word5,d=getword(d);%分割出第六个字符word6,d=getword(d);%分割出第七个字符word7,d=getword(d);figure(9);subplot(2,7,1),imshow(word1),title('1');subplot(2,7,2),imshow(word2),title('2');subplot(2,7,3),imshow(word

40、3),title('3');subplot(2,7,4),imshow(word4),title('4');subplot(2,7,5),imshow(word5),title('5');subplot(2,7,6),imshow(word6),title('6');subplot(2,7,7),imshow(word7),title('7');m,n=size(word1);word1=imresize(word1,40 20);word2=imresize(word2,40 20);word3=imresize(word3,40 20);word4=imresize(word4,40 20);word5=imresize(word5,40 20);word6=imresize(word6,40 20);word7=imresize(word7,40 20);subplot(2,7,8),imshow(word1),title('1');subplot(2,7,9),imshow(word2),title('2');subp