车牌识别中的车牌定位

该程序参考的是一篇硕士论文：王璐《基于MATLAB的车牌识别系统研究》，并参考了百百度文库中的作者成果，链接为http://wenku.baidu.com/view/23dde718fc4ffe473368abb8.html 

以下是程序，每一步都标识出操作目的，做起来很流畅，但是经过多方尝试，发觉这段代码只能识别出画质清晰，车牌面积比重较大的情况下，所以目前的工作还有待改进。

close all
clear all
clc
%%%3.1 车辆图像预处理
I=imread('C:\Users\Administrator\Desktop\car1.jpg');
gray=rgb2gray(I);
gray=histeq(gray);
gray=medfilt2(gray);
s=strel('disk',5);
Bgray=imopen(gray,s);
Egray=imsubtract(gray,Bgray);

%%%3.2 车牌边缘提取(小波变换)
Egray=double(Egray);
[c,s]=wavedec2(Egray,2,'db1');
cv2=detcoef2('v',c,s,2);

%%%3.3.1结构元素提取
a1=imclearborder(cv2,8);
st=ones(1,14);
bg1=imclose(a1,st);
bg2=imopen(bg1,st);
bg3=imopen(bg2,[1,1,1,1,1,1,1]');
bg=uint8(bg3);
level = graythresh(bg);
bw2=im2bw(bg,level);
se=strel('disk',3);
bw=imopen(bw2,se);
figure(1)
subplot(1,3,1);imshow(Egray);
subplot(1,3,2);imshow(bg3,[]);
subplot(1,3,3);imshow(bw,[]);

%%%3.3.2提取候选区域
[L,num]=bwlabel(bw,8);
feature=regionprops(L,'basic');
area=[feature.Area];
boundbox=[feature.BoundingBox];
RGB=label2rgb(L,'spring','k','shuffle');
lx=0;
for l=1:1:num
col=boundbox((l-1)*4+1);     
row=boundbox((l-1)*4+2); 
width=boundbox((l-1)*4+3);
hight=boundbox((l-1)*4+4);
rato=width/hight;    
if (row>5)&(col>5)&(rato>2)&(rato<4)
lx=lx+1;break;
end
end

%%显示车牌区域
bwD=bw((row-1):(row+hight-1),(col-1):(col+width-1));
grayD=gray((4*row-20):(4*row+4*hight+10),(4*col-20):(4*col+4*width+10));
figure(2);
subplot(1,2,1);imshow(bwD);
subplot(1,2,2);imshow(grayD);

%%%3.4.1 车牌水平方向定位算法
%%水平方向一阶差分
[m,n]=size(grayD);
f=zeros(m,n-1);
for i=1:1:m
line=grayD(i,:);
f(i,:)=diff(line);
f(i,:)=abs(f(i,:));
end
%水平投影
F=sum(f,2);
w=fspecial('average',[3,1]);
Th=imfilter(F,w);
Th=uint8(Th);
Th=im2bw(Th,0.6);

for i=1:1:m
if (Th(i+1,1)-Th(i,1))==1
strow=i+1;
else if (Th(i+1,1)-Th(i,1))==-1
endrow=i;
break;
end
end
end

%%竖直方向一阶差分
g=zeros(m-1,n);
for j=1:1:n
line=grayD(:,j);
g(:,j)=diff(line);
g(:,j)=abs(g(:,j));
end
%竖直投影
G=sum(g);
w=fspecial('average',[1,3]);
Tv=imfilter(G,w);
Tv=uint8(Tv);
Tv=im2bw(Tv,0.4);

for j=1:1:n
if (Tv(1,j+1)-Tv(1,j))==1
strcol=j+1;
else if (Tv(1,j+1)-Tv(1,j))==-1
endcol=j;
break;
end
end
end

%%%车牌精确定位后的图像
Des=grayD(strow:endrow,strcol:endcol);
figure(3);
subplot(2,1,1);imshow(grayD);
subplot(2,1,2);imshow(Des);

原图

figure(1)

figure(2)

figure(3)