车牌识别matlab程序.doc

filename, pathname = uigetfile(*.jpg, JPEG文件(*.jpg);*.bmp,BMP文件(*.bmp););if(filename = 0), return, endglobal FILENAME %声明全局变量FILENAME = pathname filename;I=imread(FILENAME);imshow(I); %显示图像I I1=rgb2gray(I);%RGB图转化为灰度图 figure,imshow(I1); w1=medfilt2(I1);figure,imshow(w1); s1=histeq(w1,256);figure,imshow(s1); t1=imadjust(s1); figure,imshow(t1); I2=edge(t1,robert,0.15,both); %用ROBERT算子提取图像边缘 figure,imshow(I2); se=1;1;1;I3=imerode(I2,se); %弱化二进制图像I2的边缘figure,imshow(I3);%为定位车牌，将白色区域膨胀，腐蚀去无关的小物件，包括车牌字符(下面两句)se=strel(rectangle,25,20); %用来腐蚀的形状为矩形，面积20*25I4=imclose(I3,se);figure,imshow(I4); I5=bwareaopen(I4,2000); %去除图像中面积过小的，可以肯定不是车牌的区域。 figure,imshow(I5); y,x,z=size(I5); %z=1。 y，x分别为I5图像的高和宽 myI=double(I5)%myI=I5;tic %begin横向扫描 white_y=zeros(y,1); %white_y为y行1列的零矩阵 for i=1:y for j=1:x if(myI(i,j,1)=1) %如果myI(i,j,1)即myI图像中坐标为(i,j)的点为白色 %则white_y的相应行的元素white_y(i,1)值加1 white_y(i,1)= white_y(i,1)+1; end end end temp MaxY=max(white_y); %temp为向量white_y的元素中的最大值，MaxY为该值的索引（在向量中的位置） PY1=MaxY; while (white_y(PY1,1)=80)&(PY11) PY1=PY1-1; end PY2=MaxY; while (white_y(PY2,1)=80)&(PY2y) PY2=PY2+1; end IY=I(PY1:PY2,:,:); %IY为原始图像I中截取的纵坐标在PY1：PY2之间的部分 %end横向扫描 %begin纵向扫描 white_x=zeros(1,x); for j=1:x for i=PY1:PY2 if(myI(i,j,1)=1) white_x(1,j)= white_x(1,j)+1; end end end PX1=1; while (white_x(1,PX1)3)&(PX1x) PX1=PX1+1; end PX2=x; while (white_x(1,PX2)PX1) PX2=PX2-1; end %end纵向扫描 PX1=PX1-2; PX2=PX2+2; Plate=I(PY1:PY2,:); t=toc figure,plot(white_y);grid figure,plot(white_x);grid figure,imshow(IY); figure,imshow(Plate);% if isrgd(Plate); Plate3=rgb2gray(Plate); %else % Plate2=Plate; %end plate4=medfilt2(Plate3); plate5=histeq(plate4,256);Plate2=imadjust(plate5); g_max=double(max(max(Plate2); g_min=double(min(min(Plate2); t=round(g_max-(g_max-g_min)/2); m,n=size(Plate2); Plate2=im2bw(Plate2,t/256); figure,imshow(Plate2); plate=bwareaopen(Plate2,20); figure,imshow(plate); y1,x1,z1=size(plate); plate1=double(plate);% tt=1; Y1=zeros(y1,1); for i=1:y1 for j=1:x1 if(plate1(i,j,1)=1) Y1(i,1)=Y1(i,1)+1; end end end py1=1; py0=1; while(Y1(py0,1)20)&(py0=20)&(py1y1) py1=py1+1; end plate=plate(py0:py1,:,:); figure,imshow(plate); X1=zeros(1,x1); for j=1:x1 for i=1:y1 if(plate1(i,j,1)=1) X1(1,j)=X1(1,j)+1; end end end figure plot(0:x1-1,X1) px0=1; px1=1; for i=1:7 while(X1(1,px0)3)&(px0=3)&(px1x1)|(px1-px0)10) px1=px1+1; end Z=plate(:,px0:px1,:); switch strcat(Z,num2str(i) case Z1 PIN0=Z; case Z2 PIN1=Z; case Z3