图像处理车牌识别系统设计实验报告书Word格式.docx
《图像处理车牌识别系统设计实验报告书Word格式.docx》由会员分享,可在线阅读,更多相关《图像处理车牌识别系统设计实验报告书Word格式.docx(31页珍藏版)》请在冰点文库上搜索。
hsi=rgb2hsi(img);
fori=1:
x
forj=1:
y
if(hsi(i,j,1)<
=0.65&
&
hsi(i,j,1)>
=0.55&
hsi(i,j,2)>
0.4);
%0.40.
hsi(i,j,2)=0;
hsi(i,j,3)=1;
;
else
hsi(i,j,3)=0;
end
end
origonImg=hsi2rgb(hsi);
rgbnew=origonImg;
(2)灰度校正:
由于牌照图象在拍摄时受到种种条件的限制和干扰,图象的灰度值往往与实际景物不完全匹配,这将直接影响到图象的后续处理。
如果造成这种影响的原因主要是由于被摄物体的远近不同,使得图象中央区域和边缘区域的灰度失衡,或是由于摄像头在扫描时各点的灵敏度有较大的差异而产生图象灰度失真,或是由于曝光不足而使得图像的灰度变化围很窄。
这时就可以采用灰度校正的方法来处理,增强灰度的变化围、丰富灰度层次,以达到增强图象的对比度和分辨率。
灰度校正图graynew=rgb2gray(rgbnew);
(3)边缘检测:
边缘是指图像局部亮度变化显著的部分,是图像风、纹理特征提取和形状特征提取等图像分析的重要基础。
所以在此我们要对图像进行边缘检测。
图象增强处理对图象牌照的可辩认度的改善和简化后续的牌照字符定位和分割的难度都是很有必要
实现代码:
graynew=rgb2gray(rgbnew);
graynew=imfill(graynew,'
holes'
fori=1:
3
se1=[111;
111;
111];
graynew=imdilate(graynew,se1);
se=[1;
1;
1];
graynew=imerode(graynew,se);
%腐蚀Imerode(X,SE).其中X是待处理的像,%SE是结构元素对象
graynew=getcenter(graynew);
graynew=double(graynew);
[rowcol]=size(graynew);
1.2车牌定位
牌照的定位是在经图象预处理后的原始灰度图象中确定牌照的具体位置,并将包含牌照字符的一块子图象从整个图象中分割出来,供字符识别子系统识别之用。
图为定位出的车牌
车牌裁剪:
colArray=zeros(row,1);
row
col
if(graynew(i,j,1)==1)
colArray(i,1)=colArray(i,1)+1;
[maxcountmaxIndex]=max(colArray);
roughY=maxIndex;
while((colArray(roughY,1)>
=5)&
(roughY>
1))
roughY=roughY-1;
y1=roughY;
(roughY<
row))
roughY=roughY+1;
y2=roughY;
rowArray=zeros(1,col);
forj=y1:
y2
ifgraynew(j,i)==1
rowArray(1,i)=rowArray(1,i)+1;
[maxcountmaxIndex]=max(rowArray);
roughX=maxIndex;
while((rowArray(1,roughX)>
(roughX>
roughX=roughX-1;
x1=roughX;
(roughX<
col))
roughX=roughX+1;
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('
%读取图像
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('
开运算后图像'
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,[00500100]);
%对图像进行剪切,保留图像的%一条直线,减小运算量
theta=0:
179;
%定义theta角度围
r=radon(bw1,theta);
%对图像进行Radon变换
figure,imshow(r)
[m,n]=size(r);
c=90;
m
forj=1:
n
ifr(1,1)<
r(i,j)
r(1,1)=r(i,j);
c=j;
end%检测Radon变换矩阵中的峰值所对应的列坐标
rot=90-c;
%确定旋转角度
A=imrotate(I1,rot,'
crop'
%对图像进行旋转矫正
2.3计算行方向的像素,用来去除行方向上的边框
计算行像素:
Y1=zeros(y,1);
if(I4(i,j,1)==1)
Y1(i,1)=Y1(i,1)+1;
%计算I3第j列有几个一
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;
%找到去除边框后上边的位置
while((Y1(Py1,1)>
(Py1<
y))
Py1=Py1+1;
%找到去除边框后下边的位置
Z1=I4(Py0:
Py1,:
:
%将二值图像上下边框去除
figure(3);
imshow(Z1),title('
将二值图像上下边框去除后图像'
2.4计算列方向的像素,用来去除垂直方向上的边框
[y,x,z]=size(Z1);
%计算此时图像的大小
计算列像素
X1=zeros(1,x);
if(Z1(i,j,1)==1)
X1(1,j)=X1(1,j)+1;
figure(4);
plot(0:
x-1,X1),title('
列方向像素点灰度值累计和'
列值'
累计像素'
去除垂直边框
x1=fix(x/2)+1;
5
while(i~=5)
while((X1(1,x1)>
=15)&
(x1<
x))
x1=x1+1;
%找到去除边框后右边的位置
i=i+1;
while((X1(1,x1)<
15)&
x)&
(i<
4))
%从车牌中间开始寻找字符间隙,直到找到第四个间隙%为止,即找到去除边框后车牌右边的位置
x0=fix(x*45/440);
%找到第一个字符的位置
2
while(i~=2)
while((X1(1,x0)>
(x0>
x0=x0-1;
%从车牌的第一个字符开始寻找第一个字符间隙,找到%去除边框左边的位置
while((X1(1,x0)<
5)&
1)&
i~=2)
%找到去除边框后左边的位置
Z2=Z1(:
x0:
x1,:
%将二值图像左右边框去除
figure(5);
imshow(Z2),title('
将二值图像垂直边框去除后图像'
2.5再次计算列方向的像素,用来去除车牌上的圆点
[y,x,z]=size(Z2);
if(Z2(i,j,1)==1)
End
去除车牌中的圆点:
Z2=double(Z2);
x1=fix(x*105/409);
x2=fix(x*128/409);
forj=x1:
x2
Z2(i,j)=0;
figure(6);
去除圆点后图像'
2.6用投影法分割车牌,把分割出来的子图转成标准子图
用投影法分割
Px0=1;
Px1=1;
y=fix(x*90/409);
7
while((X1(1,Px0)<
10)&
(Px0<
Px0=Px0+1;
%找到待分割字符左边的位置
Px1=Px0;
a=1;
while(((X1(1,Px1)>
=10)&
(Px1<
x))||((Px1<
((a/y)<
=0.5)))
Px1=Px1+1;
%找到待分割字符右边的位置
a=Px1-Px0;
Z3=Z2(:
Px0:
Px1,:
%二值化图像分割后
figure(7);
subplot(1,7,i);
imshow(Z3);
%将二值化图像分割后显示出来
2.7变换为标准子图
Z4=imresize(Z3,[4020]);
%将分割后的二值图变换为标准子图
figure
(1);
imshow(Z4),title('
标准子图'
%将标准子图显示出来
Px0=Px1;
imwrite(Z4,strcat('
imgAfterSplit/'
num2str(i),'
.jpg'
));
%把标准子图存到imagAfterSplit中
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;
forI=1:
ii=int2str(I);
%将整型数据转换为字符串型数据
t=imread(['
ii,'
%依次读入七位车牌字符
SegBw2=imresize(t,[4020],'
nearest'
%对读入的字符进行缩放
SegBw2=im2bw(SegBw2);
ifI==1%第一位汉字识别
kmin=37;
kmax=67;
elseifI==2%第二位A~Z字母识别
kmin=11;
kmax=36;
else%第三位以后是字母或数字识别;
即I>
=3
kmin=1;
index=getFinalTemp(kmin,kmax,codes,SegBw2,I);
CarNum(l)=codes(index);
l=l+1;
%xlabel(['
第三步:
识别结果为:
CarNum],'
Color'
'
b'
getFinalTemp.m程序代码(返回最终识别的单个字符index下标):
function[index]=getFinalTemp(minIndex,maxIndex,codes,SegBw2,numIndex)
ifnumIndex<
3%后五位不能出现i、o字母;
fori=minIndex:
maxIndex
path=strcat('
template/'
codes(i));
%调用getFromLib遍历文件夹找出同一字符最佳的模板
maxSame=getFromLib(path,SegBw2);
near(i)=maxSame;
%记录下字符与模板i同的点个数
else
%ifi~=find(codes='
l'
)&
i~=find(codes='
O
ifcodes(i)~='
codes(i)~='
O'
nearest=max(near);
%差别最小的,不同字符模板
index=find(near==nearest);
%找出差别最小的模板的下标
getFromLib.m程序代码(遍历单个字符模板库,返回最佳的模板):
function[maxSame]=getFromLib(path,carNum)
files=dir([path'
/*.jpg'
length(files)
file=imread([path'
/'
num2str(i)'
]);
ifislogical(file)==0
filebw=im2bw(file);
else
filebw=file;
subImg=carNum-filebw;
%以上相当于两幅图相减得到第三幅图进行匹配
sameDots=subImg==0;
%记录下字符与模板一样的点个数
sumDots=sum(sameDots(:
arr(i)=sumDots;
maxSame=max(arr);
4、系统界面的实现:
界面我们小组一共分有了两个模块:
一个是父模块车牌识别系统,进行显图,定位,分割和识别的操作;
另一个是子模块是字符模块,是为了增加字符的模块进库,提高车牌的识别效率。
父模块——车牌识别系统
图一
界面分别有载入图片,车牌定位,车牌分割,车牌识别和增加字符模块五个按钮,每个按钮对应不同的功能。
载入车牌:
选择你需要识别的带有车牌的图片并显示出来在原始图片的位置。
车牌定位:
部进行定位函数对车牌进行定位,并在原图上用黄线圈出车牌相应的位置。
车牌分割:
将定位好的车牌进行分割,分割出7个独立的字符模块。
车牌识别:
对每个分割出来的字符进行识别,并文本输出来。
增加字符模块:
点击进入子界面,进行字符模块的增加,提高车牌识别的效率。
效果如下图:
图二
子模块——字符模块
图三
字符模块分为5大部分,分别是车牌显示,选择添加模板的位置,选择要添加的模板,添加字符模块和显示字符模块。
车牌显示:
将刚在父模块里车牌分割出来的图片显示在这个对应的位置上。
选择添加模板的位置:
选择你需要添加的字符在图片对应的位置,如“粤”在第一位就选择“1”。
选择要添加的模板:
选择你需要添加进库的字符对应的数字,字母或者文字,如要添加第一位“粤”字就选择粤“粤”字。
添加字符模板:
当你选择好“选择添加模板的位置”和“选择要添加的模板”后,就可以按这个按钮确认了。
如图四
显示字符模块:
添加成功后你就可以点击这个按钮查看对应的某字符的字符模板库里的所有“粤”字模板了。
如图五
图四
图五
四、设计结果分析
从上面结果看车牌识别成功,但有一些车牌