【正文】 imwrite(b,39。39。)。figure(8)。subplot(3,2,1),imshow(b),title(39。39。)g_max=double(max(max(b)))。g_min=double(min(min(b)))。T=round(g_max(g_maxg_min)/3)。 % T 為二值化的閾值[m,n]=size(b)。d=(double(b)=T)。 % d:二值圖像imwrite(d,39。39。)。figure(8)。subplot(3,2,2),imshow(d),title(39。39。)figure(8),subplot(3,2,3),imshow(d),title(39。39。)% 濾波h=fspecial(39。average39。,3)。d=im2bw(round(filter2(h,d)))。imwrite(d,39。39。)。figure(8),subplot(3,2,4),imshow(d),title(39。39。)% 某些圖像進行操作% 膨脹或腐蝕% se=strel(39。square39。,3)。 % 使用一個3X3的正方形結(jié)果元素對象對創(chuàng)建的圖像膨脹% 39。line39。/39。diamond39。/39。ball39。...se=eye(2)。 % eye(n) returns the nbyn identity matrix 單位矩陣[m,n]=size(d)。if bwarea(d)/m/n= d=imerode(d,se)。elseif bwarea(d)/m/n= d=imdilate(d,se)。endimwrite(d,39。39。)。figure(8),subplot(3,2,5),imshow(d),title(39。39。)% 尋找連續(xù)有文字的塊，若長度大于某閾值，則認為該塊有兩個字符組成，需要分割d=qiege(d)。[m,n]=size(d)。figure,subplot(2,1,1),imshow(d),title(n)k1=1。k2=1。s=sum(d)。j=1。while j~=n while s(j)==0 j=j+1。 end k1=j。 while s(j)~=0 amp。amp。 j=n1 j=j+1。 end k2=j1。 if k2k1=round(n/) [val,num]=min(sum(d(:,[k1+5:k25])))。 d(:,k1+num+5)=0。 % 分割 endend% 再切割d=qiege(d)。% 切割出 7 個字符y1=10。y2=。flag=0。word1=[]。while flag==0 [m,n]=size(d)。 left=1。wide=0。 while sum(d(:,wide+1))~=0 wide=wide+1。 end if widey1 % 認為是左側(cè)干擾 d(:,[1:wide])=0。 d=qiege(d)。 else temp=qiege(imcrop(d,[1 1 wide m]))。 [m,n]=size(temp)。 all=sum(sum(temp))。 two_thirds=sum(sum(temp([round(m/3):2*round(m/3)],:)))。 if two_thirds/ally2 flag=1。word1=temp。 % WORD 1 end d(:,[1:wide])=0。d=qiege(d)。 endend% 分割出第二個字符[word2,d]=getword(d)。% 分割出第三個字符[word3,d]=getword(d)。% 分割出第四個字符[word4,d]=getword(d)。% 分割出第五個字符[word5,d]=getword(d)。% 分割出第六個字符[word6,d]=getword(d)。% 分割出第七個字符[word7,d]=getword(d)。figure(9),imshow(word1),title(39。139。)。figure(10),imshow(word2),title(39。239。)。figure(11),imshow(word3),title(39。339。)。figure(12),imshow(word4),title(39。439。)。figure(13),imshow(word5),title(39。539。)。figure(14),imshow(word6),title(39。639。)。figure(15),imshow(word7),title(39。739。)。[m,n]=size(word1)。% 商用系統(tǒng)程序中歸一化大小為 40*20,此處演示word1=imresize(word1,[40 20])。word2=imresize(word2,[40 20])。word3=imresize(word3,[40 20])。word4=imresize(word4,[40 20])。word5=imresize(word5,[40 20])。word6=imresize(word6,[40 20])。word7=imresize(word7,[40 20])。figure(16),subplot(3,7,8),imshow(word1),title(39。139。)。subplot(3,7,9),imshow(word2),title(39。239。)。subplot(3,7,10),imshow(word3),title(39。339。)。subplot(3,7,11),imshow(word4),title(39。439。)。subplot(3,7,12),imshow(word5),title(39。539。)。subplot(3,7,13),imshow(word6),title(39。639。)。subplot(3,7,14),imshow(word7),title(39。739。)。imwrite(word1,39。39。)。imwrite(word2,39。39。)。imwrite(word3,39。39。)。imwrite(word4,39。39。)。imwrite(word5,39。39。)。imwrite(word6,39。39。)。imwrite(word7,39。39。)。liccode=char([39。039。:39。939。 39。A39。:39。Z39。 39。蘇豫陜魯39。])。 %建立自動識別字符代碼表 SubBw2=zeros(40,20)。l=1。for I=1:7 ii=int2str(I)。 t=imread([ii,39。.jpg39。])。 SegBw2=imresize(t,[40 20],39。nearest39。)。 if l==1 %第一位漢字識別 kmin=37。 kmax=40。 elseif l==2 %第二位 A~Z 字母識別 kmin=11。 kmax=36。 else l=3 %第三位以后是字母或數(shù)字識別 kmin=1。 kmax=36。 end for k2=kmin:kmax fname=strcat(39。字符模板\39。,liccode(k2),39。.jpg39。)。 SamBw2 = imread(fname)。 for i=1:40 for j=1:20 SubBw2(i,j)=SegBw2(i,j)SamBw2(i,j)。 end end % 以上相當于兩幅圖相減得到第三幅圖 Dmax=0。 for k1=1:40 for l1=1:20 if ( SubBw2(k1,l1) 0 | SubBw2(k1,l1) 0 ) Dmax=Dmax+1。 end end end Error(k2)=Dmax。 end Error1=Error(kmin:kmax)。 MinError=min(Error1)。 findc=find(Error1==MinError)。 l=l+1。endfigure(17),subplot(3,1,2),imshow(dw),title (39。車牌號碼: 陜 A B A 2 2 339。)。%子程序：（getword子程序）function [word,result]=getword(d)word=[]。flag=0。y1=8。y2=。 while flag==0 [m,n]=size(d)。 wide=0。 while sum(d(:,wide+1))~=0 amp。amp。 wide=n2 wide=wide+1。 end temp=qiege(imcrop(d,[1 1 wide m]))。 [m1,n1]=size(temp)。 if widey1 amp。amp。 n1/m1y2 d(:,[1:wide])=0。 if sum(sum(d))~=0 d=qiege(d)。 % 切割出最小范圍 else word=[]。flag=1。 end else word=qiege(imcrop(d,[1 1 wide m]))。 d(:,[1:wide])=0。 if sum(sum(d))~=0。 d=qiege(d)。flag=1。 else d=[]。 end end end result=d。 % （qiege子程序）function e=qiege(d)[m,n]=size(d)。top=1。bottom=m。left=1。right=n。 % initwhile sum(d(top,:))==0 amp。amp。 top=m top=top+1。endwhile sum(d(bottom,:))==0 amp。amp。 bottom1 bottom=bottom1。endwhile sum(d(:,left))==0 amp。amp。 leftn left=left+1。endwhile sum(d(:,right))==0 amp。amp。 right=1right=right1。enddd=rightleft。hh=bottomtop。e=imcrop(d,[left top dd hh])