车牌识别---数字图像课设(共19页)

1、数字图像处理课程设计报告(bogo)课设题目：车牌识别学 院：信息科学与工程学院专 业：电子与信息工程班 级：0802503姓 名：周亚楠学 号：080250303指导教师：于海雁哈尔滨工业大学（威海）2011年11月23日 不要删除行尾的分节(fn ji)符，此行不会被打印哈尔滨工业大学（威海）课程设计报告PAGE II- - PAGE 22 -TOC o 1-3 h u HYPERLINK l _Toc27700 一. 课程设计任务(rn wu) PAGEREF _Toc27700 3 HYPERLINK l _Toc17437 二. 课程设计原理(yunl)及设计方案 PAGEREF _

2、Toc17437 4 HYPERLINK l _Toc2855 三课程设计的步骤(bzhu)和结果 PAGEREF _Toc2855 5 HYPERLINK l _Toc26988 三. 设计体会 PAGEREF _Toc26988 19 HYPERLINK l _Toc5474 四. 参考文献 PAGEREF _Toc5474 20 课程设计任务(rn wu)在交通管理过程中，通常采用视频监控方式对闯红灯和超速等违章车辆进行监督。对违章车辆，需要自动检测车牌信息，提取(tq)车牌号码，以便查找车主信息和监督管理。国内常用的一般车牌通常是是蓝底白字，长宽比3:1。车牌定位系统的目的在于正确获取

3、整个图像中车牌的区域，并识别出车牌号。通过设计实现车牌识别系统，能够提高学生分析问题和解决问题的能力，还能培养(piyng)一定的科研能力。要求：1、对车牌图像进行预处理，然后进行车牌定位；2、进行字符分割；2、对车牌中的数字和字母进行提取和识别（对汉字不作要求）；（提高部分）3、要求自行设计方案、编写代码实现上述功能。课程设计原理(yunl)及设计方案 1. 基本原理 车牌识别系统作为数字摄像、计算机信息管理、图像分割(fng)和图形识别技术在智能交通领域的应用，是智能交通管理系统中重要的组成部分。车牌识别(shbi)技术可应用于道路交通监控、交通事故现场勘察、交通违章自动记录、高速公路超速

4、管理系统、小区智能化管理等方面，为智能交通管理提供了高效、实用的手段。 目前世界各国都在进行适用于本国汽车牌照的自动识别研究，美、日、韩等国已有相关系统（基于传感器）问世。引进这些系统费用比较高、而且由于各国车牌和实际的交通环境不同，引进的系统往往无法满足我国城市的需求，而国内市场上虽然已有产品投入使用，但是在后续处理时很大程度上仍然需要人工识别，所以对车牌识别技术的研究依然是目前高科技领域的热门课题之一。 车牌识别系统的成功设计、开发和应用具有相当大的社会效益、经济效益和学术意义。 图像拍摄预处理字符分割字符识别2.系统框图三课程设计的步骤(bzhu)和结果1.总体设计方案(fng n):

