【车牌识别】基于matlab GUI模板匹配新能源车牌识别【含Matlab源码 865期】

举报
海神之光 发表于 2022/05/29 03:36:07 2022/05/29
【摘要】 一、简介 1 在进行图像处理任务之前,首先要明确处理对象和处理流程 处理对象: 新能源车牌和蓝色车牌 处理流程: 2 由于要对蓝色和新能源车牌进行识别,而且新能源车牌是渐变色,所以直接通过设定阈值的方...

一、简介

1 在进行图像处理任务之前,首先要明确处理对象和处理流程
处理对象: 新能源车牌和蓝色车牌
处理流程:
在这里插入图片描述
2 由于要对蓝色和新能源车牌进行识别,而且新能源车牌是渐变色,所以直接通过设定阈值的方法去识别两种车牌并没有单独识别蓝色车牌的可靠性高。这里,我根据车牌字符的特征进行车牌位置定位:
2.1 根据字符特征初步找出车牌位置
在这里插入图片描述
2.2 根据车牌长宽比进行初步筛选
我们可以用最小外接矩形算法框出上一步已选出的区域,然后根据车牌长宽比设定阈值,将满足阈值的矩形筛选出来,当然但还是会存在一部分干扰区域。

2.3 对矩形内像素进行颜色统计
通过对初步筛选后剩余区域分别进行像素的颜色统计,我们即可准确找出车牌位置,与此同时车牌的类型(新能源或蓝牌)也可以在这里判断出来

3 倾斜校正
这一步是为后面垂直投影法分割字符和模板匹配做准备

简单来讲,使用霍夫线算法检测车牌上下边框所在直线,计算这条直线倾斜角度,之后对车牌进行旋转即可。当然复杂点的方法可以具体看我的源程序,具体效果如下图:
在这里插入图片描述

4 二值化、形态学处理、字符识别
4.1 二值化和形态学处理这里就不在赘述,很多博客已经讲的非常清晰
4.2 字符分割
在这里需要对字符进行准确的分割,在分割之前还是可能会有干扰,需要再次分割出车牌字符的准确位置。这里我的方法是:先向左投影,设定阈值,切割出上下边界,同理在切割出左右边界。但是这只适用于车牌边框明显的情况下,所以我不建议你使用这个方法,慎用!!!
在这里插入图片描述
到这里就可以使用垂直投影法进行字符分割,这里需要注意的是新能源车牌有八位,蓝牌有七位
分割结果如下图所示:
在这里插入图片描述
4.3 字符识别
这里使用模板匹配,因为新能源车牌字符和蓝牌字符不一致,所以要制作两套字符模板。另外需要注意的是,根据《中华人民共和国机动车号牌》序号编码规则,26个英文字符中“I“、”O“不能使用,所以不存在"1"和"I"、"0"和"O"字符混淆的问题

二、源代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% resize的目标尺寸650*X
% 调参
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;


se1 = strel('rectangle',[15 15]);                                 %创建一个平坦的矩形结构
se2 = strel('rectangle',[12 20]);                                 %创建一个平坦的矩形结构
se3 = strel('rectangle',[12 15]);
imgraw=imread('C:\Users\lenovo\Desktop\licenseplate\0.png');                %读取图

% se1 = strel('rectangle',[15 25]);                                 %创建一个平坦的矩形结构
% se2 = strel('rectangle',[10 25]);                                 %创建一个平坦的矩形结构
% se3 = strel('rectangle',[14 19]);
% imgraw=imread('E:\matlab workspace\licenseplate\车牌2.jpg');                %读取图像


[m,~]=size(imgraw);
times = 650/m;
img=imresize(imgraw,times,'nearest');
g= fspecial('gaussian',[5,5],1);                                  %高斯滤波器
g_img = imfilter(img, g, 'replicate');                            %图像滤波
gray_img=rgb2gray(g_img);                                         %转为灰色图像        
open = imopen(gray_img,se1);                                      %图像开操作
sub_img = imsubtract(gray_img,open);                              %图像相减
BW_img = imbinarize(sub_img);                                     %使用Otsu方法确定全局阈值进行二值化
edge_img = edge(BW_img,'canny');
edge_img1 = imclose(edge_img,se2);                                %图像闭操作
edge_img2 = imopen(edge_img1,se3);                                %图像开操作

%%%%%%%%%%%%%%%%%%%%%%%%%%%使用车牌长宽比进行判断%%%%%%%%%%%%%%%%%%%%%%%%%%%

[B,L] = bwboundaries(edge_img2,'noholes');  
num_rect = 0;                                                               %符合要求的矩形数量
matrix_rect = {};                                                           %符合要求的矩形的坐标元胞矩阵
for k = 1:length(B)
    boundary = B{k};
    c = boundary(:,2);
    r = boundary(:,1);
    [rectx,recty,area,perimeter] = minboundrect(c,r,'a');                    %rectx,recty为逆时针排列
    width = round(sqrt((recty(1)-recty(2))^2+(rectx(1)-rectx(2))^2));
    leth = round(sqrt((recty(2)-recty(3))^2+(rectx(2)-rectx(3))^2));
    if width > leth
        temp = width;
        width = leth;
        leth = temp;
    end
    ratio = leth / width;
    if ratio >2.5 && ratio <4.5
        num_rect = num_rect + 1;
        matrix_rect{num_rect} = [rectx,recty];
    end      
end    

