matlab车牌识别代码

【核心代码】

car=imread('car1.jpg');
imshow(car);
%color image to gray
car1=rgb2gray(car);
figure,imshow(car1);
car2=histeq(car1);
figure,imshow(car2);
%edge detection
temp=edge(car1,'sobel');
figure,imshow(temp);
[x,y]=size(temp);
d=2;
temp1=temp;
%continuous peoperty test
for i=1:x
    for j=d 1:y-d
        if(2*temp(i,j)-temp(i,j-d)-temp(i,j d)==0)
            temp1(i,j)=0;
        else
            temp1(i,j)=temp1(i,j);
        end
    end
end
figure,imshow(temp1);
%character compensation
se=ones(3);
temp2=imclose(temp1,se);
figure,imshow(temp2);
temp3=temp2;
th=20;
th1=70;
%property transformation
for i=1:x
    bp=[];
    bp1=[];
    pos=[];
    len=[];
    ind=find(temp2(i,:)==1);
    [m,n]=size(ind);
    index=0;
    if n~=0
        for j=1:n-1
            if ind(1,j 1)-ind(1,j)>th
                index=index 1;
                bp(1,index)=ind(j 1);
                bp1(1,index)=ind(j);
            end
        end
        [x1,y1]=size(bp);
        if y1~=0
            pos(1,1)=ind(1,1);
            pos(2:y1 1)=bp;
            len(1,1)=bp1(1,1)-pos(1,1);
            for k=1:y1-1
                len(1,k 1)=bp1(1,k 1)-bp(1,k);
            end
            len(1,y1 1)=ind(1,n)-pos(1,y1 1);
            for k=1:y1 1
                if len(k)>=th1
                    temp2(i,pos(k):pos(k) len(k))=1;
                else
                    temp2(i,pos(k):pos(k) len(k))=0;
                end
            end
        elseif ind(n)-ind(1)>=th1
            temp2(i,ind(1):ind(n))=1;
        else
            temp2(i,ind(1):ind(n))=0;
        end
    end
end
figure,imshow(temp2);
%final locate
connComp=bwconncomp(temp2);
features=regionprops(connComp,'basic');
borders=[features.BoundingBox];
area=[features.Area];
for i=1:connComp.NumObjects
    leftx=borders((i-1)*4 1);
    lefty=borders((i-1)*4 2);
    width=borders((i-i)*4 3);
    height=borders((i-1)*4 4);
    if area(i)<1000 || area(i)>5000
        temp2(connComp.PixelIdxList{i})=0;
    elseif width/height>6 || width/height<3
        temp2(connComp.PixelIdxList{i})=0;
    end
end
figure,imshow(temp2);
%calculate plate coordinates
[x,y]=size(temp2);
xcounts=zeros(x,1);
for i=1:x
    for j=1:y
        if(temp2(i,j)==1)
            xcounts(i,1)=xcounts(i,1) 1;
        end
    end
end
[mv,maxx]=max(xcounts);
px1=maxx;
while((xcounts(px1,1)>=5)&&(px1>1))
    px1=px1-1;
end
px2=maxx;
while((xcounts(px2,1)>=5)&&(px2<x))
    px2=px2 1;
end
ycounts=zeros(1,y);
for j=1:y
    for i=px1:px2
        if(temp2(i,j)==1)
            ycounts(1,j)=ycounts(1,j) 1;
        end
    end
end
py1=1;
while((ycounts(1,py1)<3)&&(py1<y))
    py1=py1 1;
