图像处理车牌识别系统设计实验报告书.doc

图像处理

班别：11医学应用（2）班
组长：
组员：
指导老师：

目录
一、摘要 2
二、 设计原理 3
1、车牌的定位研究 3
2、字符分割的研究 3
3、字符识别的研究 3
三、 详细设计步骤 3
1、车牌定位 4
1.1 图像的预处理 4
1.2车牌定位 7
2、字符分割 10
2.1对读入图像进行预处理操作 11
2.2图像校正 12
2.3去除水平方向上的边框 13
2.4去除垂直方向上的边框 15
2.5去除车牌上的圆点 17
3、字符识别 20
3.1建立字符模板数据库 20
3.2对分割字符进行匹配 22
4、系统界面的实现 25
四、 设计结果分析 29
五、 设计体会 29

车牌识别系统的设计

一、摘要
车牌是一辆汽车独一无二的信息，因此，对车辆牌照的识别技术可以作为辨识一辆车最为有效的方法。随着ITS(智能交通系统)的高速发展，对车牌识别技术的研究也随之发展。从根本上讲，牌照识别应用了先进的图像处理，模式识别，人工智能技术来获取，处理，解释，记录拍照的图像。目前， 汽车牌照的自动识别技术已经得到了广泛应用。汽车牌照自动识别整个处理过程分为预处理、边缘提取、车牌定位、字符分割、字符识别五大模块，其中字符识别过程主要由以下3个部分组成：①正确地分割文字图像区域；②正确的分离单个文字；③正确识别单个字符。用MATLAB软件编程来实现每一个部分,最后识别出汽车牌照。在研究的同时对其中出现的问题进行了具体分析、处理。
二、 设计原理
车牌自动识别是一项利用车辆的动态视频或静态图像进行车牌号码、车牌颜色自动识别的模式识别技术。其核心包括车牌定位算法、车牌字符分割算法和字符识别算法等。
输入要处理图像
车牌定位
预处理

单字符识别
字符分割
字符特征提取

图1 牌照识别系统原理图
主要研究内容如下：
1、车牌的定位研究。
先进行图像的预处理，包括RGB彩色图像的灰度化、图像灰度拉伸、图像边缘检测、灰度图的二值化等;车牌定位采用基于水平和垂直投影分布特征的方法。
2、字符分割的研究。
先对定位后的车牌图像进行预处理，然后按照车牌的先验信息, 用区域增长算法来确定候选车牌的字符区域。
3、字符识别的研究。
对于提取出的单个字符，先进行归一化操作，再与给定的模板做对比，识别出字符。
三、 详细设计步骤
1. 车牌定位：车辆牌照识别整个系统主要是由车牌定位和字符识别两部分组成，其中车牌定位又可以分为图像预处理及边缘提取模块和牌照的定位；
导入图像
图像灰度校正
基于颜色的基础移除大对象
对图像进行压缩处理
流程图：

对灰度校正后图像进行形态学处理提取边缘

计算X,Y方向车牌区域寻找到车牌，完成车牌定位
移除图像中小对象
平滑处理

1.1 图像的预处理
为了用于牌照的分割和牌照字符的识别，原始图象应具有适当的亮度，较大的对比度和清晰可辩的牌照图象。但由于该系统的摄像部分工作于开放的户外环境，加之车辆牌照的整洁度、自然光照条件、拍摄时摄像机与牌照的矩离和角度以及车辆行驶速度等因素的影响，牌照图象可能出现模糊、歪斜和缺损等严重缺陷，因此需要对原始图象进行识别前的预处理。