subplot(451); imshow(img); title('原始图像');
subplot(452); imshow(g_img);    title('高通滤波后图像');
subplot(453); imshow(gray_img);    title('灰度图像');
subplot(454); imshow(open);    title('图像开操作');
subplot(455); imshow(sub_img);    title('图像相减');
subplot(456); imshow(BW_img);    title('图像二值化');
subplot(457); imshow(edge_img);    title('图像边缘');
subplot(458); imshow(edge_img2);    title('图像开闭操作');
subplot(459); imshow(img);    title('原图框出车牌');
hold on
for k = 1:num_rect
    rect = matrix_rect{k};
    rx = rect(:,1);
    ry = rect(:,2);
    line(rx,ry,'Color','red');
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%判断颜色、二次判断车牌%%%%%%%%%%%%%%%%%%%%%%%%%%%

p=0;  %计算满足颜色条件的次数
color = 'no';
for o=1:num_rect
    z=matrix_rect{o};
    mid1x=(z(1,1)+z(2,1))/2;
    mid1y=(z(1,2)+z(2,2))/2;
    mid2x=(z(3,1)+z(4,1))/2;
    mid2y=(z(3,2)+z(4,2))/2;
    green=0;blue=0;
    hsv= rgb2hsv(img);  %RGB转换为HSV
    img_count = (abs(ceil(mid2x+0.601)-floor(mid1x+0.601)))*(abs(ceil(mid2y+0.001)-floor(mid1y+0.001)));
    subplot(4,5,10);  imshow(hsv);

    for i=ceil(mid1y+1):floor(mid2y-1)
        if mid2x==mid1x
            for j=floor(mid2x+1):ceil(mid1x+1)
                h = hsv(i,j,1);
                s = hsv(i,j,2);
                v = hsv(i,j,3);
                if 0.40< h && h <= 0.46 &&  s>0.4 && s<=1 && v>0.3 && v<1
                    green = green+1;
                elseif 0.56<h && h<=0.72 && s>0.4 && s<=1 && v>0.4 && v<1
                    blue =blue +1;                
                end
            end       
            
        else
            for j=ceil(mid2x+1):floor(mid1x-1)
                h = hsv(i,j,1);
                s = hsv(i,j,2);
                v = hsv(i,j,3);

                if 0.40< h && h <= 0.46 &&  s>0.4 && s<=1 && v>0.3 && v<1
                    green = green+1;
                elseif 0.56<h && h<=0.72 && s>0.4 && s<=1 && v>0.4 && v<1
                    blue =blue +1;                
               end
            end
        end
    end
    if green*2 >= img_count && green>blue
        color = 'green';
        p=p+1;
        b=z;   %满足条件的元胞
    elseif blue*2 >= img_count && blue>green
        color = 'blue';
        p=p+1;
        b=z;   %满足条件的元胞
    end
end

if (color == 'blue')&&(p == 1)
    flag = 0;
    style = '识别到的是蓝色';
elseif (color == 'green')&&(p == 1)
    flag = 1;
    style = '识别到的是新能源';
end
function image=cutting(img,flag) 
    %获取图片大小
    [row,col] = size(img); 
    %设定边界
    img(:,col)=1;
    img(:,1)=1; 
    %图像取反,方便后续处理,根据黑色判断边界
    all = sum(~img);  %行向量
    %figure('NumberTitle','off','Name','投影'),bar(a),title('投影');
    left = 1;
    right = 1;
    count1 =0; %记录字符个数
    for i = 1:col-1  %从左到右扫描
        if all(i)==0&&all(i+1)~=0 %左侧边界
            left = i;
        end
        if all(i)~=0&&all(i+1)==0 %右侧边界
            right=i;
        else
            right =0; %左右两侧均为黑色边框
        end
        if right~=0
            count1 = count1+1;
            left_bound(count1) = left; %存放左边界
            right_bound(count1) = right; %存放右边界
        end    
    end
    for i = 1:count1 %???
        if left_bound(i)<fix(col/8)
            left_bound(i)=left_bound(1);
        end
    end
    count2 =1; %k记录满足条件的字符个数
    right_boundary = zeros(1,count1); %存放右边界
    left_boundary = zeros(1,count1); %存放左边界
    for i = 1:count1
        if (right_bound(i) - left_bound(i))>(fix(col/15)) %判断字符宽度是否满足要求 问题:字符宽度怎么求
            right_boundary(count2) = right_bound(i);
            left_boundary(count2) = left_bound(i);
            count2 = count2+1;
        end
    end


  
 
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
  • 115
  • 116
  • 117
  • 118
  • 119
  • 120
  • 121
  • 122
  • 123
  • 124
  • 125
  • 126
  • 127
  • 128
  • 129
  • 130
  • 131
  • 132
  • 133
  • 134
  • 135
  • 136
  • 137
  • 138
  • 139
  • 140
  • 141
  • 142
  • 143
  • 144
  • 145
  • 146
  • 147
  • 148
  • 149
  • 150
  • 151
  • 152
  • 153
  • 154
  • 155
  • 156
  • 157
  • 158
  • 159
  • 160
  • 161
  • 162
  • 163
  • 164
  • 165
  • 166
  • 167
  • 168
  • 169
  • 170
  • 171
  • 172
  • 173
  • 174
  • 175

三、运行结果

在这里插入图片描述

四、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]乔有田.基于数学形态学和灰度跳变特征的车牌定位方法[J].工业控制计算,2015(2):77-78

文章来源: qq912100926.blog.csdn.net,作者:海神之光,版权归原作者所有,如需转载,请联系作者。

原文链接:qq912100926.blog.csdn.net/article/details/116301231

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。