基于matlab的物料大小分级算法的实现①.docVIP

基于matlab的物料大小分级算法的实现① 　　[摘 要]本文以摄取的苹果图像为例,利用matlab对苹果的大小分级算法进行了仿真,为后续开发农产品物料图像分级系统的研究做铺垫。 　　[关键词]物料 大小分级 matlab 视觉检测 　　[中图分类号]P272[文献标识码]A[文章编号]1007-9416(2009)11-0071-02 　　 　　分级是把物料按其特征分为不同等级。筛分是根据物料颗粒大小将混合物料进行分级,采用机械式的分级筛是常用的高效的方法,但非坚果类农产品物料多半不能筛分,因为机械震动造成的碰撞和摩擦易造成被分级物料的损坏。近年来基于图像例如物料的形状、颜色、大小等的新型分级方法不断涌现,拓展了传统局限于物料大小的分级概念,使分级的领域、精度、效率、智能化程度以及环保得到深化和提高。基于图像的物料分级技术是我国的薄弱环节,不少产品靠从国外引进,深入研究和提高我国基于图像的物料分级水平,是十分有意义的工作。 　　采用嵌入式的视觉检测系统对物料大小进行自动分级可以提高效率。相对于人工分级具有很多优点,例如可以一次性地实现多因素分级,无损分级,分级标准统一,降低机械噪音和能耗,提供各种分类指标的统计信息,为物料的日后的品质改进,产品的销售提供可靠数据。因此,研究更快更准确的智能化的图像分级系统对现代化农业有着举足轻重的作用。 　　表面缺陷和形状大小是影响苹果品质的重要指标,本文从判断苹果图像大小出发,对大小分级的几种算法在PC机的平台上用MATLAB仿真软件进行分析和比较。 　　 　　1 图像预处理过程介绍(见图1) 　　 　　1.1 图像采集及背景去除 　　目前,在试验中苹果手工放置在一个光照匣子中,里面放进一个背景光源,cmos摄像头对其进行图像采集,拍摄的图片为 80*64的.jpg的图像。一幅图像可以表示为一个像素为m*n的矩阵。 　　1.2 灰度处理 　　灰度处理就是将真彩图变为只有黑白颜色的灰度图,但也并不是纯黑或者纯白。一幅灰度图像就是一个数据矩阵,它的数据代表了在一定范围内的颜色灰度值。 　　1.3 滤波消噪 　　滤波消噪有很多种方法,小波分析,平滑法,中值滤波法,数字滤波器等。本文采用中值滤波medfilt2函数对??度图像滤波消噪,它在一定程度上消除干扰噪声和点状噪声。 　　1.4 二值化并取反(最佳阈值) 　　图像的二值化就是将图像的像素仅以0和1来表示,每个像素只取2个灰度值。利用函数graythresh取得最佳的阈值。大于阀值的像素以黑色表示,小于阈值的像素以白色表示。对图像求反就是将灰度图的灰度值翻转,就是让黑变白,让白变黑。 　　1.5 边缘检测 　　两个灰度不同的区域存在着边缘,二值化后图像的矩阵只存在着0,1这两个灰度值,本文利用bwperim函数进行边缘检测。检测后边缘的像素为1,其它像素为0。 　　1.6 特征提取 　　对于图像边缘的面积,周长,偏心率,似圆度等一些需要关注并获得的参数进行提取获得,本文对图像边缘的面积参数进行了计算和分析。 　　 　　2 苹果大小分级算法及matlab仿真试验 　　 　　对于农产品物料分级中的尺寸分级,Paulus[2]等提出了三种分级依据的尺寸:(1)面积,(2)直径,(3)体积。本文以半径为分级的尺寸。 　　2.1 苹果的截面积 　　利用截面积求像素点的方法,像素点越多,尺寸越大。仿真程序: 　　v=imread(d:\A.jpg);% 读入图像。 　　imshow(v) % 显示真彩图。 　　w=rgb2gray(v); %化为灰度图像。 　　imshow(w)%显示灰度图。 　　L=medfilt2(w); %对灰度图像中值滤波。 　　level=graythresh(L); %获得灰度图像的最佳阈值。 　　bw=im2bw(L,level);%二值化图像。 　　imshow(~bw)%显示二值化的图像。 　　bwarea(~bw)%计算二值化后的像素点,得到苹果的像素点(见图2)。 　　结果:苹果A,B的像素点数分别为S1A=1241.1,S1B=1026.5。由s=3.14*r*r, r1A=19.88,r1B=18.08(单位:像素)。 　　2.2 通过求取圆心的坐标来求得平均半径 　　图像进行处理后,得到苹果边缘的像素,将边缘像素点提取出来,并存为直角坐标的数组,得到m个数组。像素的横坐标方向是从左向右,纵坐标方向是从上向下。第一步,求取苹果边缘像素的中心点O,O点坐标(x0,y0)。其中,。第二步,求边缘的M个点与中心点O(x0,y0)的距离,记为d1,d2……dM。第三步,求的距离的平均值r2作为苹果的半径,r2=(d1+d2+…+dM)/M