《数字图像处理第09章数学形态学及其应用》由会员分享，可在线阅读，更多相关《数字图像处理第09章数学形态学及其应用（40页珍藏版）》请在金锄头文库上搜索。
1、上上上上一一一一页页页页下下下下一一一一页页页页home第第9章章 数学形态学及其应用数学形态学及其应用上上上上一一一一页页页页下下下下一一一一页页页页home9.1 形态学简介（1）运算由集合运算（如并、交、补等）来定义；（2）所有的图像都必须以合理的方式转换为集合。形态学算子的性能主要以几何方式进行刻画，更适合视觉信息的处理和分析。基本思想：基本思想：利利用用结结构构元元素素作作为为“探探针针”在在图图像像中中不不断断移移动动，在在此此过过程程中中收收集集图图像像的的信信息息、分分析析图图像像各各部部分分间间的的相相互互关关系，从而了解图像的结构特征。系，从而了解图像的结构特征。 上上上上
2、一一一一页页页页下下下下一一一一页页页页home图图9.1 9.1 数学形态学的方法数学形态学的方法上上上上一一一一页页页页下下下下一一一一页页页页home根据探测研究图像的不同结构特点，结构元素可携带形态、大小、灰度、色度等信息。不同点的集合形成具有不同性质的结构元素。由于不同的结构元素可以用来检测图像不同侧面的特征，因此设计符合人的视觉特性的结构元素是分析图像的重要步骤。最基本的形态学运算有：最基本的形态学运算有：腐蚀，膨胀，开启，闭合。用这些算子及其组合来进行图像形状和结构的分析及处理，可以解决抑制噪声、特征提取、边缘检测、形状识别、纹理分析、图像恢复与重建等方面的问题。上上上上一一一一
3、页页页页下下下下一一一一页页页页home数学形态学的特性：（1）反映的是一幅图像中像素点间的逻辑关系，而不是简单的数值关系。（2）是一种非线性的图像处理方法，并且具有不可逆性(后边具体解释)。（3）可以并行实现(串行和并行的区别)。（4）可以用来描述和定义图像的各种集合参数和特征。上上上上一一一一页页页页下下下下一一一一页页页页home9.1.2 几个基本概念几个基本概念1.击中与击不中击中与击不中设有两幅图像A和B，如果AB（空集），那么称B击中（hit）A，记为BA，；否则，如果AB=，称B击不中（miss）A。ABB上上上上一一一一页页页页下下下下一一一一页页页页home2.移动和反射移
4、动和反射设A是一幅数字图像，a是A的元素；b是一个点，那么定义A移动到b后的结果为 (A)b(a)b| aA （9.1）一幅数字图像A关于原点的反射定义为 AVa| -aA （9.2）上上上上一一一一页页页页下下下下一一一一页页页页home3. 结构元素结构元素根据不同的图像分析目的，常用的结构元素有方形、扁平形、圆形等。在多尺度形态学分析中，结构元素的大小可以变化，但结构元素的尺寸一般地要明显小于目标图像的尺寸。 上上上上一一一一页页页页下下下下一一一一页页页页home9.2 9.2 二值形态学二值形态学二值图像是数字图像的重要子集，指灰度值只取两种值的图像。两个灰度值可取为0和1。二值形态
5、学处理算法都是以膨胀,腐蚀这两种最基本的运算为基础的。一般设集合A为图像集合，集合B为结构元素，数学形态学运算是用B对A进行操作。二值二值分割：分割：上上上上一一一一页页页页下下下下一一一一页页页页home9.2.1 二值腐蚀二值腐蚀集合A（输入图像）被集合B（结构元素）腐蚀: (9.3) 图9.2 腐蚀示意图dd/8d/8dd/4d/4A AB B A AB Bd上上上上一一一一页页页页下下下下一一一一页页页页home9.2.2 二值膨胀二值膨胀腐蚀运算的对偶运算，可以直接定义，也可通过对补集的腐蚀来定义，即以AC表示集合A的补集， 表示B关于坐标原点的反射。集合集合A被集合被集合B膨胀膨胀
6、表示为： 也可表示为也可表示为上上上上一一一一页页页页下下下下一一一一页页页页home图9.3 膨胀示意图上上上上一一一一页页页页下下下下一一一一页页页页home腐蚀和膨胀操作的直观解释腐蚀和膨胀操作的直观解释 腐蚀是对图像内部作滤波处理，而膨胀是利用结构元素对图像补集进行填充，因而它是对图像外部作滤波处理。腐蚀具有收缩图像的作用，膨胀具有扩大图像的作用。上上上上一一一一页页页页下下下下一一一一页页页页home9.2.3 二值开运算二值开运算有两种二次运算起着非常重要的作用开运算闭运算（开运算的对偶运算） 。从结构元素填充的角度看，它们具有更为直观的几何形式。上上上上一一一一页页页页下下下下一
7、一一一页页页页home开运算的定义假设A仍为输入图像，B为结构元素，利用B对A作开运算，用符号AB表示，其定义为：上上上上一一一一页页页页下下下下一一一一页页页页home图图9.4 9.4 利用圆盘作开运算利用圆盘作开运算 上上上上一一一一页页页页下下下下一一一一页页页页home9.2.4 二值闭运算二值闭运算n闭运算是开运算的对偶运算，定义为先作膨胀后作腐蚀。u开、闭运算也互为对偶运算上上上上一一一一页页页页下下下下一一一一页页页页home图图9.5 9.5 利用圆盘作闭运算利用圆盘作闭运算上上上上一一一一页页页页下下下下一一一一页页页页home开运算具有磨光图像外边界的作用即去除毛刺；闭运
8、算具有磨光图像内边界的作用即填充破损。上上上上一一一一页页页页下下下下一一一一页页页页home二值形态学应用一、形态学滤波器二、图像骨架提取VC+MATLAB上上上上一一一一页页页页下下下下一一一一页页页页home9.3 灰值形态学灰值形态学在灰度图像形态处理中，输入和输出的图像都是灰度级形式的输入和输出像素值是在最低灰度值到最高灰度值之间。上上上上一一一一页页页页下下下下一一一一页页页页home9.3.1 灰值腐蚀灰值腐蚀形态学源于填充的概念灰值形态学处理的对象是图像信号波形的拓扑特性，结构元素也是一个信号。二值形态学中，集合的交、并运算起到关键作用在灰值形态学中这两种运算对应于极大和极小运
9、算。可利用填充、极大/极小概念直接定义灰值运算。上上上上一一一一页页页页下下下下一一一一页页页页home用结构元素g对输入图像f (x, y)进行灰值腐蚀记为 l用一维函数对其进行简化，定义为l由于结构元素必须在信号的下方，故空间平移结构元素的定义域必为信号定义域的子集，否则腐蚀在该点没有定义。l结构元素从信号的下面对信号产生滤波作用，这与结构元素从内部对二值图像滤波的情况是相似的。 上上上上一一一一页页页页下下下下一一一一页页页页home图9.6 灰值腐蚀运算上上上上一一一一页页页页下下下下一一一一页页页页home9.3.2 灰值膨胀灰值膨胀灰值膨胀可用灰值腐蚀的对偶运算来定义。采用求极大值
10、的方法，即在位于信号下方的条件下，求上推结构所能达到的最大值。利用结构元素的反射，求将信号限制在结构元素的定义域内时，上推结构元素使其超过信号时的最小值来定义。 上上上上一一一一页页页页下下下下一一一一页页页页home用结构元素g对输入图像f(x, y)进行灰值膨胀 (fg)(s, t)=maxf(s-x, t-y)+g(x, y)|s-x, t-y Df, x+yDg 用一维函数对其进行简化，定义为 (fg)(s)=maxf(s-x)+g(x)| s-x Df, xDg 步骤：步骤：对结构元素g的定义域Dg 中的每一个点x将信号f平移x，然后，再对每次平移信号的值加上g(x)，这样对于结构元
11、素定义域中的每个点都得到一个信号，对所有这些信号逐点取其最大值，便可得到膨胀结果。上上上上一一一一页页页页下下下下一一一一页页页页home图9.7 灰值膨胀运算 (a)灰值膨胀过程 (b)灰值膨胀结果上上上上一一一一页页页页下下下下一一一一页页页页home9.3.3 灰值开运算灰值开运算灰值开和灰值闭运算是腐蚀和膨胀的组合运算。先作腐蚀再作膨胀的迭代运算：这两种运算也为对偶运算，并且都可用填充概念来说明。 上上上上一一一一页页页页下下下下一一一一页页页页home图图9.8 9.8 灰值开运算灰值开运算 上上上上一一一一页页页页下下下下一一一一页页页页home9.3.4 灰值闭运算灰值闭运算根据
12、对偶性定义，灰值闭运算定义为： fg = (fg)g (9.14)灰值闭运算具有扩展性滤波结果总位于原始图像的上方。它从图像的上方磨光图像灰值表面向下突出的尖峰（即波谷）。 上上上上一一一一页页页页下下下下一一一一页页页页home图9.9 灰值闭运算上上上上一一一一页页页页下下下下一一一一页页页页home图9.10 细胞组织图像的灰值形态运算上上上上一一一一页页页页下下下下一一一一页页页页home9.3.59.3.5 灰值形态学梯度灰值形态学梯度利用扁平结构元素g对f作腐蚀和膨胀可得到f的局部极大和极小值，与数字差分定义的梯度相应。形态学梯度的定义为： u为了更好地获得边缘检测的效果，可以将形
13、态学梯度与阈值结合起来使用。 上上上上一一一一页页页页下下下下一一一一页页页页home图9.11 核磁共振图像的形态学梯度上上上上一一一一页页页页下下下下一一一一页页页页home9.3.69.3.6 高帽变换和低帽变换高帽变换和低帽变换通过这两种变换可以得到灰度图像中一些重要的标记点。在较亮的背景中求暗的像素点或在较暗的背景中求亮的像素点；检测受到噪声污染图像中的边缘等。为了使上述效果更明显，对变换后的图像也可以作阈值处理。上上上上一一一一页页页页下下下下一一一一页页页页home高帽变换从一幅原始图像f中减去对其作开运算后得到图像WHT(f)，其定义为： WHT(f) = f (fg) (9.
14、16)其中，g为结构元素。高帽变换是一种波峰检测器它在较暗的背景中求亮的像素点很有效。 上上上上一一一一页页页页下下下下一一一一页页页页home低帽变换与高帽变换相对偶的算子，定义为： BHT(f) = (fg) f (9.17)低帽变换是一种波谷检测器适合于在较亮的背景中求暗的像素点。上上上上一一一一页页页页下下下下一一一一页页页页home图9.12 利用Top-Hat和Bot-hat变换检测图像峰值和谷值上上上上一一一一页页页页下下下下一一一一页页页页home9.3.7 开闭运算和闭开运算开闭运算和闭开运算形态开、闭运算作为最基本的形态滤波运算。在实际的图像处理中，仅仅采用形态开和闭的滤波效果往往不能令人满意。在基本的形态开、闭运算的基础上设计出形态开闭和形态闭开组合滤波器，发挥其更好的滤波性能。上上上上一一一一页页页页下下下下一一一一页页页页home开闭运算和闭开运算互为对偶操作形态开闭运算能够滤除图像中形状小于结构元素的亮噪声。形态闭开运算可以滤除图像中形状小于结构元素的暗噪声。开闭运算： fg= fgg (9.18)闭开运算： fg = fgg (9.19)
