第七章图像分析基础

第七章图像分析基础第1页，课件共164页，创作于2023年2月9/7/20231本章主要内容：1.图像分割2.图像描述3.图像分析第2页，课件共164页，创作于2023年2月9/7/202321.图像分割图像分割的概念把图像分解成构成它的部件和对象的过程有选择性地定位感兴趣对象在图像中的位置和范围第3页，课件共164页，创作于2023年2月9/7/202331.图像分割图像分割的基本思路从简到难，逐级分割控制背景环境，降低分割难度把焦点放在增强感兴趣对象，缩小不相干图像成分的干扰上第4页，课件共164页，创作于2023年2月9/7/202341.图像分割图像分割的基本策略分割算法基于灰度值的两个基本特性：不连续性和相似性检测图像像素灰度级的不连续性，找到点、线（宽度为1）、边（不定宽度）。先找边，后确定区域。第5页，课件共164页，创作于2023年2月9/7/202351.图像分割图像分割的基本策略检测图像像素的灰度值的相似性，通过选择阈值，找到灰度值相似的区域，区域的外轮廓就是对象的边第6页，课件共164页，创作于2023年2月9/7/202361.图像分割1.1边界分割法点的检测线的检测边的检测第7页，课件共164页，创作于2023年2月9/7/202371.图像分割1.1边界分割法点的检测用空域的高通滤波器来检测孤立点例：R=(-1*8*8+128*8)/9=(120*8)/9=960/9=106设：阈值：T=6488881288888-1-1-1-18-1-1-1-1图像模板第8页，课件共164页，创作于2023年2月9/7/202381.图像分割1.1边界分割法点的检测如果R的值等于0，说明当前检测点与的灰度值与周围点的相同当R的值足够大时，说明该点的值与周围的点非常不同，是孤立点。通过阈值T来控制如T=32、64、128等

|R|>

T便检测到一个孤立点第9页，课件共164页，创作于2023年2月9/7/202391.图像分割1.1边界分割法-1-1-1222-1-1-1-1-12-12-12-1-1-12-1-12-1-12-12-1-1-12-1-1-12线的检测通过比较典型模板的计算值，确定一个点是否在某个方向的线上第10页，课件共164页，创作于2023年2月9/7/2023101.图像分割1.1边界分割法111555111111555111111555111线的检测R1=-6+30=24R2=-14+14=0R3=-14+14=0

R4=-14+14=0第11页，课件共164页，创作于2023年2月9/7/2023111.图像分割1.1边界分割法边的检测边界的定义：一段边是两个具有相对不同灰度值特性的区域的边界线适用于：假定问题中的区域是非常类似的，两个区域之间的过渡，仅仅根据灰度的不连续性便可确定不适用于：当假定不成立时，阈值分割技术一般来说比边缘检测更加实用第12页，课件共164页，创作于2023年2月9/7/2023121.图像分割1.1边界分割法边的检测分割对象分割对象第13页，课件共164页，创作于2023年2月9/7/2023131.图像分割1.1边界分割法边的检测基本思想：计算局部微分算子一阶微分二阶微分边界图像截面图第14页，课件共164页，创作于2023年2月9/7/2023141.图像分割1.1边界分割法边的检测一阶微分：用梯度算子来计算特点：对于亮的边，边的变化起点是正的，结束是负的。对于暗边，结论相反。常数部分为零。用途：用于检测图像中边的存在第15页，课件共164页，创作于2023年2月9/7/2023151.图像分割1.1边界分割法边的检测二阶微分：通过拉普拉斯来计算特点：二阶微分在亮的一边是正的，在暗的一边是负的。常数部分为零。用途：1）二次导数的符号，用于确定边上的像素是在亮的一边，还是暗的一边。2）0跨越，确定边的准确位置第16页，课件共164页，创作于2023年2月9/7/2023161.图像分割1.1边界分割法边的检测梯度算子函数f(x,y)在(x,y)处的梯度为一个向量： f=[f/x,f/y]计算这个向量的大小为：f=mag(f)=[(f/x)2+(f/y)2]1/2近似为:

f|x|+|y|z2z8z5z3z9z6z1z7z4第17页，课件共164页，创作于2023年2月9/7/2023171.图像分割1.1边界分割法边的检测梯度算子梯度的方向角为： (x,y)=tan(y/x)Sobel算子为：

x=(z7+2z8+z9)

-(z1+2z2+z3)y=(z3+2z6+z9)

-(z1+2z4+z7)-220-110-110000-1-1-2112xy第18页，课件共164页，创作于2023年2月9/7/2023181.图像分割1.1边界分割法边的检测Sobel梯度算子的使用与分析 1.直接计算y、x可以检测到边的存在， 以及从暗到亮，从亮到暗的变化2.仅计算|x|，产生最强的响应是正交 于x轴的边；|y|则是正交于y轴的边。

3.由于微分增强了噪音，平滑效果是Sobel 算子特别引人注意的特性第19页，课件共164页，创作于2023年2月9/7/2023191.图像分割1.1边界分割法边的检测拉普拉斯二维函数f(x,y)的拉普拉斯是一个二阶的微分定义为：

2f=[2f/x2,2f/y2]可以用多种方式被表示为数字形式。对于一个3x3的区域，经验上被推荐最多的形式是：

2f

=4z5–(z2+z4+z6+z8)z2z8z5z3z9z6z1z7z4第20页，课件共164页，创作于2023年2月9/7/2023201.图像分割1.1边界分割法边的检测拉普拉斯定义数字形式的拉普拉斯的基本要求是，作用于中心像素的系数是一个正数，而且其周围像素的系数为负数，系数之和必为0。-1-1400-100-1第21页，课件共164页，创作于2023年2月9/7/2023211.图像分割1.1边界分割法边的检测拉普拉斯算子的分析：缺点：对噪音的敏感；会产生双边效果； 不能检测出边的方向应用：拉普拉斯算子不直接用于边的检测， 通常只起第二位的角色；检测一个像素是在边的亮的一边还是暗的一边利用零跨越，确定边的位置第22页，课件共164页，创作于2023年2月9/7/2023221.图像分割1.2边缘连接法局部处理法Hough变换第23页，课件共164页，创作于2023年2月9/7/2023231.图像分割1.2边缘连接法边缘连接法由于噪音的原因，边界的特征很少能够被完整地描述，在亮度不一致的地方会中断。因此典型的边检测算法后面总要跟随着连接过程和其它边界检测过程，用来归整边像素，成为有意义的边。第24页，课件共164页，创作于2023年2月9/7/2023241.图像分割1.2边缘连接法局部连接处理连接处理的时机和目的连接处理的原理局部连接算法描述第25页，课件共164页，创作于2023年2月9/7/2023251.图像分割1.2边缘连接法局部连接处理连接处理的时机和目的：

