实现车牌自动定位-Read

1、实现车牌自动定位刚才在 url=/url付费网站上找到的一段程序代码. 具体怎么样没研究过 , 是否可行也有待考证. 在说下我现在在做毕业设计 , 寻找各种各样的这方面的资料 .( 现在在看模式识别和智能系统方面的书 ). 要是有一起研究的请加QQ:178043048.大家一起讨论讨论 ,O0j N6_ y6a+c 这个程序我没有经过修改. 原来就是这样子的 . $d L3o S0| 7Z%L%程序功能：实现车牌自动定位%作者：重庆大学 田建国 r-M2S:e,Y (O%编写时间 :2007.04.10016 021%修改时间 :2007.07.06;修改内容 : 增加合并区域操作，第88 行

2、开始-A.U%修改时间： 2007.12.21 ；修改内容： 77，84 行 对车牌长宽比阈值进行修改，%同时对候选车牌区域增加长宽限制，即该区域不能占整个图像的比例不能太小和太大%对检测出来无车牌的情况进行提示B C6p6N8a% 可修改阈值如下 :threshold1, threshold2,功能 : 区域合并的条件 函数 mycombine的参数%threshold3,功能 : 候选车牌区域长宽比阈值a)c7PU/r%threshold4功能 : 车牌区域占整个图像比列阈值9s4(Zclear;clc;close all%阈值设置 , 可修改查看阈值具体最用 Ctrl+F查找相关阈值名

3、%可根据实际情况修改threshold1=0.5;threshold2=1/20;n m/Be3tthreshold3=2.5;threshold4=1/7,1/3,1/12,1/6;%顺序为宽的最小最大 , 高的最小最大 i.F K$Mfilename=0003.bmp;I=im2gray(filename);%调用自编函数读取图像，并转化为灰度图象；-+q A,r,F Vtic%计时开始height,width=size(I);%预处理I_edge=zeros(height,width);r-?9gi5R3d6afor i=1:width-1I_edge(:,i)=abs(I(:,i+1)

4、-I(:,i);end l6n B K J N 4RI_edge=(255/(max(max(I_edge)-min(min(I_edge)*(I_edge-min(min(I_edge);I_edge,y1=select(I_edge,height,width);%调用select函数BW2 = I_edge;%一些形态学处理8h2n:Gc/ASE=strel(rectangle,6,6);IM2=bwareaopen(BW2,20);%开运算，消除细小物体_)m6YAIM2=imerode(IM2,SE);% 腐蚀IM3=imdilate(IM2,SE);%膨胀-W&e6o-*m.oSE=

5、strel(rectangle,5,5);IM3=imdilate(IM3,SE);%膨胀%计算区域的外接矩形boundaries = bwboundaries(IM3);%调用函数 bwboundariesk=length(boundaries);hkkk=k;rectange=cell(1,k);for i=1:ku&GQL+BI/mQ9Sb = boundariesi;,v%P%-ez8| m;Crectangei(1)=min(b(:,1);rectangei(2)=max(b(:,1);%矩形的行起点与终点rectangei(3)=min(b(:,2);rectangei(4)=ma

6、x(b(:,2);%矩形的列起点与终点 -B r7O+O r K7Xrec_centeri=fix(rectangei(1)+rectangei(2)/2),fix(rectangei(3)+rectangei(4)/2);%矩形中心rec_ratio(i)=(rectangei(4)-rectangei(3)/(rectangei(2)-rectangei(1);%矩形长宽比 _;? 5endbp=rectange;%wSg+H,b%对上面参数和变量的说明：rectange 为一胞元，用于存放每个图像块的左上和右下两个点的坐标；%存放格式为：rectangek=x1,x2,y1,y2；x1,

7、x2分别为行坐标，y1,y2为列坐标1E+f,m4O%rec_center 为一胞元 , 用于存放每个图像块的中心坐标,p_centerk=x,y;x,y分别为行 , 列坐标 37?1Y&U L4Z%rec_ratio为一矩阵，用来存放图像块的长宽比例%合并临近区域 %h-q+U0RT6V(K&m4Mrectange,rec_center,rec_ratio=mycombine(rectange,height,width,threshold1,threshold2);%调用自编函数 mycombine，其中最后两个为阈值，需要调整k=length(rectange); %更新区域个数bp=re

8、ctange;1f-D$Z:dZ d5dyXa2_kkk=k;%合并结束%1q?!B.P6T-b:iD.q3p m%根据区域比例判断是否为车牌区域%p+LU(jk!m=1;T=0.6*max(rec_ratio);%0.8参数需要调整 pn,Y ?T=min(T,threshold3);%阈值选取for i=1:k9Ax2g;O6b%imshow(I(rectangei(1):rectangei(2),rectangei(3):rectangei(4)if(rec_ratio(i)=T&rec_ratio(i)=width*threshold4(1) &disw=height*threshol

9、d4(3) &dish1)rectange=plate;$_;1O*C0#;N-nk=l-1; *!G F-y.rend V a 2k$ r:g +H%判定结束 %7o7m-ttoc%计时结束 2V!.P5c +clear edge_IM3;clear x; clear y;% 清空部分变量T*i4|g2u/p(?#vn8Dif(k=1)V$Bh#e%imwrite(I(index(1)-2:index(2),index(3):index(4),cp.jpg);e%n Lx0Xend 4B%显示figure;z g%eQ#Hm;U%I(rectange1(1)-2:rectange1(1),r

10、ectange1(3)-2:rectange1(4)+2)=255;a1q!t_O s* G g3jP.X*g%I(rectange1(2):rectange1(2)+2,rectange1(3)-2:rectange1(4)+2)=255;n;Mb1u-U 8Z P%I(rectange1(1)-2:rectange1(2)+2,rectange1(3)-2:rectange1(3)=255;%I(rectange1(1)-2:rectange1(2)+2,rectange1(4):rectange1(4)+2)=255;J#X eimshow(I);:C c8|+K!K -x f8simw

11、rite(I,I.jpg)a!t7qsubplot(221);imshow(I);I(rectange1(1),rectange1(3):rectange1(4)=255;h7P9FI(rectange1(2),rectange1(3):rectange1(4)=255;subplot(222);imshow(BW2);subplot(223);imshow(IM2);K&L)R;M0wsubplot(224);imshow(IM3);hold onfor i=1:kkk,W95R_6B 3bplot(bpi(3):bpi(4),bpi(1)*ones(1,bpi(4)+1-bpi(3),g); Iw)Hplot(bpi(3):bpi(4),bpi(2)*ones(1,bpi(4)+1-bpi(3),g);-up1D-Q*? y/W;j _6Zc z0Cplot(bpi(3)*ones(1,bpi(2)+1-bpi(1),bpi(1):bpi(2),g);plot(bpi(4)*ones(1,bpi(2)+1-bpi(1),bpi(1):bpi(2),g)end%显示figure,A9f #D8oh Ofor i=1:k