车牌识别matlab程序1

[, pathname] = uigetfile({'*.jpg', 'JPEG文件(*.jpg)';'*.bmp','BMP文件(*.bmp)';});

if( == 0), return, end

global %声明全局变量

= [pathname ];

I=imread();

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*25

I4=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)&&(PY1>1))

PY1=PY1-1;

end

PY2=MaxY;

while ((white_y(PY2,1)>=80)&&(PY2

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)&&(PX1

PX1=PX1+1;

end

PX2=x;

while ((white_x(1,PX2)<3)&&(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

end

py1=py0;

while((Y1(py1,1)>=20)&&(py1

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)