第八章数学形态学及其应用.doc

摘要： 第八章数学形态学及其应用8.1 引言8.2 二值形态学8.3 灰值形态学8.4 形态学的应用8.5 应用实例细化 8.1 引言8.1.1 数学形态学数学形态学(Mathematical Morphology)诞生于1964年，是由法国巴黎矿业学院博士生赛拉(J. Serra)和导师马瑟荣，在从事铁矿核的定量岩石学分析及预测其开采价值的研究中提出“击中/击不中变换”， 并在理论层面上第一次引入了形态学的表达式，建立了颗粒分析方法。他们的工作奠定了这门学科的理论基础， 如击中/击不中变换、开闭运算、布尔模型及纹理分析器的原型等。数学形态学的基本思想是用具有一定形态的结构元素去量度和提取图像中的对应形状以达到对图像分析和识别的目的。数学形态学的数学基础和所用语言是集合论，因此它具有完备的数学基础，这为形态学用于图像分析和处理、形态滤波器的特性分析和系统设计奠定了坚实的基础。数学形态学的应用可以简化图像数据，保持它们基本的形状特性，并除去不相干的结构。数学形态学的算法具有天然的并行实现的结构， 实现了形态学分析和处理算法的并行，大大提高了图像分析和处理的速度。数学形态学是由一组形态学的代数运算子组成的，它的基本运算有4个： 膨胀（或扩张）、腐蚀（或侵蚀）、开启和闭合， 它们在二值图像和灰度图像中各有特点。基于这些基本运算还可推导和组合成各种数学形态学实用算法，用它们可以进行图像形状和结构的分析及处理，包括图像分割、特征抽取、边界检测、图像滤波、图像增强和恢复等。数学形态学方法利用一个称作结构元素的“探针”收集图像的信息，当探针在图像中不断移动时便可考察图像各个部分之间的相互关系，从而了解图像的结构特征。数学形态学基于探测的思想，与人的FOA(Focus Of Attention) 的视觉特点有类似之处。作为探针的结构元素，可直接携带知识（形态、大小、甚至加入灰度和色度信息）来探测、研究图像的结构特点。数学形态学的基本思想及方法适用于与图像处理有关的各个方面，如基于击中/击不中变换的目标识别，基于流域概念的图像分割， 基于腐蚀和开运算的骨架抽取及图像编码压缩，基于测地距离的图像重建，基于形态学滤波器的颗粒分析等。迄今为止， 还没有一种方法能像数学形态学那样既有坚实的理论基础，简洁、朴素、统一的基本思想，又有如此广泛的实用价值。有人称数学形态学在理论上是严谨的，在基本观念上却是简单和优美的。数学形态学是一门建立在严格数学理论基础上的学科，其基本思想和方法对图像处理的理论和技术产生了重大影响。事实上数学形态学已经构成一种新的图像处理方法和理论，成为计算机数字图像处理的一个重要研究领域， 并且已经应用在多门学科的数字图像分析和处理的过程中。这门学科在计算机文字识别， 计算机显微图像分析(如定量金相分析， 颗粒分析)， 医学图像处理（例如细胞检测、心脏的运动过程研究、脊椎骨癌图像自动数量描述），图像编码压缩， 工业检测(如食品检验和印刷电路自动检测)，材料科学， 机器人视觉，汽车运动情况监测等方面都取得了非常成功的应用。另外，数学形态学在指纹检测、经济地理、合成音乐和断层X光照像等领域也有良好的应用前景。形态学方法已成为图像应用领域工程技术人员的必备工具。目前，有关数学形态学的技术和应用正在不断地研究和发展。 8.1.2 基本符号和术语1. 元素和集合在数字图像处理的数学形态学运算中，把一幅图像称为一个集合。对于二值图像而言，习惯上认为取值为1的点对应于景物中心，用阴影表示，而取值为0的点构成背景，用白色表示， 这类图像的集合是直接表示的。考虑所有值为1的点的集合为A， 则A与图像是一一对应的。对于一幅图像A，如果点a在A的区域以内， 那么就说a是A的元素，记为aA，否则，记作aA， 如图81（a）所示。图8-1 元素与集合间的关系a A b B A (a) (b) 2. 交集、并集和补集两个图像集合A和B的公共点组成的集合称为两个集合的交集记为AB，即AB=aaA且aB。两个集合A和B的所有元素组成的集合称为两个集合的并集，记为AB，即AB=a aA或aB。对一幅图像A，在图像A区域以外的所有点构成的集合称为A的补集，记为AC，即AC=aa A。交集、并集和补集运算是集合的最基本的运算，如图8-2所示。图8-2 集合的交集、并集和补集B A B A B A AB AB AC 3.击中（Hit）与击不中（Miss）设有两幅图像A和B，如果AB，那么称B击中A，记为BA， 其中是空集合的符号；否则，如果AB=， 那么称B击不中A， 如图8-3所示。 图8-3 击中与击不中（a） B击中A； （b） B击不中A (a) (b) A B B A 4.平移和反射设A是一幅数字图像（见图8-4(a)），b是一个点（见图8-4(b)），那么定义A被b平移后的结果为Abab| aA，即取出A中的每个点a的坐标值，将其与点b的坐标值相加，得到一个新的点的坐标值a+b，所有这些新点所构成的图像就是A被b平移的结果，记为A+b，如图8-4（c）所示。图8-4 平移与反射0 1 2 3 4 x y 12345 (a) 0 1 2 3 4 x 123 y (b) b 0 1 2 3 4 x y 12345 y 1234x 4 3 2 1 0 (c) (d) 5. 目标和结构元素被处理的图像称为目标图像，一般用大写英文字母表示。为了确定目标图像的结构，必须逐个考察图像各部分之间的关系，并且进行检验，最后得到一个各部分之间关系的集合。在考察目标图像各部分之间的关系时，需要设计一种收集信息的“探针”， 称为“结构元素”。“结构元素”一般用大写英文字母表示，例如用S表示。在图像中不断移动结构元素， 就可以考察图像之间各部分的关系。一般，结构元素的尺寸要明显小于目标图像的尺寸。 8.2 二值形态学二值形态学中的运算对象是集合。设A为图像集合，S为结构元素，数学形态学运算是用S对A进行操作。需要指出，实际上结构元素本身也是一个图像集合。对每个结构元素可以指定一个原点，它是结构元素参与形态学运算的参考点。应注意， 原点可以包含在结构元素中，也可以不包含在结构元素中，但运算的结果常不相同。以下用阴影代表值为1的区域，白色代表值为0的区域，运算是对值为1的区域进行的。二值形态学中两个最基本的运算腐蚀与膨胀，如图8-5所示。图8-5 腐蚀与膨胀示意图二值图像腐蚀膨胀 8.2.1 腐蚀腐蚀是最基本的一种数学形态学运算。对一个给定的目标图像X和一个结构元素S， 想象一下将S在图像上移动。在每一个当前位置x， S+x只有三种可能的状态（见图8-6）： (1) S+xX； (2) S+xXC； (3) S+xX与S+xXC均不为空。? 图8-6 S+x的三种可能的状态x Sx1 Sx2 Sx3 第一种情形说明S+x与X相关最大，第二种情形说明S+x与X 不相关，而第三种情形说明S+x与X只是部分相关。因而满足式(8-1)的点x的全体构成结构元素与图像最大相关点集，这个点集称为S对X的腐蚀(简称腐蚀，有时也称X用S腐蚀)，记为XS。腐蚀也可以用集合的方式定义，即式(8-2)表明，X用S腐蚀的结果是所有使S平移x后仍在X中的x的集合。换句话说，用S来腐蚀X得到的集合是S完全包括在X中时S 的原点位置的集合。上式也可以帮助我们借助相关概念来理解腐蚀操作。XS = x | S + x ? X （8-2）式(8-2)表明，X用S腐蚀的结果是所有使S平移x后仍在X中的x 的集合。换句话说，用S来腐蚀X得到的集合是S完全包括在X中时S的原点位置的集合。上式也可以帮助我们借助相关概念来理解腐蚀操作。腐蚀在数学形态学运算中的作用是消除物体边界点。如果结构元素取33的像素块，腐蚀将使物体的边界沿周边减少一个像素。腐蚀可以把小于结构元素的物体(毛刺、小凸起)去除，这样选取不同大小的结构元素，就可以在原图像中去掉不同大小的物体。如果两个物体之间有细小的连通，那么当结构元素足够大时通过腐蚀运算可以将两个物体分开。例8-1 腐蚀运算图解。图8-7给出腐蚀运算的一个简单示例。其中， 图8-7(a)中的阴影部分为集合X， 图8-7 (b)中的阴影部分为结构元素S， 而图(c)中黑色部分给出了XS 的结果。由图可见，腐蚀将图像（区域）收缩小了。图8-7 腐蚀运算示例(a) (b) (c) XS = S + X | ?sS (8-3) 如果S包含了原点，即OS， 那么XS将是X的一个收缩， 即XSX（当OS时）；如果S不包含原点，那么XSX未必成立。如果结构元素S关于原点O是对称的， 那么S=SV，因此X S=XSV, 但是，如果S关于原点O不是对称的，那么X被S腐蚀的结果与X被SV腐蚀的结果是不同的。利用腐蚀运算的定义式可以直接设计腐蚀变换的算法。但有时为了更方便，常使用腐蚀的另一种表达式，即? ? 式(8-3)可从定义式(8-2)中推出，它把腐蚀表示为图像平移的交，这在某些并行处理环境中特别有用。图8-8给出了一个例子。图8-8 腐蚀表示为图像平移的交O x S y X X(S) 根据上述理论，利用VC+可以编写一个实现腐蚀运算的函数MakeErosion()。/*/函数名称： BOOL MakeErosion() /基本功能： 本函数对图像数据执行腐蚀操作/参数说明： /int *nMask 结构元素数组指针/int nMaskLen 结构元素长度（以点数为计数单位） /unsigned char*pOut输出图像数据指针/unsigned char*pIn输入图像数据指针/intnWidthBytes图像宽度（以字节表示）/intnWidth图像宽度（以像素表示） /intnHeight图像高度（以像素表示）/返回值： BOOL 成功返回TRUE， 失败返回FALSE/*BOOL CMorphPro: : MakeErosion(int *nMask, int nMaskLen, unsigned char *pOut, unsigned char *pIn, int nWidthBytes, int nWidth, int nHeight) /若传入的图像数据为空， 将无法完成操作， 直接返回if(pOut = = NULL | pIn = = NULL) return FALSE; /定义变量int x, y, k; unsigned char Mark; /执行腐蚀操作 for( y = 0; y nHeight; y+ )unsigned char *pOutTemp = pOut; pOutT