骨架提取与重建.docVIP

形态学骨架提取与重建摘要：本文首先介绍了三种骨架抽取算法, 并选取了最大圆盘的形态学骨架抽取算法作为研究对象,在c+编程的实验基础上，展示了基于形态学的图像骨架提取和图像重建的过程。关键词：数学形态学；骨架提取；骨架重建骨架是图象几何形态的一种重要拓扑描述。利用骨架表示原始图像，可以在保持图象重要拓扑特征的前提下, 减少图象中的冗余信息。因此, 被广泛应用于生物形状描述, 模式识别, 工业检测, 定量金 相以及图象压缩编码等领域。1骨架抽取的几种模型随着研究的深入, 人们提出了多种骨架化算法, 下面首先讨论一 下三种最常用的算法。 1.1中轴变换法: 中轴可以用下面的例子来形象的说明, 设想 t = 0 时刻, 将目标边界各处同时点燃, 火的前沿以匀速向目标蔓延, 当前沿相交时火焰熄灭, 火焰熄灭点的集合就构成了中轴。通常我们对目标图像的细化处理, 就是将图像上的文字、曲线、直线等几何元素的线条沿着中心轴线将其细化成一个象素宽的线条的处理过程。 1.2 影响区骨架法: 首先介绍一下测地距离和测地影响区的概念。测地距离以x , y为端点, 包含在A 中的最短距离。 在图1 中为欧氏距离, 为测地距离。图1 测地距离 图2 最大圆盘方式定义的中轴 测地影响区 ：设B 为A 中一个集合, 它由多个连通子集组成, 记为。连通成分 的测地影响区 , 由A 中距 的测地距离比距B 的其他连通成分距离都要小的所有点集组成。 (1)其中, k 为B 中连通成分的个数。由此得到影响区骨架( SK I Z ) 的定义。影响区骨架( S K I Z ) ：在A 中不属于任何一个B 的测地影响区的点集, 构成影响区骨架, 用表示: (2)1.3 基于最大圆盘的形态学骨架抽取算法B l um 与N ag l e 等人对区域的骨架作了很形象的描述。他们设想采用一个直径大小可任意改变的圆盘, 连接成分可以由一系列的最大圆盘来描述。最大圆盘与连接成分的轮廓相切, 连接成分的骨架则是这一系列最大内切圆盘圆心的连线。如图2 所示。以骨架上这些点为圆心的圆盘满足两个条件:(1)这些圆位于图像内部( 即被图像完全包容) ;(2)这些圆都是最大的,即在图像内部不可能再找到另一个圆能完全包容这个圆。通过以上的说明可以看到在求图像的骨架时应满足两个条件:在求骨架的过程中,图像应该有规律的缩小;在图像逐步缩小的过程中应使图像的欧拉数, 即连通性质保持不变。下面用数学形态学来描述求骨架的过程。对于n = 0, 1, 2, 定义骨架子集Skel( S ; n ) 为图像S 内所有最大圆盘nB的圆心 构成的集合。从骨架的定义可知, 骨架是所有骨架子集的并: (3) 定义1设为欧氏空间上的集合,的骨架记为 , 为骨架的子集, 即)对应于的最大内切圆半径为 的骨架, 则数学形态学对骨架的描述为 (4) 式中指半径为r的圆, 为具有微小半径值的圆。欧氏网格空间中, 圆可近似认为是一个小的结构元素 。圆则近似地认为是有离散半径的圆:( n 次) 。由此可以得出数学形态学对骨架的另一种定义形式。定义2设X为欧氏网格空间上的集合, X的骨架记为S ( X ) , 则数学形态学对骨架可描述为: (5) 式中为的第个骨架子集, 为满足和的n 值, 即至少被腐蚀次才会变为空图像。为个的膨胀, 即。由定义1可以得出计算骨架的迭代过程。从本质上讲, 定义1 与定义2 其实是一致的, 即定义2 中的每一个子集 中的所有点都是用最大圆概念描述骨架定义的方法中所指的大小为 的最大圆圆心。2数学形态学骨架抽取算法定义为包括图像及背景的积累集。 (1)(2)(3)If , Then and go to 7.(4)(5)(6) and go to 2.(7)(8) Return用数学形态学方法提取的骨架具有下列性质：(1) 若则(2) 若选时则(3) 骨架的各个子集之间是不连通的，即所得到的骨架为间断的.(4) 若为的第个子集中的元素，即，那么结构元中至少包含二个的边界点.3 骨架重建设二值图像的集合为,中边界的集合为,为骨架集合中的点,为到集合中最近点的距离,若中存在至少二点、;满足:;即P中至少有二个点x、y,m到、的距离相等且就是到集合中最近点的距离,则所有的的集合就是 的骨架.按照此定义,通常的用骨架重建原图像的方法,是以骨架上每一个点为圆心,以半径为作圆盘,所有圆盘的并集就是重建的原始图像.用数学形态方求出骨架后,也可用形态学变换的方法重建原始图像,这实际上是求骨架的逆运算.若设骨架函数为:对任一,若,则.这样,对每一骨架上的元素都对应一个骨架函数,且记录了该元素所在的是第个骨架子集,即:若则从重建原图像: (I)或 (I)中的骨架表达式表明,图像可由各个骨架子集做一次膨胀运算获得.也可用骨架上的各个元素与