植物叶片脉络图像提取方法设计

1、植物叶片脉络图像提取方法设计摘要图像分割是根据某种图像特征的一致要求， 将一幅图像分割为多个区域。 它 是图像理解和物体识别的第一步， 也是一经典难题。 随着科学技术的发展， 它广 泛的应用于各个领域。植物叶片是植物的重要器官， 是研究植物生长状况的重要依据， 也是比较容 易获得的研究素材。 叶片脉络和形态是研究植物的重要出发点， 叶片脉络的提取 是实现完整的植物叶特征提取及无损植物叶特征提取中重要的一步。 传统的植物 形态数据测量方法主要是用尺测量叶片长度等手工方法，这种方法由于受人为、 自然等因素的影响较多， 导致测量精度差、测量结果不理想， 基本上已经被淘汰。 利用图像处理技术进行植物脉

2、络提取， 是随着计算机软、 硬件技术及信息技术的 进步而迅速发展起来的新兴研究领域。本文介绍了两种利用图象分割技术提取叶片脉络的基本方法：分水岭法， 二值化法。通过 MATLAB 进行了仿真实验，并将几种不同的分割算法结果进行了 比较关键词 ： 图像分割，特征提取，分水岭算法 , 二值化一 分水岭算法1 分水岭算法综述分水岭算法是新近发展起来的数学形态学图像分割方法。将这种方法用于 图像处理是由S.Beuchre , F.Meyer最早提出的。Beucher和Vincent等人的研 究工作使得分水岭算法的理论得以建立，并且大量用于灰度图像的分割。分水岭算法的一个实例就是从太平洋到大西洋的一条水

3、流把美国分为两个 区域。分水岭界限分开了这样的两个区域,被称为太平洋和大西洋的“集水盆 地"。由于这些“集水盆地”的存在,相对应的,这两个大洋之间就有一些最小 连接。现在,在图像处理领域和数学形态学领域, 数字图像经常被认为是拓扑框架。 在给定图像的拓扑表示中, 每个点的数字化值代表这个点的梯度值。 这种表示非 常有用，首先，这样能够使研究者对于被研究图像的转换有一个更清晰的认识。 另外，利用这种表示，一些概念比如极小区、集水盆地、分水岭等都可以在数字 图像中更好的定义。近年来分水岭算法在灰度图像分割算法的研究中受到重视。分水岭算法的 主要应用是从背景中提取近乎一致的目标对象，那些在

4、灰度级上变化较小的区 域的梯度值也较小，而不同目标对象在相邻接触的灰度变化往往比区域内部像 素间的差异要大得多，对应的梯度值也较大。因此，实际上分水岭算法与图像 的梯度有更大的关系，而不是图像本身，一般以梯度图像作为参考图像进行分 割。分水岭算法有多种实现方法，较典型的有基于标记的分水岭变换和欧氏距 离应设法。传统的分水岭方法模拟一个水浸入的过程。首先将图像像素的灰度看做高 度，图像中不同灰度值的区域就对应于山峰和山谷盆地，盆地的边缘就是所谓 的分水岭。分水岭的来历就是将一幅图像逐渐浸入一个湖中，图像的最低点先 进水，然后水逐渐浸入整个山谷盆地。当水位到达盆地的边缘高度就将溢出， 这时在水溢出

5、处建立堤坝，如此直到整个图像沉入水中，所建立的堤坝就成为 分开盆地的分水岭。图1演示了这种方法的工作原理。图1分水岭算法假定图中的物体灰度值低，而背景的灰度值高。该图显示了沿扫描线的灰 度分布，该线穿过两个靠的很近的物体。图像最初在一个灰度值上基于分水岭算法的红外图像分割方法二值化。灰度值把图像分割成正确数目的物体，但它 们的边界偏向物体内部。随后阈值逐渐增加，每一次增加 1个灰度级，物体的边 界将随着阈值增加而扩展。当边界相互接触时，这些物体并没有合并。因此， 初次接触的点变成了相邻物体间的最终边界。这个过程在阈值达到背景灰度级之前终止，也就是说，在被恰当分割的物体的边界正确地确定时终止2基

