車牌識別的matlab程序

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

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

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

4、9;）;%輸出圖像邊緣bg1=imclose（grd,strel（'rectangle',5,19）;%取矩形框的閉運算figure,imshow（bg1）;title（'圖像閉運算5,19'）;%輸出閉運算的圖像bg3=imopen（bg1,strel（'rectangle',5,19）;%取矩形框的開運算figure,imshow（bg3）;title（'圖像開運算5,19'）;%輸出開運算的圖像bg2=imopen（bg3,strel（'rectangle',19,1）;%取矩形框的開運算figure,ims

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

6、39;k', 'shuffle'）; %標志圖像向RGB圖像轉換figure,imshow（RGB）;title（'圖像彩色標記'）;%輸出框架的彩色圖像lx=0;for l=1:num width=BoundingBox（l-1）*4+3）;%框架寬度的計算 hight=BoundingBox（l-1）*4+4）;%框架高度的計算 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;%開始列 startrow=BoundingBox（l-1）*4+2）-2;%開始行 width=BoundingBox（l-1）*4+3）+8;%車牌寬 hight=BoundingBox（l-1）*4+4）+2;%車牌高 rato=width/hight;%計算車牌長寬比 if rato>2 & rato<4 break; endendsbw1=bw2（startrow:startrow+hight,startcol:startcol+width-1）; %

8、獲取車牌二值子圖subcol1=Sgray（startrow:startrow+hight,startcol:startcol+width-1）;%獲取車牌灰度子圖figure,subplot（2,1,1）,imshow（subcol1）;title（'車牌灰度子圖'）;%輸出灰度圖subplot（2,1,2）,imshow（sbw1）;title（'車牌二值子圖'）;%輸出車牌的二值圖%Step6 計算車牌水平投影，并對水平投影進行峰谷分析histcol1=sum（sbw1）; %計算垂直投影histrow=sum（sbw1'）; %計算水平投影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（'車牌二值子圖'）;%輸出二值圖%對水平投影進行峰谷分析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;%上升點 markrow1（l）=count1;%谷寬度（下降點至下一個上升點） l=l+1; end count1=0; endendmarkrow2=diff（markrow）;%峰距離（上升點至下一個上升點）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）;%下降點 markrow4（k）=markrow3（k）-markrow（k）;%峰寬度（上升點至下降點） markrow5（k）=markrow3（k）-double（uint16（markrow4（k）/2）;%峰中心位置end %Step7 計算車牌旋轉角度%（1）在上升點至下降點找第一個為1的點m2,n2=size

12、（sbw1）;%sbw1的圖像大小m1,n1=size（markrow4）;%markrow4的大小maxw=max（markrow4）;%最大寬度為字符if markrow4（1） = maxw%檢測上邊 ysite=1; k1=1; for l=1:n2 for k=1:markrow3（ysite）%從頂邊至第一個峰下降點掃描 if sbw1（k,l）=1 xdata（k1）=l; ydata（k1）=k; k1=k1+1; break; end end endelse %檢測下邊 ysite=n1; if markrow4（n1） =0 if markrow4（n1-1） =maxw

13、ysite= 0; %無下邊 else ysite= n1-1; end end if ysite =0 k1=1; for l=1:n2 k=m2; while k>=markrow（ysite） %從底邊至最后一個峰的上升點掃描 if sbw1（k,l）=1 xdata（k1）=l; ydata（k1）=k; k1=k1+1; break; end k=k-1; end end endend %（2）線性擬合，計算與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）旋轉車牌圖象subcol = imrotate（subcol1,angle,'bilinear','crop'）; %旋轉車牌圖象sbw = imrotate（sbw1,angle,'bilinear','crop'）;%旋轉圖像figure,subplot（2,1,1）,imshow（subcol）;title（'車牌灰度子圖'）;%輸出車牌旋轉后的灰度圖像標題顯示車牌灰度子圖subplot

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

16、ubplot（2,1,2）,bar（histrow）; title（'水平投影（旋轉后）'）;figure,subplot（2,1,1）,bar（histrow）; title（'水平投影（旋轉后）'）;subplot（2,1,2）,imshow（sbw）;title（'車牌二值子圖（旋轉后）'）;%去水平（上下）邊框,獲取字符高度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 計算車牌垂直投影，去掉車牌垂直邊框，獲取車牌及字符平均寬度histcol=sum（sbw2）; %計算垂直投影figure,subplot（2,1,1）,bar（histcol）;title（'垂直投影（去水平邊框后）'）;%輸出車牌的垂直投影圖像subplot（2,1,2）,imshow（sbw2）; %輸出垂直投影圖像title（'車牌字符高度：

18、 ',int2str（maxhight）,'Color','r'）;%輸出車牌字符高度%對垂直投影進行峰谷分析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; %字符上升點 markcol1（l）=

19、count1; %谷寬度（下降點至下一個上升點） l=l+1; end count1=0; endendmarkcol2=diff（markcol）;%字符距離（上升點至下一個上升點）m1,n1=size（markcol2）;n1=n1+1;markcol（l）=width;markcol1（l）=count1;markcol2（n1）=markcol（l）-markcol（l-1）;%Step10 計算車牌上每個字符中心位置，計算最大字符寬度maxwidthl=0;for k=1:n1 markcol3（k）=markcol（k+1）-markcol1（k+1）;%字符下降點 markcol

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

21、變換為22行´14列標準子圖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列標準子圖 subplot（2,n1,l）,imshow（SegGray）; if l=7 title（'車牌字符寬度： ',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 將計算計算獲取的字符圖像與樣本庫進行匹配，自動識別出字符代碼。liccode=char（'0':'9' 'A':'Z' '粵桂海云貴川京津滬'）; %建立自動識別字符代碼表 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