【身份牌识别】基于matlab GUI模板匹配身份牌识别【含Matlab源码 1354期】

一、模板匹配算法简介

1 概述
模式识别就是通过计算机，用数学模型求解的方法研究模式的自动处理和判读。在模式识别的各种方法中，模板匹配是最容易的一种，其数学模型易于建立，通过模板匹配对数字图像模式识别有助于我们了解数学模型在数字图像中的应用。

2 模板匹配算法
2.1 相似性测度求匹配
模板匹配的实际操作思路很简单：拿已知的模板，和原图像中同样大小的一块区域去对。最开始时，模板的左上角点和图像的左上角点是重合的，拿模板和原图像中同样大小的一块区域去对比，然后平移到下一个像素，仍然进行同样的操作， ……所有的位置都对完后,差别最小的那块就是我们要找的物体。

以上所描述的是相似性测度法求匹配的求解思路，其在计算机中操作的如图2所示。设模板T叠放在搜索图上平移，被模板覆盖搜索图下的那个图像叫做子图Si , j，i , j 为这块子图的左上角像素点在S图的坐标，称为参考点，从图2可知，i , j 的取值范围是：1<i ，j <N- M+1. 现在可以比较T和Si , j的内容。若两者一致,则T和S之差为零. 因此,可用下列公式(1) 和公式(2) 来衡量T和Si , j的相似程度。

在(2) 式中第3项表示模板总能量，是一个与(i , j) 无关的常数；第1项是模板覆盖下子图的能量，它随着(i , j) 的位置缓慢地改变；第2项表示的子图与模板的互相关系，随着(i , j) 的改变而改变，当T和Si , j匹配时这项取值最大。因此可用下列相关函数(3) 作相似性测度。

当矢量t 和S1之间的夹角为0时,即当S1(i , j) =kt 时(k为常量) ,有R(i , j) =1,否则R(i , j) <1. 显然R(i , j) 越大,模板T和Si , j就越相似,点(i , j) 就是我们要寻找的匹配点。

2.2 序贯相似性检测的算法
用相关法求匹配的计算量很大，因为模板要在(N- M+1)2个参考位置上作相关计算，除了在匹配点外，其它点作的都是无用功。因此，人们提出一种叫序贯相似性检测的算法，简称SSDA(Sequential SimiliarityDetectionAlgorithm) 其要点是：
在数字图像中，SSDA法用公式(6) 计算图像f ( x, y) 在点(i , j) 的非相似度m(i , j) 作为匹配尺度。式中(i , j) 表示的不是模板中心坐标，而是它左上角坐标。模板的大小为n ×m。

如果在(i , j) 处图像中有和模板一致的图案时，则m(i , j) 值很小，反之则很大。特别是模板和搜索图下的子图部分，完全不一致的场合下，如果在模板内的各像素与图像重合部分对应的像素灰度差的绝对值依次增加，其和会急剧增大。因此，在作加法时，如果灰度差的绝对值部分和超过某一阈值时，就认为这个位置不存在和模板一致的图案，从而转移到下一个位置上计算m(i , j)。并且在这模板下的各像素点计算中止，因此能大幅度地缩短计算时间,提高匹配速度。
根据上述思路，我们可以进一步改进SSDA算法。就是把在图像上的模板移动分为粗检索和细检索2个阶段进行。首先进行粗检索，它不是让模板每次移动1个像素，而是每隔若干的像素把模板和图像重叠，并且计算匹配的尺度，从而求出待寻找的图案大致存在的范围。然后在这个范围内，让模板每隔1个像素移动1次，根据匹配的尺度确定寻找图案的所在位置。这样，整体上计算模板匹配的次数减少了，计算时间缩短，匹配的速度提高了。但用这种方法具有漏掉图像中最恰当位置的危险性。

2.3 相关算法
2个函数的相关性定义,可用公式(7)表示：

f*表示f 的复共轭。我们知道相关理论类似卷积理论，F( u, v) 和 H( u, v) 分别表示f ( x, y) 和h( x, y)的傅立叶变换. 根据卷积理论有

可知卷积是空间域过滤和频率域过滤之间的纽带。相关的重要用途在于匹配。在匹配中，f ( x, y) 是一幅包含物体或区域的图像。如果想要确定f 是否包含有感兴趣的物体或区域，让h( x, y) 作为那个物体的区域(通常称该图像为模板)。如果匹配成功，2个函数的相关值会在h找到f 中相应点的位置上达到最大。从上面分析可知，相关算法可以有2种方法：可在空间域进行，也可在频率域进行。