（1） 一般的车牌识别只对小对象进行移除，但是有时候因为拍照原因，即使对图像进行了很好的预处理，还是不能排除一些比较大的又和车牌比较相像的地方，例如图（1）中，车后面的玻璃窗仅仅进行灰度处理和形态学处理，平滑处理是无法排除其对车牌定位的影响的，因此需要对大对象移除（移除大对象前需对图像进行统一化--压缩）。
图（1）
实现的代码：
img=image;
I=img;
[x,y,z]=size(img);
if x>2000||y>2000
img=img(1:2:end,1:2:end,:);
[x,y,z]=size(img);
end
hsi=rgb2hsi(img);
for i=1:x
for j=1:y
if (hsi(i, j ,1)<=0.65&&hsi(i, j,1)>=0.55&&hsi(i, j, 2)>0.4);%0.4 0.
hsi(i, j,2)=0;
hsi(i, j,3)=1;
;
else
hsi(i, j,2)=0;
hsi(i, j,3)=0;
end
end
end
origonImg=hsi2rgb(hsi);
rgbnew=origonImg;
（2） 灰度校正：
由于牌照图象在拍摄时受到种种条件的限制和干扰，图象的灰度值往往与实际景物不完全匹配，这将直接影响到图象的后续处理。如果造成这种影响的原因主要是由于被摄物体的远近不同，使得图象中央区域和边缘区域的灰度失衡，或是由于摄像头在扫描时各点的灵敏度有较大的差异而产生图象灰度失真，或是由于曝光不足而使得图像的灰度变化范围很窄。这时就可以采用灰度校正的方法来处理，增强灰度的变化范围、丰富灰度层次，以达到增强图象的对比度和分辨率。

灰度校正图 graynew=rgb2gray(rgbnew);

（3） 边缘检测：
边缘是指图像局部亮度变化显著的部分，是图像风、纹理特征提取和形状特征提取等图像分析的重要基础。所以在此我们要对图像进行边缘检测。图象增强处理对图象牌照的可辩认度的改善和简化后续的牌照字符定位和分割的难度都是很有必要

实现代码：origonImg=hsi2rgb(hsi);
rgbnew=origonImg;
graynew=rgb2gray(rgbnew);
graynew=imfill(graynew,'holes');
for i=1:3
se1=[1 1 1;1 1 1;1 1 1];
graynew=imdilate( graynew,se1);
se=[1;1;1];
graynew=imerode(graynew,se);% 腐蚀Imerode(X,SE).其中X是待处理的像，%SE是结构元素对象
graynew=getcenter(graynew);
graynew=double(graynew);
[row col]=size(graynew);
1.2车牌定位
牌照的定位是在经图象预处理后的原始灰度图象中确定牌照的具体位置，并将包含牌照字符的一块子图象从整个图象中分割出来，供字符识别子系统识别之用。

图为定位出的车牌
实现代码：
车牌裁剪：
[row col]=size(graynew);
colArray=zeros(row,1);
for i=1:row
for j=1:col
if (graynew(i,j,1)==1)
colArray(i,1)=colArray(i,1)+1;
end
end
end
[maxcount maxIndex]=max(colArray);
roughY=maxIndex;
while ((colArray(roughY,1)>=5)&&(roughY>1))
roughY=roughY-1;
end
y1=roughY;
roughY=maxIndex;
while ((colArray(roughY,1)>=5)&&(roughY<row))
roughY=roughY+1;
end
y2=roughY;
rowArray=zeros(1,col);
for i=1:col
for j=y1:y2
if graynew(j,i)==1
rowArray(1,i)=rowArray(1,i)+1;
end
end
end
[maxcount maxIndex]=max(rowArray);
roughX=maxIndex;
while ((rowArray(1,roughX)>=5)&&(roughX>1))
roughX=roughX-1;
end
x1=roughX;
roughX=maxIndex;
while ((rowArray(1,roughX)>=5)&&(roughX<col))
roughX=roughX+1;
end
x2=roughX;
dw=img(y1:y2,x1:x2,:);
greenframe=drawframe(img,y1,y2,x1,x2);
imwrite(greenframe,'imgAfterLocation/greenframe.jpg');
imwrite(dw,'imgAfterLocation/dw.jpg');

2、字符分割：
传统的字符分割方法有投影法、模板匹配法、聚类分析法。由于投影法比较准确、编程较简单，且易于实际操作，能满足在复杂环境下，所以我们采用的是投影法分割车牌字符。

车牌分割总流程图：

各部分代码实现和运行效果：

2.1、对读入图像进行预处理操作
I=imread('imgAfterLocation/dw.jpg');%读取图像
I1 = rgb2gray(I);%彩色图像转化为灰度图像
T=graythresh(I1);%找到灰度图像的阈值
[y,x,z]=size(I1);%计算I1各维的大小
Se=strel('disk',fix(y/45));%创建一个平坦的圆盘形结构元素
I2=imopen(I1,Se);%取圆盘形的开运算
subplot(2,2,3);%图一的第三幅图
imshow(I2),title('开运算后图像');
Se=strel('diamond',fix(y/140));%创建一个平坦的菱形结构元素
I3=imclose(I2,Se);%取菱形结构的闭运算
I4=double(I3);%变为双精度

