車(chē)牌識(shí)別地matlab程序

1、精選優(yōu)質(zhì)文檔-傾情為你奉上附錄車(chē)牌識(shí)別程序clear ;close all;%Step1 獲取圖像 裝入待處理彩色圖像并顯示原始圖像Scolor = imread('3.jpg');%imread函數(shù)讀取圖像文件%將彩色圖像轉(zhuǎn)換為黑白并顯示Sgray = rgb2gray(Scolor);%rgb2gray轉(zhuǎn)換成灰度圖figure,imshow(Scolor),title('原始彩色圖像');%figure命令同時(shí)顯示兩幅圖figure,imshow(Sgray),title('原始黑白圖像');%Step2 圖像預(yù)處理 對(duì)Sgray 原始黑白圖

2、像進(jìn)行開(kāi)操作得到圖像背景s=strel('disk',13);%strel函數(shù)Bgray=imopen(Sgray,s);%打開(kāi)sgray s圖像figure,imshow(Bgray);title('背景圖像');%輸出背景圖像%用原始圖像與背景圖像作減法，增強(qiáng)圖像Egray=imsubtract(Sgray,Bgray);%兩幅圖相減figure,imshow(Egray);title('增強(qiáng)黑白圖像');%輸出黑白圖像%Step3 取得最佳閾值，將圖像二值化fmax1=double(max(max(Egray);%egray的最大值并輸出雙

3、精度型fmin1=double(min(min(Egray);%egray的最小值并輸出雙精度型level=(fmax1-(fmax1-fmin1)/3)/255;%獲得最佳閾值bw22=im2bw(Egray,level);%轉(zhuǎn)換圖像為二進(jìn)制圖像bw2=double(bw22);%Step4 對(duì)得到二值圖像作開(kāi)閉操作進(jìn)行濾波figure,imshow(bw2);title('圖像二值化');%得到二值圖像grd=edge(bw2,'canny')%用canny算子識(shí)別強(qiáng)度圖像中的邊界figure,imshow(grd);title('圖像邊緣提取

4、9;);%輸出圖像邊緣bg1=imclose(grd,strel('rectangle',5,19);%取矩形框的閉運(yùn)算figure,imshow(bg1);title('圖像閉運(yùn)算5,19');%輸出閉運(yùn)算的圖像bg3=imopen(bg1,strel('rectangle',5,19);%取矩形框的開(kāi)運(yùn)算figure,imshow(bg3);title('圖像開(kāi)運(yùn)算5,19');%輸出開(kāi)運(yùn)算的圖像bg2=imopen(bg3,strel('rectangle',19,1);%取矩形框的開(kāi)運(yùn)算figure,ims

5、how(bg2);title('圖像開(kāi)運(yùn)算19,1');%輸出開(kāi)運(yùn)算的圖像%Step5 對(duì)二值圖像進(jìn)行區(qū)域提取，并計(jì)算區(qū)域特征參數(shù)。進(jìn)行區(qū)域特征參數(shù)比較，提取車(chē)牌區(qū)域L,num = bwlabel(bg2,8);%標(biāo)注二進(jìn)制圖像中已連接的部分Feastats = imfeature(L,'basic');%計(jì)算圖像區(qū)域的特征尺寸Area=Feastats.Area;%區(qū)域面積BoundingBox=Feastats.BoundingBox;%x y width height車(chē)牌框架大小RGB = label2rgb(L, 'spring', &#

6、39;k', 'shuffle'); %標(biāo)志圖像向RGB圖像轉(zhuǎn)換figure,imshow(RGB);title('圖像彩色標(biāo)記');%輸出框架的彩色圖像lx=0;for l=1:num width=BoundingBox(l-1)*4+3);%框架寬度的計(jì)算 hight=BoundingBox(l-1)*4+4);%框架高度的計(jì)算 if (width>98 & width<160 & hight>25 & hight<50)%框架的寬度和高度的范圍 lx=lx+1; Getok(lx)=l; endend

7、for k= 1:lx l=Getok(k); startcol=BoundingBox(l-1)*4+1)-2;%開(kāi)始列 startrow=BoundingBox(l-1)*4+2)-2;%開(kāi)始行 width=BoundingBox(l-1)*4+3)+8;%車(chē)牌寬 hight=BoundingBox(l-1)*4+4)+2;%車(chē)牌高 rato=width/hight;%計(jì)算車(chē)牌長(zhǎng)寬比 if rato>2 & rato<4 break; endendsbw1=bw2(startrow:startrow+hight,startcol:startcol+width-1); %

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

9、ure,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('車(chē)牌二值子圖');%輸出二值圖%對(duì)水平投影進(jìn)行峰谷分析meanrow=mean(histrow);%求水平投影的平

10、均值minrow=min(histrow);%求水平投影的最小值levelrow=(meanrow+minrow)/2;%求水平投影的平均值count1=0;l=1;for k=1:hight if histrow(k)<=levelrow count1=count1+1; else if count1>=1 markrow(l)=k;%上升點(diǎn) markrow1(l)=count1;%谷寬度(下降點(diǎn)至下一個(gè)上升點(diǎn)) l=l+1; end count1=0; endendmarkrow2=diff(markrow);%峰距離(上升點(diǎn)至下一個(gè)上升點(diǎn))m1,n1=size(markrow

11、2);n1=n1+1;markrow(l)=hight;markrow1(l)=count1;markrow2(n1)=markrow(l)-markrow(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 %Step7 計(jì)算車(chē)牌旋轉(zhuǎn)角度%(1)在上升點(diǎn)至下降點(diǎn)找第一個(gè)為1的點(diǎn)m2,n2=size

12、(sbw1);%sbw1的圖像大小m1,n1=size(markrow4);%markrow4的大小maxw=max(markrow4);%最大寬度為字符if markrow4(1) = maxw%檢測(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

13、ysite= 0; %無(wú)下邊 else ysite= n1-1; end end if ysite =0 k1=1; for l=1:n2 k=m2; while k>=markrow(ysite) %從底邊至最后一個(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ì)算與x夾角fresult = fit(xdata',ydata','poly1'); %poly1 Y = p1*x+p2p1=fresult

14、.p1;angle=atan(fresult.p1)*180/pi; %弧度換為度，360/2pi, pi=3.14%(3)旋轉(zhuǎn)車(chē)牌圖象subcol = imrotate(subcol1,angle,'bilinear','crop'); %旋轉(zhuǎn)車(chē)牌圖象sbw = imrotate(sbw1,angle,'bilinear','crop');%旋轉(zhuǎn)圖像figure,subplot(2,1,1),imshow(subcol);title('車(chē)牌灰度子圖');%輸出車(chē)牌旋轉(zhuǎn)后的灰度圖像標(biāo)題顯示車(chē)牌灰度子圖subplot

15、(2,1,2),imshow(sbw);title('');%輸出車(chē)牌旋轉(zhuǎn)后的灰度圖像title('車(chē)牌旋轉(zhuǎn)角: ',num2str(angle),'度' ,'Color','r');%顯示車(chē)牌的旋轉(zhuǎn)角度%Step8 旋轉(zhuǎn)車(chē)牌后重新計(jì)算車(chē)牌水平投影，去掉車(chē)牌水平邊框，獲取字符高度histcol1=sum(sbw); %計(jì)算垂直投影histrow=sum(sbw'); %計(jì)算水平投影figure,subplot(2,1,1),bar(histcol1);title('垂直投影(旋轉(zhuǎn)后)');s

16、ubplot(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('車(chē)牌二值子圖(旋轉(zhuǎn)后)');%去水平(上下)邊框,獲取字符高度maxhight=max(markrow2);findc=find(markrow2=maxhight);rowtop=markrow(findc);rowbot=markrow(findc+1)-markrow1(f

17、indc+1);sbw2=sbw(rowtop:rowbot,:); %子圖為(rowbot-rowtop+1)行maxhight=rowbot-rowtop+1; %字符高度(rowbot-rowtop+1)%Step9 計(jì)算車(chē)牌垂直投影，去掉車(chē)牌垂直邊框，獲取車(chē)牌及字符平均寬度histcol=sum(sbw2); %計(jì)算垂直投影figure,subplot(2,1,1),bar(histcol);title('垂直投影(去水平邊框后)');%輸出車(chē)牌的垂直投影圖像subplot(2,1,2),imshow(sbw2); %輸出垂直投影圖像title('車(chē)牌字符高度：

18、 ',int2str(maxhight),'Color','r');%輸出車(chē)牌字符高度%對(duì)垂直投影進(jìn)行峰谷分析meancol=mean(histcol);%求垂直投影的平均值mincol=min(histcol);%求垂直投影的平均值levelcol=(meancol+mincol)/4;%求垂直投影的1/4count1=0;l=1;for k=1:width if histcol(k)<=levelcol count1=count1+1; else if count1>=1 markcol(l)=k; %字符上升點(diǎn) markcol1(l)=

19、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;markcol(l)=width;markcol1(l)=count1;markcol2(n1)=markcol(l)-markcol(l-1);%Step10 計(jì)算車(chē)牌上每個(gè)字符中心位置，計(jì)算最大字符寬度maxwidthl=0;for k=1:n1 markcol3(k)=markcol(k+1)-markcol1(k+1);%字符下降點(diǎn) markcol

20、4(k)=markcol3(k)-markcol(k); %字符寬度(上升點(diǎn)至下降點(diǎn)) markcol5(k)=markcol3(k)-double(uint16(markcol4(k)/2);%字符中心位置end markcol6=diff(markcol5); %字符中心距離(字符中心點(diǎn)至下一個(gè)字符中心點(diǎn))maxs=max(markcol6); %查找最大值，即為第二字符與第三字符中心距離findmax=find(markcol6=maxs);markcol6(findmax)=0;maxwidth=max(markcol6);%查找最大值，即為最大字符寬度%Step11 提取分割字符,并

21、變換為22行´14列標(biāo)準(zhǔn)子圖l=1;m2,n2=size(subcol);figure;for k=findmax-1:findmax+5 cleft=markcol5(k)-maxwidth/2; cright=markcol5(k)+maxwidth/2-2; if cleft<1 cleft=1; cright=maxwidth; end if cright>n2 cright=n2; cleft=n2-maxwidth; end SegGray=sbw(rowtop:rowbot,cleft:cright); SegBw1=sbw(rowtop:rowbot,cl

22、eft:cright); SegBw2 = imresize(SegBw1,22 14);%變換為22行´14列標(biāo)準(zhǔn)子圖 subplot(2,n1,l),imshow(SegGray); if l=7 title('車(chē)牌字符寬度： ',int2str(maxwidth),'Color','r'); end subplot(2,n1,n1+l),imshow(SegBw2); fname=strcat('c:worksamimage',int2str(k),'.jpg'); imwrite(SegBw2,f

23、name,'jpg') l=l+1;end%Step12 將計(jì)算計(jì)算獲取的字符圖像與樣本庫(kù)進(jìn)行匹配，自動(dòng)識(shí)別出字符代碼。liccode=char('0':'9' 'A':'Z' '粵桂海云貴川京津滬'); %建立自動(dòng)識(shí)別字符代碼表 SubBw2=zeros(22,14);l=1;m2,n2=size(sbw);for k=findmax-1:findmax+5 cleft=markcol5(k)-maxwidth/2; cright=markcol5(k)+maxwidth/2-2; if cleft<1 cleft=1; cright=maxwidth; end if cright>n2 cright=n2; cleft=n2-maxwidth; end SegBw1=sbw(rowtop:rowbot,cleft:cright); SegBw2 = imresize(SegBw1,22 14