时机：对做过边界检测的图像进行目的：连接间断的边连接处理的原理：

用比较梯度算子的响应强度和梯度方向确定两个点是否同属一条边第26页，课件共164页，创作于2023年2月9/7/2023261.图像分割1.2边缘连接法局部连接处理连接处理的原理：通过比较梯度，确定两个点的连接性：

对于点(x’,y’)，判断其是否与邻域内的边界点(x,y)相似，当：

|f

(x,y)–f

(x’,y’)|

T

其中T是一个非负的阈值第27页，课件共164页，创作于2023年2月9/7/2023271.图像分割1.2边缘连接法局部连接处理连接处理的原理：比较梯度向量的方向角

对于点(x’,y’)，判断其是否与邻域内的边界点(x,y)的方向角相似，当：

|

(x,y)–

(x’,y’)|<A

其中A是一个角度阈值第28页，课件共164页，创作于2023年2月9/7/2023281.图像分割1.2边缘连接法局部连接处理连接处理的原理：当梯度值和方向角都是相似的，则点(x’,y’)，与边点界(x,y)是连接的

第29页，课件共164页，创作于2023年2月9/7/2023291.图像分割1.2边缘连接法局部连接处理局部连接算法描述：1）设定A、T的阈值大小，确定邻域的大小2）对图像上每一个像素的邻域点进行分析， 判断是否需要连接。3）记录像素连接的情况，另开一个空间， 给不同的边以不同的标记。4）最后，删除孤立线段，连接断开的线段。第30页，课件共164页，创作于2023年2月9/7/2023301.图像分割1.2边缘连接法Hough变换问题的提出Hough变换的基本思想算法实现Hough变换的扩展第31页，课件共164页，创作于2023年2月9/7/2023311.图像分割1.2边缘连接法Hough变换问题的提出在找出边界点集之后，需要连接，形成完整的边界图形描述第32页，课件共164页，创作于2023年2月9/7/2023321.图像分割1.2边缘连接法Hough变换的基本思想对于边界上的n个点的点集，找出共线的点集和直线方程。对于任意两点的直线方程：y=ax+b，构造一个参数ab的平面，从而有如下结论：ab第33页，课件共164页，创作于2023年2月9/7/2023331.图像分割1.2边缘连接法Hough变换的基本思想xy平面上的任意一条直线，对应在参数ab平面上都有一个点。过xy平面一个点(x,y)的所有直线，构成参数ab平面上的一条直线。如果点(x1,y1)与点(x2,y2)共线，那么这两点在参数ab平面上的直线将有一个交点。在参数ab平面上相交直线最多的点，对应的xy平面上的直线就是我们的解。第34页，课件共164页，创作于2023年2月9/7/2023341.图像分割1.2边缘连接法Hough变换的基本思想abA第35页，课件共164页，创作于2023年2月9/7/2023351.图像分割1.2边缘连接法Hough变换算法实现由于垂直直线a,为无穷大，我们改用极坐标形式：

xcos

+ysin

=

参数平面为

，对应不是直线而是正弦曲线使用交点累加器，或交点统计直方图，找出相交线段最多的参数空间的点，然后找出该点对应的xy平面的直线线段。第36页，课件共164页，创作于2023年2月9/7/2023361.图像分割1.2边缘连接法Hough变换的扩展Hough变换不只对直线，也可以用于圆：（x–c1）2+(y-c2)2=c32

这时需要三个参数的参数空间。第37页，课件共164页，创作于2023年2月9/7/2023371.图像分割1.3阈值分割法阈值分割法通过交互方式得到阈值通过直方图得到阈值通过边界特性选择阈值简单全局阈值分割分割连通区域基于多个变量的阈值第38页，课件共164页，创作于2023年2月9/7/2023381.图像分割1.3阈值分割法阈值分割法阈值分割法的基本思想：确定一个合适的阈值T（阈值选定的好坏是此方法成败的关键）。将大于等于阈值的像素作为物体或背景，生成一个二值图像。Iff(x,y)

Tset255Elseset0在四邻域中有背景的像素，既是边界像素。025525502550255255255第39页，课件共164页，创作于2023年2月9/7/2023391.图像分割1.3阈值分割法阈值分割法阈值分割法的特点：适用于物体与背景有较强对比的情况，重要的是背景或物体的灰度比较单一。（可通过先求背景，然后求反得到物体）这种方法总可以得到封闭且连通区域的边界。灰度值f(x0,y0)T第40页，课件共164页，创作于2023年2月9/7/2023401.图像分割1.3阈值分割法通过交互方式得到阈值基本思想：在通过交互方式下，得到对象（或背景）的灰度值，比得到阈值T容易得多。假设：对象的灰度值(也称样点值)为f(x0,y0),且： T=f(x0,y0)–R有： f(x,y)

T f(x,y)

f(x0,y0)–R |f(x,y)–f(x0,y0)|

R其中R是容忍度，可通过试探获得。第41页，课件共164页，创作于2023年2月9/7/2023411.图像分割1.3阈值分割法通过交互方式得到阈值实施方法：（1）通过光标获得样点值f(x0,y0)（2）选取容忍度R（3）if|f(x,y)–f(x0,y0)|

