车牌分割matlab开发程序

suffer

来源于人民邮电出版社Visual C++/MATLAB 图像处理与识别实用案例精选，具体原理自己看书吧

1. I=imread('Car.jpg');
   
2. [y,x,z]=size(I);
   
3. myI=double(I);
   
4. %%%%%%%%%%% RGB to HSI %%%%%%%%
   
5. tic % 测定算法执行的时间，开始计时
   
6. 
   
7. %%%%%%%%%%% 统计分析 %%%%%%%%%%%%%%%
   
8. %%%%%%%% Y 方向 %%%%%%%%%%
   
9. Blue_y=zeros(y,1);
   
10. for i=1:y
    
11. for j=1:x
    
12. if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3)<=142)&&(myI(i,j,3)>=119)))
    
13. % 蓝色RGB的灰度范围
    
14. Blue_y(i,1)= Blue_y(i,1)+1; % 蓝色象素点统计
    
15. end
    
16. end
    
17. end
    
18. [temp MaxY]=max(Blue_y); % Y方向车牌区域确定
    
19. PY1=MaxY;
    
20. while ((Blue_y(PY1,1)>=5)&&(PY1>1))
    
21. PY1=PY1-1;
    
22. end
    
23. PY2=MaxY;
    
24. while ((Blue_y(PY2,1)>=5)&&(PY2<y))
    
25. PY2=PY2+1;
    
26. end
    
27. IY=I(PY1:PY2,:,:);
    
28. 
    
29. %%%%%%%% X 方向 %%%%%%%%%%
    
30. Blue_x=zeros(1,x); % 进一步确定X方向的车牌区域
    
31. for j=1:x
    
32. for i=PY1:PY2
    
33. if((myI(i,j,1)<=30)&&((myI(i,j,2)<=62)&&(myI(i,j,2)>=51))&&((myI(i,j,3)<=142)&&(myI(i,j,3)>=119)))
    
34. Blue_x(1,j)= Blue_x(1,j)+1;
    
35. end
    
36. end
    
37. end
    
38. 
    
39. PX1=1;
    
40. while ((Blue_x(1,PX1)<3)&&(PX1<x))
    
41. PX1=PX1+1;
    
42. end
    
43. PX2=x;
    
44. while ((Blue_x(1,PX2)<3)&&(PX2>PX1))
    
45. PX2=PX2-1;
    
46. end
    
47. PX1=PX1-2; % 对车牌区域的修正
    
48. PX2=PX2+2;
    
49. 
    
50. Plate=I(PY1:PY2,PX1-2:PX2,:);
    
51. t=toc % 读取计时
    
52. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
53. figure,imshow(I);
    
54. figure,plot(Blue_y);grid
    
55. figure,plot(Blue_x);grid
    
56. figure,imshow(IY);
    
57. figure,imshow(Plate);
    
58. 
    
59. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
60. %%文字识别子函数
    
61. 
    
62. function [Numeral]=Recognition(StrokeTop,StrokeLeft,StrokeRight,StrokeBottom,Comp)
    
63. % 采用四边的轮廓结构特征和笔划统计（仅针对 0 和 8）识别残缺数字
    
64. % Comp 是用于识别 0和8 的底部补充信息
    
65. StrT='T';
    
66. StrL='T';
    
67. StrR='T';
    
68. StrB='T';
    
69. RStr='T'; % 用于保存识别出的数字
    
70. [temp XT]=size(StrokeTop);
    
71. [temp XL]=size(StrokeLeft);
    
72. [temp XR]=size(StrokeRight);
    
73. %[temp XB]=size(StrokeBottom);
    
74. 
    
75. for Ti=2:XT
    
76. if (StrokeTop(Ti)=='C')
    
77. if ((XL==2)&&(XR==2))
    
78. if ((Comp>=3)||((StrokeBottom(2)~='C')&&(StrokeLeft(2)=='C')&&(StrokeRight(2)=='C')))
    
79. RStr='8';
    
80. else
    
81. RStr='0';
    
82. end
    
83. 
    
84. else if ((StrokeLeft(XL)=='L')&&(StrokeLeft(XL-1)=='P')&&(StrokeLeft(2)~='C'))
    
85. RStr='2';
    
86. else if ((StrokeLeft(2)=='C')&&(XL>=3)&&(StrokeLeft(3)=='P'))
    
87. RStr='9';
    
88. else if (XL>2)
    
89. for Li=2:XL
    
90. if (StrokeLeft(Li)=='P')
    
91. RStr='3';
    
92. end
    
93. end
    
94. else if (XL==2)
    
95. for Ri=2:XR-1
    
96. if (StrokeRight(Ri)=='P')
    
97. RStr='6';
    
98. end
    
99. end
    
100. end
     
101. end
     
102. end
     
103. end
     
104. end
     
105. else if (StrokeTop(Ti)=='V') % Top
     
106. if ((XR==2)&&(StrokeRight(2)=='C')) % 数字 3 右端只有一个结构
     
107. RStr='3';
     
108. else if ((XR==2)&&((StrokeLeft(2)=='P')||(StrokeLeft(3)=='P')||(StrokeLeft(XL)=='V')))
     
109. RStr='7';
     
110. else if (XR>2)
     
111. for Ri=2:XR
     
112. if (StrokeRight(Ri)=='P')
     
113. RStr='5';
     
114. end
     
115. end
     
116. end
     
117. end
     
118. end
     
119. else if ((StrokeTop(Ti)=='L')) % Top
     
120. if ((StrokeRight(XR)=='V')&&(XR<=2))
     
121. RStr='4';
     
122. else
     
123. for Ri=2:XR
     
124. if ((StrokeRight(Ri)=='L')||(StrokeRight(Ri)=='C'))
     
125. RStr='6';
     
126. end
     
127. end
     
128. end
     
129. end
     
130. 
     
131. end %对应于'V'
     
132. end % 对应于'C'
     
133. end % 对应于第一个 for 循环
     
134. Numeral=RStr; 来源研学论坛
     
135. Matlab实现
     
136. ① 装入待处理图象并显示原始图象，如①图所示：
     
137. clear;
     
138. close all;
     
139. i=imread('原始图象');
     
140. figure,imshow(i);
     
141. ......
     
142. Matlab程序：
     
143. clear;
     
144. close all;
     
145. i=imread('原始图象');
     
146. figure,imshow(i);
     
147. %装入待处理图象并显示原始图象
     
148. s=strel('disk',10);
     
149. bg=imopen(i,s);
     
150. figure,imshow(bg);
     
151. %对‘原始图象’进行开操作得到图象背景
     
152. bw=imsubtract(i,bg);
     
153. figure,imshow(bw);
     
154. %用原始图象与背景图象作减法
     
155. %figure,imhist(i);
     
156. %求图象直方图
     
157. level=0.5;
     
158. bw2=im2bw(bw,level);
     
159. figure,imshow(bw2);
     
160. %通过图象的直方图取得阈值，将图象二值化

复制代码

ginger

好的，学习中！

mengmengmeng

:victory: ding ha

ChaChing

个人看法:若能给齐Car.jpg此文件, 可能较容易学习!