2.4 幅度排序相关算法
这种算法有2个步骤组成：
第1步，把实时图中的各个灰度值按幅度的大小排成列的形式，然后在对它进行二进制(或三进制) 编码，根据二进制排序的序列，把实时图变换为二进制阵列的一个有序的集合{ Cn, n =1,2, …, N}。这一过程称之为幅度排序预处理。
第2步，序贯地将这些二进制序列与基准图进行由粗到细的相关，直到确定出匹配点为止。由于篇幅的限制，这里就不列出例子了。

2.5 分层搜索的序惯判决算法
这种分层搜索算法是直接基于人们先粗后细寻找事物的惯例而形成的，例如，在中国地图上找肇庆的位置时，可以先找广东省这个地域，这过程称为粗相关。然后在这个地域中，再仔细确定肇庆的位置，这叫做细相关。很明显，利用这种方法，可以很快找出肇庆的位置。因为在这过程中省略了寻找广东省以外区域所需的时间，这种方法称为分层搜索的序贯判决法，利用这种思想形成的分层算法具有相当高的搜索速度。限于篇幅,这里只给出这种操作的一般思路。

3 总结
模式匹配本质就是应用数学。模板匹配过程如下：①将图像数字化,按顺序取出每个点的像素值；②代入事先建立的数学模型进行预处理；③选择一种合适的算法进行模式匹配；④将匹配后图像的坐标列出或直接在原图显示.。在模板匹配操作中最为关键的是如何建立数学模型,这是正确匹配的核心。

二、部分源代码