Rset255 elseset0第42页，课件共164页，创作于2023年2月9/7/2023421.图像分割1.3阈值分割法通过直方图得到阈值基本思想边界上的点的灰度值出现次数较少T第43页，课件共164页，创作于2023年2月9/7/2023431.图像分割1.3阈值分割法通过直方图得到阈值取值的方法：取直方图谷底，为最小值的灰度值为阈值T缺点：会受到噪音的干扰，最小值不是预 期的阈值，而偏离期望的值；改进：取两个峰值之间某个固定位置，如中间位置上。由于峰值代表的是区域内外的典型值，一般情况下，比选谷底更可靠，可排除噪音的干扰第44页，课件共164页，创作于2023年2月9/7/2023441.图像分割1.3阈值分割法通过直方图得到阈值T第45页，课件共164页，创作于2023年2月9/7/2023451.图像分割1.3阈值分割法通过直方图得到阈值对噪音的处理 对直方图进行平滑处理，如最小二乘法，等不过点插值。第46页，课件共164页，创作于2023年2月9/7/2023461.图像分割1.3阈值分割法通过边界特性选择阈值基本思想：如果直方图的各个波峰很高、很窄、对称，且被很深的波谷分开时，有利于选择阈值。为了改善直方图的波峰形状，我们只把区域边缘的像素绘入直方图，而不考虑区域中间的像素。用微分算子，处理图像，使图像只剩下边界中心两边的值。第47页，课件共164页，创作于2023年2月9/7/2023471.图像分割1.3阈值分割法通过边界特性选择阈值基本思想：这种方法有以下优点：1)在前景和背景所占区域面积差别很大时，不会造一个灰度级的波峰过高，而另一个过低2)边缘上的点在区域内还是区域外的概率是相等的，因此可以增加波峰的对称性3)基于梯度和拉普拉斯算子选择的像素，可以增加波峰的高度第48页，课件共164页，创作于2023年2月9/7/2023481.图像分割1.3阈值分割法通过边界特性选择阈值算法的实现：1）对图像进行梯度计算，得到梯度图像。2）得到梯度值最大的那一部分（比如10%）的像素直方图3）通过直方图的谷底，得到阈值T。如果用拉普拉斯算子，不通过直方图，直接得到阈值，方法是使用拉普拉斯算子过滤图像，将0跨越点对应的灰度值为阈值T.第49页，课件共164页，创作于2023年2月9/7/2023491.图像分割1.3阈值分割法简单全局阈值分割基本思想：用前述方法获得阈值T，并产生一个二值图，区分出前景对象和背景算法实现：规定一个阈值T，逐行扫描图像。凡灰度级大于T的，颜色置为255；凡灰度级小于T的，颜色置为0。

适用场合：明度图像是可以控制的情况，例如用于工业监测系统中。第50页，课件共164页，创作于2023年2月9/7/2023501.图像分割1.3阈值分割法分割连通区域基本思想：用前述方法获得阈值T，并产生一个二值图，区分出单独的连通前景对象和背景区域算法实现：规定一个阈值T，上下左右4个方向进行逐行扫描图像凡灰度级大于T的，颜色置为255；凡灰度级小于T的，颜色置为0。第51页，课件共164页，创作于2023年2月9/7/2023511.图像分割1.3阈值分割法分割连通区域适用场合：印前等。先左后右，先上半部分、后下半部分第52页，课件共164页，创作于2023年2月9/7/2023521.图像分割1.3阈值分割法基于多个变量的阈值基本思想：把前面的方法扩展到多维空间，则寻找波谷的过程，变为寻找点簇的过程。算法实现：各维分量波谷之间进行逻辑与运算，从波谷重合的点，得到实际的阈值T。应用场合：有多个分量的颜色模型，如RGB模型、CMYK模型、HSI模型第53页，课件共164页，创作于2023年2月9/7/2023531.图像分割1.4面向区域的分割基本概念通过像素集合的区域增长区域分裂与合并第54页，课件共164页，创作于2023年2月9/7/2023541.图像分割1.4面向区域的分割基本概念目标：将区域R划分为若干个子区域R1,R2,…,Rn，这些子区域满足5个条件：1)完备性：2)连通性：每个Ri都是一个连通区域3)独立性：对于任意i≠j，Ri∩Rj=Ф第55页，课件共164页，创作于2023年2月9/7/2023551.图像分割1.4面向区域的分割基本概念4)单一性：每个区域内的灰度级相等， P（Ri）=TRUE，i=1,2,…,n5)互斥性：任两个区域的灰度级不等， P（Ri∪Rj）=FALSE，i≠j第56页，课件共164页，创作于2023年2月9/7/2023561.图像分割1.4面向区域的分割通过像素集合的区域增长算法实现：1）根据图像的不同应用选择一个或一组种子，它或者是最亮或最暗的点，或者是位于点簇中心的点。2）选择一个描述符（条件）3）从该种子开始向外扩张，首先把种子像素加入集合，然后不断将与集合中各个像素连通、且满足描述符的像素加入集合4）上一过程进行到不再有满足条件的新结点加入集合为止。第57页，课件共164页，创作于2023年2月9/7/2023571.图像分割1.4面向区域的分割通过像素集合的区域增长算法实现：区域A区域B

种子像素第58页，课件共164页，创作于2023年2月9/7/2023581.图像分割1.4面向区域的分割区域分裂与合并算法实现：1）对于图像中灰度级不同的区域，均分为四个子区域。2）如果相邻的子区域所有像素的灰度级相同，则将其合并。3）反复进行上两步操作，直至不再有新的分裂与合并为止。第59页，课件共164页，创作于2023年2月9/7/2023591.图像分割1.4面向区域的分割区域分裂与合并算法实现：第60页，课件共164页，创作于2023年2月9/7/2023601.图像分割1.4面向区域的分割区域分裂与合并算法实现：实际应用中还可作以下修改：P(Ri)的定义为：1）区域内多于80%的像素满足不等式 |zj-mi|<=2σi， 其中：zj是区域Ri中第j个点的灰度级， mi是该区域的平均灰度级， σi是区域的灰度级的标准方差。2）当P(Ri)=TRUE时，将区域内所有像素的灰度级置为mi。第61页，课件共164页，创作于2023年2月9/7/2023611.图像分割1.5数学形态学形态学（Morphology）原是对于动植物调查时采取的某种形式的研究。数学形态学（MathematicalMorphology）是分析几何形状和结构的数学方法，它建立在集合代数的基础上，是用集合论方法定量描述集合结构的学科。1985年之后，数学形态学逐渐成为分析图像几何特征的工具。第62页，课件共164页，创作于2023年2月9/7/2023621.图像分割1.5数学形态学数学形态学包括一组基本的形态学运算子：腐蚀（Erosion）、膨胀（Dilation）、开（Opening）、闭（Closing）等。运用这些算子及其组合来进行图像形状和结构的分析及处理。形态学的理论基础是集合论。在图像处理中形态学的集合代表着黑白和灰度图像的形状，如黑白图像中的所以黑像素点组成了此图像的完全描述。通常我们选择图像中感兴趣的目标图像区域像素集合来进行形态学变换。第63页，课件共164页，创作于2023年2月9/7/2023631.图像分割1.5数学形态学基本运算