end
py2=y;
while((ycounts(1,py2)<3)&&(py2>py1))
py2=py2-1;
end
%car plate image to gray
temp3=car(px1:px2,py1:py2,:);
figure,imshow(temp3);
temp4=rgb2gray(temp3);
figure,imshow(temp4);
%car plate image binarize
g_max=double(max(max(temp4)));
g_min=double(min(min(temp4)));
T=round(g_max-(g_max-g_min)/2);
[m,n]=size(temp4);
temp5=double(temp4)>=T;
figure,imshow(temp5);title('temp5');
%median filter
h=fspecial('average',3);
temp6=im2bw(round(filter2(h,temp5)));
figure,imshow(temp6);
%Character diviSion
[m,n]=size(temp5);
char_width=[100];
char_position=[100];
%vertical diviSion
s1=sum(temp5);
content1=find(s1>=10 & s1<=n*0.8);
[m2,n2]=size(content1);
height1=content1(1);
height2=content1(n2);
%hotiZontai diviSion
s=sum(temp5(height1:height2,0:1));
j=1;
i=1;
flag=1;
content=find(s~=0);
[m1,n1]=size(content);
for i=2:n1
    if(content(i)-content(i-1)>1)
        flag=flag 1;
        char_position(flag)=content(i);
        char_width(flag-1)=content(i-1)-char_position(flag-1);
    end
end
char_position(1)=1;
char_width(1)=char_width(1)-1;
char_width(flag)=content(n1)-char_position(flag);
%write the chatacters IntO p1CS
word1=temp5(height1:height2,char_position(1):char_position(1) char_width(1));
figure,imshow(word1);
word2=temp5(height1:height2,char_position(2):char_position(2) char_width(2));
figure,imshow(word2);
word3=temp5(height1:height2,char_position(3):char_position(3) char_width(3));
figure,imshow(word3);
word4=temp5(height1:height2,char_position(4):char_position(4) char_width(4));
figure,imshow(word4);
word5=temp5(height1:height2,char_position(5):char_position(5) char_width(5));
figure,imshow(word5);
word6=temp5(height1:height2,char_position(6):char_position(6) char-width(6));
figure,imshow(word6);
word7=temp5(height1:height2,char_position(7):char_position(7) char-width(7));
figure,imshow(word7);
%character normialization
word1=imresize(word1,[32,16]);
word2=imresize(word2,[32,16]);
word3=imresize(word3,[32,16]);
word4=imresize(word4,[32,16]);
word5=imresize(word5,[32,16]);
word6=imresize(word6,[32,16]);
word7=imresize(word7,[32,16]);
imshow(word1);
imshow(word2);
imshow(word3);
imshow(word4);
imshow(word5);
imshow(word6);
imshow(word7);
imwrite(word1,'1.jpg');
imwrite(word2,'2.jpg');
imwrite(word3,'3.jpg');
imwrite(word4,'4.jpg');
imwrite(word5,'5.jpg');
imwrite(word6,'6.jpg');
imwrite(word7,'7.jpg');

%file name list
liccode=char(['0':'9' 'A':'Z' '京津沪渝冀晋辽吉黑苏浙皖闽赣鲁豫鄂湘粤琼川贵云陕甘青藏桂蒙新宁港']);

subbw=zeros(32,16);
l=1;
%compare with the standard models
for i=1:7
    ii=int2str(i);
    t=imread([ii,'.jpg']);
    seg=im2bw(t);
    if i==1
        kmin=37;
        kmax=68;
    elseif i==2
        kmin=11;
        kmax:36;
    elseif i>=3
        kmin=1;
        kmax=36;
    end

    for k2=kmin:kmax
        fname=strcat('',liccode(k2),'.bmp');
        t1=imread(fname);
        sam=im2bw(t1);
        for i=1:32
            for j=1:16
                subbw(i,j)=seg(i,j)-sam(i,j);
            end
        end
        dmax=0;
        for i=l:32
            for J=1:16
                if(subbw(i,j)~=0)
                end
            end
            dmax=dmax 1;
        end
        error(k2)=dmax;
    end
%find the minimum number of 1s
    Error1=error(kmin:kmax);
    MinError=min(Error1);
    findc=find(Error1==MinError);
    Code(1*2-1)=liccode(findc(1) kmin-1);
    Code(1*2)=' ';
    l=l 1;
end