License Plate Recognition.rar_license plate_车牌_车牌识别

close all clc % 清空命令窗口的所有输入和输出，类似于清屏 %I=uigetfile('.jpg','输入一个需要识别的车牌图像'); I=imread('D:\Program Files\MATALAB\R2013b\img6.jpg');%读取图像? figure();subplot(3,2,1),imshow(I),title('原始图像'); I1=rgb2gray(I);%转化为灰度图像?subplot(3,2,2),imshow(I1),title('灰度图像');? %I2=edge(I1,'roberts',0.09,'both');%采用robert算子进行边缘检测? I2=edge(I1,'roberts',0.08,'both');%采用robert算子进行边缘检测? subplot(3,2,3),imshow(I2),title('边缘检测后图像'); se=[1;1;1];%线型结构元素? I3=imerode(I2,se);%腐蚀图像? subplot(3,2,4),imshow(I3),title('腐蚀后边缘图像'); se=strel('rectangle',[25,25]);%矩阵结构元素 I4=imclose(I3,se);%图像聚类、填充图像? subplot(3,2,5),imshow(I4),title('填充后图像'); I5=bwareaopen(I4,2000);%去除聚团灰度值小于2000的部分? subplot(3,2,6),imshow(I5),title('形态滤波后图像'); [y,x,z]=size(I5); I6=double(I5); Y1=zeros(y,1); for i=1:y for j=1:x if(I6(i,j,1)==1) Y1(i,1)=Y1(i,1)+1; end end end [temp MaxY]=max(Y1); figure(); subplot(3,2,1),plot(0:y-1,Y1),title('行方向像素点灰度值累计和'),xlabel('行值'),ylabel('像素'); PY1=MaxY; %while ((Y1(PY1,1)>=50)&&(PY1>1)) while ((Y1(PY1,1)>=50)&&(PY1>1)) PY1=PY1-1; end PY2=MaxY; %while ((Y1(PY2,1)>=50)&&(PY2<y)) while ((Y1(PY2,1)>=50)&&(PY2<y)) PY2=PY2+1; end IY=I(PY1:PY2,:,:); X1=zeros(1,x); for j=1:x for i=PY1:PY2 if(I6(i,j,1)==1) X1(1,j)=X1(1,j)+1; end end end subplot(3,2,2),plot(0:x-1,X1),title('列方向像素点灰度值累计和'),xlabel('列值'),ylabel('像数'); PX1=1; while ((X1(1,PX1)<3)&&(PX1<x)) PX1=PX1+1; end PX2=x; while ((X1(1,PX2)<3)&&(PX2>PX1)) PX2=PX2-1 end PX1=PX1-1; PX2=PX2+1; %PX2=PX2-1; dw=I(PY1:PY2,PX1:PX2,:); %dw=I(PY1:PY2-5,PX1:PX2,:); subplot(3,2,3),imshow(dw),title('定位剪切后的彩色车牌图像') if isrgb(dw) %Im_original=imread('len.ppm'); %if (ndims(Im_original)==3) I1 = rgb2gray(dw);%将RGB图像转化为灰度图像? else I1=dw; end g_max=double(max(max(I1))); g_min=double(min(min(I1))); T=round(g_max-(g_max-g_min)/3);%?T?为二值化的阈值? [m,n]=size(I1);%?d:二值图像? %h=graythresh(I1);? I1=im2bw(I1,T/256); subplot(3,2,4); imshow(I1),title('二值化车牌图像'); I2=bwareaopen(I1,20); subplot(3,2,5); imshow(I2),title('形态学滤波后的二值化图像'); [y1,x1,z1]=size(I2); I3=double(I2); TT=1; %%%%%%%去除图像顶端和底端的不感兴趣区域%%%%% Y1=zeros(y1,1); for i=1:y1 for j=1:x1 if(I3(i,j,1)==1) Y1(i,1)= Y1(i,1)+1 end end end Py1=1; Py0=1; while ((Y1(Py0,1)<20)&&(Py0<y1)) Py0=Py0+1; end Py1=Py0; while((Y1(Py1,1)>=20)&&(Py1<y1)) Py1=Py1+1; end I2=I2(Py0:Py1,:,:); subplot(3,2,6); imshow(I2),title('目标车牌区域'); X1=zeros(1,x1); for j=1:x1 for i=1:y1 if(I3(i,j,1)==1) X1(1,j)= X1(1,j)+1; end end end figure(5); plot(0:x1-1,X1),title('列方向像素点灰度值累计和'),xlabel('列值'),ylabel('累计像素量'); Px0=1; Px1=1; %%%%%%%%%%%%分割字符%%%%%%%%%%%%%%%%%%? for i=1:7 while ((X1(1,Px0)<3)&&(Px0<x1)) Px0=Px0+1; end Px1=Px0; while(((X1(1,Px1)>=3)&&(Px1<x1))||((Px1-Px0)<10)) Px1=Px1+1; end Z=I2(:,Px0:Px1,:); switch strcat('Z',num2str(i)) case 'Z1' PIN0=Z; case 'Z2' PIN1=Z; case 'Z3' PIN2=Z; case 'Z4' PIN3=Z; case 'Z5' PIN4=Z; case 'Z6' PIN5=Z; otherwise PIN6=Z; end figure(3); subplot(1,7,i); imshow(Z); Px0=Px1; end