机器视觉的道路识别结题报告.docx

1、机器视觉的道路识别结题报告机器视觉的道路识别结题报告篇一：机器视觉的车辆牌照识别深圳稻草人自动化培训 机器视觉的车辆牌照识别，更具自动化与智能化 之前在提起自动化或是智能化时，人们会不自觉的想到工业生产，这是因为自动化这个字眼进入中国，确实是以工业领域为登陆平台的。但是，随着时代的进步，自动化已经不再局限于工业领域，而是广泛的存在于各种领域，甚至与我们的生活息息相关。智能化交通系统是交通领域科学技术发展的产物，不管是对道路行驶车辆的监控拍照，还是公路收费领域的电子自动化收费系统，都是智能化的表现。与此同时，我们也不难发现，不管智能化交通系统应用在哪里，其都是通过对车辆牌照的识别来实现的，而目前

2、基于机器视觉的车辆牌照识别系统，是最为自动化与智能化的。 就像公民的身份证件一样，车辆牌照就是车辆的身份证明，是最真实、清晰的标志。车辆牌照识别系统能够自动获取车辆行进过程中的动态数据，准确识别车辆牌照字符，并进行快速的判断与处理。基于机器视觉的车辆牌照识别系统主要包括图像采集与处理、车牌定位、字符分割与识别等几个部分。 其一：图像采集与处理 机器视觉系统的主要功能就是将被摄目标转换成图像信号，再由图像处理系统转变成数字信号，最终根据这些信号进行各种的运算来抽取目标的特征进行比对。那么，基于机器视觉的车辆牌照识别系统就可以将车牌图像由模拟信号转换为数字信号，通过图像处理对车牌图像进行增强或是变

3、换，突出车牌的主要特征，为下一步的车牌定位打好基础。 其二：车牌定位车牌定位是基于机器视觉的车辆牌照识别系统的关键点，也是难点，噪声等都会影响到定位的精确度。其主要包括图像预处理，车辆牌照粗定位，车辆牌照精确定位等几个过程。图像预处理主要是区分已经变成电信号的信息、去除信号中的噪声与杂质信号、对文字格式进行规范，使获取的信息变得条理清晰。车牌照粗定位部分的作用就是大体选出几个待判断识别的候选车牌区域，为精确定位缩小判断范围。最后，车牌的精确定位便可以从候选区域判断出真正的区域并给出相应坐标。 (转自： 小 草 范文网:机器视觉的道路识别结题报告) 其三：字符分割与识别 车辆牌照的字符分割就是将

4、定位后的车牌区域中的每个字符从整个图像中分割出来，使其成为单个字符，以便于字符的识别。而字符识别就是基于机器视觉的车辆牌照识别系统最重要、最终的目的，它将对前面所有环节共同所得到的结果进行识别处理，最终得到需求的清晰准确的车辆牌照字符。 篇二：机器视觉 实验三报告 模板匹配法实现车牌识别实验三报告模板匹配法实现车牌识别一、实验目的结合印刷体字符识别方法，用模板匹配法实现车牌识别。要求：能实现车牌定位、字符分割和车牌中数字0-9的识别。二、实验设备微机三、实验内容及步骤1 上机编写程序。2 调试程序。3 根据实验结果，撰写实验报告。四、实验报告(一)对汽车图像进行图像转换、图像增强和边缘检测等:

5、1.载入车牌图像：I=imread(car1.jpg);figure(1),imshow(I);title(original image)；%将车牌的原图显示出来，结果如下： 2.将彩图转换为灰度图并绘制直方图： I1=rgb2gray(I);%将彩图转换为灰度图figure(2),subplot(1,2,1),imshow(I1);title(gray image);figure(2),subplot(1,2,2),imhist(I1);title(灰度图直方图);%绘制灰度图的直方图 结果如下所示： 3.用roberts算子进行边缘检测： I2=edge(I1,roberts,0.18,b

6、oth);%选择阈值0.18，用roberts算子进行边缘检测figure(3),imshow(I2);title(roberts operator edge detection image)；结果如下：4.图像实施腐蚀操作：se=1;1;1;I3=imerode(I2,se);%对图像实施腐蚀操作，即膨胀的反操作figure(4),imshow(I3);title(corrosion image); 5.平滑图像 se=strel(rectangle,25,25);%构造结构元素以正方形构造一个seI4=imclose(I3,se);%图像聚类、填充图像figure(5),imshow(I4

7、);title(smothing image); 结果如下所示： 6.删除二值图像的小对象 I5=bwareaopen(I4,XX);%去除聚团灰度值小于XX的部分figure(6),imshow(I5);title(remove the small objects);%用imshow函数显示滤波后图像 结果如下所示： （二）车牌定位 y,x,z=size(I5);%返回I5各维的尺寸，存储在x,y,z中 myI=double(I5);%将I5转换成双精度tic%tic表示计时的开始，toc表示计时的结束Blue_y=zeros(y,1);%产生一个y*1的零阵for i=1:yfor j=1