%运行界面部分
warning('off');
close all;
chos=0;%初始界面按钮参数
possibility=3;%退出界面按钮参数
while chos~=possibility
chos=menu('身份牌信息识别系统','选择身份牌建立模板','  选择身份牌开始进行识别  ','》》结束《《');%各个按钮名称，第一个为该图形化窗口的名称
if chos==1 %如果获取的按键为'选择身份牌建立模板'按键
clc;
clear all;
close all;
warning off
[fn,pn,fi]=uigetfile('*.*','选择图片');%选择图片
bmp=imread([pn fn]);%显示原始图像
figure(1)
imshow(bmp);title('原始名牌图像');
%%%%%%%%%%图像HSV处理%%%%%%%%%%
mysize=size(bmp);%判断图像色度
if numel(mysize)>2 %判断是否是rgb图像
myhsv=rgb2hsv(bmp);%对图像进行HSV处理
h=myhsv(:,:,1);
s=myhsv(:,:,2);
v=myhsv(:,:,3);
%%%%%%%%%%图像中值滤波%%%%%%%%%%
bmp=medfilt2(v,[3,3]);
%%%%%%%%%%图像二值处理%%%%%%%%%%
A=im2bw(v,0.3);%对图像进行二值化处理
elseA=bmp;
end
%%%%%%%%图像做形态学处理%%%%%%%%%
zfwz=[100 320 150 100];
B=bwmorph(A,'spur',inf);%对图像进行形态学处理
B=im2bw(B,0.5);
rectangle('position',zfwz,'edgecolor','r');%字符画框
C=imcrop(B,zfwz);%画框部分图像截取
C=imcomplement(C);%对图像进行黑白反转
figure(2)
imshow(C);
pics=dir('*.bmp');MBS=length(pics);
[pathstr,name,ext]= fileparts(mfilename('fullpath'));
imwrite(C,[pathstr,'\',num2str(MBS),'.bmp'],'bmp')
chos=0;
possibility=3;
endif chos==2  %如果获取的按键为'选择身份牌开始进行识别'按键
clc;
clear all;
close all;
warning off
[fn,pn,fi]=uigetfile('*.*','选择图片');%选择图片
bmp=imread([pn fn]);%显示原始图像
figure(1)
imshow(bmp);title('原始名牌图像');
%%%%%%%%%%图像HSV处理%%%%%%%%%%
myhsv=rgb2hsv(bmp);%对图像进行HSV处理
h=myhsv(:,:,1);
s=myhsv(:,:,2);
v=myhsv(:,:,3);
figure(2)
subplot(131);imshow(h);title('HSV转化后的H分量');
subplot(132);imshow(s);title('HSV转化后的S分量');
subplot(133);imshow(v);title('HSV转化后的V分量');
%%%%%%%%%%图像中值滤波%%%%%%%%%%
bmp=medfilt2(v,[3,3]);
% figure(2);subplot(221);imshow(bmp);title('中值滤波图像');
%%%%%%%%%%图像二值处理%%%%%%%%%%
A=im2bw(v,0.5);%对图像进行二值化处理
%%%%%%%%%%%%%%%%%%%%%倾斜校正%%%%%%%%%%%%%%%%%%%%%%
qingxiejiao=rando_bianhuan(A);
A=imrotate(A,qingxiejiao,'bilinear','crop');
figure,imshow(A);title('倾斜校正');%取值为负值向右旋转
[pathstr,name,ext]= fileparts(mfilename('fullpath'));
imwrite(A,[pathstr,'\xxx.bmp'],'bmp')
%%%%%%%%%%%%%%%%%%%%%裁去黑边%%%%%%%%%%%%%%%%%%%%%%
a = imread([pathstr,'\xxx.bmp']);        %读取图片
delete([pathstr,'\xxx.bmp']);
A=caiquhei(a);
figure(3),imshow(A);
%%%%%%%%图像做形态学处理%%%%%%%%%
zfwz=[100 320 150 100];
B=bwmorph(A,'spur',inf);%figure,imshow(bw);title('擦除之前');
B=im2bw(B,0.5);
figure(4)
imshow(B);
rectangle('position',zfwz,'edgecolor','r');%字符画框
C=imcrop(B,zfwz);%画框部分图像截取
C=imcomplement(C);%对图像进行黑白反转
figure(5)
imshow(C);
mobanpun=maopaopp(imread('0.bmp'),imread('1.bmp'),imread('2.bmp'),imread('3.bmp'));
jieguohanzi  = mobanpp(mobanpun,C);
hs=msgbox(['名牌信息识别结果为识别结果为：',jieguohanzi],'result');
pause(1);
close(hs);
%disp(jieguohanzi)
chos=0;
possibility=3;
end end
%用于定位
clc;
clear all;
close all;
warning off
m=1;%图片的亮度系数对与识别成功率影响很大，如果图片本身对比度很高值设为1即可，如果光线较暗可适当提高数值反之亦然。
[fn,pn,fi]=uigetfile('*.*','选择图片');%选择图片
bmp=imread([pn fn]);%显示原始图像
figure(1)
imshow(bmp);
title('原始名牌图像');
%%%%%%%%%%图像HSV处理%%%%%%%%%%
myhsv=rgb2hsv(bmp);%对图像进行HSV处理
h=myhsv(:,:,1);
s=myhsv(:,:,2);
v=myhsv(:,:,3);
%%%%%%%%%%图像中值滤波%%%%%%%%%%
bmp=medfilt2(v,[3,3]);
% figure(2);subplot(221);imshow(bmp);title('中值滤波图像');
%%%%%%%%%%图像二值处理%%%%%%%%%%
A=im2bw(v,0.3);%对图像进行二值化处理
%%%%%%%%图像做形态学处理%%%%%%%%%
zfwz=[100 320 150 100];
B=bwmorph(A,'spur',inf);%figure,imshow(bw);title('擦除之前');
B=im2bw(B,0.5);
imshow(B);
B=imcrop(B,[100 300 200 150]);%将含有名牌欣喜的部分剪裁下来
grayimg = B;BWimg = grayimg;[width,height]=size(grayimg);%灰度图数据传递，导入宽高
A=im2bw(B,0.5);     % thresh=0.5 表示将灰度等级在128以下的像素全部变为黑色，将灰度等级在128以上的像素全部变为白色。
bw = edge(A,'sobel','vertical'); figure(2); imshow(bw);title("边缘图像"); % 垂直边缘检测
Z = strel('rectangle',[30,18]);%连通域处理
bw_close=imclose(bw,Z);%连通域闭处理
bw_open = imopen(bw,Z);%连通域开处理
showImg = grayimg;%灰度图数据传递
%图像数据二值化处理,如果图像尺度过大会使处理时间几何倍增加，所以尽量使用小尺度图片，（你的图片完全OK）我提供个测试图片就属于过大的了
for i=1:width  for j=1:heightif(BWimg(i,j) == 255)showImg(i,j)= grayimg(i,j);else showImg(i,j)= 0;endend
end
figure(3);
imshow(showImg);%可视化显示二值化图像
[l,m] = bwlabel(bw_close);  %连续域标签
status=regionprops(l,'BoundingBox');%图像区域度量
centroid = regionprops(l,'Centroid');%字符区域度量
imshow(B);hold on;
a=[0 0 0 0];
if m>1
for i=1:m
if status(i).BoundingBox(1,3) > status(i).BoundingBox(1,4)%归一化预处理status(i).BoundingBox(1,4)=status(i).BoundingBox(1,3);
elsestatus(i).BoundingBox(1,3)=status(i).BoundingBox(1,4);
endstatus(i).BoundingBox=status(i).BoundingBox+a;rectangle('position',status(i).BoundingBox,'edgecolor','g');%字符画框text(centroid(i,1).Centroid(1,1)-25,centroid(i,1).Centroid(1,2)-25, num2str(i),'Color', 'r') %字符框标号
end
end

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社，2020.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社，2013.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社，2013.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社，2015.