第六章图像表示与描述.ppt

1、第6章 表示与描述,图像描述实例：指纹图像描述,（a）Gabor滤波器； （b）小波变换； （c）细节点 （分叉点、端点）,主要内容 6.1 表示方法 6.2 边界描绘子 6.3 区域描绘子 6.4 运用主分量进行描述（自学）,6.1 表示方法 1 链码 （1）概念 链码是对边界点的一种编码表示方法，其基本思想是利用一系列具有特定长度和方向的相连的直线段来表示目标的边界。,举例：若设起始点O的坐标为（5，5），则分别用如下4方向和8方向链码按逆时针顺序表示区域边界： 4方向链码： （5， 5）1 1 1 2 3 2 3 2 3 0 0 0； 8方向链码： （5， 5）2 2 2 4 5 5 6

2、 0 0 0。,（2）链码表示的特点 只有边界的起点需用绝对坐标表示，其余点都可只用接续方向来代表偏移量。与用坐标值相比，链码表达可大大减少边界表示所需的数据量。,（3）存在的问题 直接对分割所得的目标边界编码，有可能出现如下问题： A、产生的码串通常很长； B、噪声等干扰会导致小的边界变化而使链码发生与目标整体形状无关的较大变动。 （4）改进措施 对原边界以较大的网格重新采样，并把与原边界点最接近的大网格点定为新的边界点。,总结重采样的规则。,（5）链码的起点归一化 A、问题的引出 链码起点的选择常是很关键的。对同一个边界，如用不同的边界点作为链码起点，得到的链码是不同的。 B、归一化 给定

3、1个从任意点开始而产生的链码，可把它看作1个由各方向数构成的自然数。将这些方向数依1个方向循环以使它们所构成的自然数的值最小。我们将这样转换后所对应的链码起点作为这个边界的归一化链码的起点。,另一起点：33221010,（6） 链码的旋转归一化 A、问题 用链码表示给定目标的边界时，如果目标旋转，则链码会发生变化。 B、解决方法 利用链码的一阶差分来重新构造1个序列(1个表示原链码各段之间方向变化的新序列)。这相当于把链码进行旋转归一化。,原码 旋转90度码 链码 10103322 21210033 4方向差分： 33133030 33133030,练习： 分别写出下面目标的4连接和8连接链码

4、、最小值链码和差分链码。分别以（1，2）、（2，1 ）点为起点，顺时针进行。,2 多边形近似,（1）问题的引出 实际应用中的数字边界常由于噪声、采样等的影响而有许多较小的不规则处，这些不规则处常对链码和边界段表达产生较明显的干扰影响。 （2）多边形方法的基本思想 多边形是一系列线段的封闭集合，它可用来逼近大多数使用的曲线到任意的精度。 在实际中多边形表达的目的是要用尽量少的线段来代表边界并保持边界的基本形状，从而用较简单的形式来表达和描述边界。,（3）基于收缩的最小周长多边形法 将边界看成是有弹性的线，将组成边界的像素系列的内外边各看成一堵墙，如将线拉紧则可到最小周长多边形。,（4）聚合技术,

5、A、算法步骤： 1）沿着边界选两个相邻的点对，计算首尾连接直线段与原始折线段的误差R。 2）如果误差R小于预先设置的阈值T。去掉中间点，选新点对与下一相邻点对，重复1）；否则，存储线段的参数，置误差为0，选被存储线段的终点为起点，重复1）2）。 3）当程序的第一个起点被遇到，程序结束。,R,R T,B、聚合算法存在的问题： 顶点一般不对应于边界的拐点（如拐角）。因为新的线段直到超过误差的阈值才开始。,（4）拆分技术,算法步骤： 1）连接边界线段的两个端点（如果是封闭边界，连接最远点）； 2）如果最大正交距离大于阈值，将边界分为两段，最大值点定位一个顶点。重复1）； 3）如果没有超过阈值的正交距

