图像表示与描述

图像表示与描述第一页，共九十五页，编辑于2023年，星期五图像描述实例：指纹图像描述（a）Gabor滤波器；（b）小波变换；（c）细节点（分叉点、端点）第二页，共九十五页，编辑于2023年，星期五主要内容6.1表示方法6.2边界描绘子6.3区域描绘子6.4运用主分量进行描述（自学）第三页，共九十五页，编辑于2023年，星期五6.1表示方法1链码

（1）概念链码是对边界点的一种编码表示方法，其基本思想是利用一系列具有特定长度和方向的相连的直线段来表示目标的边界。第四页，共九十五页，编辑于2023年，星期五

举例：若设起始点O的坐标为（5，5），则分别用如下4方向和8方向链码按逆时针顺序表示区域边界：4方向链码：（5，5）111232323000；8方向链码：（5，5）2224556000。第五页，共九十五页，编辑于2023年，星期五（2）链码表示的特点只有边界的起点需用绝对坐标表示，其余点都可只用接续方向来代表偏移量。与用坐标值相比，链码表达可大大减少边界表示所需的数据量。

第六页，共九十五页，编辑于2023年，星期五（3）存在的问题直接对分割所得的目标边界编码，有可能出现如下问题：A、产生的码串通常很长；

B、噪声等干扰会导致小的边界变化而使链码发生与目标整体形状无关的较大变动。（4）改进措施对原边界以较大的网格重新采样，并把与原边界点最接近的大网格点定为新的边界点。第七页，共九十五页，编辑于2023年，星期五总结重采样的规则。第八页，共九十五页，编辑于2023年，星期五（5）链码的起点归一化

A、问题的引出

链码起点的选择常是很关键的。对同一个边界，如用不同的边界点作为链码起点，得到的链码是不同的。

B、归一化给定1个从任意点开始而产生的链码，可把它看作1个由各方向数构成的自然数。将这些方向数依1个方向循环以使它们所构成的自然数的值最小。我们将这样转换后所对应的链码起点作为这个边界的归一化链码的起点。

第九页，共九十五页，编辑于2023年，星期五另一起点：33221010第十页，共九十五页，编辑于2023年，星期五（6）链码的旋转归一化A、问题用链码表示给定目标的边界时，如果目标旋转，则链码会发生变化。B、解决方法利用链码的一阶差分来重新构造1个序列(1个表示原链码各段之间方向变化的新序列)。这相当于把链码进行旋转归一化。第十一页，共九十五页，编辑于2023年，星期五原码旋转90度码链码10103322212100334方向差分：3313303033133030

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

A、算法步骤：1）沿着边界选两个相邻的点对，计算首尾连接直线段与原始折线段的误差R。2）如果误差R小于预先设置的阈值T。去掉中间点，选新点对与下一相邻点对，重复1）；否则，存储线段的参数，置误差为0，选被存储线段的终点为起点，重复1）2）。3）当程序的第一个起点被遇到，程序结束。第十六页，共九十五页，编辑于2023年，星期五RR<T第十七页，共九十五页，编辑于2023年，星期五第十八页，共九十五页，编辑于2023年，星期五B、聚合算法存在的问题：

顶点一般不对应于边界的拐点（如拐角）。因为新的线段直到超过误差的阈值才开始。第十九页，共九十五页，编辑于2023年，星期五（4）拆分技术

算法步骤：1）连接边界线段的两个端点（如果是封闭边界，连接最远点）；2）如果最大正交距离大于阈值，将边界分为两段，最大值点定位一个顶点。重复1）；3）如果没有超过阈值的正交距离，结束。第二十页，共九十五页，编辑于2023年，星期五第二十一页，共九十五页，编辑于2023年，星期五3标记

