MATLAB大作业

更新时间:2024-05-04 15:11:01 阅读量: 综合文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

基于MATLAB识别图片中的文字 知识点:图像识别 图像处理 计算机视觉

实验目的:以MATLAB为工具,识别图片中的文字。 问题描述:

交通是一个十分重要的问题,红绿灯上的摄像头可以辨别车主的身份。现在,我们是否能够用matlab设计一个图像识别的系统,通过车牌号的来识别相应的数字。 问题分析及模型建立: 第一步:确定车牌号的区域 对于一张图像

在matlab中,在处理图像元素时用(x,y,z)三坐标的形式表达的灰度集时,其中,x和y是空间坐标,z是f 在坐标(x,y)处的值,就可以表达图像在该坐标轴上的点。再将,将f(x,y)的数值简单地显示为一个矩阵,就可以定量地表达了一幅数值图像。

矩阵中的每个元素称为像素。

所以,如果我们想要截取一段车牌号的图像,我们就可以扫描图像每一个点的像素。然后我们进行恢复处理,算子边缘检测,腐蚀,闭操作,删除小对象得到

我们就可以利用RGB的值找到白色区域的边界,然后我们利用边界,截取我们的原图像就可以得到我们所需要的车牌号的区域了。

第二步:进行文字的分割

在进行分割之前,就需要将我们所获得的车牌号的区域的图像进行进一步的处理,突出我们需要的文字部分。 灰度处理

彩色的图片占用的空间比较大,处理会花费很长的时间,先将图片进行灰度处理

二值化

再将图片二值化,将图片的256个灰度级强制减少到只有{0, 255}2个灰度级,更加方便了后

续的文字识别识别。

均值滤波

均值滤波是对是对图像进行局部平均, 以平均值来代表该像素点的灰度值。矩形滤波器对这个二维矢量的每一个分量进行独立的平滑处理。获得一张比较干净清晰的图像。

腐蚀和膨胀

腐蚀可以分割独立的图像元素,膨胀用于连接相邻的元素,更加利于图像的分割 进行完图像的预处理,就可以开始图像的切割了

首先进行边缘的切割,处理完后的图像为黑底白字。其中黑色为0,白色1为。这样我们就可以从四个边缘进行扫描,若行或者列像素之和0的话,就为背景,这样就可以继续向内推移,直到和不为0的位置为止。这样就进行了边缘的切割。

接下来,先进行文字的切割,文字和字母、数字在长度和宽度上有较大的差别,所以需要单独的拿出来。因为号码是平铺过来的,所以我们进行列的扫描。先从最左边进行扫描,若列像素之和不为0的话,则说明有文字,直到和为0的时候停止。但是如果扫描出来的宽度太小的话就说明背景可能有白色区域干扰。这时,应该将刚刚扫描出来的区域置黑,再次进行扫描,截取图像。同理,数字和字母可同样按着上面的方式扫描,截取。

截取后的图像

模板的配对

所有的文字数字和字母都截取下来后,我们需要将所得的图片统一的量化为同样大小的图像方便比对。

统一量化后的图像

完成了这部的操作之后,就需要建立一个小型的模板库。将车牌号常见的文字、十个数字和26个英文字母存储为上面同样大小的二值化图片,并且以对应字符的名称命名。

将所有的模板放在一起,按照一定的顺序排列。然后,我们就可以调用我们已经切割好的图片和模板内部的图片进行相减,寻找差别最小的图片进行匹配,这样就可以识别出车牌号上

的字符了。

问题求解: matlab代码:

主函数

function [d]=main(jpg) close all clc

I=imread('Car.jpg');%读入车牌图像Car.jpg I1=rgb2gray(I);

I2=edge(I1,'roberts',0.15,'both'); se=[1;1;1]; I3=imerode(I2,se);

se=strel('rectangle',[25,25]); I4=imclose(I3,se); I5=bwareaopen(I4,2000);

figure(1),imshow(I5);title('从对象中移除小对象');

[y,x,z]=size(I5);%图像以(y,x,z)表达灰度集,x,y为坐标,z为对应位置的值 myI=double(I5);%转化为双精度,便于确定范围

tic % 测定算法执行的时间,开始计时 Blue_y=zeros(y,1); for i=1:y for j=1:x

if(myI(i,j,1)==1) % 寻找y轴方向的白色区域