6、离，结束。,3 标记,（1）基本思想 标记是边界的1-D泛函表达，其基本思想是把2-D的边界用1-D的较易描述的函数形式来表达。 （2）最简单的标记方法 先对给定的物体求出质心，然后把边界点与质心的距离作为角度的函数就得到一种标记。,（3）存在问题 函数过分依赖于旋转和比例的变化。 （4）改进措施-旋转不变 A、选择离质心最远的点作为起点； B、选择从质心到本征轴最远的点作为起点； C、使用差分链码的方法。,（5）改进措施-比例不变 对函数进行正则化，使函数值总是分布在相同的值域里，比如说0，1。 A、利用长短轴进行正则化； B、利用所有边界样本进行正则化。,4 边界分段 （1）基本概念 A、

7、一个任意集合S（区域）的凸起外缘H是：包含S 的最小凸起的集合。 B、H-S的差的集合被称为集合S的凸起补集D。,S,S,D,S + D = H,（2）分段算法： 给进入和离开凸起补集D的变换点打标记来划分边界段。 （3）优点：不依赖于方向和比例的变化,S,（4）存在问题 噪音的影响，导致出现零碎的划分。 （5）改进措施 先平滑边界，或用多边形逼近边界，然后再分段。,5 骨架,（1）基本思想 表示一个平面区域结构形状的一种重要方法是把它削减成图形。这种削减可以通过细化（也称为抽骨架）算法，获取区域的骨架来实现。 （2）Blum的中轴变换方法（MAT） 设：R是一个区域，B为R的边界点，对于R中

8、的点p，找p在B上“最近”的邻居。如果p有多于一个的邻居，称它属于R的中轴（骨架）,p,R,B,（3）存在问题：计算量大,（4）算法改进思想 在保证产生正确的骨架的同时，改进算法的效率。比较典型的是一类细化算法，它们不断删去边缘，但保证删除满足： A、不移去端点 B、不破坏连通性 C、不引起区域的过度腐蚀,（5）一种细化二值区域的算法 假设区域内的点值为1，背景值为0。由两个基本操作组成。,A、基本操作1 对于满足以下四个条件的边界点打标记准备删除： (a) 2N(p1)6 （ N(p1)=p2+p3+p9，是点p1邻域中1的个数） (b) S(p1) = 1 (S(p1)是按p2,p3,p9

9、，p2顺序，0-1转换的个数) (c) p2 * p4 * p6 = 0 （p2 、p4 、p6 至少有一个0） (d) p4 * p6 * p8 = 0 （p4 、p6 、p8 至少有一个0）,p9,p2,p1,p8,p3,p4,p7,p6,p5,p9,p2,p1,p8,p3,p4,p7,p6,p5,p9,p2,p1,p8,p3,p4,p7,p6,p5,所有条件都满足，才打删除标记。删除并不立即进行，而是等到对所有边界点都打完标记后，再把作了标记的点一起删除 举例: N(p1) = 4 S(p1) = 3 p2*p4*p6 = 0 p4*p6*p8 = 0 第2个条件没满足不打标记,0,0,

10、p1,1,1,0,1,0,1,p9,p2,p1,p8,p3,p4,p7,p6,p5,p9,p2,p1,p8,p3,p4,p7,p6,p5,B、基本操作2 条件(a)、(b)与操作1相同，条件(c)、(d)改为： c) p2* p4* p8= 0 d) p2* p6* p8= 0,p9,p2,p1,p8,p3,p4,p7,p6,p5,p9,p2,p1,p8,p3,p4,p7,p6,p5,随堂练习：（骨架抽取）,6.2 边界描绘子 1 一些简单的描绘子 （1）边界的长度,A、定义：区域的边界长度。 B、计算方法 1)周长用边界所占面积表示， 也即边界点数之和， 每个点占面积为1的一个小方块。,2)

11、 当把像素看作一个个点时，则周长用链码表示。此时，当链码值为奇数时，其长度记作 ； 当链码值为偶数时，其长度记作1。即周长p表示为,随堂练习：（周长计算）,（2）边界的直径,Diam(B) = maxD(pi, pj),（3）边界线的离心率：长轴和短轴的比率。,A、边界最大轴a：是连接距离最远的两个点的线段。 B、边界最小轴b：与最大轴垂直，且其长度确定的包 围盒刚好包围边界。 C、基本矩形: 包围边界的矩形。,（4）曲率 定义为斜率的改变率，描述了边界上各点沿边界方向的变化量。用相邻边界线段（描述为直线）的斜率差作为在边界线交点处的曲率描述子。,a,k1,k2,交点a处的曲率为 dk = k