（1）基本思想标记是边界的1-D泛函表达，其基本思想是把2-D的边界用1-D的较易描述的函数形式来表达。（2）最简单的标记方法先对给定的物体求出质心，然后把边界点与质心的距离作为角度的函数就得到一种标记。第二十二页，共九十五页，编辑于2023年，星期五第二十三页，共九十五页，编辑于2023年，星期五（3）存在问题函数过分依赖于旋转和比例的变化。（4）改进措施-----旋转不变A、选择离质心最远的点作为起点；B、选择从质心到本征轴最远的点作为起点；C、使用差分链码的方法。（5）改进措施-----比例不变对函数进行正则化，使函数值总是分布在相同的值域里，比如说[0，1]。A、利用长短轴进行正则化；B、利用所有边界样本进行正则化。第二十四页，共九十五页，编辑于2023年，星期五4边界分段（1）基本概念A、一个任意集合S（区域）的凸起外缘H是：包含S的最小凸起的集合。B、H-S的差的集合被称为集合S的凸起补集D。SSDS+D=H第二十五页，共九十五页，编辑于2023年，星期五（2）分段算法：给进入和离开凸起补集D的变换点打标记来划分边界段。（3）优点：不依赖于方向和比例的变化S第二十六页，共九十五页，编辑于2023年，星期五（4）存在问题噪音的影响，导致出现零碎的划分。（5）改进措施先平滑边界，或用多边形逼近边界，然后再分段。第二十七页，共九十五页，编辑于2023年，星期五4骨架（1）基本思想表示一个平面区域结构形状的一种重要方法是把它削减成图形。这种削减可以通过细化（也称为抽骨架）算法，获取区域的骨架来实现。（2）Blum的中轴变换方法（MAT）设：R是一个区域，B为R的边界点，对于R中的点p，找p在B上“最近”的邻居。如果p有多于一个的邻居，称它属于R的中轴（骨架）第二十八页，共九十五页，编辑于2023年，星期五pRB（3）存在问题：计算量大第二十九页，共九十五页，编辑于2023年，星期五（4）算法改进思想在保证产生正确的骨架的同时，改进算法的效率。比较典型的是一类细化算法，它们不断删去边缘，但保证删除满足：A、不移去端点B、不破坏连通性C、不引起区域的过度腐蚀第三十页，共九十五页，编辑于2023年，星期五（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顺序，0-1转换的个数)(c)p2*p4*p6=0

（p2、p4、p6至少有一个0）(d)p4*p6*p8=0

（p4、p6、p8至少有一个0）p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第三十一页，共九十五页，编辑于2023年，星期五所有条件都满足，才打删除标记。删除并不立即进行，而是等到对所有边界点都打完标记后，再把作了标记的点一起删除举例: N(p1)=4

S(p1)=3 p2*p4*p6=0 p4*p6*p8=0第2个条件没满足不打标记00p1110101p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第三十二页，共九十五页，编辑于2023年，星期五B、基本操作2条件(a)、(b)与操作1相同，条件(c)、(d)改为：c’)p2*p4*p8=0d’)p2*p6*p8=0p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第三十三页，共九十五页，编辑于2023年，星期五随堂练习：（骨架抽取）第三十四页，共九十五页，编辑于2023年，星期五第三十五页，共九十五页，编辑于2023年，星期五第三十六页，共九十五页，编辑于2023年，星期五6.2边界描绘子1一些简单的描绘子

（1）边界的长度A、定义：区域的边界长度。

B、计算方法1)周长用边界所占面积表示，也即边界点数之和，每个点占面积为1的一个小方块。2)当把像素看作一个个点时，则周长用链码表示。此时，当链码值为奇数时，其长度记作 ；当链码值为偶数时，其长度记作1。即周长p表示为第三十七页，共九十五页，编辑于2023年，星期五随堂练习：（周长计算）第三十八页，共九十五页，编辑于2023年，星期五（2）边界的直径Diam(B)=max[D(pi,pj)]（3）边界线的离心率：长轴和短轴的比率。A、边界最大轴a：是连接距离最远的两个点的线段。B、边界最小轴b：与最大轴垂直，且其长度确定的包围盒刚好包围边界。C、基本矩形:包围边界的矩形。第三十九页，共九十五页，编辑于2023年，星期五边界最大轴a边界最小轴b基本矩形第四十页，共九十五页，编辑于2023年，星期五（4）曲率

定义为斜率的改变率，描述了边界上各点沿边界方向的变化量。用相邻边界线段（描述为直线）的斜率差作为在边界线交点处的曲率描述子。

ak1k2交点a处的曲率为

dk=k1–k2

第四十一页，共九十五页，编辑于2023年，星期五在一个边界点的曲率的符号描述了边界在该点的凹凸性。

P1P2第四十二页，共九十五页，编辑于2023年，星期五第四十三页，共九十五页，编辑于2023年，星期五2形状数（1）定义：形状数是链码的最小值的差分码。例如：基于4—方向的链码为:10103322,差分码为：33133030,形状数为：03033133。（2）形状数的阶(order)形状数序列的长度(即码的个数)。对闭合曲线，阶总是偶数。对凸形区域，阶对应边界外包矩形的周长。

