matlab车牌识别代码 出现错误怎么弄？

1.                                 
   
2. function [d]=main(jpg)
   
3. close all
   
4. clc
   
5. I=imread('car1.jpg');
   
6. figure(1),imshow(I);title('原图')
   
7. I1=rgb2gray(I);
   
8. figure(2),subplot(1,2,1),imshow(I1);title('灰度图');
   
9. figure(2),subplot(1,2,2),imhist(I1);title('灰度图直方图');
   
10. I2=edge(I1,'roberts',0.15,'both');
    
11. figure(3),imshow(I2);title('robert算子边缘检测')
    
12. se=[1;1;1];
    
13. I3=imerode(I2,se);
    
14. figure(4),imshow(I3);title('腐蚀后图像');
    
15. se=strel('rectangle',[25,25]);
    
16. I4=imclose(I3,se);
    
17. figure(5),imshow(I4);title('平滑图像的轮廓');
    
18. I5=bwareaopen(I4,2000);
    
19. figure(6),imshow(I5);title('从对象中移除小对象');
    
20. [y,x,z]=size(I5);
    
21. myI=double(I5);
    
22. tic
    
23. Blue_y=zeros(y,1);
    
24. for i=1:y
    
25.     for j=1:x
    
26.              if(myI(i,j,1)==1)
    
27.   
    
28.                 Blue_y(i,1)= Blue_y(i,1)+1;%蓝色像素点统计
    
29.             end  
    
30.      end      
    
31. end
    
32. [temp MaxY]=max(Blue_y);%Y方向车牌区域确定
    
33. PY1=MaxY;
    
34. while ((Blue_y(PY1,1)>=5)&&(PY1>1))
    
35.         PY1=PY1-1;
    
36. end   
    
37. PY2=MaxY;
    
38. while ((Blue_y(PY2,1)>=5)&&(PY2<y))
    
39.         PY2=PY2+1;
    
40. end
    
41. IY=I(PY1:PY2,:,:);
    
42. %%%%%% X方向 %%%%%%%%%
    
43. Blue_x=zeros(1,x);%进一步确定x方向的车牌区域
    
44. for j=1:x
    
45.      for i=PY1:PY2
    
46.             if(myI(i,j,1)==1)
    
47.                 Blue_x(1,j)= Blue_x(1,j)+1;               
    
48.             end  
    
49.      end      
    
50. end
    
51.   
    
52. PX1=1;
    
53. while ((Blue_x(1,PX1)<3)&&(PX1<x))
    
54.        PX1=PX1+1;
    
55. end   
    
56. PX2=x;
    
57. while ((Blue_x(1,PX2)<3)&&(PX2>PX1))
    
58.         PX2=PX2-1;
    
59. end
    
60. PX1=PX1-1;%对车牌区域的校正
    
61. PX2=PX2+1;
    
62.   dw=I(PY1:PY2-8,PX1:PX2,:);
    
63. t=toc;
    
64. figure(7),subplot(1,2,1),imshow(IY),title('行方向合理区域');
    
65. figure(7),subplot(1,2,2),imshow(dw),title('定位剪切后的彩色车牌图像')
    
66. imwrite(dw,'dw.jpg');
    
67. [filename,filepath]=uigetfile('dw.jpg','输入一个定位裁剪后的车牌图像');
    
68. jpg=strcat(filepath,filename);
    
69. a=imread('dw.jpg');
    
70. b=rgb2gray(a);
    
71. imwrite(b,'1.车牌灰度图像.jpg');
    
72. figure(8);subplot(3,2,1),imshow(b),title('1.车牌灰度图像')
    
73. g_max=double(max(max(b)));
    
74. g_min=double(min(min(b)));
    
75. T=round(g_max-(g_max-g_min)/3); % T 为二值化的阈值
    
76. [m,n]=size(b);
    
77. d=(double(b)>=T);  % d:二值图像
    
78. imwrite(d,'2.车牌二值图像.jpg');
    
79. figure(8);subplot(3,2,2),imshow(d),title('2.车牌二值图像')
    
80. figure(8),subplot(3,2,3),imshow(d),title('3.均值滤波前')
    
81. 
    
82. % 滤波
    
83. h=fspecial('average',3);
    
84. d=im2bw(round(filter2(h,d)));
    
85. imwrite(d,'4.均值滤波后.jpg');
    
86. figure(8),subplot(3,2,4),imshow(d),title('4.均值滤波后')
    
87. 
    
88. % 某些图像进行操作
    