集合关系

设A和S为R2的子集，A为为物体区域，B为某种结构单元，则B结构单元对A的关系有三类：S包含于A,S包含于A,S击不中(MISS)A，

平移，记A平移x为Ax，定义为

第64页，课件共164页，创作于2023年2月9/7/2023641.图像分割1.5数学形态学膨胀S为结构单元，广义的膨胀定义为当S为3×3结构元时，广义膨胀就为一般意义上的膨胀。一般意思上的膨胀是将与物体边界接触的背景像素合并到物体中的过程。如果物体是个圆，进行一次膨胀后，它的直径会增大两个像素。如果两个物体在某处用少于三个像素分开，膨胀后这两个物体就合并成为一个物体了。第65页，课件共164页，创作于2023年2月9/7/2023651.图像分割1.5数学形态学腐蚀S为结构单元，广义的腐蚀运算定义为当S为3×3结构元时，广义腐蚀就为一般意义上的腐蚀。简单的腐蚀运算是将一个物体沿边界减小的过程，在物体的周边较少一个像素。如果物体是一个圆，则进行一次腐蚀运算后，它的直径减少2。第66页，课件共164页，创作于2023年2月9/7/2023661.图像分割1.5数学形态学腐蚀和膨胀示意图

二值图

腐蚀

膨胀

腐蚀二值图膨胀第67页，课件共164页，创作于2023年2月9/7/2023671.图像分割1.5数学形态学开和闭运算

腐蚀运算后再进行膨胀运算的组合运算称为开运算（Opening）。开运算的效果：删除小物体；将物体拆分为小物体；平滑大物体边界而不明显改变它们的面积；膨胀运算后再进行腐蚀运算的组合运算称为闭运算（Closing）。闭运算的效果：填充物体的小洞；连接相近的物体；平滑物体的边界而不明显改变它们的面积。第68页，课件共164页，创作于2023年2月9/7/2023681.图像分割1.5数学形态学

腐蚀和膨胀的衍生运算

腐蚀的反复进行会导致物体消失，而膨胀的反复进行的结果是所以物体都合并到一起了。我们可以改变这些过程来产生一些别的效果以适应实际的应用。

收缩保持单个像素的物体不变的腐蚀运算称为收缩(Shrinking)。细化我们可以修改腐蚀计算过程来保持物体不被分开。首先我们进行有条件的常规的腐蚀过程，我们只是将要删除的像素打上标记而并不真正删除；然后逐步访问打上标记的像素，如果删除该标记像素不会分开物体，就删除它，否则就保留它。以上过程就是细化(Thinning)。细化的结果是把曲线型物体变成一个像素宽的线型图。第69页，课件共164页，创作于2023年2月9/7/2023691.图像分割1.5数学形态学骨骼化

和细化相关的一个算子是骨骼化，也被称为中轴变换或火烧草场算法。中轴是和边界上至少两点相切的圆的圆心的轨迹。可以用火烧草场来说明，设物体区域上铺满了草料，火从物体边界同时均匀地烧起，最后草场全部烧光火熄灭的地方就是它的骨骼或骨架。中轴变换的火烧草场算法示意图修剪在很多情况下，细化或骨骼化过程会留下很多短刺，这些是有两三个像素点的分支。这些短刺是由于边界上的单个像素的摆动引起的。短刺可以用3×3的算子来移去端点，然后重新建立删去的分支。第70页，课件共164页，创作于2023年2月9/7/2023701.图像分割1.5数学形态学加厚不把相近物体合并的膨胀过程称为加厚(Thickening)。和细化过程一样，它也可以分两步完成。和它互补的操作是对背景进行细化，任何一种腐蚀类的操作都伴随着膨胀类的操作作用与互补的图像区域上。一些分割技术使用非常紧凑的边界来包围物体来保证不出现物体的错误合并。通常，用来分割物体最好的边界总是太紧，给后续的测量带来困难。加厚操作可以对此进行修正，它增大边界而不合并物体。第71页，课件共164页，创作于2023年2月9/7/2023712.图像描述描述的基本概念表示法设计边界描述关系描述第72页，课件共164页，创作于2023年2月9/7/2023722.图像描述2.1描述的基本概念基本概念在用前一章的方法，把图像分割后，为了进一步的处理，分割后的图像一般要进行形式化的描述解决形式化表达问题一般有两种选择：1）根据区域的外部特征来进行形式化表示2）根据区域的内部特征（比较区域内部的象素值）来来进行形式化表示第73页，课件共164页，创作于2023年2月9/7/2023732.图像描述2.1描述的基本概念基本概念选择表达方式，要本着使数据变得更有利于下一步的计算工作。下一步工作是基于所选的表达方式描述这个区域，一般情况下：1）如果关注的焦点是形状特性，选择外部表示方式2）如果关注的焦点是反射率特性，如颜色、文理时，选择内部表示方式。3）所选表示方式，应该对尺寸、变换、旋转等变量尽可能的不敏感第74页，课件共164页，创作于2023年2月9/7/2023742.图像描述2.2表示法设计链码多边形逼近外形特征边界分段区域骨架第75页，课件共164页，创作于2023年2月9/7/2023752.图像描述2.2表示法设计链码定义：1）链码是一种边界的编码表示法。 2）用边界的方向作为编码依据。为简化边界的描述。一般描述的是边界点集。012301472354-链码8-链码第76页，课件共164页，创作于2023年2月9/7/2023762.图像描述2.2表示法设计链码算法：给每一个线段一个方向编码。有4-链码和8-链码两种编码方法。从起点开始，沿边界编码，至起点被重新碰到，结束一个对象的编码。第77页，课件共164页，创作于2023年2月9/7/2023772.图像描述2.2表示法设计链码举例：4-链码：000033333322222211110011第78页，课件共164页，创作于2023年2月9/7/2023782.图像描述2.2表示法设计链码问题1：1）链码相当长。2）噪音会产生不必要的链码。改进1：1）加大网格空间。2）依据原始边界与结果的接近程度，来确定新点的位置。第79页，课件共164页，创作于2023年2月9/7/2023792.图像描述2.2表示法设计链码举例：4-链码：003332221101第80页，课件共164页，创作于2023年2月9/7/2023802.图像描述2.2表示法设计链码问题2：1）由于起点的不同，造成编码的不同2）由于角度的不同，造成编码的不同改进2：1）从固定位置作为起点(最左最上)开始编码2）通过使用链码的首差代替码子本身的方式第81页，课件共164页，创作于2023年2月9/7/2023812.图像描述2.2表示法设计链码循环首差链码：用相邻链码的差代替链码例如：4-链码10103322循环首差为：33133030循环首差：1-2=-1(3) 3-0=3