第四十四页，共九十五页，编辑于2023年，星期五序号为4、6、8的形状数举例：

序号4链码：0321首差：3333形状：3333序号6链码：003221首差：303303形状：033033序号8链码：00032221首差：30033003形状：00330033第四十五页，共九十五页，编辑于2023年，星期五序号为6的形状数举例：

序号6链码：033211首差：330330形状：033033序号6链码：003221首差：303303形状：033033形状数与方向无关第四十六页，共九十五页，编辑于2023年，星期五序号为8的形状数举例：

序号8链码：03032211首差：33133030形状：03033133序号8链码：00332211首差：30303030形状：03030303序号8链码：00323211首差：30331330形状：03033133第四十七页，共九十五页，编辑于2023年，星期五（3）存在问题 虽然链码的首差是不依赖于旋转的，但一般情况下边界的编码依赖于网格的方向。（4）改进措施（规整化网格方向） 大多数情况下，将链码网格与基本矩形对齐，即可得到一个唯一的形状数。规整化网格方向的一种算法如下：A、首先确定形状数的序号n；B、在序号为n的矩形形状数中，找出一个与给定形状的基本矩形的离心率最接近的形状数的矩形。第四十八页，共九十五页，编辑于2023年，星期五C、然后再用这个矩形与基本矩形对齐，构造网格。D、用获得链码的方法得到链码；E、再得到循环首差；F、首差中的最小循环数即为形状数。例如:如果n=12，所有序号为12的矩形（即周长为12）为2*4，3*3，1*5。如果2*4矩形的离心率最接近于给定边界的基本矩形的离心率，我们建立一个2*4的网格。第四十九页，共九十五页，编辑于2023年，星期五链码：000033222121首差：300030300313形状：0003030031330123第五十页，共九十五页，编辑于2023年，星期五第五十一页，共九十五页，编辑于2023年，星期五3傅里叶描绘子（1）基本方法：A、将XY平面中的曲线段转化为复平面上的1个序列，从而用复数的形式来表示给定边界上每个点（x,y）。对1个由N个点组成的封闭边界，从任一点开始绕边界1周就得到1个复数序列：

s(k)=x(k)+jy(k)

y0y1xx(k)=xky(k)=yk第五十二页，共九十五页，编辑于2023年，星期五B、进行离散傅立叶变换

N-1a(u)=1/N∑s(k)exp(-j2uk/N)u=0,1,…,N-1

u=0

系数a(u)被称为边界的傅立叶描述子第五十三页，共九十五页，编辑于2023年，星期五C、选取整数MN-1，由a(0)，a(1)，a(M-1)描述二维封闭边界。这时，对应于边界的点数没有改变，但在重构每一个点所需要的计算项大大减少了。如果边界点数很大，M一般选为2的指数次方的整数。第五十四页，共九十五页，编辑于2023年，星期五（2）M的选取与描述符的关系在上述方法中，相当于对于u>M-1的部分舍去不予计算。由于傅立叶变换中高频部分对应于图像的细节描述，因此M取得越小，细节部分丢失得越多。进行逆傅立叶变换（重构）

M-1s’(k)=∑a(u)exp(j2uk/N) k=0,1,…,N-1

u=0M=4M=61M=62N=64思考题：如何进行空间曲线平滑？第五十五页，共九十五页，编辑于2023年，星期五（3）使用价值A、较少的傅立叶描述子（如4个），就可以获取边界本质的整体轮廓；B、这些带有边界信息的描述子，可以用来区分明显不同的边界。（4）优点A、使用复数作为描述符，对于旋转、平移、放缩等操作和起始点的选取不十分敏感。B、几何变换的描述子可通过对函数作简单变换来获得。几何变换傅立叶描述子原形a(u)旋转a(u)=a(u)ej平移a(u)=a(u)+xy(u)放缩a(u)=a(u)起点a(u)=a(u)e-j2k0u/N第五十六页，共九十五页，编辑于2023年，星期五4统计矩（1）基本思想：将描述形状的任务减少至描述一个一维函数，边界段和特征的形状可以用矩量来量化地描述。（2）统计矩的定义A、把边界当作直方图函数：g(r)第五十七页，共九十五页，编辑于2023年，星期五

B、定义：

K-1

n(r)=∑(ri-m)ng(ri)

i=0K-1

其中m=∑rig(ri)