2.2图像校正
代码：
bw1=edge(I1,'sobel', 'horizontal');%用Sobel水平算子对图像边化
figure,imshow(bw1)
bw1=imcrop(bw1,[0 0 500 100]);%对图像进行剪切，保留图像的 %一条直线，减小运算量
figure,imshow(bw1)
theta=0:179;%定义theta角度范围
r=radon(bw1,theta);%对图像进行Radon变换
figure,imshow(r)
[m,n]=size(r);
c=90;
for i=1:m
for j=1:n
if r(1,1)<r(i,j)
r(1,1)=r(i,j);
c=j;
end
end
end %检测Radon变换矩阵中的峰值所对应的列坐标

rot=90-c;%确定旋转角度
A=imrotate(I1,rot,'crop');%对图像进行旋转矫正

2.3计算行方向的像素，用来去除行方向上的边框
代码：
计算行像素：
Y1=zeros(y,1);
for i=1:y
for j=1:x
if(I4(i,j,1)==1)
Y1(i,1)= Y1(i,1)+1;%计算I3第j列有几个一
end
end
end
figure(2);
plot(Y1,0:y-1),title('行方向像素点灰度值累计和'),xlabel('累计像素量'),ylabel('行');

去除行方向边框：
Py0=fix(y/2);
Py1=fix(y/2)+1;
while ((Y1(Py0,1)>=30)&&(Py0>2))
Py0=Py0-1;%找到去除边框后上边的位置
end
while ((Y1(Py1,1)>=30)&&(Py1<y))
Py1=Py1+1;%找到去除边框后下边的位置
end
Z1=I4(Py0:Py1,:,:);%将二值图像上下边框去除
figure(3);
imshow(Z1),title('将二值图像上下边框去除后图像');

2.4计算列方向的像素，用来去除垂直方向上的边框
代码：
[y,x,z]=size(Z1);%计算此时图像的大小
计算列像素
X1=zeros(1,x);
for j=1:x
for i=1:y
if(Z1(i,j,1)==1)
X1(1,j)= X1(1,j)+1;%计算I3第j列有几个一
end
end
end
figure(4);
plot(0:x-1,X1),title('列方向像素点灰度值累计和'),xlabel('列值'),ylabel('累计像素');

去除垂直边框
x1=fix(x/2)+1;
for i=1:5
while (i~=5)
while ((X1(1,x1)>=15)&&(x1<x))
x1=x1+1;%找到去除边框后右边的位置
end
i=i+1;
while ((X1(1,x1)<15)&&(x1<x)&&(i<4))
x1=x1+1;%从车牌中间开始寻找字符间隙，直到找到第四个间隙%为止，即找到去除边框后车牌右边的位置
end
end
end
x0=fix(x*45/440);%找到第一个字符的位置
for i=1:2
while (i~=2)
while ((X1(1,x0)>=5)&&(x0>2))
x0=x0-1;%从车牌的第一个字符开始寻找第一个字符间隙，找到%去除边框左边的位置
end
i=i+1;
while ((X1(1,x0)<5)&&(x0>1)&&i~=2)
x0=x0-1;%找到去除边框后左边的位置
end
end
end
Z2=Z1(:,x0:x1,:);%将二值图像左右边框去除
figure(5);
imshow(Z2),title('将二值图像垂直边框去除后图像');

2.5再次计算列方向的像素，用来去除车牌上的圆点
代码：
[y,x,z]=size(Z2);%计算此时图像的大小
X1=zeros(1,x);
for j=1:x
for i=1:y
if(Z2(i,j,1)==1)
X1(1,j)= X1(1,j)+1;%计算I3第j列有几个一
end
end
End
去除车牌中的圆点：
[y,x,z]=size(Z2);
Z2=double(Z2);
x1=fix(x*105/409);
x2=fix(x*128/409);
for i=1:y
for j=x1:x2
Z2(i,j)=0;
end
end
figure(6);
imshow(Z2),title('去除圆点后图像');