0-1=-1(3) 3-3=0 1-0=1 2-3=-1(3) 0-1=-1(3) 2-2=0第82页，课件共164页，创作于2023年2月9/7/2023822.图像描述2.2表示法设计链码应用背景：如果边界的本身对于旋转和比例修改来说是无变化的，使用链码才是正确的。一般来说这是不可能的，实际应用时还需要改进。用链码后，对象只要用1)起点坐标，2)周长（边界点数）3)链码，4)对象编号，就可以描述。链码一般用于一幅图像中有多个对象的情况，对单个对象不适用。第83页，课件共164页，创作于2023年2月9/7/2023832.图像描述2.2表示法设计多边形逼近基本思想：用最少的多边形线段，获取边界形状的本质。寻找最小基本多边形的方法一般有两种：点合成法和边分裂法第84页，课件共164页，创作于2023年2月9/7/2023842.图像描述2.2表示法设计多边形逼近点合成算法思想举例：第85页，课件共164页，创作于2023年2月9/7/2023852.图像描述2.2表示法设计多边形逼近合成点算法：1）沿着边界选两个相邻的点对，计算首尾连接直线段与原始折线段的误差。2）如果误差小于预先设置的阈值。去掉中间点，选新点对与下一相邻点对，重复1）；否则，存储线段的参数，置误差为0，选被存储线段的终点为起点，重复1）2）。3）当程序的第一个起点被遇到，程序结束。第86页，课件共164页，创作于2023年2月9/7/2023862.图像描述2.2表示法设计多边形逼近合成点算法的问题：

顶点一般不对应于边界的拐点（如拐角）。因为新的线段直到超过误差的阈值才开始。下面讲到的分裂法可用于缓解这个问题第87页，课件共164页，创作于2023年2月9/7/2023872.图像描述2.2表示法设计多边形逼近边分裂算法思想举例：

第88页，课件共164页，创作于2023年2月9/7/2023882.图像描述2.2表示法设计多边形逼近分裂边算法：（1）连接边界线段的两个端点（如果是封闭边界，连接最远点）；（2）如果最大正交距离大于阈值，将边界分为两段，最大值点定位一个顶点。重复（1）；（3）如果没有超过阈值的正交距离，结束。第89页，课件共164页，创作于2023年2月9/7/2023892.图像描述2.2表示法设计外形特征基本思想：外形特征是一种用一维函数表达边界的方法。基本思想是把边界的表示降到一维函数第90页，课件共164页，创作于2023年2月9/7/2023902.图像描述2.2表示法设计外形特征函数定义——质心角函数：边上的点到质心的距离r，作为夹角的

的函数。A

r

r()2A第91页，课件共164页，创作于2023年2月9/7/2023912.图像描述2.2表示法设计外形特征举例：A

r

r()A第92页，课件共164页，创作于2023年2月9/7/2023922.图像描述2.2表示法设计外形特征问题：函数过分依赖于旋转和比例的变化改进：对于旋转——两种改进：a.选择离质心最远的点作为起点b.选择从质心到主轴最远的点作为起点对于比例变换：对函数进行正则化，使函数值总是分布在相同的值域里，比如说[0，1]第93页，课件共164页，创作于2023年2月9/7/2023932.图像描述2.2表示法设计边界分段基本概念：一个任意集合S（区域）的凸起外缘H是：包含S的最小凸起的集合H-S的差的集合被称为集合S的凸起补集DSSDS+D=H第94页，课件共164页，创作于2023年2月9/7/2023942.图像描述2.2表示法设计边界分段分段算法：给进入和离开凸起补集的变换点打标记来划分边界段。优点：不依赖于方向和比例的变化S第95页，课件共164页，创作于2023年2月9/7/2023952.图像描述2.2表示法设计边界分段问题：噪音的影响，导致出现零碎的划分。解决的方法：先平滑边界，或用多边形逼近边界，然后再分段第96页，课件共164页，创作于2023年2月9/7/2023962.图像描述2.2表示法设计区域骨架基本思想表示一个平面区域结构形状的重要方法是把它削减成图形。这种削减可以通过细化（也称为抽骨架）算法，获取区域的骨架来实现Blum的中轴变换方法（MAT）设:R是一个区域，B为R的边界点，对于R中的点p，找p在B上“最近”的邻居。如果p有多于一个的邻居，称它属于R的中轴（骨架）第97页，课件共164页，创作于2023年2月9/7/2023972.图像描述2.2表示法设计区域骨架问题：计算量大pRB第98页，课件共164页，创作于2023年2月9/7/2023982.图像描述2.2表示法设计区域骨架算法改进思想在保证产生正确的骨架的同时，改进算法的效率。比较典型的是一类细化算法，它们不断删去边缘，但保证删除满足：（1）不移去端点（2）不破坏连通性（3）不引起区域的过度腐蚀第99页，课件共164页，创作于2023年2月9/7/2023992.图像描述2.2表示法设计区域骨架一种细化二值区域的算法假设区域内的点值为1，背景值为0这个方法由对给定区域的边界点连续进行两个基本操作构成这里边界点是指任何值为1且至少有一个8邻域上的点为0的象素第100页，课件共164页，创作于2023年2月9/7/20231002.图像描述2.2表示法设计区域骨架基本操作1对于满足以下四个条件的边界点打标记准备删除：(a)2