8、:xif(myI(i,j,1)=1)%如果myI(i,j,1)即myI的图像中坐标为(i,j)的点值为1，即该点为车牌背景颜色蓝色%则Blue_y(i,1)的值加1Blue_y(i,1)= Blue_y(i,1)+1;%蓝色像素点统计endendendtemp MaxY=max(Blue_y);%Y方向车牌区域确定%temp为向量yellow_y的元素中的最大值，MaxY为该值的索引PY1=MaxY;while (Blue_y(PY1,1)=5)&(PY11)PY1=PY1-1;endPY2=MaxY;while (Blue_y(PY2,1)=5)&(PY2PY2=PY2+1;endIY=I(

9、PY1:PY2,:,:);%行方向车牌区域确定% X方向 %Blue_x=zeros(1,x);%进一步确定x方向的车牌区域for j=1:xfor i=PY1:PY2if(myI(i,j,1)=1)篇三：北航机器视觉第一次作业八边形识别八边形识别本文将四个八边形的图像称为源图像（test.jpg），将单个多边形的图像称为待匹配图像（test1.jpg、test2.jpg）。步骤1：角点检测函数octagon_corner.m函数是利用课本上提到的Harris角点探测器，对整个图像提取角点。 步骤2：图像的预处理及角点匹配octagon_search.m和octagon_source.m函数，

10、分别对源图像和待匹配图像处理，用于将检测出的角点与各自归属的八边形匹配。步骤3：交比计算octagon_cal_much.m和octagon_cal_one.m函数，分别用于对源图像和待匹配图像计算交比。 步骤4：找到最像的八边形octagon_main.m主函数，将待匹配图像与源图像匹配。 匹配结果1： 匹配结果2： 代码 function theta=angle1(x1,y1,x2,y2,x3,y3)theta=acosd(dot(x1-x2,y1-y2,x3-x2,y3-y2)/(norm(x1-x2,y1-y2)*norm(x3-x2,y3-y2); functionoctagon_C

11、R=octagon_cal_much(point_devide) number_octagon=size(point_devide,3); octagon_angle=zeros(8,4,number_octagon);octagon_CR=zeros(8,number_octagon);x=zeros(4,1);y=zeros(4,1);l=zeros(4,1);for i=1:number_octagonfor j=1:8xx=point_devide(j,1,i);yy=point_devide(j,2,i);for k=1:4if(j+k=8)x(k)=point_devide(8,1

12、,i);y(k)=point_devide(8,2,i);elsex(k)=point_devide(rem(j+k,8),1,i);y(k)=point_devide(rem(j+k,8),2,i);endendoctagon_angle(j,1,i)=sind(angle1(x(1),y(1),xx,yy,x(3),y(3);octagon_angle(j,2,i)=sind(angle1(x(2),y(2),xx,yy,x(3),y(3);octagon_angle(j,3,i)=sind(angle1(x(1),y(1),xx,yy,x(4),y(4);octagon_angle(j,

13、4,i)=sind(angle1(x(2),y(2),xx,yy,x(4),y(4); octagon_CR(j,i)=(octagon_angle(j,1,i)/octagon_angle(j,2,i)/(octagon_angle(j,3,i)/octagon_angle(j,4,i); end endfunctionoctagon_CR=octagon_cal_one(point_devide)number_octagon=size(point_devide,3);octagon_angle=zeros(8,4);octagon_CR=zeros(8,1);x=zeros(4,1);y=

14、zeros(4,1);l=zeros(4,1);for j=1:8xx=point_devide(j,1);yy=point_devide(j,2);for k=1:4if(j+k=8)x(k)=point_devide(8,1);y(k)=point_devide(8,2);elsex(k)=point_devide(rem(j+k,8),1);y(k)=point_devide(rem(j+k,8),2);endendoctagon_angle(j,1)=sind(angle1(x(1),y(1),xx,yy,x(3),y(3);octagon_angle(j,2)=sind(angle1

15、(x(2),y(2),xx,yy,x(3),y(3);octagon_angle(j,3)=sind(angle1(x(1),y(1),xx,yy,x(4),y(4);octagon_angle(j,4)=sind(angle1(x(2),y(2),xx,yy,x(4),y(4); octagon_CR(j)=(octagon_angle(j,1)/octagon_angle(j,2)/(octagon_angle(j,3)/octagon_angle(j,4); endfunction ROW COLUMN=octagon_corner(Image)%Image = im2uint8(rgb

16、2gray(Image); dx = -1 0 1;-1 0 1;-1 0 1; %dx：横向Prewitt差分模版 Ix2 = filter2(dx,Image).2;Iy2 = filter2(dx,Image).2;Ixy = filter2(dx,Image).*filter2(dx,Image); %生成 9*9高斯窗口。窗口越大，探测到的角点越少。 h= fspecial(gaussian,9,2);A = filter2(h,Ix2); % 用高斯窗口差分Ix2得到AB = filter2(h,Iy2);C = filter2(h,Ixy);ow = size(Image,1);

17、ncol = size(Image,2);Corner = zeros(ow,ncol); %矩阵Corner用来保存候选角点位置,初值全零，值为1的点是角点%真正的角点在137和138行由(row_ave,column_ave)得到%参数t:点(i,j)八邻域的“相似度”参数，只有中心点与邻域其他八个点的像素值之差在%（-t,+t）之间，才确认它们为相似点，相似点不在候选角点之列t=20;%我并没有全部检测图像每个点，而是除去了边界上boundary个像素，%因为我们感兴趣的角点并不出现在边界上boundary=8;for i=boundary:ow-boundary+1for j=boun

18、dary:ncol-boundary+1nlike=0; %相似点个数if Image(i-1,j-1)Image(i,j)-t & Image(i-1,j-1)nlike=nlike+1;endif Image(i-1,j)Image(i,j)-t & Image(i-1,j)nlike=nlike+1;endif Image(i-1,j+1)Image(i,j)-t & Image(i-1,j+1)nlike=nlike+1;endif Image(i,j-1)Image(i,j)-t & Image(i,j-1)nlike=nlike+1;endif Image(i,j+1)Image(

19、i,j)-t & Image(i,j+1)nlike=nlike+1;endif Image(i+1,j-1)Image(i,j)-t & Image(i+1,j-1)nlike=nlike+1;endif Image(i+1,j)Image(i,j)-t & Image(i+1,j)nlike=nlike+1;endif Image(i+1,j+1)Image(i,j)-t & Image(i+1,j+1)nlike=nlike+1;endifnlike=2 &nlikeCorner(i,j)=1;%如果周围有0，1，7，8个相似与中心的（i,j） %那(i,j)就不是角点，所以，直接忽略e

20、nd;end;end;CRF = zeros(ow,ncol); % CRF用来保存角点响应函数值,初值全零CRFmax = 0; % 图像中角点响应函数的最大值，作阈值之用 t=0.04; %建议取0.04（机器视觉）% 计算CRF%工程上常用CRF(i,j) =det(M)/trace(M)计算CRF，那么此时应该将下面第105行的 %比例系数t设置大一些，t=0.1对采集的这几幅图像来说是一个比较合理的经验值 for i = boundary:ow-boundary+1for j = boundary:ncol-boundary+1if Corner(i,j)=1 %只关注候选点M =

21、A(i,j) C(i,j);C(i,j) B(i,j);CRF(i,j) = det(M)-t*(trace(M)2;if CRF(i,j) CRFmaxCRFmax = CRF(i,j);end;endend;end;%CRFmaxcount = 0; % 用来记录角点的个数t=0.1;% 下面通过一个3*3的窗口来判断当前位置是否为角点for i = boundary:ow-boundary+1for j = boundary:ncol-boundary+1if Corner(i,j)=1 %只关注候选点的八邻域if CRF(i,j) t*CRFmax& CRF(i,j) CRF(i-1,

22、j-1) .& CRF(i,j) CRF(i-1,j) & CRF(i,j) CRF(i-1,j+1) .& CRF(i,j) CRF(i,j-1) & CRF(i,j) CRF(i,j+1) .& CRF(i,j) CRF(i+1,j-1) & CRF(i,j) CRF(i+1,j).&CRF(i,j) CRF(i+1,j+1)count=count+1;%这个是角点，count加1else % 如果当前位置（i,j）不是角点，则在Corner(i,j)中删除对该候选角点的记录Corner(i,j) = 0;end;end;end;end;% toc(t1)p=9;pp=(p-1)/2;RO

23、W=zeros(1,count);COLUNM=zeros(1,count);number=1;for i=boundary:ow-boundary+1for j=boundary:ncol-boundary+1column_ave=0;row_ave=0;k=0;if Corner(i,j)=1for x=i-pp:i+pp %7*7邻域for y=j-pp:j+ppif Corner(x,y)=1% 用算数平均数作为角点坐标，如果改用几何平均数求点的平均坐标，对角点的提取意义不大row_ave=row_ave+x;column_ave=column_ave+y;k=k+1;endendendend