一、车牌识别简介

车牌识别技术起源于20世纪80年代初期，图像模式识别和计算机视觉在其中起到至关重要的作用。随着计算机技术的迅速普及和进步，80年代中后期，车牌识别系统逐步投入市场并使用，但识别精度和速度都不理想。目前，发达国家的车牌识别系统已广泛用于市场，其中以色列和新加坡公司的车牌识别系统较为领先。由于车牌设计的不同，不存在一种通用的车牌检测技术。因中国车牌包括汉字，使得国外车牌识别系统不能直接用于国内车牌识别，需要中国自主研发。中国的车牌识别技术研究起步略晚于西方国家。目前国内技术领先的是中科院的“汉王眼”和香港的视觉科技公司。但是，这些投入市场的车牌识别系统的使用都有一定的局限性，在车牌图像质量差、车牌倾斜、光照条件不理想等情况时，识别率会大幅度降低。
为了克服上述问题，本文对不同光照条件下获得的车牌图像进行了识别和分析。由于车牌涉及保密信息，暂时没有公开的大规模数据集可供使用和实验对比。本文先对车牌数据进行采集，统一处理320×240分辨率的图片，建立测试用数据库。数据库分为两类：测试库1为光照条件较好，无阴影车牌；测试库2为光照条件较差，有阴影车牌。采用传统的模板匹配算法建立车牌识别系统，模板图片采用中值滤波算法进行平均处理，统一建立分辨率为20×40的识别模板库。为方便用户使用， 利用Matlab建立了图形用户交互界面(GUI) 。通过对测试库车牌实验， 可以看到车牌成像质量对识别率
有很大影响。因此在建立车牌识别系统时，适当角度的补光可有效提高识别的准确率。
1 图像的预处理
由于车牌多是通过交通监控等条件获得，因周边交通环境、拍摄角度、光照和实时性要求等多种原因，得到的照片直接用于车牌的识别准确率难以得到保证。因此，预处理图像可以提高识别的准确性。
1.1图像的灰度化
由相机拍摄的牌照通常是彩色图像， 即RGB图像为了提高车牌识别系统的速度、节省内存， 先将RGB图像转换为灰度图.红、蓝、绿三种基本颜色按不同的比例可以组成任意一种颜色，所以车牌照片的每一个颜色都可由不同比例的红、蓝、绿组成。灰度图像可以用数组I表示，数组且I的数据类型一般有整数和双精度两种。通常0代表黑色，255代表白色。采用适合于人类视觉系统的转换方法，如公式(1)所示。

其中，系数a，b，c要大于0，且a+b+c=1，Y表示灰度值，R、G、B分别表示红色、绿色和蓝色，a=0.299、b=0.584、c=0.117.图1显示了带有车牌照片的原始车牌图像，图2显示了在公式(1)的灰度处理之后车牌图像的灰度图像。

图1 车牌原图

图2 车牌灰度图
1.2 图像二值化
图像二值化不仅可以大大减少数据量，还可以突出图像的目标轮廓，有利于后续的图像定位和分割处理。在车牌图像的二值化处理中，被确定为目标区域的像素具有大于或等于阈值的灰度值，计算如公式(2)所示。经过测试，当阈值th=0.76时，效果较好，结果如图3所示。

图3 车牌的二值化图
1.3 图像的边缘检测
本文采用Roberts算子来进行边缘检测。采用该算子的计算量小， 速度快， 便于后续的实时处理， 其模板如表1所示。若梯度幅度G(x，y)大于设定的阈值，则判断为边缘。G(x，y)计算式如公式(3)所示，此处阈值选择为0.15。
表1 Roberts算子模板

其中，f(x，y)是图像空间(x，y)处的灰度值，G(x，y)是f(x，y)的梯度幅度。图4是边缘检测后的结果。

图4 车牌的边缘检测图

2车牌识别系统的设计
2.1车牌定位
车牌定位采用颜色特征提取，一般的车牌区域都具有很明显的特点，中国的车牌以蓝底白字居多。传统的方法一般根据车牌的色彩特征，彩色像素点统计的方法分割出合理的车牌区域假设经相机采集包含车牌的RGB图像， 水平方向记为y，垂直方向记为x。首先，确定车牌各分量分别对应的颜色范围。其次，计算水平方向上对应的像素数量和车牌的合理面积。然后，计算分割的水平方向区域中垂直方向上该颜色范围内白的像素数量，并为定位设置合理的阈值。最后，根据对应方向的范围确定车牌区域。但是这种方法的准确率较低，本文采用文献[6]提出的颜色特征提取及二值化归类方法进行算法改进，计算方法如公式(4)所示。图5为定位后的车牌。
其中的蓝色特征可用式(4)表示：

图5 车牌定位图
2.2车牌字符识别
在字符识别前，先对车牌字符分割。字符分割采用阈值分割，主要包含两个步骤【7]。
(1)确定所需要进行分割的阈值。
(2)将阈值与每个点的灰度值对照，以达到分割目的。
将分割后的图像进行归一化处理，可以有效地将字符图像的大小进行缩放以得到大小一致的字符图像，便于后续的字符识别。
模板匹配一般是数字图像处理中最常使用的识别方法之一，先建立模板库，再将字符输入到模板中寻找与之最佳匹配的模板字符。
模板图像大小为20×40，模板库由数字0-9，32个省份的简称汉字和大写英文字母(其中O一般不用作车牌字母)三部分组成。模板的部分图像如图6所示。

图6 模板库部分图片

二、部分源代码

function varargout = gui(varargin)
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @gui_OpeningFcn, ...
                   'gui_OutputFcn',  @gui_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