2.6用投影法分割车牌，把分割出来的子图转成标准子图
代码：
用投影法分割
Px0=1;
Px1=1;
y=fix(x*90/409);
for i=1:7
while ((X1(1,Px0)<10)&&(Px0<x))
Px0=Px0+1;%找到待分割字符左边的位置
end
Px1=Px0;
a=1;
while (((X1(1,Px1)>=10)&&(Px1<x))||((Px1<x)&&((a/y)<=0.5)))
Px1=Px1+1;%找到待分割字符右边的位置
a=Px1-Px0;
end
Z3=Z2(:,Px0:Px1,:);%二值化图像分割后
figure(7);
subplot(1,7,i);
imshow(Z3);%将二值化图像分割后显示出来

2.7变换为标准子图
Z4=imresize(Z3,[40 20]);%将分割后的二值图变换为标准子图
figure(1);
subplot(1,7,i);
imshow(Z4),title('标准子图');%将标准子图显示出来
Px0=Px1;
imwrite(Z4,strcat('imgAfterSplit/',num2str(i),'.jpg'));%把标准子图存到imagAfterSplit中
End

3、字符识别：
牌照字符识别 :字符识别方法目前主要有基于模板匹配算法和基于人工神经网络算法。基于模板匹配算法首先将分割后的字符二值化,并将其尺寸大小缩放为字符数据库中模板的大小，然后与所有的模板进行匹配，最后选最佳匹配作为结果。基于人工神经元网络的算法有两种：一种是先对待识别字符进行特征提取，然后用所获得特征来训练神经网络分配器；另一种方法是直接把待处理图像输入网络，由网络自动实现特征提取直至识别出结果。实际应用中，牌照识别系统的识别率与牌照质量和拍摄质量密切相关。牌照质量会受到各种因素的影响，如生锈、污损、油漆剥落、字体褪色、牌照被遮挡、牌照倾斜、高亮反光、多牌照、假牌照等等；实际拍摄过程也会受到环境亮度、拍摄亮度、车辆速度等等因素的影响。这些影响因素不同程度上降低了牌照识别的识别率，也正是牌照识别系统的困难和挑战所在。为了提高识别率，除了不断的完善识别算法，还应该想办法克服各种光照条件，使采集到的图像最利于识别。
本系统字符识别方法采用基于模板匹配算法。
3.1建立字符模板数据库
模板库的合理建造是字符识别准确的关键之一，所以在字符识别之前必须把模板库设置好。汽车牌照的字符一般有 7 个，大部分车牌第一位是汉字，通常代表车辆所属省份，或是军种 警别等有特定含义的字符简称；紧接其后的为字母与数字。车牌字符识别与一般文字识别在于它的字符数有限， 十个阿拉伯数字0～9， 26 个大写英文字母 A～Z 以及相关的车牌用汉字：京、沪、苏、台、港、澳、甲、乙、丙、使、领、学、试、境、消、边、警等，以及新式军牌中的汉字南、兰、广、北、沈、济、空、海等；车牌颜色：蓝、白、黑、黄等。所以建立字符模板库也极为方便。
为了提高识别的准确率，本系统的为每个字符都建立一个模板库，并且可以动态添加字符到模板库，随着字符模板的增加，识别率也随之提高。
建立模板数据库时必须对这些图片进行统一处理，因为对前面处理分割后的车牌图像的测量得知单个字符的最佳宽高比是1：2，所以将这些图片归一化为40×20大小；因为之后的字符识别考虑使用模板匹配算法进行字符识别，所以再将上面归一化后的模板图像的样本排列在一起构成40×20的矩阵样本。以下分别是粤和A字符模板：