5、为了进行牌照识别，需要以下(yxi)几个基本的步骤：a.牌照定位，定位图片中的牌照位置；b.牌照字符分割，把牌照中的字符分割出来；c.牌照字符识别，把分割好的字符进行识别，最终组成牌照号码。牌照识别过程中，牌照颜色的识别依据算法不同，可能在上述不同步骤实现，通常与牌照识别互相配合、互相验证。(1）牌照定位:自然环境下，汽车图像背景复杂、光照不均匀，如何在自然背景中准确地确定牌照区域是整个识别过程的关键。首先对采集到的视频图像进行大范围相关搜索，找到符合汽车牌照特征的若干区域作为候选区，然后对这些侯选区域做进一步分析、评判，最后选定一个最佳的区域作为牌照区域，并将其从图象中分割出来。流程图:导入

6、原始图像图像预处理增强效果图像边缘提取对图像开闭运算车牌定位 （2）牌照字符分割 :计算水平投影进行车牌水平校正去掉车牌的框架分析垂直投影找到每个字符中心位置按左右宽度切割出字符完成牌照区域的定位后，再将牌照区域分割成单个字符，然后进行识别。字符分割一般采用垂直投影法。由于字符在垂直方向上的投影必然在字符间或字符内的间隙处取得局部最小值的附近，并且这个位置应满足牌照的字符书写格式、字符、尺寸限制和一些其他条件。利用垂直投影法对复杂环境下的汽车图像中的字符分割有较好的效果。 （3）牌照(pizho)字符识别 :字符识别方法目前主要有基于模板匹配算法和基于人工神经网络算法。基于模板匹配算法首先将分

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

8、牌照识别系统的困难和挑战所在。为了提高识别率，除了不断的完善识别算法，还应该想办法克服各种光照条件，使采集到的图像最利于识别。2.各模块(m kui)的实现:（1）输入待处理的原始图像:clear ;close all;%Step1 获取图像 装入待处理彩色图像并显示原始图像Scolor = imread(3.jpg);%imread函数读取图像文件图2.1原始图像（2）图像的灰度化:彩色图像包含着大量的颜色信息，不但在存储上开销很大，而且在处理上也会降低系统的执行速度，因此在对图像进行识别等处理中经常将彩色图像转变为灰度图像，以加快处理速度。由彩色转换为灰度的过程叫做灰度化处理。选择的标准是

9、经过灰度变换后，像素的动态范围增加，图像的对比度扩展，使图像变得更加清晰、细腻、容易识别。%将彩色图像转换为黑白并显示Sgray = rgb2gray(Scolor);%rgb2gray转换成灰度图figure,imshow(Sgray),title(原始黑白图像);图2.2原始黑白(hibi)图像（3）对原始图像进行开操作(cozu)得到图像背景图像:s=strel(disk,13);%strei函数(hnsh)Bgray=imopen(Sgray,s);%打开sgray s图像figure,imshow(Bgray);title(背景图像);%输出背景图像图2.3背景图像（4）灰度图像与背

10、景图像作减法，对图像进行增强处理:Egray=imsubtract(Sgray,Bgray);%两幅图相减figure,imshow(Egray);title(增强黑白图像);%输出黑白图像图2.4黑白图像（5）取得最佳(zu ji)阈值，将图像二值化:二值图像是指整幅图像画面内仅黑、白二值的图像。在实际的车牌处理系统中，进行图像二值变换的关键是要确定合适的阀值，使得字符与背景能够分割开来，二值变换的结果图像必须要具备良好的保形性，不丢掉有用的形状信息，不会产生额外的空缺等等。车牌识别系统要求处理的速度(sd)高、成本低、信息量大，采用二值图像进行处理，能大大地提高处理效率。阈值处理(chl)

11、的操作过程是先由用户指定或通过算法生成一个阈值，如果图像中某中像素的灰度值小于该阈值，则将该像素的灰度值设置为0或255，否则灰度值设置为255或0。fmax1=double(max(max(Egray);%egray的最大值并输出双精度型fmin1=double(min(min(Egray);%egray的最小值并输出双精度型level=(fmax1-(fmax1-fmin1)/3)/255;%获得最佳阈值bw22=im2bw(Egray,level);%转换图像为二进制图像bw2=double(bw22);figure,imshow(bw2);title(图像二值化);%得到二值图像图2.

12、5二值图像（6）边缘检测:两个具有不同灰度值的相邻区域之间总存在边缘，边缘就是灰度值不连续的结果，是图像分割、纹理特征提取和形状特征提取等图像分析的基础。为了对有意义的边缘点进行分类，与这个点相联系的灰度级必须比在这一点的背景上变换更有效，我们通过门限方法来决定一个值是否有效。所以，如果一个点的二维一阶导数比指定的门限大，我们就定义图像中的次点是一个边缘点，一组这样的依据事先定好的连接准则相连的边缘点就定义为一条边缘。经过一阶的导数的边缘检测，所求的一阶导数高于某个阈值，则确定该点为边缘点，这样会导致检测的边缘点太多。可以通过求梯度局部最大值对应的点，并认定为边缘点，去除非局部最大值，可以检测

13、出精确的边缘。一阶导数的局部最大值对应二阶导数的零交叉点，这样通过找图像强度的二阶导数的零交叉点就能找到精确边缘点。grd=edge(bw2,canny)%用canny算子识别强度图像中的边界figure,imshow(grd);title(图像边缘提取);%输出图像边缘图2.6像边缘(binyun)提取（7）对得到图像(t xin)作开操作进行滤波:数学形态非线性滤波，可以用于抑制噪声，进行特征提取、边缘检测、图像分割等图像处理问题。腐蚀是一种消除边界点的过程，结果是使目标缩小，孔洞增大，因而可有效的消除孤立噪声点；膨胀是将与目标物体接触的所有背景点合并到物体中的过程，结果是使目标增大，孔洞

14、缩小，可填补目标物体中的空洞，形成连通域。先腐蚀后膨胀的过程称为(chn wi)开运算，它具有消除细小物体，并在纤细处分离物体和平滑较大物体边界的作用；先膨胀后腐蚀的过程称为闭运算，具有填充物体内细小空洞，连接邻近物体和平滑边界的作用。对图像做了开运算和闭运算，闭运算可以使图像的轮廓线更为光滑，它通常用来消掉狭窄的间断和长细的鸿沟，消除小的孔洞，并弥补轮廓线中的断裂。bg1=imclose(grd,strel(rectangle,5,19);%取矩形框的闭运算figure,imshow(bg1);title(图像闭运算5,19);%输出闭运算的图像bg3=imopen(bg1,strel(re

15、ctangle,5,19);%取矩形框的开运算figure,imshow(bg3);title(图像开运算5,19);%输出开运算的图像bg2=imopen(bg3,strel(rectangle,19,1);%取矩形框的开运算figure,imshow(bg2);title(图像开运算19,1);%输出开运算的图像图2.7.1闭运算的图像 图2.7.2开运算的图像图2.7.3开运算(yn sun)的图像（8）对二值图像(t xin)进行区域提取，并计算区域特征参数。进行区域(qy)特征参数比较，提取车牌区域:a.对图像每个区域进行标记，然后计算每个区域的图像特征参数：区域中心位置、最小包含矩

16、形、面积。L,num = bwlabel(bg2,8);%标注二进制图像中已连接的部分Feastats = imfeature(L,basic);%计算图像区域的特征尺寸Area=Feastats.Area;%区域面积BoundingBox=Feastats.BoundingBox;%x y width height车牌的框架大小RGB = label2rgb(L, spring, k, shuffle); %标志图像向RGB图像转换figure,imshow(RGB);title(图像彩色标记);%输出框架的彩色图像图2.8.1彩色图像b. 计算出包含所标记的区域的最小宽和高，并根据先验知识

17、，比较谁的宽高比更接近实际车牌宽高比，将更接近的提取并显示出来。计算矩形的高度框架的宽度和高度的范围车牌的开始列车牌的开始行计算车牌长宽比获取车牌二值子图计算矩形的宽度 2.8.2灰度子图和二值子图（9）对水平投影(tuyng)进行峰谷分析:对水平投影进行峰谷分析(fnx),计算出车牌上边框、车牌字符投影、车牌下边框的波形峰上升点、峰下降点、峰宽、谷宽、峰间距离、峰中心位置参数。histcol1=sum(sbw1); %计算垂直(chuzh)投影histrow=sum(sbw1); %计算水平投影figure,subplot(2,1,1),bar(histcol1);title(垂直投影（含边

18、框）);%输出垂直投影subplot(2,1,2),bar(histrow); title(水平投影（含边框）);%输出水平投影图2.9.1垂直投影(tuyng)和水平投影figure,subplot(2,1,1),bar(histrow); title(水平投影(tuyng)（含边框）);%输出水平投影subplot(2,1,2),imshow(sbw1);title(车牌二值子图);%输出(shch)二值图对水平投影进行峰谷分析:求水平投影的最小值取阈值计算谷宽度计算峰距离计算下降点找到峰中心位置求水平投影的平均值图2.9.2水平投影和二值图 程序流程图（10）计算车牌旋转角度:a.车牌倾

19、斜的原因导致投影效果峰股谷不明显，在这里需要做车牌矫正处理。这里采取的线性拟合的方法,计算出车牌上边或下边图像值为1的点拟合直线与水平X轴的夹角。求最大宽度为字符检测上边从顶边至第一个峰下降点扫描从底边至最后一个峰的上升点扫描找第一个为1的点标示出图像大小程序(chngx)流程图 (2)线性拟合，计算(j sun)与x夹角fresult = fit(xdata,ydata,poly1); %poly1表示(biosh)一介拟合 Y = p1*x+p2p1=fresult.p1;angle=atan(fresult.p1)*180/pi; %弧度换为度，360/2pi, pi=3.14 (3)旋

20、转车牌图象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);%显示车牌的旋转角度图2.10.1旋转后的灰度图像和旋转角度b

21、.旋转车牌后重新计算车牌水平投影，去掉车牌水平边框，获取字符高度:histcol1=sum(sbw); %计算(j sun)垂直投影histrow=sum(sbw); %计算水平(shupng)投影figure,subplot(2,1,1),bar(histcol1);title(垂直投影(tuyng)（旋转后）);subplot(2,1,2),bar(histrow); title(水平投影（旋转后）);图2.10.2垂直投影（旋转后）和水平投影（旋转后）figure,subplot(2,1,1),bar(histrow); title(水平投影（旋转后）);subplot(2,1,2),i

22、mshow(sbw);title(车牌二值子图（旋转后）);图2.10.3水平投影（旋转后）和车牌二值子图（旋转后）（11）去水平（上下）边框,获取字符高度: a.通过以上水平投影、垂直投影分析计算，获得了车牌字符高度、字符顶行与尾行、字符宽度、每个字符的中心位置，为提取分割字符具备了条件。maxhight=max(markrow2);findc=find(markrow2=maxhight);rowtop=markrow(findc);rowbot=markrow(findc+1)-markrow1(findc+1);sbw2=sbw(rowtop:rowbot,:); %子图为(rowbo

23、t-rowtop+1)行maxhight=rowbot-rowtop+1; %字符高度(rowbot-rowtop+1)b.计算车牌垂直投影，去掉车牌垂直边框，获取车牌及字符平均宽度histcol=sum(sbw2); %计算垂直投影figure,subplot(2,1,1),bar(histcol);title(垂直投影（去水平边框后）);%输出车牌的垂直投影图像subplot(2,1,2),imshow(sbw2); %输出垂直投影图像title(车牌字符高度： ,int2str(maxhight),Color,r);%输出车牌字符高度%对垂直投影(tuyng)进行峰谷分析求垂直投影的最小

24、值取阈值计算字符上升点计算谷宽度计算字符距离找到字符中心位置求垂直投影的平均值图2.11垂直投影图像和车牌字符(z f)高度 程序流程图c.计算车牌上每个字符(z f)中心位置，计算最大字符宽度maxwidthl=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(markc

25、ol5); %字符中心距离（字符中心点至下一个字符中心点）maxs=max(markcol6); %查找最大值，即为第二字符与第三字符中心距离findmax=find(markcol6=maxs);markcol6(findmax)=0;maxwidth=max(markcol6);%查找最大值，即为最大字符宽度d.提取分割字符,并变换为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;

26、if cleftn2 cright=n2; cleft=n2-maxwidth; end SegGray=sbw(rowtop:rowbot,cleft:cright); SegBw1=sbw(rowtop:rowbot,cleft:cright); SegBw2 = imresize(SegBw1,22 14); %变换(binhun)为22行*14列标准子图 subplot(2,n1,l),imshow(SegGray); if l=7 title(车牌字符(z f)宽度： ,int2str(maxwidth),Color,r); end subplot(2,n1,n1+l),imshow(SegBw2); fname=strcat(F:MATLABworksamimage,int2str(k),.jpg);%保存子图备选(bi xun)入样本库,并建立样本库 imwrite(SegBw2,fname,jpg) l=l+1;EndGUI界面设计(shj)体会在设计的过程(guchng)中，遇到一些问题，比如说在字符分隔显示时，由于最初的原代码是显示在一个figure上，而在GUI界面上要显示在7个不同的小图上，这