%结束初始化
function gui_OpeningFcn(hObject, ~, handles, varargin)
handles.output = hObject;
guidata(hObject, handles);
function varargout = gui_OutputFcn(~, ~, handles) 
varargout{1} = handles.output;

% ======================输入图像===============================
function pushbutton1_Callback(hObject, ~, handles)
[filename, pathname]=uigetfile({'*.jpg';'*.bmp'}, 'File Selector');
I=imread([pathname '\' filename]);
handles.I=I;
guidata(hObject, handles);
axes(handles.axes1);
imshow(I);title('原图');

% ======================图像处理===============================
function pushbutton2_Callback(~, ~, handles)
I=handles.I;
I1=rgb2gray(I);
I2=edge(I1,'roberts',0.18,'both');
axes(handles.axes2);
imshow(I1);title('灰度图');
axes(handles.axes3);
imshow(I2);title('边缘检测');
se=[1;1;1];
I3=imerode(I2,se);%腐蚀操作
se=strel('rectangle',[25,25]);
I4=imclose(I3,se);%图像聚类，填充图像
I5=bwareaopen(I4,2000);%去除聚团灰度值小于2000的部分
[y,x,z]=size(I5);%返回15各维的尺寸，存储在x,y,z中
myI=double(I5);
tic      %tic计时开始，toc结束
Blue_y=zeros(y,1);%产生一个y*1的零针
for i=1:y
    for j=1:x
        if(myI(i,j,1)==1)%如果myI图像坐标为（i，j）点值为1，即背景颜色为蓝色，blue加一
            Blue_y(i,1)=Blue_y(i,1)+1;%蓝色像素点统计
        end
    end
end
[temp MaxY]=max(Blue_y);
%Y方向车牌区域确定
%temp为向量yellow_y的元素中的最大值，MaxY为该值得索引
PY1=MaxY;
while((Blue_y(PY1,1)>=5)&&(PY1>1))
    PY1=PY1-1;
end
PY2=MaxY;
while((Blue_y(PY2,1)>=5)&&(PY2<y))
    PY2=PY2+1;
end
IY=I(PY1:PY2,:,:);
%X方向车牌区域确定
Blue_x=zeros(1,x);%进一步确认x方向的车牌区域
for j=1:x
    for i=PY1:PY2
        if(myI(i,j,1)==1)
            Blue_x(1,j)=Blue_x(1,j)+1;
        end
    end
end
PX1=1;
while((Blue_x(1,PX1)<3)&&(PX1<x))
    PX1=PX1+1;
end
PX2=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;
axes(handles.axes4);imshow(dw),title('定位车牌');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
imwrite(dw,'dw.jpg');%将彩色车牌写入dw文件中
a=imread('dw.jpg');%读取车牌
b=rgb2gray(a);%将车牌图像转换为灰度图
imwrite(b,'灰度车牌.jpg');%将灰度图写入文件
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,'二值化.jpg');
%均值滤波前
%滤波
h=fspecial('average',3);
%建立预定义的滤波算子，average为均值滤波，模板尺寸为3*3
d=im2bw(round(filter2(h,d)));%使用指定的滤波器h对h进行d即均值滤波
imwrite(d,'均值滤波.jpg');
%某些图像进行操作
%膨胀或腐蚀
se=eye(2);%单位矩阵
[m,n]=size(d);%返回信息矩阵
if bwarea(d)/m/n>=0.365%计算二值图像中对象的总面积与整个面积的比是否大于0.365
    d=imerode(d,se);%如果大于0.365则进行腐蚀
elseif bwarea(d)/m/n<=0.235%计算二值图像中对象的总面积与整个面积的比值是否小于0.235
    d=imdilate(d,se);%%如果小于则实现膨胀操作
end
imwrite(d,'膨胀.jpg');
%寻找连续有文字的块，若长度大于某阈值，则认为该块有两个字符组成，需要分割
d=qiege(d);
[~,n]=size(d);
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;
    end
    k2=j-1;
    if k2-k1>=round(n/6.5)
        [val,num]=min(sum(d(:,[k1+5:k2-5])));
        d(:,k1+num+5)=0;%分割
    end
end
%再切割
d=qiege(d);
%切割出7个字符
y1=10;
y2=0.25;
flag=0;
word1=[];
while flag==0
    [m,n]=size(d);
    left=1;
    wide=0;
    while sum(d(:,wide+1))~=0
        wide=wide+1;
    end
    if wide<y1 %认为是左干扰 
        d(:,[1:wide])=0;
        d=qiege(d);
    else
        temp=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>y2
            flag=1;word1=temp;%word1
        end
        d(:,[1:wide])=0;d=qiege(d);
    end
end

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社，2020.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社，2013.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社，2013.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社，2015.
[5]呙润华，苏婷婷，马晓伟.BP神经网络联合模板匹配的车牌识别系统[J].清华大学学报(自然科学版)，2013，53(9)：1221-1226.
[6]鲁扬.基于BP神经网络的车牌识别算法研究[D].大庆：东北石油大学，2018.
[7]李强,张娟.一种改进的基于模板匹配的污损车牌识别方法[J].智能计算机与应用. 2019,9(03).
[8] 梁凯.基于MATLAB的汽车车牌识别系统的设计与实现[D] .哈尔滨：黑龙江大学， 2018.
[9]刘雄飞，朱盛春.车牌字符多特征提取与BP神经网络的识别算法[J].计算机仿真，2014，31(10)：161-164，290.
[10] 曾泉， 谭北海.基于SVM和BP神经网络的车牌识别系统[J] .电子科技， 2016， 29(1) ：98-101.