Blue_y(i,1)= Blue_y(i,1)+1; % 是蓝色区域的则进行计数 end end end

[temp MaxY]=max(Blue_y);%Y方向车牌区域确定 temp(最多点数):所有行中,最多的累积 PY1=MaxY;

while ((Blue_y(PY1,1)>=5)&&(PY1>1)) PY1=PY1-1; end %Y轴方向的上限 PY2=MaxY;

while ((Blue_y(PY2,1)>=5)&&(PY2

%%%%%% X方向 %%%%%%%%%

Blue_x=zeros(1,x);%进一步确定x方向的车牌区域

for j=1:x

for i=PY1:PY2 %只需要扫描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

end %确定x轴的右边界 PX2=x;

while ((Blue_x(1,PX2)<3)&&(PX2>PX1)) PX2=PX2-1; end %确定x轴的左边界 PX1=PX1-1;%对车牌区域的校正 PX2=PX2+1;

dw=I(PY1:PY2-8,PX1:PX2,:); %确定图片的截取区域 t=toc; % 读取程序的运行时间

a=imread('dw.jpg'); %读入已经截取好的图像 A=size(a);

if length(A)==3 b=rgb2gray(a); else b=a; end

figure(8);subplot(3,2,1),imshow(b),title('1.车牌灰度图像') 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,'2.车牌二值图像.jpg');

figure(8);subplot(3,2,2),imshow(d),title('2.车牌二值图像') figure(8),subplot(3,2,3),imshow(d),title('3.均值滤波前') % 滤波

h=fspecial('average',3); %进行均值滤波 imwrite(d,'4.均值滤波后.jpg');

figure(8),subplot(3,2,4),imshow(d),title('4.均值滤波后')

% 某些图像进行操作 % 膨胀或腐蚀

se=eye(2); % eye(n)可返回一个2*2单位矩阵 [m,n]=size(d);%图像为二值图像,x,y为坐标

if bwarea(d)/m/n>=0.365 %计算二值图像中对象的总面积 d=imerode(d,se); %实现图形的腐蚀 elseifbwarea(d)/m/n<=0.235

d=imdilate(d,se);%实现图像的膨胀 end

imwrite(d,'5.膨胀或腐蚀处理后.jpg');

figure(8),subplot(3,2,5),imshow(d),title('5.膨胀或腐蚀处理后')

% 寻找连续有文字的块,若长度大于某阈值,则认为该块有两个字符组成,需要分割 d=qiege(d); %截掉边界黑色的部分 [m,n]=size(d); %图像已经切割的图像的坐标 figure,subplot(2,1,1),imshow(d),title(n) 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; % k1+num+5列的位置赋值为黑色,分割 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

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; % WORD 1 end

d(:,[1:wide])=0;d=qiege(d); end end

% 分割出第二个字符

[word2,d]=getword(d);%getword的函数与上面类似,字母、数字和文字在一些细节上需要去区分

% 分割出第三个字符 [word3,d]=getword(d); % 分割出第四个字符 [word4,d]=getword(d); % 分割出第五个字符 [word5,d]=getword(d); % 分割出第六个字符

[word6,d]=getword(d); % 分割出第七个字符 [word7,d]=getword(d);

subplot(5,7,1),imshow(word1),title('1'); subplot(5,7,2),imshow(word2),title('2'); subplot(5,7,3),imshow(word3),title('3'); subplot(5,7,4),imshow(word4),title('4'); subplot(5,7,5),imshow(word5),title('5'); subplot(5,7,6),imshow(word6),title('6'); subplot(5,7,7),imshow(word7),title('7'); [m,n]=size(word1);

% 将生成的图片归一化大小为 40*20 word1=imresize(word1,[40 20]); word2=imresize(word2,[40 20]); word3=imresize(word3,[40 20]); word4=imresize(word4,[40 20]); word5=imresize(word5,[40 20]); word6=imresize(word6,[40 20]); word7=imresize(word7,[40 20]);

subplot(5,7,15),imshow(word1),title('1');

subplot(5,7,16),imshow(word2),title('2'); subplot(5,7,17),imshow(word3),title('3'); subplot(5,7,18),imshow(word4),title('4'); subplot(5,7,19),imshow(word5),title('5'); subplot(5,7,20),imshow(word6),title('6'); subplot(5,7,21),imshow(word7),title('7'); imwrite(word1,'1.jpg'); imwrite(word2,'2.jpg'); imwrite(word3,'3.jpg'); imwrite(word4,'4.jpg'); imwrite(word5,'5.jpg'); imwrite(word6,'6.jpg'); imwrite(word7,'7.jpg');

liccode=char(['0':'9' 'A':'Z' '苏豫陕鲁晋']); %建立自动识别字符代码表 SubBw2=zeros(40,20); @*20的零矩阵 l=1; for I=1:7

ii=int2str(I); %int2str整型转换字符串 t=imread([ii,'.jpg']);

SegBw2=imresize(t,[40 20],'nearest'); %改变图片的大小 if l==1 %第一位汉字识别 kmin=37;

kmax=41; %模板中汉字所在的位置 elseif l==2 %第二位 A~Z 字母识别 kmin=11;

kmax=36; %模板中数字所在的位置

else l>=3 %第三位以后是字母或数字识别 kmin=1;

kmax=36; %模板中字母所在的位置

end

for k2=kmin:kmax

fname=strcat('字符模板\\',liccode(k2),'.jpg'); SamBw2=imread(fname);

Error(k2)=sum(sum(abs(SegBw2-SamBw2)));

%将现有字符逐个与模板字符相减,认为相减误差最小的现有字符与该模板字符匹配 end

Error1=Error(kmin:kmax); MinError=min(Error1);

findc=find(Error1==MinError);%寻找最匹配的字符

Code(l*2-1)=liccode(findc(1)+kmin-1);%将上述找到的字符位置与liccode中的字符链接上

Code(l*2)=' ';%在字符之间加空格 l=l+1; end

figure(10),imshow(a),title (['车牌号码:', Code],'Color','b');

函数getword

function [word,result]=getword(d) word=[];flag=0;y1=8;y2=0.5; while flag==0 [m,n]=size(d); wide=0;

while sum(d(:,wide+1))~=0 && wide<=n-2 wide=wide+1; end

temp=qiege(imcrop(d,[1 1 wide m])); [m1,n1]=size(temp); if widey2 d(:,[1:wide])=0; if sum(sum(d))~=0

d=qiege(d); % 切割出最小范围 else word=[];flag=1; end else

word=qiege(imcrop(d,[1 1 wide m])); d(:,[1:wide])=0; if sum(sum(d))~=0; d=qiege(d);flag=1; else d=[]; end end end %end

result=d;

qiege函数

function e=qiege(d) [m,n]=size(d); top=1; %下端 bottom=m; %下端 left=1; %左端 right=n; % 右端

while sum(d(top,:))==0 && top<=m %如果矩阵d的第top行的和为0的话 top=top+1; end

while sum(d(bottom,:))==0 && bottom>=1 %如果矩阵d的第bottom行的和为0的话 bottom=bottom-1;

end

while sum(d(:,left))==0 && left<=n %如果矩阵d的第left列的和为0的话 left=left+1; end

while sum(d(:,right))==0 && right>=1 %如果矩阵d的第right列的和为0的话 right=right-1; end

dd=right-left; %切割的宽度 hh=bottom-top; %切割的高度

e=imcrop(d,[left top ddhh]); %返回该截取的区域 结果分析:

1, 这种垂直扫描的方式,图像的中车的牌照需要接近水平,如果倾斜的比较厉害,很难只

截取车牌号的区域。

2, 扫描字符的时候,如果扫描左右结构的字符时,比如“沪”,则会截出‘氵’和‘户’两个部分,

就无法匹配字符。

3, 如果遇到受污染很严重的图片,就难切割图片的和模板中的图片匹配。

参考文献

[1]艾冬梅,李艳晴,张丽静,刘琳MATLAB与数学实验[M]机械工业出版社 [2]Rafael C.Gonzalez Richard E.Woods数字图像处理[M]电子工业出版社 [3]樊昀, 王润生.从图像中提取文字[J] .国防科技大学学报,2002,24 (1) :59-62. [4]邹浩,余龙,邹勇博,刘宇童,和振乔,李少梅.基于MATLAB的图片中文字的提取

及识别[J].西安电子科技大学 电子工程学院,西安,710126

本文来源:https://www.bwwdw.com/article/oghg.html

Top