12、1 k2,在一个边界点的曲率的符号描述了边界在该点的凹凸性。,P1,P2,2 形状数,（1）定义：形状数是链码的最小值的差分码。,例如：基于4方向的链码为:10103322,差分码为：33133030,形状数为：03033133。 （2）形状数的阶(order) 形状数序列的长度(即码的个数)。对闭合曲线，阶总是偶数。对凸形区域，阶对应边界外包矩形的周长。,序号为4、6、8的形状数举例：,序号4,链码：0321 首差：3333 形状：3333,序号6,链码：003221 首差：303303 形状：033033,序号8,链码：00032221 首差：30033003 形状：00330033,序号

13、为6的形状数举例：,序号6,链码：033211 首差：330330 形状：033033,序号6,链码：003221 首差：303303 形状：033033,形状数与方向无关,序号为8的形状数举例：,序号8,链码：03032211 首差：33133030 形状：03033133,序号8,链码：00332211 首差：30303030 形状：03030303,序号8,链码：00323211 首差：30331330 形状：03033133,（3）存在问题 虽然链码的首差是不依赖于旋转的，但一般情况下边界的编码依赖于网格的方向。 （4）改进措施（规整化网格方向）,大多数情况下，将链码网格与基本矩形对齐

14、，即可得到一个唯一的形状数。规整化网格方向的一种算法如下 ： A、首先确定形状数的序号n； B、在序号为n的矩形形状数中，找出一个与给定形状的基本矩形的离心率最接近的形状数的矩形。,C、然后再用这个矩形与基本矩形对齐，构造网格。 D、用获得链码的方法得到链码； E、再得到循环首差； F、首差中的最小循环数即为形状数。 例如: 如果n=12，所有序号为12的矩形（即周长为12）为2*4，3*3，1*5。如果2*4矩形的离心率最接近于给定边界的基本矩形的离心率，我们建立一个2*4的网格。,链码：000033222121 首差：300030300313 形状：000303003133,0,1,2,3

15、,3 傅里叶描绘子,（1）基本方法： A、 将XY平面中的曲线段转化为复平面上的1个序列，从而用复数的形式来表示给定边界上每个点（x,y）。对1个由N个点组成的封闭边界，从任一点开始绕边界1周就得到1个复数序列： s(k) = x(k) + jy(k),y0,y1,x,x(k) = xk y(k) = yk,B、进行离散傅立叶变换 N-1 a(u) =1/N s(k)exp(-j2uk/N) u=0,1,N-1 u=0 系数a(u)被称为边界的傅立叶描述子,C、选取整数 MN-1，由a(0)，a(1)，a(M-1)描述二维封闭边界。 这时，对应于边界的点数没有改变，但在重构每一个点所需要的计算

16、项大大减少了。如果边界点数很大，M一般选为2的指数次方的整数。,（2）M的选取与描述符的关系 在上述方法中，相当于对于u M-1的部分舍去不予计算。由于傅立叶变换中高频部分对应于图像的细节描述，因此M取得越小，细节部分丢失得越多。 进行逆傅立叶变换（重构） M-1 s(k) = a(u)exp(j2uk/N) k=0,1,N-1 u=0,M=4,M=61,M=62,N=64,思考题：如何进行空间曲线平滑？,（3）使用价值 A、较少的傅立叶描述子（如4个），就可以获取边界本质的整体轮廓； B、这些带有边界信息的描述子，可以用来区分明显不同的边界。,（4）优点 A、使用复数作为描述符，对于旋转、平