6、本概念设f(P)为灰度级函数，p为像素坐标，门z2的相邻像素pl N(p),距离函数为dist ( p, p')(在4连通或8连通域中，对所有的邻域点，距离为 1)定义1：较低斜坡LS(p)二maxP' N ( p)仃(P) - f(P') dist (p, p')f(P')空 f (P)(1.1)定义2:斜率最大的相邻点NLS ( p)='f (P) -f (P')、pS(p) = LS(p), f(p') £ f(p) >dist (p, p')( 1.2)定义3:分水岭分割如果每个区域最低处mi有唯一

7、的标记L(mi),称为分水岭分割。对于每NLS ( p) = '的 p 二匚：，则(1.3)(1.4)p- NLS (p), L(p)二 L(p)定义4:集水盆C (mJ = p L( p) = L(mj3分水岭算法的实现通常分水岭算法的实现过程采用模拟淹没过程。首先，检测出局部的最低点，并用数字逐一标识出；接着，利用等级队列模拟淹没过程。这样的序列由H个先进先出(FIFO)序列组成，每个序列有H灰度级。第h个FIFO序列的图像中灰度 级为h的像素数。这个数据结构是用于存储进行处理的像素顺序。初始时1等级队列中为淹没的初始点。灰度级为 h的点进入队列的第h个FIFO队列。等级队 列的灰

8、度级由低到高。当这个点被淹没后，将它从队列中删除，并将它的标记传 递给相邻的点；后一个点进入队列，直到所有的点被淹没，从而使各个区域正确划分另外，A.Bieniek提出了基于连接像素的有效的分水岭算法.对于像素点p, 它的相邻点NLS(p)的标识将作为此点处理的前提条件。然而，局部的最低点并 没有斜率最大的相邻点。因此对于这些点，应采用另一种连接关系：最低点 p 的所有邻域点中和p点灰度值相同的点属于同一区域。因此，局部最小的位置应有这一区域中位置编码最小的点的位置号来代表。在处理过程中，除了输入图像之外，还应有图像来存储每个像素的标识。这里，用像素点的位置来标识各个区域。这个算法只需对图像进

9、行 次扫描，就可完成分割。具体算法为输入图像f，输出图像|。1) 第 1 次扫描图像像素 pfq- p ； for each(pN(p)and f p' ： f p)N(p) and fp' ： f (p)if (fp'H： f(p) qp'；if (q = p) l【P- q'；else l p - PLATEAU ； 2)第二次扫描图像像素 pif (l p = PLATEAU )LP P ；for each(P' N p® ( p) and f p' = f p)r 一 FIND (l, p) ； r'一 FIND

10、 (l, p') ； l门一 lr' 一 min( r,r') ； 3)第3次扫描图像像素plp FIND (L, P) ； FIND (l,u) for (r :一 u ； lr = r', r :一 lr)；Return r； 在第3次扫描中，有较小的相邻点的像素p，它的标记为q NLS (P)。如果 像素没有较小的相邻点，那么，这一点为局部最小点，记为PLATEAU在第2轮扫描中，计算每3个区域最小的标记。利用函数 FIND(l,p)，将区域最小的标记 用这一区域最小点中最小的坐标值表示。在第3轮的扫描中，将所有的像素的标记用最小坐标值替代。这样，图像中的

11、每一个点均满足定义3中有关分水岭分割的定义，从而完成了分水岭图像分割。分水岭算法原理流程图如图2所示:图2分水岭原理流程图分水岭算法程序：M=imread（'叶子 4.bmp'）;figure,imshow(M) M1 = watershed(M);figure,imshow(M1)分水岭分割结果如图3所示:图3分水岭分割结果二值化法1二值化方法基本概念在数字图像的处理中，二值化图像占有很重要的地位。特别是在实用的图像 处理系统中，以二值图像处理为核心构成的系统是很多的，要进行二值图像的处理与分析，这是因为在实用的图像处理系统中，要求处理的速度高、成本低，直 接对信息量大的灰度

12、图像进行处理花销太大，不是上策，而二值化后的图像能够用几何学中的概念进行分析和特征描述，比起灰度图像来说方便得多，因而二值图像处理目前己成为图像处理中的一个独立的重要分支而获得广泛的应用。所谓图像的二值化是指通过设定阈值(threshold)把灰度图像变成仅用两个 值分别表示图像的目标和背景的二值图像。对于叶片图像，可以将其看作是叶脉 线和叶体线两种区域，叶片图像的二值化就是将这两种区域分割开来，并以像素 值0来量化叶脉线区域，以像素值255量化叶体线区域，由此得到的图像就是叶 片的二值化图像。假设一幅灰度图的像素值为f(x，y)，设阈值T,贝卩f (x, y)=丿0f (x,y) _Tf (

13、x,y)乞T(2.1)通常，二值图中的1来表示目标子图，0来表示背景子图。叶片的二值化方 法很多，关键在于阈值T的选取，而T的取值方法又取决于所采用的二值化技术, 以阈值的设置是二值化操作的关键所在。 而且二值化方法的优劣直接关系二值化效果的好坏，只有得到较好的值化效果图，才能提高最后的特征匹配率，从而提 高处理后图像的质量5. 3. 1基于灰度空间分布的二值化方法实现由于叶片图像是由叶脉线和叶体线相间组成的， 在叶片灰度图中，叶脉线的 灰度值比叶体线的灰度值要大，因此可以利用灰度值空间分布特点进行叶片图像 的二值化。基于灰度空间分布的二值化方法不考虑叶片图像特有的方向信息，只 从叶片图像灰度

14、值出发进行阈值T的计算。阈值可以分为两类：全局阈值和局部 阈值。全局阈法是对整个图像采用一个阈值进行划分； 局部阈值法是将图像分成 一些子块，对于一块选定一个阈值。全局阈值二值化是给出一个针对整幅图像的阈值 T，所有大于T的像素点都 被二化为I，小于等于T的像素点都被二值化为0。对于质量较低的灰度图像，常用局部动态阈值法26，该方法由像素的灰度值 以及该像素周围点的局部灰度特性来确定阈值，从而进行二值化。本文所采用的全局阈值算法是一种简单有效的二值化方法旳。全局值法是与灰度直方图紧密联系的。灰度直方图的横坐标表示图像的灰度， 纵坐标表示该灰 度在图像全体像素中出现的额度。 图4是一个简单图像的态度直方图，灰度直方 图明显地区分了对象与背景。我们选取适当的阈值就可以将图像中感兴趣的域分 割出来。二值化方法设计程序为：I = imread('叶子 4.bmp'); %读入图像imshow(l)%显示原图像l=l>T;%T为选取阈值figure, imshow(l);%显示处理后图像叶片的灰度直方图如图4所示。图4叶片灰度直方图二值化方法处理结果如图5所示:图5不同阈值对应的二值化处理图像由处理图像可看出，当阈值为125左右时分割效果