89. % 膨胀或腐蚀
    
90. % se=strel('square',3);  % 使用一个3X3的正方形结果元素对象对创建的图像进行膨胀
    
91. % 'line'/'diamond'/'ball'...
    
92. se=eye(2); % eye(n) returns the n-by-n identity matrix 单位矩阵
    
93. [m,n]=size(d);
    
94. if bwarea(d)/m/n>=0.365
    
95.     d=imerode(d,se);
    
96. elseif bwarea(d)/m/n<=0.235
    
97.     d=imdilate(d,se);
    
98. end
    
99. imwrite(d,'5.膨胀或腐蚀处理后.jpg');
    
100. figure(8),subplot(3,2,5),imshow(d),title('5.膨胀或腐蚀处理后')
     
101. 
     
102. % 寻找连续有文字的块，若长度大于某阈值，则认为该块有两个字符组成，需要分割
     
103. d=qiege(d);
     
104. [m,n]=size(d);
     
105. figure,subplot(2,1,1),imshow(d),title(n)
     
106. k1=1;k2=1;s=sum(d);j=1;
     
107. while j~=n
     
108.     while s(j)==0
     
109.         j=j+1;
     
110.     end
     
111.     k1=j;
     
112.     while s(j)~=0 && j<=n-1
     
113.         j=j+1;
     
114.     end
     
115.     k2=j-1;
     
116.     if k2-k1>=round(n/6.5)
     
117.         [val,num]=min(sum(d(:,[k1+5:k2-5])));
     
118.         d(:,k1+num+5)=0;  % 分割
     
119.     end
     
120. end
     
121. % 再切割
     
122. d=qiege(d);
     
123. % 切割出 7 个字符
     
124. y1=10;y2=0.25;flag=0;word1=[];
     
125. while flag==0
     
126.     [m,n]=size(d);
     
127.     left=1;wide=0;
     
128.     while sum(d(:,wide+1))~=0
     
129.         wide=wide+1;
     
130.     end
     
131.     if wide<y1   % 认为是左侧干扰
     
132.         d(:,[1:wide])=0;
     
133.         d=qiege(d);
     
134.     else
     
135.         temp=qiege(imcrop(d,[1 1 wide m]));
     
136.         [m,n]=size(temp);
     
137.         all=sum(sum(temp));
     
138.         two_thirds=sum(sum(temp([round(m/3):2*round(m/3)],:)));
     
139.         if two_thirds/all>y2
     
140.             flag=1;word1=temp;   % WORD 1
     
141.         end
     
142.         d(:,[1:wide])=0;d=qiege(d);
     
143.     end
     
144. end
     
145. % 分割出第二个字符
     
146. [word2,d]=getword(d);
     
147. % 分割出第三个字符
     
148. [word3,d]=getword(d);
     
149. % 分割出第四个字符
     
150. [word4,d]=getword(d);
     
151. % 分割出第五个字符
     
152. [word5,d]=getword(d);
     
153. % 分割出第六个字符
     
154. [word6,d]=getword(d);
     
155. % 分割出第七个字符
     
156. [word7,d]=getword(d);
     
157. subplot(5,7,1),imshow(word1),title('1');
     
158. subplot(5,7,2),imshow(word2),title('2');
     
159. subplot(5,7,3),imshow(word3),title('3');
     
160. subplot(5,7,4),imshow(word4),title('4');
     
161. subplot(5,7,5),imshow(word5),title('5');
     
162. subplot(5,7,6),imshow(word6),title('6');
     
163. subplot(5,7,7),imshow(word7),title('7');
     
164. [m,n]=size(word1);
     
165. 
     
166. % 商用系统程序中归一化大小为 40*20,此处演示
     
167. word1=imresize(word1,[40 20]);
     
168. word2=imresize(word2,[40 20]);
     
169. word3=imresize(word3,[40 20]);
     
170. word4=imresize(word4,[40 20]);
     
171. word5=imresize(word5,[40 20]);
     
172. word6=imresize(word6,[40 20]);
     
173. word7=imresize(word7,[40 20]);
     
174. 
     
175. subplot(5,7,15),imshow(word1),title('1');
     
176. subplot(5,7,16),imshow(word2),title('2');
     
177. subplot(5,7,17),imshow(word3),title('3');
     
178. subplot(5,7,18),imshow(word4),title('4');
     
179. subplot(5,7,19),imshow(word5),title('5');
     
180. subplot(5,7,20),imshow(word6),title('6');
     