17、移、放缩等操作和起始点的选取不十分敏感。 B、几何变换的描述子可通过对函数作简单变换来获得。,4 统计矩,（1）基本思想： 将描述形状的任务减少至描述一个一维函数，边界段和特征的形状可以用矩量来量化地描述。 （2）统计矩的定义 A、把边界当作直方图函数：g(r),B、定义： K-1 n(r) = (ri- m)n g(ri) i=0 K-1 其中 m = ri g(ri) i=1 这里K是边界上点的数目, n(r)是边界的矩量,（3）矩量的优点 A、实现是直接的； B、附带了一种关于边界形状的“物理”解释 C、对于旋转的不敏感性 D、为了使大小比例不敏感，可以通过伸缩r的范围来将大小正则化。,

18、6.3 区域描绘子 1 一些简单的描绘子 （1）区域面积：区域中像素的数目。,（2）区域周长：区域边界的长度。,（3）致密度：(周长)2/面积 。,（4）其他简单描绘子 如最大值、最小值、中值、均值、重心、方差等。,实例：利用面积描绘子从 图像中提取信息,2 拓扑描绘子 （1）拓扑性质 研究一种图像在没有撕裂和连接的情况下（橡皮伸展变形），不受任何变形影响的性质。 （2）孔洞数H、连通分量的数目C、欧拉数E E=C-H,（3）拓扑网络与欧拉数的关系 V-Q+F=C-H=E （V顶点数、Q边数、F面数）,在分割后的图像中提取最大连通分量,3 纹理 （1）基础 A、定义 纹理是由许多相互接近的、

19、互相编织的元素构成， 它们富有周期性。,B、纹理描绘子分类,a) 统计方法。（采用统计方法对纹理进行分析。） b) 结构方法。（从图像结构的观点出发，则认为纹理是结构，由纹理基元按一定规律排列而成。 采用句法结构方法。 ） c) 频谱法。（分析纹理的频域特征）,（2）统计方法 A、基于直方图的统计矩,B、其他的基于直方图的纹理量度,实例,课堂练习： 计算图像的均值、方差、平滑度、三次矩、一致度、平均熵。,存在问题： 没有利用像素之间的相对位置关系。,C、灰度共生矩阵,基本方法：取图像中任意一点(x，y)及偏离它的另一点(x+a, y+b)，设该点对的灰度值为(g1，g2)。对于整幅图像，统计出

20、每种(g1，g2)值出现的概率p(g1，g2)，并排列成方阵，称为联合概率矩阵，也叫做共生矩阵。再由共生矩阵计算五个统计量。 具体步骤： a) 由原始图像生成满足位置算子P的点对矩阵；,0 0 0 1 2 1 1 0 1 1 2 2 1 0 0 1 1 0 2 0 0 0 1 0 1,P：在右下方的一个像素,b) 计算发生的概率（点对数/图像中满足P的点对总数）； 灰度共生矩阵,c) 计算下列统计量，以描述灰度共生矩阵的“内容”特性。,0 0 0 1 2 1 1 0 1 1 2 2 1 0 0 1 1 0 2 0 0 0 1 0 0,随堂练习：P定义为“在右边的一个像素”,分析：由于上图象中的

21、灰度沿水平方向和垂直方向均有较高频率的变化，所以其共生矩阵图中大部分项均不为零。下图中的灰度在较大范围内变化缓慢，所以其共生矩阵图中仅有主对角线上的元素取较大的值。,纹理特征匹配举例：从1万张图片中检索的结果,（2）结构性方法,A、基本思想 图像中各个部分间的结构关系是二维的，而串是一维的，期望找到一种方法把二维关系转化为一维的串。 B、通过产生规则来生成结构,阶梯结构关系 1) S-aA 2) A-bS 3) A-b 其中S、A是变量,(1,3),(1,2,1,3),(1,2,1,2,1,3),a,a,a,b,b,b,a,a,b,b,a,b,（3）频谱方法,A、纹理描述中常用的三个性质 a) 傅立叶频谱中突起的峰值对应纹理模式的显现度; b) 这些峰在频域平面的位置对应模式的基本周期和方向; c) 如果利用滤波把周期性成分除去， 剩下的非周期性部分可用统计方法描述。,B、频谱的整体性描述,