車牌識(shí)別-MATLAB代碼(共10頁(yè))

1、clear ;close all;%Step1 獲取圖像 裝入待處理彩色圖像并顯示(xinsh)原始圖像Scolor = imread(C:UserslenovoDesktop程序(chngx)chepail2.jpg);%將彩色圖像轉(zhuǎn)換為黑白(hibi)并顯示Sgray = rgb2gray(Scolor);figure,imshow(Scolor),title(原始彩色圖像);figure,imshow(Sgray),title(原始灰度圖像);%bw3=medfilt2(Sgray,5 5);%figure,imshow(bw3);title(圖像中值濾波);%bw4=imfilter(

2、bw3,fspecial(laplacian);%figure,imshow(bw4);title(圖像中值濾波);%Egray=imsubtract(bw3,bw4);%figure,imshow(Egray);%Step2 圖像預(yù)處理 對(duì)Sgray 原始黑白圖像進(jìn)行開操作得到圖像背景s=strel(disk,10);Bgray=imopen(Sgray,s);figure,imshow(Bgray);title(背景圖像);%用原始圖像與背景圖像作減法，增強(qiáng)圖像Egray=imsubtract(Sgray,Bgray);figure,imshow(Egray);title(增強(qiáng)黑白圖像);

3、%Step3 取得最佳閾值，將圖像二值化%fmax1=double(max(max(Egray);%fmin1=double(min(min(Egray);%level=(fmax1-(fmax1-fmin1)/3)/255;%level=graythresh(Egray);%bw3=filter2(fspecial(average,3),Egray);%Egray=medfilt2(Egray,3 3);%figure,imshow(Egray);title(圖像中值濾波);bw2=im2bw(Egray,level);figure,imshow(bw2);title(圖像二值化);%bw4

4、=im2bw(bw3,level);%figure,imshow(bw4);title(圖像二值化);%Step4 對(duì)得到二值圖像作開閉操作進(jìn)行(jnxng)濾波grd=edge(bw2,canny);figure,imshow(grd);title(圖像(t xin)邊緣提取);%bg1=imclose(grd,strel(rectangle,5,19);%figure,imshow(bg1);title(圖像(t xin)閉運(yùn)算5,19);bg2=bwareaopen(bg1,2000);figure,imshow(bg2);title(小對(duì)象);bg2=imclose(bg2,strel

5、(rectangle,20,50);%bg3=imclose(bg2,strel(rectangle,20,50);bg2=imopen(bg2,strel(rectangle,5,20);figure,imshow(bg2);title(圖像開運(yùn)算5,19);bg2=imopen(bg2,strel(rectangle,19,1);figure,imshow(bg2);title(圖像開運(yùn)算19,1);%Step5 對(duì)二值圖像進(jìn)行區(qū)域提取，并計(jì)算區(qū)域特征參數(shù)。進(jìn)行區(qū)域特征參數(shù)比較，提取車牌區(qū)域L,num =bwlabel(bg2,8);Feastats = regionprops(L,bas

6、ic);Area=Feastats.Area;BoundingBox=Feastats.BoundingBox; %x y width heightRGB = label2rgb(L, spring, k, shuffle); figure,imshow(RGB);title(圖像彩色標(biāo)記);lx=0;%Getok=zeros(1,5);for l=1:num width=BoundingBox(l-1)*4+3); hight=BoundingBox(l-1)*4+4); if (width250 & width70 & hight2 & rato4 break; endendsbw1=bw

7、2(startrow:startrow+hight,startcol:startcol+width-1); %獲取(huq)車牌二值子圖subcol1=Sgray(startrow:startrow+hight,startcol:startcol+width-1); %獲取(huq)車牌灰度子圖figure,subplot(2,1,1),imshow(subcol1);title(車牌灰度子圖);subplot(2,1,2),imshow(sbw1);title(車牌二值子圖);%Step6 計(jì)算車牌水平投影，并對(duì)水平投影進(jìn)行峰谷分析histcol1=sum(sbw1); %計(jì)算垂直投影his

8、trow=sum(sbw1); %計(jì)算水平投影figure,subplot(2,1,1),bar(histcol1);title(垂直投影（含邊框）);subplot(2,1,2),bar(histrow); title(水平投影（含邊框）);figure,subplot(2,1,1),bar(histrow); title(水平投影（含邊框）);subplot(2,1,2),imshow(sbw1);title(車牌二值子圖);%對(duì)水平投影進(jìn)行峰谷分析meanrow=mean(histrow);minrow=min(histrow);levelrow=(meanrow+minrow)/2;c

9、ount1=0;l=1;%global markrow;for k=1:hight if histrow(k)=1 markrow(l)=k; %上升點(diǎn) markrow1(l)=count1; %谷寬度（下降點(diǎn)至下一個(gè)上升點(diǎn)） l=l+1; end count1=0; endendmarkrow2=diff(markrow); %峰距離（上升(shngshng)點(diǎn)至下一個(gè)上升點(diǎn)）m1,n1=size(markrow2);n1=n1+1;markrow(l)=hight;markrow1(l)=count1;markrow2(n1)=markrow(l)-markrow(l-1);l=0;for

10、 k=1:n1 markrow3(k)=markrow(k+1)-markrow1(k+1); %下降(xijing)點(diǎn) markrow4(k)=markrow3(k)-markrow(k); %峰寬度(kund)（上升點(diǎn)至下降點(diǎn)） markrow5(k)=markrow3(k)-double(uint16(markrow4(k)/2); %峰中心位置end %Step7 計(jì)算車牌旋轉(zhuǎn)角度%(1)在上升點(diǎn)至下降點(diǎn)找第一個(gè)為1的點(diǎn)m2,n2=size(sbw1);m1,n1=size(markrow4);maxw=max(markrow4); %最大寬度為字符if markrow4(1) = m

11、axw %檢測(cè)上邊 ysite=1; k1=1; for l=1:n2 for k=1:markrow3(ysite) %從頂邊至第一個(gè)峰下降點(diǎn)掃描 if sbw1(k,l)=1 xdata(k1)=l; ydata(k1)=k; k1=k1+1; break; end end endelse %檢測(cè)下邊 ysite=n1; if markrow4(n1) =0 if markrow4(n1-1) =maxw ysite= 0; %無(wú)下邊 else ysite= n1-1; end end if ysite =0 k1=1; for l=1:n2 k=m2; while k=markrow(y

12、site) %從底邊(d bin)至最后一個(gè)峰的上升點(diǎn)掃描 if sbw1(k,l)=1 xdata(k1)=l; ydata(k1)=k; k1=k1+1; break; end k=k-1; end end endend %(2)線性擬合，計(jì)算(j sun)與x夾角fresult = fit(xdata,ydata,poly1); %poly1 Y = p1*x+p2p1=fresult.p1;angle=atan(fresult.p1)*180/pi; %弧度(hd)換為度，360/2pi, pi=3.14%(3)旋轉(zhuǎn)車牌圖象subcol = imrotate(subcol1,angle

13、,bilinear,crop); %旋轉(zhuǎn)車牌圖象sbw = imrotate(sbw1,angle,bilinear,crop);figure,subplot(2,1,1),imshow(subcol);title(車牌灰度子圖);subplot(2,1,2),imshow(sbw);title();title(車牌旋轉(zhuǎn)角: ,num2str(angle),度 ,Color,r);%Step8 旋轉(zhuǎn)車牌后重新計(jì)算車牌水平投影，去掉車牌水平邊框，獲取字符高度histcol1=sum(sbw); %計(jì)算垂直投影histrow=sum(sbw); %計(jì)算水平投影figure,subplot(2,1,

14、1),bar(histcol1);title(垂直投影（旋轉(zhuǎn)后）);subplot(2,1,2),bar(histrow); title(水平投影（旋轉(zhuǎn)后）);figure,subplot(2,1,1),bar(histrow); title(水平投影（旋轉(zhuǎn)后）);subplot(2,1,2),imshow(sbw);title(車牌二值子圖（旋轉(zhuǎn)后）);%對(duì)水平投影進(jìn)行峰谷分析clear markrow,markrow1,markrow2,markrow3,markrow4,markrow5meanrow=mean(histrow);minrow=min(histrow);levelrow=

15、(meanrow+minrow)/2;count1=0;l=1;for k=1:hight if histrow(k)=1 markrow(l)=k; %上升(shngshng)點(diǎn) markrow1(l)=count1; %谷寬度（下降(xijing)點(diǎn)至下一個(gè)上升點(diǎn)） l=l+1; end count1=0; endendmarkrow2=diff(markrow); %峰距離(jl)（上升點(diǎn)至下一個(gè)上升點(diǎn)）m1,n1=size(markrow2);n1=n1+1;markrow(l)=hight;markrow1(l)=count1;markrow2(n1)=markrow(l)-mark

16、row(l-1);l=0;for k=1:n1 markrow3(k)=markrow(k+1)-markrow1(k+1); %下降點(diǎn) markrow4(k)=markrow3(k)-markrow(k); %峰寬度（上升點(diǎn)至下降點(diǎn)） markrow5(k)=markrow3(k)-double(uint16(markrow4(k)/2); %峰中心位置end %去水平（上下）邊框,獲取字符高度maxhight=max(markrow2);findc=find(markrow2=maxhight);rowtop=markrow(findc);rowbot=markrow(findc+1)-m

17、arkrow1(findc+1);sbw2=sbw(rowtop:rowbot,:); %子圖為(rowbot-rowtop+1)行maxhight=rowbot-rowtop+1; %字符高度(rowbot-rowtop+1)%Step9 計(jì)算車牌垂直投影，去掉車牌垂直邊框，獲取車牌及字符平均寬度histcol=sum(sbw2); %計(jì)算垂直投影figure,subplot(2,1,1),bar(histcol);title(垂直投影（去水平邊框后）);subplot(2,1,2),imshow(sbw2); title(車牌字符(z f)高度： ,int2str(maxhight),Co

18、lor,r);%對(duì)垂直投影進(jìn)行(jnxng)峰谷分析meancol=mean(histcol);mincol=min(histcol);levelcol=(meancol+mincol)/4;count1=0;l=1;for k=1:width if histcol(k)=1 markcol(l)=k; %字符(z f)上升點(diǎn) markcol1(l)=count1; %谷寬度（下降點(diǎn)至下一個(gè)上升點(diǎn)） l=l+1; end count1=0; endendmarkcol2=diff(markcol); %字符距離（上升點(diǎn)至下一個(gè)上升點(diǎn)）m1,n1=size(markcol2);n1=n1+1;m

19、arkcol(l)=width;markcol1(l)=count1;markcol2(n1)=markcol(l)-markcol(l-1);%Step10 計(jì)算車牌上每個(gè)字符中心位置，計(jì)算最大字符寬度maxwidthl=0;for k=1:n1 markcol3(k)=markcol(k+1)-markcol1(k+1); %字符下降點(diǎn) markcol4(k)=markcol3(k)-markcol(k); %字符寬度（上升點(diǎn)至下降點(diǎn)） markcol5(k)=markcol3(k)-double(uint16(markcol4(k)/2); %字符中心位置end markcol6=dif

20、f(markcol5); %字符中心距離（字符中心點(diǎn)至下一個(gè)字符中心點(diǎn)）maxs=max(markcol6); %查找最大值，即為第二字符與第三字符中心距離findmax=find(markcol6=maxs);markcol6(findmax)=0;maxwidth=max(markcol6); %查找(ch zho)最大值，即為最大字符寬度%Step11 提取分割字符(z f),并變換為40行*20列標(biāo)準(zhǔn)子圖l=1;m2,n2=size(sbw2);%subcol換為sbw2figure;for k=findmax-1:findmax+5 cleft=markcol5(k)-maxwidt

21、h/2; cright=markcol5(k)+maxwidth/2-2; if cleftn2 cright=n2; cleft=n2-maxwidth; end SegGray=subcol(rowtop:rowbot,cleft+2:cright); SegBw1=sbw(rowtop:rowbot,cleft+2:cright); SegBw2 = imresize(SegBw1,40 20); %變換(binhun)為40行*20列標(biāo)準(zhǔn)子圖 subplot(2,n1,l),imshow(SegGray); if l=7 title(車牌字符寬度： ,int2str(maxwidth)

22、,Color,r); end subplot(2,n1,n1+l),imshow(SegBw2); %fname=strcat(G:MatlablibS,Filename,int2str(k),.bmp);%保存子圖備選入樣本庫(kù),并建立樣本庫(kù) fname=strcat(int2str(k),.BMP); imwrite(SegBw2,fname) ; l=l+1;end%Step12 將計(jì)算計(jì)算獲取的字符圖像與樣本庫(kù)進(jìn)行匹配，自動(dòng)識(shí)別出字符代碼。%在進(jìn)行匹配前需將存入的子圖名按以下編碼規(guī)則更改：%Sam0 Sam9 09 數(shù)字 110%SamASamZ AZ 字母 1136%Sam粵Sam滬

23、省別 漢字 3745liccode=char(0:9 A:Z 蘇豫陜魯); %建立自動(dòng)識(shí)別字符代碼表 l=1;m2,n2=size(sbw2);for k=findmax-1:findmax+5 cleft=markcol5(k)-maxwidth/2; cright=markcol5(k)+maxwidth/2-2; if cleftn2 cright=n2; cleft=n2-maxwidth; end SegBw1=sbw(rowtop:rowbot,cleft+2:cright); SegBw2 = imresize(SegBw1,40 20 ,bilinear); %變換(binhun)為32行*16列標(biāo)準(zhǔn)子圖 if l