181. subplot(5,7,21),imshow(word7),title('7');
     
182. imwrite(word1,'1.jpg');
     
183. imwrite(word2,'2.jpg');
     
184. imwrite(word3,'3.jpg');
     
185. imwrite(word4,'4.jpg');
     
186. imwrite(word5,'5.jpg');
     
187. imwrite(word6,'6.jpg');
     
188. imwrite(word7,'7.jpg');
     
189. liccode=char(['0':'9' 'A':'Z' '苏豫陕鲁京辽浙']);  %建立自动识别字符代码表  
     
190. SubBw2=zeros(40,20);
     
191. l=1;
     
192. for I=1:7
     
193.       ii=int2str(I);
     
194.      t=imread([ii,'.jpg']);
     
195.       SegBw2=imresize(t,[40 20],'nearest');
     
196.         if l==1                 %第一位汉字识别
     
197.             kmin=37;
     
198.             kmax=43;
     
199.         elseif l==2             %第二位 A~Z 字母识别
     
200.             kmin=11;
     
201.             kmax=36;
     
202.         else l>=3               %第三位以后是字母或数字识别
     
203.             kmin=1;
     
204.             kmax=36;
     
205.         
     
206.         end
     
207.         
     
208.         for k2=kmin:kmax
     
209.             fname=strcat('F:\edge 下载\车牌识别3 第二版\车牌识别\程序与图像\字符模板\',liccode(k2),'.jpg');
     
210.             SamBw2 = imread(fname);
     
211.             for  i=1:40
     
212.                 for j=1:20
     
213.                     SubBw2(i,j)=SegBw2(i,j)-SamBw2(i,j);
     
214.                 end
     
215.             end
     
216.            % 以上相当于两幅图相减得到第三幅图
     
217.             Dmax=0;
     
218.             for k1=1:40
     
219.                 for l1=1:20
     
220.                     if  ( SubBw2(k1,l1) > 0 | SubBw2(k1,l1) <0 )
     
221.                         Dmax=Dmax+1;
     
222.                     end
     
223.                 end
     
224.             end
     
225.             Error(k2)=Dmax;
     
226.         end
     
227.         Error1=Error(kmin:kmax);
     
228.         MinError=min(Error1);
     
229.         findc=find(Error1==MinError);
     
230.         Code(l*2-1)=liccode(findc(1)+kmin-1);
     
231.         Code(l*2)=' ';
     
232.         l=l+1;
     
233. end
     
234. figure(10),imshow(dw),title (['车牌号码:', Code],'Color','b');
     
235. 
     
236. function [word,result]=getword(d)
     
237. word=[];flag=0;y1=8;y2=0.5;
     
238.     while flag==0
     
239.         [m,n]=size(d);
     
240.         wide=0;
     
241.         while sum(d(:,wide+1))~=0 && wide<=n-2
     
242.             wide=wide+1;
     
243.         end
     
244.         temp=qiege(imcrop(d,[1 1 wide m]));
     
245.         [m1,n1]=size(temp);
     
246.         if wide<y1 && n1/m1>y2
     
247.             d(:,[1:wide])=0;
     
248.             if sum(sum(d))~=0
     
249.                 d=qiege(d);  % 切割出最小范围
     
250.             else word=[];flag=1;
     
251.             end
     
252.         else
     
253.             word=qiege(imcrop(d,[1 1 wide m]));
     
254.             d(:,[1:wide])=0;
     
255.             if sum(sum(d))~=0;
     
256.                 d=qiege(d);flag=1;
     
257.             else d=[];
     
258.             end
     
259.         end
     
260.     end
     
261. %end
     
262.           result=d;
     
263.          
     
264.          
     
265.           function e=qiege(d)
     
266. [m,n]=size(d);
     
267. top=1;bottom=m;left=1;right=n;   % init
     
268. while sum(d(top,:))==0 && top<=m
     
269.     top=top+1;
     
270. end
     
271. while sum(d(bottom,:))==0 && bottom>=1
     
272.     bottom=bottom-1;
     
273. end
     
274. while sum(d(:,left))==0 && left<=n
     
275.     left=left+1;
     
276. end
     
277. while sum(d(:,right))==0 && right>=1
     
278.     right=right-1;
     
279. end
     
280. dd=right-left;
     
281. hh=bottom-top;
     
282. e=imcrop(d,[left top dd hh]);
     

复制代码

经过各位大神的知道问题解决了，就是在根目录下放一张名为car1的图片就可以了