i=1这里K是边界上点的数目,n(r)是边界的矩量第五十八页，共九十五页，编辑于2023年，星期五（3）矩量的优点A、实现是直接的；B、附带了一种关于边界形状的“物理”解释C、对于旋转的不敏感性D、为了使大小比例不敏感，可以通过伸缩r的范围来将大小正则化。

第五十九页，共九十五页，编辑于2023年，星期五6.3区域描绘子1一些简单的描绘子

（1）区域面积：区域中像素的数目。

（2）区域周长：区域边界的长度。第六十页，共九十五页，编辑于2023年，星期五

（3）致密度：(周长)2/面积。

（4）其他简单描绘子如最大值、最小值、中值、均值、重心、方差等。第六十一页，共九十五页，编辑于2023年，星期五实例：利用面积描绘子从图像中提取信息第六十二页，共九十五页，编辑于2023年，星期五2拓扑描绘子

（1）拓扑性质研究一种图像在没有撕裂和连接的情况下（橡皮伸展变形），不受任何变形影响的性质。（2）孔洞数H、连通分量的数目C、欧拉数EE=C-H第六十三页，共九十五页，编辑于2023年，星期五（3）拓扑网络与欧拉数的关系V-Q+F=C-H=E（V顶点数、Q边数、F面数）第六十四页，共九十五页，编辑于2023年，星期五在分割后的图像中提取最大连通分量第六十五页，共九十五页，编辑于2023年，星期五3纹理

（1）基础A、定义纹理是由许多相互接近的、互相编织的元素构成，它们富有周期性。第六十六页，共九十五页，编辑于2023年，星期五第六十七页，共九十五页，编辑于2023年，星期五B、纹理描绘子分类a)统计方法。（采用统计方法对纹理进行分析。）b)结构方法。（从图像结构的观点出发，则认为纹理是结构，由纹理基元按一定规律排列而成。采用句法结构方法。

）c)频谱法。（分析纹理的频域特征）第六十八页，共九十五页，编辑于2023年，星期五（2）统计方法A、基于直方图的统计矩第六十九页，共九十五页，编辑于2023年，星期五B、其他的基于直方图的纹理量度第七十页，共九十五页，编辑于2023年，星期五实例第七十一页，共九十五页，编辑于2023年，星期五课堂练习：计算图像的均值、方差、平滑度、三次矩、一致度、平均熵。第七十二页，共九十五页，编辑于2023年，星期五

存在问题：没有利用像素之间的相对位置关系。第七十三页，共九十五页，编辑于2023年，星期五C、灰度共生矩阵

基本方法：取图像中任意一点(x，y)及偏离它的另一点(x+a,y+b)，设该点对的灰度值为(g1，g2)。对于整幅图像，统计出每—种(g1，g2)值出现的概率p(g1，g2)，并排列成方阵，称为联合概率矩阵，也叫做共生矩阵。再由共生矩阵计算五个统计量。

具体步骤：a)由原始图像生成满足位置算子P的点对矩阵；0001211011221001102000100P：在右下方的一个像素第七十四页，共九十五页，编辑于2023年，星期五b)计算发生的概率（点对数/图像中满足P的点对总数）；第七十五页，共九十五页，编辑于2023年，星期五c)计算下列统计量，以描述灰度共生矩阵的“内容”特性。第七十六页，共九十五页，编辑于2023年，星期五0001211011221001102000100随堂练习：P定义为“在右边的一个像素”第七十七页，共九十五页，编辑于2023年，星期五分析：由于上图象中的灰度沿水平方向和垂直方向均有较高频率的变化，所以其共生矩阵图中大部分项均不为零。下图中的灰度在较大范围内变化缓慢，所以其共生矩阵图中仅有主对角线上的元素取较大的值。

第七十八页，共九十五页，编辑于2023年，星期五纹理特征匹配举例：从1万张图片中检索的结果第七十九页，共九十五页，编辑于2023年，星期五（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)aaabbbaabbab第八十页，共九十五页，编辑于2023年，星期五第八十一页，共九十五页，编辑于2023年，星期五（3）频谱方法A、纹理描述中常用的三个性质a)傅立叶频谱中突起的峰值对应纹理模式的显现度;b)这些峰在频域平面的位置对应模式的基本周期和方向;c)如果利用滤波把周期性成分除去，剩下的非周期性部分可用统计方法描述。第八十二页，共九十五页，编辑于2023年，星期五