N(p1)

6其中N(p1)是点p1的邻域中1的个数，即： N(p1)=p2+p3+…+p9(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第101页，课件共164页，创作于2023年2月9/7/20231012.图像描述2.2表示法设计区域骨架所有条件都满足，才打删除标记。删除并不立即进行，而是等到对所有边界点都打完标记后，再把作了标记的点一起删除。举例: N(p1)=4 S(p1)=3 p2*p4*p6=0 p4*p6*p8=0第2个条件没满足不打标记00p1110101p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第102页，课件共164页，创作于2023年2月9/7/20231022.图像描述2.2表示法设计区域骨架基本操作2条件(a)、(b)与操作1相同条件(c)、(d)改为：c’)

p2*p4*p8=0d’)

p2*p6*p8=0p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第103页，课件共164页，创作于2023年2月9/7/20231032.图像描述2.2表示法设计区域骨架细化算法细化算法的一轮操作包括：按操作1，给边界点打标记——删除点按操作2，给边界点打标记——删除点这个基本过程反复进行，直至没有点可以删除为止。此时算法终止。第104页，课件共164页，创作于2023年2月9/7/20231042.图像描述2.2表示法设计区域骨架算法分析：1）条件a)的分析：当轮廓点p1的8邻域上有1个或7个值为1的点时，不满足条件a。 有1个点说明：p1是骨架上的终点，显然不能删除有7个点说明：如果删除p1会引起区域的腐蚀2）条件b)的分析：当p1在宽度为1的笔划上时，不满足条件b。因而该条件保证了骨架的连续性。第105页，课件共164页，创作于2023年2月9/7/20231052.图像描述2.2表示法设计区域骨架算法分析：

（3）当(p4=0orp6=0)or(p2=0andp8=0)时，条件c,d同时满足。满足这个条件的点可能是右边、下边、左上角的边界点。任何一种情况下，p1都不是骨架的一部分，应被删除。当(p4=0andp6=0)or(p2=0orp8=0)时，条件c’,d’同时满足。满足这个条件的点可能是左边、上边、右下角的边界点，应被删除。p9p2p1p8p3p4p7p6p5p9p2p1p8p3p4p7p6p5第106页，课件共164页，创作于2023年2月9/7/20231062.图像描述2.2表示法设计区域骨架例:

第107页，课件共164页，创作于2023年2月9/7/20231072.图像描述2.3边界描述简单描述子形状数傅立叶描述子矩量第108页，课件共164页，创作于2023年2月9/7/20231082.图像描述2.3边界描述简单描述子边界的周长：是最简单的描述符之一。沿轮廓线计算象素的个数，给出了一个长度的近似估计边界的直径：边界B的直径是：

Diam(B)=max[D(pi,pj)]

D是欧氏距离或几何距离，pi,pj是边界上的点。直径的长度和直径的两个端点连线（这条线被称为边界的主轴）的方向，是关于边界的有用的描述符。第109页，课件共164页，创作于2023年2月9/7/20231092.图像描述2.3边界描述简单描述子边界的直径举例第110页，课件共164页，创作于2023年2月9/7/20231102.图像描述2.3边界描述简单描述子边界的曲率： 曲率被描述为斜率的变化率。近似：用相邻边界线段（描述为直线）的斜率差作为在边界线交点处的曲率描述子。

交点a处的曲率为

dk=k1–k2

其中k1、k2为相邻线段的斜率ak1k2第111页，课件共164页，创作于2023年2月9/7/20231112.图像描述2.3边界描述简单描述子边界的凸线段点： 当顶点p上的斜率是非负时，称其为凸线段上的点边界的凹线段点： 当顶点p上的斜率为负时，称其为凹线段上的点第112页，课件共164页，创作于2023年2月9/7/20231122.图像描述2.3边界描述简单描述子边界的凸线段点P1：边界的凹线段点P2：P1P2第113页，课件共164页，创作于2023年2月9/7/20231132.图像描述2.3边界描述形状数形状数定义：最小循环首差链码。循环首差链码：用相邻链码的差代替链码例如：4-链码10103322循环首差为：33133030循环首差:1-2=-1(3) 3-0=3 0-1=-1(3) 3-3=0 1-0=1 2-3=-1(3) 0-1=-1(3) 2-2=0第114页，课件共164页，创作于2023年2月9/7/20231142.图像描述2.3边界描述形状数形状数定义：例如：4-链码：10103322 循环首差：33133|030 形状数：03033133形状数序号n的定义： 形状数中阿拉伯数字的个数，对于封闭边界序号是偶数。如order4、6、8。第115页，课件共164页，创作于2023年2月9/7/20231152.图像描述2.3边界描述形状数形状数例如：

序号4链码：0321首差：3333形状：3333序号6链码：003221首差：303303形状：033033序号8链码：00032221首差：30033003形状：00330033第116页，课件共164页，创作于2023年2月9/7/20231162.图像描述2.3边界描述形状数形状数例如：

序号6链码：033211首差：330330形状：033033序号6链码：003221首差：303303形状：033033第117页，课件共164页，创作于2023年2月9/7/20231172.图像描述2.3边界描述形状数问题： 虽然链码的首差是不依赖于旋转的，但一般情况下边界的编码依赖于网格的方向。改进： 规整化网格方向，具体方法如下：第118页，课件共164页，创作于2023年2月9/7/20231182.图像描述2.3边界描述形状数几个基本概念：边界最大轴a:是连接距离最远的两个点的线段边界最小轴b:与最大轴垂直，且其长度确定的包围盒刚好包围边界。边界离心率c：最大轴长度与最小轴长度的比