3.2对分割字符进行匹配
Recognition.m程序代码（最终返回识别的车牌号）：
function [CarNum ] = recognition ()
codes=char(['0':'9' 'A':'Z' '京津沪渝吉辽鲁豫冀鄂湘晋青皖苏赣浙闽粤琼陕甘云川贵黑藏蒙桂新宁']);
%建立自动识别字符代码表；'京津沪渝吉辽鲁豫冀鄂湘晋青皖苏赣浙闽粤琼陕甘云川贵黑藏蒙桂新宁'
% 编号：0-9分别为 1-10;A-Z分别为 11-36;
% 京 津 沪 渝 吉 辽 鲁 豫 冀 鄂 湘 晋 青 皖 苏
% 赣 浙 闽 粤 琼 陕 甘 云 川 贵 黑 藏 蒙 桂 新 宁
% 37 --67
SubBw2=zeros(40,20);
l=1;
for I=1:7
ii=int2str(I); % 将整型数据转换为字符串型数据
t=imread(['imgAfterSplit/',ii,'.jpg']);% 依次读入七位车牌字符
SegBw2=imresize(t,[40 20],'nearest'); % 对读入的字符进行缩放
SegBw2=im2bw(SegBw2);
if I==1 % 第一位汉字识别
kmin=37; kmax=67;
elseif I==2 % 第二位 A~Z 字母识别
kmin=11; kmax=36;
else % 第三位以后是字母或数字识别 ; 即I>=3
kmin=1; kmax=36;
end
index=getFinalTemp(kmin,kmax,codes,SegBw2,I);
CarNum(l)=codes(index);
l=l+1;
end
% xlabel(['第三步：识别结果为: ', CarNum],'Color','b');

getFinalTemp.m 程序代码（返回最终识别的单个字符index 下标）：

