数字图像处理及应用：第七章图像分析.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第,7,章,图像分析,对图像进行增强、恢复、编码等处理时，输入是图像，所要求的输出是一幅近似于输入的图像，这是此类处理的一个特点。图像处理的另一个主要分支是图像分析或景物分析。这类处理的输入仍然是图像，但是所要求的输出是已知图像或景物的描述。这类处理基本上用于自身图像分析和模式识别一类的领域。,1,）输入是文字组成的二值图像，输出是读出该段文字；,2,）输入是血球照片，输出是血球数量；,3,）输入是细胞图像，输出是细胞类型；,4,）输入是遥感照片，输出是地貌、植被描述等。,这些都是图像分析的典型例子。,描述一般是针对图像或景物中的特定区域或目标。为了描述，首先要进行分割，有些分割运算可直接用于整个图像，而有些分割算法只适用于已被局部分割的图像。例如，分割染色体的处理，可先用设置门限的方法把染色体和背景分割开来，然后可采用尺寸大小、形状等准则进一步将其分割成单个染色体。,值得注意的一点是，没有唯一的、标准的分割方法，因此，也就没有规定成功分割的准则。,本章只讨论一些最基本的分割、描述方法。,7.1,分割,(segmentation),分割的目的是把图像空间分成一些有意义的区域。例如一幅航空照片，可民分割成工业区、住宅区、湖泊、森林等等。可以以逐个像素为基础去研究图像分割，也可以利用在规定领域中的某些图像信息去分割。分割的依据可建立在相似性和非连续性两个基本概念之上。,7.1.1,灰度阈值法分割,7.1.2,样板匹配,7.1.3,区域生长,7.1.4,区域聚合,假如有如下形状的直方图,图,7,1,图像 的直方图,由直方图可以知道图像 的大部分像素灰度值较低，其余像素较均匀地分布在其他灰度级上。由此可以推断这幅图像是由有灰度级的物体叠加在一个暗背景上形成的。可以设一个阈值 ，把直方图分成两个部分，如图所示。的选择要本着如下原则：应尽可能包含与背景相关连的灰度级，而 则应包含物体的所有灰度级。,当扫描这幅图像时，从 到 之间的灰度变化就指示出有边界存在。当然，为了找出水平方向和垂直方向上的边界，要进行两次扫描。也就是说，首先确定一个门限 ，然后执行下列步骤：,第一，对 的每一行进行检测，产生的图像的灰度将遵循如下规则,(7,1),式中 是指定的边缘灰度级，是背景灰度级。,第二，对 的每一列进行检测，产生的图像的灰度将遵循下述规则,(7,2),为了得到边缘图像，可采用下述关系,（,7,3,）,为了提高边缘抽取能力，一种方法是把图像变成二值图像。例如，图像 的灰度级范围是 ，设 是 和 之间的一个数，那么 可由式,(7,4),表示,(7,4),另一方法是把规定的灰度级范围变换为，而把范围以外的灰度级变换为,0,，例如,(7,5),(7,6),那么，在分割中如何设置最佳阈值呢？,假设一幅图像是由背景和物体组成。其中，物体像素的灰度级具有正态概率密度 ，其均值为 ，方差为 ；而背景像素的灰度级也具有正态概率密度 ，其均值为 ，方 差为 。,物体占图像总面积的比为 ，背景占总面积的比为 ，所以这幅图像总的灰度级概率密度为,(7,7),假设对图像设置一阈值 ，并且把小于 的全部点称为目标物体点，而把大于等于 的所有点称为背景点。,把背景错归为物体点的概率为 ，,把物体点错归为背景点的概率为 ，,则有,(7,8),(7,9),总的错分概率为,(7,10),要求得式,(8,10),的最小阈值，可将上式对 微分，并令其结果为,0,，则得到,(7,11),因为,(7,12),(7,13)代入式,(8,11),，并取对数,（,7,14),或者,（,7,15),由这个二次方程可以求解出 值。如果 ，那么,(7,16),这就是最佳门限,对于复杂图像，在许多情况下对整幅图像用单一阈值不能给出良好的分割结果。例如，图像是在光亮背景上的暗物体，但由于照射光的不均匀，虽然物体与背景始终有反差，但在图像的某一部分物体和背景两者都比另一部分亮。因此，在图像的一部分能把物体和背景精确地分开的阈值，对另一部分来说，可能把太多的背景也当作物体分割下来了。,克服这一缺点有如下一些方法：,如果已知在图像上的位置函数描述不均匀照射，就可以设法利用灰度级校正技术进行校正，然后采用单一阈值来分割；,另外一种方法是把图像分成小块，并对每一块设置局部阈值。但是，如果某块图像只含物体或只含背景，那么对这块图像就找不到阈值。这时，可以由附近的像块求得的局部阈值用内插法给此像块指定一个阈值。,8.1.1,灰度阈值法分割,8.1.2,样板匹配,8.1.3,区域生长,8.1.4,区域聚合,在数字图像处理中，样板是为了检测某些不变区域特性而设计的阵列。样板可根据检测目的不同而分为点样板、线样板、梯度样板、正交样板等等。,点样板的例子如图,7,2,所示。下面用一幅具有恒定强度背景的图像来讨论。,1,）、点样板,-1,-1,-1,-1,8,-1,-1,-1,-1,图,7,2,点样板,假定小 块 之 间 的 距 离 大 于,，这里 、分别是在,x,和,y,方向的取样距离，用点样板的检测步骤如下：,样板中心（标号为,8,）沿着图像从一个像素移到另一个像素，在每一个位置上，把处在样板内的图像的每一点的值乘以样板的相应方格中指示的数字，然后把结果相加。,如果在样板区域内所有图像的像素有同样的值，则其和为零。另一方面，如果样板中心位于一个小块的点上，则其和不为零。,如果小块在偏离样板中心的位置上，其和也不为零，但其响应幅度比起这个小块位于样板中心的情况时要小一些，这时，可以采用阈值法清除这类较弱的响应，如果其幅度值超过阈值，就意味着小块被检测出来了；如果低于阈值则忽略掉。,例如，设 代表,3,3,模板的权，并使 为模板内各像素的灰度值。从上述方法来看，应求两个矢量的积，即：,(7,24),式中,(7,25),(,7,26),设置一阈值 ，如果,(,7,27),我们认为小块已检测出来了。这个步骤可很容易地推广到,n,n,大小的样板，不过此时要处理,n,2,维矢量。,线检测样板如图,7,3,所示。其中，样板,(a),沿一幅图像移动，它将对水平取向的线（一个像素宽度）有最强的响应。对于恒定背景，当线通过样板中间一行时出现最大响应；,样板,(b),对,45,方向的那些线具有最好响应；样板,(c),对垂直线有最大响应；样板,(d),则对,-45,方向的那些线有最好的响应。,图,7,3,线样板,设 是图,7,3,中四个样板的权值组成的九维矢量。与点样板的操作步骤一样，在图像中的任一点上，线样板的各个响应为 ，这 里,i,=1,、,2,、,3,、,4,。,此处,X,是样板面积内九个像素形成的矢量。给定一个特定的,X,，希望能确定在讨论问题的区域与四个线样板中的哪一个有最相近的匹配。如果第,i,个样板响应最大，则可以断定,X,和第,i,个样板最相近。,换言之，如果对所有的 值，除 外，有：,(7,28),就 可 以 说 和 第 个 样 板 最 接近。如果 ，,=2,、,3,、,4,，可以断定 代表的区域有水平线的性质。,对于边缘检测来说也同样遵循上述原理。通常采用的方法是执行某种形式的二维导数。类似于离散梯度计算，考虑,大小的模板，如图,7,4,所示。,图,7,4 3 3,样板,考虑,的图像区域，及 分别用下式表示,(7,29),(7,30),采用绝对值的一种定义为,在 点的梯度为,(7,31),(7,32),梯度模板如图,7,5,所示。,图,7,5,梯度样板,把图,7,5,的区域与式,(7,29),比较，可以看出 为第一行和第三行的差，其中最靠近,e,的元素（,b,和,h,）,的加权等于角偶上权值的两倍，因此，代表在,x,方向上导数的估值。式,(7,31),和式,(7,32),可用图,7,5,中两个样板来实现。,边缘检测也可以表示成矢量，其形式与线样板检测相同。如果 代表所讨论的图像区域，则7,33),(7,34),这里 ，是图,7,5,中的两个样板矢量。分别代表它们的转置。,这样，梯度公式,(7,31),和式,(7,32),变为式,(7,35),和式,(7,36),的形式。,(7,35),(7,36),检测点、线和边缘的矢量公式可用下边的方法实现。,假定有两个只有三个元素的样板，此时，则有两个矢量 和 ，它们都是三维的。又假定 和 都是正交的和归一 化 的，因 此，它 们 都 有 单 位 幅 值。,图,7,6 X,向单位矢量 的投影,和 项分别等于在相应矢量,和 上,X,的投影。对于 来说,(,7,37),这里 是两个矢量间的夹角。,因为 ，因此有,(,7,38),这就是 在 上的投影，这种情况如图,7,6,所示。对 来说，亦然。,现在假 定 有 三 个 正 交 的 单 位 矢,量 、分别与三个三点样板相对应，,那么乘积 、代表在三个矢,量 、上的投影。其几何关系 如 图,7,7,所示。,图,7,7 X,向 确定的子空间的投影,假定样板和是检测线的，而是检测点的，代表的这个区域是更象一条线呢还是更象一个点呢？为了回答这一问题，把 投影到 、的子空间上去，和子空间的夹角可以说明 更接近于线还是更接近于点。,(,7,39),这可以从图,77,的几何关系上看出来。在,由 和 所确定的平面上投影的幅度可由式,(735),表示，而 的幅度由下式表示：,和其投影间的夹角为7,40),同理，向 子空间上投影的夹角可由下式表示：,(8,41),那么，如果 ，就说明 所代表的区域更接近于线特性而不是点特性。,如果考虑,3,3,的模板，则问题就成为维的，但前边讨论的概念仍然适用。这里，需要个维正交矢量形成一个完整的基。这个模板如图,7,8,所示。其中前四个模板,(a),、,(b),、,(c),、,(d),适合于边缘检测；,(e),、,(f),、,(g),、,(h),四块模板适合于检测线；最后一块模板,(,i,),则正比于一幅图像中模板所在区域的像素平均值。,图,7,8,正交模板（边缘）,图,7,8,正交模板,如果由 代表的,3,3,区 域，并 假 定 矢 量 1,2,3,，,9,是归一化的，从前边的讨论有：,(8,42),(8,43),(8,44),式中 分别是 向边缘、线和平均子空间投影的幅度。,同样道理，有,(8,45),(8,46),(8,47),采用这种检测方法可扩展到其它基与维数，只要基本矢量是正交的就可以。,式中 是 与它在边缘、线和平均子空间投影之间的夹角。,7.1.1,灰度阈值法分割,7.1.2,样板匹配,7.1.3,区域生长,7.1.4,区域聚合,分割的目的是要把一幅图像划分成一些区域，对于这个问题的最直接的方法是把一幅图像分成满足某种判据的区域，也就是说，把点组成区域。为了实现分组，要解决三个基本问题：,1,）、确定区域的数目；,2,）、要确定一个区域与其他区域相区别的特征；,3,）、确定一个产生有意义分割的相似性判据。,分割区域的一种方法叫区域生长或区域生成。假定区域的数目以及在每个区域中单个点的位置已知，则可推导一种算法。从一个已知点开始，加上与已知点相似的邻近点形成一个区域。这个相似性准则可以是灰度级、彩色、组织、梯度或其他特性。,相似性的测度可以由所确定的阈值来判定。它的方法是从满足检测准则的点开始，在各个方向上生长区域。当其邻近点满足检测准则就并入小块区域中，当新的点被合并后再用新的区域重复这一过程，直到没有可接受的邻近点时，生成过程终止。,一个简单的例子。,相似性准则：邻近点的灰度级与物体的平均灰,度级的差小于。,(a),是输入图像，生长点是,9,；,(b),是第一步接受的邻近点；,(c),是第二步接受的邻近点；,(d),生长点是。,图,7,9,区域生长简例,5,5,8,6,4,8,9,7,2,2,8,3,3,3,3,3,8,8,相似性准则是邻近点的灰度级与物体的平均灰度级的差小于,8,5,5,8,6,4,8,9,7,2,2,8,3,3,3,3,3,8,8,第二步：,8,7,5,5,8,6,4,8,9,7,2,2,8,3,3,3,3,3,6,另一种区域生长：,7,8,当生成任意物体时，接受准则可以以结构为基础，而不是以灰度级或对比度为基础。为了把候选的小群点包含在物体中，可以检测这些小群点，而不是检测单个点，如果它们的结构与物体的结构充分并且足够相似时就接受它们。,另外，还可以使用界线检测对生成建立“势垒”，如果在“势垒”的近邻点和物体之间有界线，则不能把这邻近点接受为物体中的点。,8.1.1,灰度阈值法分割,8.1.2,样板匹配,8.1.3,区域生长,8.1.4,区域聚合,区域聚合可直接用于图像分割。它要求聚合中的各个点必须在平面上相邻接而且特性相似。区域聚合的步骤是首先检查图像的测度集，以确定在测度空间中聚合的位置和数目，然后把这些聚合的定义用于图像，以得到区域聚合。一般区域聚合技术可以说明如下：,（,1,）,定义一个等价关系。例如，最简单的等价关系可定义为 。也就是说，如果 ，就说明 与 等价。任何在点的格子上的等价关系又可划分为等价类。例如 的取值范围为,0,到,63,，就可以产生,64,个等价类的模板。如果关系满足，它的值等于，否则为。,（,2,）确定邻接性，连接性可以用点 的邻点来定义。如连接邻点，连接邻点等等。连接邻点是四个非对角线上的个邻点，连接则是环绕的个邻点。,（,3,）通过这样的连接关系可以定义一个属于 的子集，这个子集形成一个区域。在这个区域中，任何点都与 有关。利用等价模板可分成最大的连接区域。然后，这些最大的连结区域又可以象搭积木一样形成有意义的分割。,年布赖斯和芬尼玛提出一种分割方法。这个方法如图,7,10,所示。,图中,(a),是具有灰度级的,3,3,的,G,阵列，图,(b),是对,S,的分割结果。其中图像格子为,G,，它是大格子,S,的子格子。,G,为,M,N,的 格子，,S,是,(2,M,+1),(2,N,+1),的大格子。,在大格子中，,G(i,j,),点位于,S,的,(2i+1),(2j+1),点上。,G,中的点与,S,中的点相对应，其中每一下标都是奇数，其余的点用来代表区域的边界。,以这种形式表现的区域，产生一种寻找最大连结区域的方法。,G,中的点与它上边和右边的点相比较，灰度级相同就合并，灰度级不同就插入边界线。把图像中的每个点都考虑了之后，整个图像就被分割成区域了。在这个例子中，由于采用了连接等价关系，因此，由图,7,10,可见，在对角线方向上的等灰度级产生了隔开的区域。,图,7,10,布赖斯和芬尼玛分割方法,