c=a/b基本矩形:包围边界的矩形。第119页，课件共164页，创作于2023年2月9/7/20231192.图像描述2.3边界描述形状数基本概念举例边界最大轴a边界最小轴b基本矩形第120页，课件共164页，创作于2023年2月9/7/20231202.图像描述2.3边界描述形状数规整化网格方向算法的思想： 大多数情况下，将链码网格与基本矩形对齐，即可得到一个唯一的形状数。对一个给定的形状序号，处理步骤如下：（1）我们找出一个序号为n的矩形，它的离心率最接近于给定形状的基本矩形的离心率。第121页，课件共164页，创作于2023年2月9/7/20231212.图像描述2.3边界描述形状数（2）然后再用这个矩形构造网格。例如:如果n=12，所有序号为12的矩形（即周长为12）为2*4，3*3，1*5。如果2*4矩形的离心率最接近于给定边界的基本矩形的离心率，我们建立一个2*4的网格。（3）再得到链码。（4）最后，再得到循环首差。（5）首差中的最小循环数即为形状数。第122页，课件共164页，创作于2023年2月9/7/20231222.图像描述2.3边界描述形状数规整化网格方向算法举例：链码：000033222121首差：300030300313形状：0003030031330123第123页，课件共164页，创作于2023年2月9/7/20231232.图像描述2.3边界描述傅立叶描述子基本思想：（1）对于XY平面上的每个边界点，将其坐标用复数表示为：

s(k)=x(k)+jy(k) k=0,1,…,N-1y0y1x0x1jyxx(k)=xky(k)=yk第124页，课件共164页，创作于2023年2月9/7/20231242.图像描述2.3边界描述傅立叶描述子基本思想：（2）进行离散傅立叶变换

N-1a(u)=1/N∑s(k)exp(-j2

uk/N)u=0,1,…,N-1

u=0

N-1s(k)=∑a(u)exp(j2

uk/N) k=0,1,…,N-1

u=0a(u)被称为边界的傅立叶描述子第125页，课件共164页，创作于2023年2月9/7/20231252.图像描述2.3边界描述傅立叶描述子基本思想：（3）选取整数M

N-1，进行逆傅立叶变换（重构）

M-1s’(k)=∑a(u)exp(j2

uk/N) k=0,1,…,N-1

u=0这时，对应于边界的点数没有改变，但在重构每一个点所需要的计算项大大减少了。如果边界点数很大，M一般选为2的指数次方的整数。第126页，课件共164页，创作于2023年2月9/7/20231262.图像描述2.3边界描述傅立叶描述符M的选取与描述符的关系

在上述方法中，相当于对于u>M-1的部分舍去不予计算。由于傅立叶变换中高频部分对应于图像的细节描述，因此M取得越小，细节部分丢失得越多。第127页，课件共164页，创作于2023年2月9/7/20231272.图像描述2.3边界描述傅立叶描述符优点1）使用复数作为描述符，对于旋转、平移、放缩等操作和起始点的选取不十分敏感。2）以上几何变换均可以通过对描述子函数作简单变换来获得。第128页，课件共164页，创作于2023年2月9/7/20231282.图像描述2.3边界描述矩量基本思想：将描述形状的任务减少至描述一个一维函数，边界段和特征的形状可以用矩量来量化地描述矩量的定义：把边界当作直方图函数：g(r)rg(r)第129页，课件共164页，创作于2023年2月9/7/20231292.图像描述2.3边界描述矩量矩量的定义：

L

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

i=1L其中 m=∑rig(ri)

i=1这里L是边界上点的数目,

n(r)是边界的矩量第130页，课件共164页，创作于2023年2月9/7/20231302.图像描述2.3边界描述矩量矩量的优点：实现是直接的附带了一种关于边界形状的“物理”解释对于旋转的不敏感性为了使大小比例不敏感，可以通过伸缩r的范围来将大小正则化。

第131页，课件共164页，创作于2023年2月9/7/20231312.图像描述2.4关系描述基本思想阶梯关系编码骨架关系编码方向关系编码内角关系编码树结构关系编码第132页，课件共164页，创作于2023年2月9/7/20231322.图像描述2.4关系描述基本思想：通过挖掘各个成分之间的结构关系来描述边界图像中各个部分间的结构关系是二维的，而串是一维的，期望找到一种方法把二维关系转化为一维的串主导思想是考虑物体各个部分的连接线段第133页，课件共164页，创作于2023年2月9/7/20231332.图像描述2.4关系描述阶梯关系编码对于如下阶梯形边界，定义两个基本元素a,babaaabbb第134页，课件共164页，创作于2023年2月9/7/20231342.图像描述2.4关系描述阶梯结构关系定义如下产生规则：

(1)S->aA (2)A->bS (3)A->b

其中S、A是变量举例：(1,3)(1,2,1,3)(1,2,12,1,3)aaabbbaabbab第135页，课件共164页，创作于2023年2月9/7/20231352.图像描述2.4关系描述骨架关系编码用有向线段来描述一个图像的各个部分（例如同构区域），这个线段是通过头尾连接等方法得到的。线段之间的不同运算代表了区域的不同组合。当图像的连通性可以通过首尾相接或其它连续的方式描述的时候，最适于使用这种串来描述。第136页，课件共164页，创作于2023年2月9/7/20231362.图像描述2.4关系描述骨架关系编码a+ba-ba×ba*baaaabbbb编码第137页，课件共164页，创作于2023年2月9/7/20231372.图像描述2.4关系描述方向关系编码跟踪对象的边界，将跟踪得到的线段按照方向或长度来编码a1a2a5a7a8a3a4a6a1a8a7a6a5a4a3a2第138页，课件共164页，创作于2023年2月9/7/20231382.图像描述2.4关系描述内角关系编码根据角度范围不同，编码为8个符号即:a1:0-45;a2:45-90;a3:90-135;…;