function [index]=getFinalTemp(minIndex,maxIndex,codes,SegBw2,numIndex)
if numIndex<3 %后五位不能出现 i 、o 字母;
for i=minIndex:maxIndex
path=strcat('template/', codes(i));
%调用 getFromLib 遍历文件夹找出同一字符最佳的模板
maxSame=getFromLib(path,SegBw2);
near(i)= maxSame; % 记录下字符与模板i同的点个数
end
else
for i=minIndex:maxIndex
% if i~=find(codes='l')&&i~=find(codes='O
if codes(i)~='l'&&codes(i)~='O'
path=strcat('template/', codes(i));
%调用 getFromLib 遍历文件夹找出同一字符最佳的模板
maxSame=getFromLib(path,SegBw2);
near(i)= maxSame; % 记录下字符与模板i同的点个数
end
end
end
nearest=max(near); % 差别最小的，不同字符模板
index=find(near==nearest); % 找出差别最小的模板的下标

getFromLib.m 程序代码（遍历单个字符模板库，返回最佳的模板）：

function [maxSame]=getFromLib(path,carNum)
files=dir([path '/*.jpg']);
for i=1:length(files)
file=imread([path '/' num2str(i) '.jpg' ]);
if islogical(file)==0
filebw=im2bw(file);
else
filebw=file;
end
subImg=carNum-filebw;
% 以上相当于两幅图相减得到第三幅图 进行匹配
sameDots=subImg==0;
% 记录下字符与模板相同的点个数
sumDots=sum(sameDots(:));
arr(i)=sumDots;
end
maxSame=max(arr);

4、 系统界面的实现：
界面我们小组一共分有了两个模块：
一个是父模块车牌识别系统，进行显图，定位，分割和识别的操作；
另一个是子模块是字符模块，是为了增加字符的模块进库，提高车牌的识别效率。

父模块——车牌识别系统

图一

界面分别有载入图片，车牌定位，车牌分割，车牌识别和增加字符模块五个按钮，每个按钮对应不同的功能。

载入车牌：选择你需要识别的带有车牌的图片并显示出来在原始图片的位置。
车牌定位：内部进行定位函数对车牌进行定位，并在原图上用黄线圈出车牌相应的位置。
车牌分割：将定位好的车牌进行分割，分割出7个独立的字符模块。
车牌识别：对每个分割出来的字符进行识别，并文本输出来。
增加字符模块：点击进入子界面，进行字符模块的增加，提高车牌识别的效率。

效果如下图：

图二

子模块——字符模块

图三

字符模块分为5大部分，分别是车牌显示，选择添加模板的位置，选择要添加的模板，添加字符模块和显示字符模块。

车牌显示：将刚在父模块里车牌分割出来的图片显示在这个对应的位置上。
选择添加模板的位置：选择你需要添加的字符在图片对应的位置，如“粤”在第一位就选择“1”。
选择要添加的模板：选择你需要添加进库的字符对应的数字，字母或者文字，如要添加第一位“粤”字就选择粤“粤”字。
添加字符模板：当你选择好“选择添加模板的位置”和“选择要添加的模板”后，就可以按这个按钮确认了。如图四
显示字符模块：添加成功后你就可以点击这个按钮查看对应的某字符的字符模板库里的所有“粤”字模板了。如图五

图四

图五

四、 设计结果分析

从上面结果看车牌识别成功，但有一些车牌会识别出错，所以在车牌识别过程中，数据库的建立很重要。只有数据库的准确才能检测出来的数据正确，因此，我们的系统能够在识别过程中不断丰富数据库，从而不断提高正确率。

五、 设计体会
经过几周的奋战，我们的课程设计终于完成了。课程设计不仅是对所学知识的应用和体会，也是对自己能力的提高。做课程设计之前，老是觉得自己什么都没学会，什么东西都不懂，也不知道学这门课的实际用途是什么，觉得很迷茫。通过这次课程设计，才明白学习是一个长期积累的过程，在以后的工作，生活中都应该不断学习，努力提高自己的知识和综合素质。通过这次课程设计，感悟最深的是学习不能只是埋头苦读，而是要合理利用资源，老师的建议，网上的资料等等。特别是在设计的过程中，组员之间不懂的大家一起讨论，听听不同的看法，学习更多的知识。还有就是知识必须通过应用才能实现其价值！有些东西以为学会了，但和真正应用起来是完全不同的，就如这次的车牌识别，想着利用常规车牌的大小定位车牌，但是真正实现起来是很难的，所以知识只有真正会用的时候才是真的学会了。这个课程设计虽然还有很多不完善的地方，但是是通过我们组一起努力奋斗，尽自己最大努力做出来的，整个设计过程中收获很多，受益匪浅。也更加明白，在今后的学习中，要多动手，多交流，特别是多和老师交流有关这方面知识的应用等等。

目 录
第一章 项目摘要 3
1.1项目基本情况 3
1.2建设目标 3
1.3建设内容及规模 4
1.4产品及去向 4
1.5效益分析 4
第二章 项目建设的可行性和必要性 5
2.1建设的必要性 5
2.2建设的可行性 5
2.3编制依据 6
2.4编制原则 9
第三章 项目建设的基础条件 9
3.1建设单位的基本情况 9
3.2项目的原料供应情况 10
3.3地址选择分析 10
第四章 产品 11
4.1沼气 11
4.2 沼气产量确定 12
4.3有机肥 13
4.4产品去向 13
第五章 沼气工程工艺设计 14
5.1工艺参数 14
5.2处理工艺选择 14
5.3工艺流程的组成 15
5.4厌氧处理工艺选择与比较 15
5.5沼气存储和净化工艺 16
5.6工艺流程 18
5.7沼气输配设施 19
5.8沼气计量设施 19
第六章 总体设计 19
6.1站内总体设计 19
6.2站外配套设计 19
第七章 土建设计 20
7.1建筑设计 20
7.2结构设计 20
第八章 电气设计 21
8.1设计依据 21
8.2设计规范 22
8.3 设计说明 22
8.4控制与保护 22
8.5防雷与接地 22
8.6配电系统 23
8.7防雷与接地 23
8.8 防爆设计 23
8.9供电负荷 23
第九章 安全、节能及消防 24
9.1安全生产 24
9.2防火消防 24
9.3节能 25
第十章 主要构（建）筑物、设备的设计参数 25
10.2 厌氧消化系统工艺参数设计 27
10.3 沼气净化系统工艺参数设计 28
10.4 沼气储存系统 28
10.5 沼肥储存系统 29
10.6配套设施区 29
第十一章 投资概算和资金筹措 30
11.1编制说明 30
11.2总投资估算表 31
11.3投资概算 33
11.4资金筹措 33
第十二章 项目实施进度和投招标 34
12.1进度安排 34
12.2招（投）标依据 34
12.3招（投）标范围 34
12.4招（投）方式 35
第十三章 项目组织与管理 35
13.1管理 35
13.2劳动定员和组织培训 37
第十四章 环境保护和安全生产 37
14.1污染源和污染物 37
14.2污染治理方案 38
14.3安全生产 39
第十五章 产品市场分析与预测 41
15.1沼气 41
15.2沼气发电 41
15.3沼液和沼渣 43
15.4（生态）农产品。 43
第十六章 社会、生态及经济效益分析 43
16.1社会效益 43
16.2生态效益 44
16.3经济效益 44
第十七章 结论 46
第十八章 附件 47

32
32