a8:315-360举例：a3a3a3a3a3a3a3a3a2a2a3a3第139页，课件共164页，创作于2023年2月9/7/20231392.图像描述2.4关系描述树结构关系树结构中每个结点的意义和结点之间的关系最为重要举例：abcd$abcdefef$第140页，课件共164页，创作于2023年2月9/7/20231403.图像分析图像分析的定义图像分析是一个发现、辨认和理解模式的过程，这些模式都与执行与图像相关的任务有关。图像分析的目标计算机图像分析的主要目的之一是，赋予某些具有感觉能力的机器，以类似人的大脑的能力。例如OCR第141页，课件共164页，创作于2023年2月9/7/20231413.图像分析自动图像分析系统行为能力的概念化分类获取、发现信息：从背景中提取有关信息学习、应用知识：抽象、归纳信息特征的学习过程，并应用到新的对象中。构造、推理知识：从不完整的信息中构造推论出新的知识，并加以应用。第142页，课件共164页，创作于2023年2月9/7/20231423.图像分析自动图像分析系统的现状我们可以设计出这种系统，但仍然缺乏理论依据。有待人类视觉认识理论的进一步研究。我们可以做出在某一应用上超过人的能力的系统，但缺乏扩展性。过分依赖应用。第143页，课件共164页，创作于2023年2月9/7/20231433.图像分析图像分析技术分类的三种基本范畴低级处理：图像获取、预处理，不需要智能中级处理：图像分割、表示与描述，需要智能高级处理：图像识别、解释，缺少理论，为降低难度，设计得更专用。第144页，课件共164页，创作于2023年2月9/7/20231443.图像分析图像分析技术分类的三种基本范畴知识库分割表示与描述识别与解释预处理图像获取低级处理高级处理中级处理第145页，课件共164页，创作于2023年2月9/7/20231453.图像分析3.1纹理分析提到纹理，人们自然会立刻想到木制家俱上的木纹、花布上的花纹等。木纹为天然纹理，花纹为人工纹理，它们反映了物体表面颜色和灰度的某种变化。这些变化与物体本身的属性相关。

第146页，课件共164页，创作于2023年2月9/7/20231463.图像分析3.1纹理分析有些图像在局部区域内呈现不规则性，而在整体上表现出某种规律性。习惯上，把这种局部不规则而宏观有规律的特性称之为纹理；以纹理特性为主导的图像，常称为纹理图像；以纹理特性为主导特性的区域，常称为纹理区域。纹理作为一种区域特性，在图像的一定区域上才能反映或测量出来。

为了定量描述纹理，多年来人们建立了许多纹理算法以测量纹理特性。这些方法大体可以分为两大类：统计分析法和结构分析法。前者从图像有关属性的统计分析出发；后者则着力找出纹理基元，然后从结构组成上探索纹理的规律。也有直接去探求纹理构成的结构规律的。

本节将主要论述纹理特征提取与分析的几种方法。第147页，课件共164页，创作于2023年2月9/7/20231473.图像分析3.2纹理分析--直方图分析法

纹理区域的灰度直方图作为纹理特征,为了研究灰度直方图的相似性，可以比较累积灰度直方图分布，计算灰度级的最大偏差或总偏差。如果限定对象，则采用这样简单的方法也能够识别纹理。但是灰度直方图不能得到纹理的二维灰度变化，即使作为一般性的纹理识别法，其能力是很低的。例如下图两种纹理具有相同的直方图，只靠直方图就不能区别这两种纹理。第148页，课件共164页，创作于2023年2月9/7/20231483.图像分析3.3纹理分析--Laws纹理能量测量法

Laws的纹理能量测量法是一种典型的一阶分析方法，在纹理分析领域中有一定影响。Laws纹理测量的基本思想是设置两个窗口：一个是微窗口，可为3×3、5×5或7×7，常取5×5用来测量以像元为中心的小区域的灰度的不规则性，以形成属性，称为微窗口滤波；另一个为宏窗口，为15×15或32×32，用来在更大的窗口上求属性的一阶统计量（常为均值和标准偏差），他称之为能量变换。整个纹理分析过程为

f(x,y)

微窗口滤波

F(x,y)

能量转换

E(x,y)

分量旋转

C(x,y)

分类

M(x,y)第149页，课件共164页，创作于2023年2月9/7/20231493.图像分析3.3纹理分析--Laws纹理能量测量法

Laws深入研究了滤波模板的选定。首先定义了一维滤波模板，然后通过卷积形成系列一维、二维滤波模板，用于检测和度量纹理的结构信息。第150页，课件共164页，创作于2023年2月9/7/20231503.图像分析3.3纹理分析--Laws纹理能量测量法

他选定的三组一维滤波模板是：

L3=[121]灰度（Level）E3=[-101]边缘（Edge）S3=[-12-1]点（Spot）

L5=[14641]E5=[-1–2021]S5=[-1020–1]W5=[-120–21]波（Wave）

R5=[1–46–41]涟漪(Ripple）

第151页，课件共164页，创作于2023年2月9/7/20231513.图像分析3.3纹理分析--Laws纹理能量测量法

L7=[1615201561E7=[-1–4–50541]S7=[-1-2141–2–1]W7=[-1030–301]R7=[1-2–14–1–21]O7=[-16–1520–156–1]振荡Oscillation）第152页，课件共164页，创作于2023年2月9/7/20231523.图像分析3.3纹理分析--Laws纹理能量测量法

1×3的矢量集是构成更大矢量的基础.每一个1×5的矢量可以由两个1×3矢量的卷积产生。1×7的矢量可以由1×3与1×5矢量卷积产生。垂直矢量和水平矢量可生成二维滤波模板。由滤波模板与图像卷积可以检测不同的纹理能量信息。所以，Laws一般选用12—15个5×5的模板。第153页，课件共164页，创作于2023年2月9/7/20231533.图像分析3.3纹理分析--Laws纹理能量测量法

以1×5矢量为基础，卷积同样维数的矢量，可获得25个5×5模板。其中最有用的是5×5的零和模板，即

其中aij是模板中的元素（i,j=1,2,3,4,5）。其中四个有最强性能的模板是：

E5S5L5S5E5L5R5R5第154页，课件共164页，创作于2023年2月9/7/20231543.图像分析3.3纹理分析--Laws纹理能量测量法

它们分别可以滤出水平边缘、高频点、V形状和垂直边缘。Laws将Brodatz的8种纹理图像拼在一起，对该图像作纹理能量测量，将每个像元指定为八个可能类中的一个，正确率达87％。可见这种纹理分析方法简单、有效。但所提供的模板较少，尚未更多地给出其变化性质，因此，应用受到一定的限制。第155页，课件共164页，创作于2023年2月9/7/20231553.图像分析3.4纹理分析--纹理分析的自相关函数法

若有一幅图像f(i，j)，i，j=0，1，…，N-1，则该图像的自相关函数定义为第156页，课件共164页，创作于2023年2月9/7/20231563.图像分析3.4纹理分析--纹理