图像描述与分析幻灯片.ppt

《图像描述与分析幻灯片.ppt》由会员分享，可在线阅读，更多相关《图像描述与分析幻灯片.ppt（135页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、图像描述与分析第1页，共135页，编辑于2022年，星期五第7章 图像描述T图像描述（图像理解）是图像处理的一个重要分支。T图像描述是用一组数量或符号（也称描述子）来表达图像物体的某些特征。T图像描述既可以是图像中各个组成部分的性质的描述，也可以是其各个部分彼此间关系的描述。第2页，共135页，编辑于2022年，星期五特征形成 根据待识别的图像，通过计算产生一组原始特征，称之为特征形成。特征提取原始特征的数量很大，或者说原始样本处于一个高维空间中，通过映射或变换的方法可以将高维空间中的特征描述用低维空间的特征来描述，这个过程就叫特征提取。图像特征第3页，共135页，编辑于2022年，星期五图像

2、特征：指图像的原始特性或属性。其中有些是视觉直接感受到的自然特征，如区域的亮度、边缘的轮廓、纹理或色彩等；有些是需要通过变换或测量才能得到的人为特征，如变换频谱、直方图、矩等。常用的图像特征有颜色特征、形状特征、幅度特征、纹理特征、变换系数特征、空间特征。图像特征第4页，共135页，编辑于2022年，星期五颜色特征是一种全局特征,描述了图像或图像区域所对应的景物的表面性质。一般颜色特征是基于像素点的特征，此时所有属于图像或图像区域的像素都有各自的贡献。由于颜色对图像或图像区域的方向、大小等变化不敏感，所以颜色特征不能很好地捕捉图像中对象的局部特征。颜色特征 第5页，共135页，编辑于2022年

3、，星期五人们的视觉系统对于景物认识的初级阶段是形状。图像经过边缘提取与分割等操作，就会得到景物的边缘和区域，也就获取了景物的的形状。任何一个景物形状特征均可由其几何属性（如长短、面积、距离、凹凸等），统计属性(如投影)和拓扑属性（如连通、欧拉数）来进行描述。形状特征 第6页，共135页，编辑于2022年，星期五第7章 图像描述与分析7.1 灰度描述7.2 边界描述7.3 区域描述7.4 纹理描述第7页，共135页，编辑于2022年，星期五7.1 灰度描述7.1.1 幅度特征7.1.2 直方图特征7.1.3 变换系数特征第8页，共135页，编辑于2022年，星期五7.1.1 幅度特征在所有的图像

4、特征中，最基本的是图像的幅度特征。可以在某一像素点或其邻域内做出幅度的测量，例如在区域NN内的平均幅度，即：可以直接从图像像素的灰度值，或从某些线性、非线性变换后构成新的图像幅度的空间来求得各式各样图像的幅度特征图。第9页，共135页，编辑于2022年，星期五7.1.1 幅度特征 a)原图 b)利用幅度特征将目标分割出来第10页，共135页，编辑于2022年，星期五一幅数字图像可以看作是一个二维随机过程的一个样本，可以用联合概率分布来描述。通过对图像的各像素幅度值可以设法估计出图像的概率分布，从而形成图的直方图特征。图像灰度的一阶概率分布定义为：M为围绕（x,y）点被测窗口内的像素总数；N(b

5、)为该窗口内灰度值为b的像素总数。7.1.2 直方图特征第11页，共135页，编辑于2022年，星期五图像的直方图特征可以提供图像信息的许多特征。图像对比度、动态范围、明暗程度等。例如直方图密集地分布在很窄的区域之内，说明图像的对比度很低；若直方图有两个峰值，说明存在着两种不同亮度的区域。一阶直方图的特征参数：rk量化层 均值：方差：斜度：7.1.2 直方图特征第12页，共135页，编辑于2022年，星期五 峭度：熵：能量：7.1.2 直方图特征第13页，共135页，编辑于2022年，星期五二阶直方图的特征是以像素对的联合概率分布为基础得出的。若两个像素f(i,j)及f(m,n)分别位于（i,

6、j）点和（m,n）点，两者的间距为|i-m|、|j-n|，并可用极坐标、表达，幅度值的联合分布为：a、b为量化的幅度值。直方图的二阶分布为：7.1.2 直方图特征表示在图像中，在方向上、径向间距为的像素对f(i,j)=a，f(m,n)=b出现的频数；M窗口中像素的总数。第14页，共135页，编辑于2022年，星期五假设图像的各像素对都是相互关联的，则P(a,b)将在阵列的对角线上密集起来。以下一些度量，用来描述围绕P(a,b)对角线能量扩散的情况：自相关：协方差：7.1.2 直方图特征惯性矩：第15页，共135页，编辑于2022年，星期五绝对值：能量：7.1.2 直方图特征熵：第16页，共13

7、5页，编辑于2022年，星期五由于图像的二维变换得出的系数反映了二维变换后图像在频率域的分布情况，因此常常用二维的傅里叶变换作为一种图像特征的提取方法。7.1.3 变换系数特征设M(u,v)是F(u,v)的平方值，即当f(x,y)的原点有了位移时，M(u,v)的值保持不变，因此M(u,v)与F(u,v)不是唯一对应的，这种性质称为位移不变性。第17页，共135页，编辑于2022年，星期五v v(m+1)u v(m)水平切口垂直切口环状切口扇状切口7.1.3 变换系数特征T如果把M(u,v)在某些规定区域内的累计值求出，也可以把图像的某些特征突出起来。u(m)u(m+1)(m+1)(m)第18页

8、，共135页，编辑于2022年，星期五T特征：图像中含有这些切口的频谱成分的含量。把这些特征提取出来，可以作为模式识别或分类系统的输入信息。这种方法已成功用于土地情况分类，放射照片病情诊断等。7.1.3 变换系数特征第19页，共135页，编辑于2022年，星期五为了描述目标物的二维形状，通常采用的方法是利用目标为了描述目标物的二维形状，通常采用的方法是利用目标物的边界来表示物体，即所谓的边界描述。当一个目标区物的边界来表示物体，即所谓的边界描述。当一个目标区域边界上的点已被确定时，就可以利用这些边界点来区别域边界上的点已被确定时，就可以利用这些边界点来区别不同区域的形状。这样做既可以节省存储信

9、息，又可以准不同区域的形状。这样做既可以节省存储信息，又可以准确的确定物体。确的确定物体。7.2 边界描述7.2.1 链码描述7.2.2 傅里叶描述子第20页，共135页，编辑于2022年，星期五7.2.1 链码描述T在数字图像中，边界或曲线是由一系列离散的像素点在数字图像中，边界或曲线是由一系列离散的像素点组成的，其最简单的表示方法是由美国学者组成的，其最简单的表示方法是由美国学者FreemanFreeman提提出的链码方法。出的链码方法。T链码实质上是一串指向符的序列，有链码实质上是一串指向符的序列，有4向链码、向链码、8向链向链码等。码等。第21页，共135页，编辑于2022年，星期五

10、4向链码 8向链码7.2.1 链码描述对任一像素点P，考虑它的8个邻近像素，指向符共有8个方向，分别用0、1.7表示。链码表示就是从某一起点开始沿曲线观察每一段的走向并用相应的指向符来表示，结果形成一个数列。可以用链码来描述任意曲线或者闭合的边界。第22页，共135页，编辑于2022年，星期五 a)原链码方向 b)逆时针旋转90 图a曲线的链码为：01122233100000765556706 图b曲线的链码为：23344455322222107770120 7.2.1 链码描述第23页，共135页，编辑于2022年，星期五用链码表示给定目标的边界时，如果目标平移，链码不会发生变化，而如果目标

11、旋转则链码会发生变化。为解决这个问题，可利用链码的一阶差分来重新构造一个表示原链码各段之间方向变化的新序列，相当于把链码进行旋转归一化。差分链码可用相邻两个方向数按反方向相减（后一个减去前一个），并对结果作模8运算得到。7.2.1 链码描述第24页，共135页，编辑于2022年，星期五 a)原链码方向 b)逆时针旋转90 图a曲线的链码为：01122233100000765556706其差分链码为：1010010670000777001116 图b曲线的链码为：23344455322222107770120 其差分链码为：10100106700007770011167.2.1 链码描述曲线旋转

12、到不同位置对应不同的链码，但其差分链码不变。第25页，共135页，编辑于2022年，星期五曲线的链码是：6022222021013444444454577012其差分链码是：2200006277121000000171201117.2.1 链码描述第26页，共135页，编辑于2022年，星期五曲线的链码是：024444424323566666676711234其差分链码是：220000627712100000171201117.2.1 链码描述第27页，共135页，编辑于2022年，星期五T链码的特殊性质 一个物体很容易实现旋转45。如果一个物体旋转n45，可由原链码加上 n 倍的模8得到链码

13、的微分，也称差分码，由原码的一阶差分求得链码差分是关于旋转不变的边界描述方法区域的一些其它性质，如面积和角点，可以由链码直接得7.2.1 链码描述第28页，共135页，编辑于2022年，星期五7.2.2 傅里叶描述子T对边界的离散傅里叶变换表达，可以作为定量描述边界形状的基础。采用傅里叶描述的一个优点是将二维的问题简化为一维问题。即将x-y平面中的曲线段转化为一维函数f(r)(在r-f(r)平面上)，也可将x-y平面中的曲线段转化为复平面上的一个序列。T转换方法：将x-y平面与复平面u-v重合，其中，实部u轴与x轴重合，虚部v轴与y轴重合。这样可用复数u+jv的形式来表示给定边界上的每个点(x

14、,y)。第29页，共135页，编辑于2022年，星期五7.2.2 傅里叶描述子T两种表示方法本质上是一致的，是点点对应的。两种表示方法本质上是一致的，是点点对应的。边界点的两种表示方法第30页，共135页，编辑于2022年，星期五现在考虑一个由N点组成的封闭边界，从任一点开始绕边界一周就得到一个复数序列，即：7.1.3 变换系数特征s(k)的离散傅里叶变换是：S(w)可称为边界的傅里叶描述，它的傅里叶反变换是：第31页，共135页，编辑于2022年，星期五在变换过程中信息没有任何增减，为我们有选择地描述边界提供了方便。取S(w)的前M个系数即可得到s(k)的一个近似：7.1.3 变换系数特征k

15、的范围不变，即在近似边界上的点数不变，但w的范围缩小了，即重建边界点所需的频率阶数减少了。傅里叶变换的高频分量对应一些细节而低频分量对应总体形状，因此用一些低频分量的傅里叶系数足以近似描述边界形状。第32页，共135页，编辑于2022年，星期五7.3 区域描述对一幅灰度图像或者彩色图像运用图像分割的方法进行处理，把其中感兴趣的像素分离出来作为目标像素，不感兴趣的部分作为背景像素，就可以得到一幅二值图像。理想情况下，希望该二值图像中的两个值准确地代表“目标”及“背景”。但实际中往往所检测道的目标中还有若干个假目标出现，还有可能提取的是多个目标，因此就需要对二值图像进行处理，实现对目标的分析。二值

16、图像包含目标的位置、形状、结构等很多重要信息，是图像分析和目标识别的依据。第33页，共135页，编辑于2022年，星期五7.3 区域描述7.3.1 几何特征7.3.2 不变矩第34页，共135页，编辑于2022年，星期五 1、像素与、像素与邻域邻域对于任意像素(i,j),(s,t)是一对适当的整数，则把像素的集合(i+s,j+t)叫做像素(i,j)的邻域.直观上看，这是像素(i,j)附近的像素形成的区域.最经常采用的是4-邻域和8-邻域 a)4-邻域 b)8-邻域 7.3.1 几何特征第35页，共135页，编辑于2022年，星期五区域面积是区域的一个基本特性，描述区域的大小。图像中的区域面积用

17、同一标记的区域内像素的个数总和来表示。区域的面积和周长 按上述表示法区域按上述表示法区域R R的面积的面积S=41。T2、区域面积、区域面积第36页，共135页，编辑于2022年，星期五 3.位置位置质心形心由于目标在图像中总有一定的面积大小，因此有必要定义目标在图像中的精确位置。目标的位置有形心、质心之分，形心为目标形状的中心，质心为目标质量的中心。第37页，共135页，编辑于2022年，星期五4.区域周长区域周长周长是区域形状的一个重要的整体参数。数字图像子集S的周长定义有不同概念，常用的有以下三种定义：(1)用边界所占面积表示，即边界点数之和；(3)区域边界8链码的长度（2）用区域中相邻

18、边缘点间距离之和来表示。在区域的边界像素中，设某像素与其水平或垂直方向上相邻边缘像素间的距离为1，与倾斜方向上相邻边缘像素间的距离为 。周长就是这些像素间距离的总和。第38页，共135页，编辑于2022年，星期五5.方向方向二阶矩轴:物体上的全部点到该线的距离平方和最小其中 是物体点到直线 的距离第39页，共135页，编辑于2022年，星期五6.距离距离 1)欧几里德距离（Euclidean）2)4-邻域距离(City-block城区距离)3)8-邻域距离（Chessboard棋盘距离）图像中两点P(x,y)和Q(u,v)之间的距离是重要的几何特性。常用以下三种方法测量：第40页，共135页，

19、编辑于2022年，星期五如图所示。街区距离和棋盘距离都是欧式距离的一种近似。第41页，共135页，编辑于2022年，星期五 下图中表示了以中心像素为原点的各像素的距离。从离开一个像素的等距离线可以看出，在欧氏距离中大致呈圆形，在棋盘距离中呈方形，在街区距离中呈倾斜45度的正方形。6.距离距离第42页，共135页，编辑于2022年，星期五 描述连通域与圆形相似程度的量。根据圆周长与圆面积的计算公式，定义圆形度的计算公式如下：其中，为连通域S的面积；为连通域S的周长。圆形度 值越大，表明目标与圆形的相似度越高。7.圆形度圆形度第43页，共135页，编辑于2022年，星期五8.矩形度矩形度 描述连通

20、域与矩形相似程度的量 其中，为连通域S的面积；是包含该连通域的最小矩形的面积。对于矩形目标，矩形度 取最大值1，对细长而弯曲的目标，则矩形度的值变得很小。第44页，共135页，编辑于2022年，星期五9.长宽比长宽比 其中，是包围连通域的最小矩形的宽度；是包围连通域的最小矩形的长度。是将细长目标与近似矩形或圆形目标进行区分时采用的形状度量。第45页，共135页，编辑于2022年，星期五10、欧拉数、欧拉数 图像的欧拉数是图像的拓扑特性之，它表明了图像的连通性。欧拉数定义为一个图中或一个区域中的连接部分数C和孔数H的差 E=C-H。下图(a)的图形有一个连接成分和一个孔，所以它的欧拉数为0，而下

21、图（b）有一个连接成分和两个孔，所以它的欧拉数为-1。具有欧拉数为具有欧拉数为0和和-1的图形的图形 可见通过欧拉数可见通过欧拉数 可用于目标识别。可用于目标识别。第46页，共135页，编辑于2022年，星期五该式描述了区域单位面积的周长大小，e值越大，表明单位面积的周长大，即区域离散，则为复杂形状；反之，则为简单形状。e值最小的区域为圆形。典型连续区域的计算结果为：圆形e=12.6；正方形e=16.0；正三角形e=20.8。11、形状复杂性、形状复杂性e：第47页，共135页，编辑于2022年，星期五（12）凹凸性凹凸性 凹凸性是区域的基本特征之一，区域凹凸性可通过以下方法进行判别：区域内任

22、意两像素间的连线穿过区域外的像素，则此区域为凹形。换句话说，区域中任意两个1像素之间，只要出现0像素，就称为凹形。相反，区域内任意两个像素的间的连线不穿过区域外的像素，则称为凸形。第48页，共135页，编辑于2022年，星期五 如图(a)第三行的区域内像素(3,3)与像素（3，8）间的像素值顺次为110001，出现了3个0像素，因此该区域为凹形。图(b)为凸形。（12）凹凸性凹凸性 第49页，共135页，编辑于2022年，星期五包含任一凹形的最小凸形称为该凹形的凸封闭包。将凸封闭包减去凹形得到凹形面积。为表示区域的凹特性，可采用凹性率E参数；E=Se/S Se为凹行面积，S为区域面积，图(c)

23、为图(a)的凸封闭包，(c)中图像减去(a)中的图像得到图(d)的结果。（12）凹凸性凹凸性 第50页，共135页，编辑于2022年，星期五7.3.2 不变矩由于图像区域的某些矩对于平移、旋转、尺度等几何变换具有一些不变的特性，因此，矩的表示方法在物体的分类与识别方面具有重要的意义。第51页，共135页，编辑于2022年，星期五对于二维连续函数 ，阶矩定义为：由于j和k可取所有的非负整数值，因此形成了一个矩的无限集。而且，这个集合完全可以确定函数f(x,y)本身。也就是说集合mjk对于函数f(x,y)是唯一的，也只有f(x,y)才具有这种特定的矩集,参数j+k称为矩的阶。为了描述物体的形状，假

24、设f(x,y)的目标物体取值为1，背景为0，即函数只反映了物体的形状而忽略其内部的灰度级细节。1.矩的定义第52页，共135页，编辑于2022年，星期五零阶矩是物体的面积：j=1,m=0,m10对二值图像来讲就是物体上所有的点x坐标的总和。m01是物体上所有的点y坐标的总和。令：1.矩的定义则：是二值图像中一个物体的质心的坐标。第53页，共135页，编辑于2022年，星期五中心矩定义为：如果f(x,y)是数字图像，则上式变为：1.矩的定义第54页，共135页，编辑于2022年，星期五 定义归一化的中心矩为：利用归一化的中心矩，可以获得对平移、缩放、镜像和旋转都不敏感的7个不变矩，定义如下：2.

25、不变矩第55页，共135页，编辑于2022年，星期五2、不变矩第56页，共135页，编辑于2022年，星期五2、不变矩 原图 原图缩小一半 逆时针旋转450逆时针旋转900 逆时针旋转1350 逆时针旋转1800 第57页，共135页，编辑于2022年，星期五2、不变矩以上图像的7个不变矩见表7-7 P139第58页，共135页，编辑于2022年，星期五7.4 纹理描述7.4.1 矩分析法 7.4.2 灰度差分统计法7.4.3 灰度共生矩阵法 7.4.4 纹理的结构分析第59页，共135页，编辑于2022年，星期五T纹理：是指人们所观察到的图像像素（或子区域）的灰度变化规律。灰度与颜色的二维变

26、化的图案，是区域的重要特征之一，习惯上把图像中这种局部不规则而宏观有规律的特性称为纹理。自然纹理：种子、草地人工纹理：织物、砖墙一般来说纹理图像中的灰度分布具有周期性，即使灰度变化是随机的，它也具有一定的统计特性。7.4 纹理描述第60页，共135页，编辑于2022年，星期五T纹理标志三要素1）某种局部的序列性在该序列更大的区域内不断重复2）序列基本元素是非随机排列组成的3）区域内任何地方都有大致相同的结构尺寸以上这些只是从感觉上看来是合理的，并不能得出定量的纹理测量。7.4 纹理描述第61页，共135页，编辑于2022年，星期五a)结构型纹理 b)随机型纹理7.4 纹理描述T根据纹理的局部统

27、计特征可以将纹理分为结构型纹理和随机型纹理。根据纹理的局部统计特征可以将纹理分为结构型纹理和随机型纹理。T结构纹理：具有独立基本结构与明显周期性的纹理（裂纹、砖墙）结构纹理：具有独立基本结构与明显周期性的纹理（裂纹、砖墙）反之称为随机型纹理（天空白云、气象云图）。反之称为随机型纹理（天空白云、气象云图）。第62页，共135页，编辑于2022年，星期五T描述纹理图像特征的参数有许多种，如 1）知道像素及邻近像素的灰度分布情况。2）检查小区域内灰度直方图，检查各小区域直方图的相似性，具有相似直方图的小区域同属一个大区域7.4 纹理描述第63页，共135页，编辑于2022年，星期五7.4 纹理描述第

28、64页，共135页，编辑于2022年，星期五7.4 纹理描述第65页，共135页，编辑于2022年，星期五(1)均值（Mean）7.4.1 矩分析法 T纹理分析的最简单方法之一是基于图像灰度直方图的矩分析法。令K为代表灰度级的随机变量，f(ki),i=0,1,2,N-1,(N灰度级数目)为对应的灰度直方图。均值给出了该图像区域平均灰度水平的估计值，它一般不反映什么具体纹理特征，但可以反映纹理的“光密度值”.第66页，共135页，编辑于2022年，星期五7.4.1 矩分析法 T表明区域灰度的离散程度，一般反映图像纹理的幅度。(2)方差（Variance）(3)扭曲度（Skewness）T反映直方

29、图的对称性，表示偏离平均灰度的像素的百分比。第67页，共135页，编辑于2022年，星期五(5)熵（Entropy）7.4.1 矩分析法 (4)峰度（Kurtosis）T反映直方图是倾向于聚集在均值附近还是散布在尾端。第68页，共135页，编辑于2022年，星期五T灰度差分统计法又称一阶统计法，通过计算图像中一对像素间灰度差分直方图来反映图像的纹理特征。令 为两个像素间的位移矢量，是位移量为 的灰度差分：T粗纹理时，位移相差为 的两像素通常有相近的灰度等级，因此，值较小，灰度差分直方图值集中在 附近；T细纹理时，位移相差为 的两像素的灰度有较大变化，值一般较大，灰度差分直方图值会趋于发散 7.

30、4.2 灰度差分统计法第69页，共135页，编辑于2022年，星期五T灰度差分统计法采用以下参数描述纹理图像的特征：7.4.2 灰度差分统计法对比度能量是灰度差分均匀性的度量，当p(i)值较平坦时，ASM值较小，当p(i)大小不均时，ASM值较大。第70页，共135页，编辑于2022年，星期五T灰度差分统计法采用以下参数描述纹理图像的特征：7.4.2 灰度差分统计法熵反映差分直方图的一致性，对于均匀分布的直方图，熵值较大。熵均值均值较小，说明p(i)值分布在i=0附近，纹理较粗糙，反之，均值较大，说明p(i)值分布远离原点，纹理较细。第71页，共135页，编辑于2022年，星期五T如果图像纹理

31、有方向性，则p(i)值得分布会随着方向矢量的变化而变化。可以通过比较不同方向上p(i)的统计量来分析纹理的方向性。例如，一幅图像在某一方向上灰度变化很小，则在该方向上得到的f(x,y)值较小，p(i)值多集中于i=0附近，它的均值较小，熵值也较小，能量值较大。7.4.2 灰度差分统计法可见，差分直方图分析方法不仅计算简单，而且能够反映纹理的空间组织情况，克服了基于灰度直方图的矩分析法不能表达纹理空间域特征的不足。第72页，共135页，编辑于2022年，星期五T灰度共生矩阵（Gray Level Co-occurrence Matrix）是由Haralick提出的一种用来分析图像纹理特征的重要方

32、法，是常用的纹理统计分析方法之一，它能较精确地反映纹理粗糙程度和重复方向。T灰度共生矩阵是建立在图像的二阶组合条件概率密度函数的基础上，即通过计算图像中特定方向和特定距离的两像素间从某一灰度过渡到另一灰度的概率，反映图像在方向、间隔、变化幅度及快慢的综合信息。7.4.3 灰度共生矩阵法 第73页，共135页，编辑于2022年，星期五T设f(x,y)为一幅NN的灰度图像，d=(dx,dy)是一个位移矢量，其中dx是行方向上的位移，dy是列方向上的位移，L为图像的最大灰度级数。T灰度共生矩阵定义为从f(x,y)的灰度为i的像素出发，统计与距离为=（dx2+dy2）1/2,灰度为j的像素同时出现的概

33、率P(i,j|d,)，7.4.3 灰度共生矩阵法 ijdxxdyy0f(x,y)=if(x+dx,y+dy)=j第74页，共135页，编辑于2022年，星期五T数学表达式为：7.4.3 灰度共生矩阵法 T灰度共生矩阵的第i行第j列元素表示图像上两个相距为、方向为、分别具有灰度级i和j的像素点对出现的次数。（x,y）是图像中的像素坐标，x,y的取值范围为0，N-1，i,j的取值范围为0，L-1,一般取00，450，900，1350。ijdxxdyyf(x,y)=if(x+dx,y+dy)=j第75页，共135页，编辑于2022年，星期五T说明：7.4.3 灰度共生矩阵法 P(i,j|d,)是一个

34、对称矩阵，其维数由图像中的灰度级数决定。若图像的最大灰度级数为L，则灰度共生矩阵为LL矩阵。这个矩阵是距离和方向的函数，在规定的计算窗口或图像区域内统计符合条件的像素对数。第76页，共135页，编辑于2022年，星期五 当两像素间位置关系选定后，就生成一定下的灰度共生矩阵P，灰度共生矩阵中的一个元素表示了一种灰度组合下出现的次数。如元素P（1，0）表示了图像上位置关系为的两像素灰度分别为1和0的情况出现的次数。7.4.3 灰度共生矩阵法 第77页，共135页，编辑于2022年，星期五T例:7.4.3 灰度共生矩阵法 设图像矩阵为00450第78页，共135页，编辑于2022年，星期五水平方向无

35、重复，变化较快水平方向数值大，重复多，纹理较粗1）00、900、1350方向上的灰度共生矩阵的对角线元素全为0，表明图像在该方向上灰度无重复、变化快，纹理细；2）450方向上灰度共生矩阵的对角线元素较大，表明图像在该方向上灰度变化慢、纹理较粗。7.4.3 灰度共生矩阵法 第79页，共135页，编辑于2022年，星期五T对于不同的，矩阵元素的定义如下：7.4.3 灰度共生矩阵法 水平方向：垂直方向：西北-东南方向东北-西南方向第80页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 第81页，共135页，编辑于2022年，星期五练习：求下面数字图像的四个方向的共生矩阵。练习：求下

36、面数字图像的四个方向的共生矩阵。一幅数字灰度图像0 0 0 0 1 1 1 10 0 0 0 1 1 1 10 0 0 0 1 1 1 10 0 0 0 1 1 1 12 2 2 2 2 2 3 32 2 2 2 2 2 3 32 2 2 2 2 2 3 32 2 2 2 2 2 3 3第82页，共135页，编辑于2022年，星期五灰度共生矩阵计算结果灰度共生矩阵计算结果 第83页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 (1)对比度（Contrast,CON）反映了图像的清晰度和纹理沟纹深浅的程度。纹理的沟纹深，对比度大，清晰；反之，对比度小，则沟纹浅，效果模糊。灰度

37、共生矩阵反映了图像灰度分布关于方向、邻域和变化幅度的综合信息，但它并不能直接提供区别纹理的特性。因此，要进一步从灰度共生矩阵中提取描述图像纹理的特征，用来定量描述纹理特性。设在取定d、参数下将灰度共生矩阵P(i,j|d,)归一化记为第84页，共135页，编辑于2022年，星期五（2）能量（角二阶矩Angular second moment,ASM）7.4.3 灰度共生矩阵法 ASM 是图像灰度分布均匀性和纹理粗细的程度。当灰度共生矩阵的元素分布较集中于主对角线附近，说明从局区域观察图像的灰度分布是均匀的。从图像的整体来观察，纹理较粗，该值较大，即粗纹理含有较多的能量；反之，纹理细，ASM较小，

38、含有较少的能量。第85页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 (3)熵（Entropy,ENT）熵是图像所具有信息量的度量，表示图像中纹理的非均匀程度或复杂程度。若图像没有纹理，灰度共生矩阵几乎为零矩阵，熵值接近0；若纹理复杂，灰度共生矩阵中的数值近似相等，熵值大。第86页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 (4)相关性（Correlation,COR）COR 描述了灰度共生矩阵中行或列元素间的相似程度，是灰度线性关系的度量。当矩阵元素值均相等时，相关值大，反之，相关性小。其中x,y,x,y分别是 Px，Py的均值和方差，Px，Py的定

39、义是：第87页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 在基于GLCM 的14 个纹理特征中,仅有4 个特征是不相关的,这4 个特征既便于计算又能给出较高的分类精度，一般采用这四个最常用的特征来提取图像的纹理特征。除了上面4个参数外，还有 （6)反差分矩（Inverse difference moment,IDM）(5)方差和（Sum variance,SV）第88页，共135页，编辑于2022年，星期五7.4.3 灰度共生矩阵法 （7）中值（8）协方差（9）同质性/逆差距:反映图像纹理的同质性，度量图像纹理局部变化的多少。其值大则说明图像纹理的不同区域间缺少变化，局部

40、非常均匀。(10)差异性（11）自相关等。第89页，共135页，编辑于2022年，星期五7.4.4 纹理的结构分析纹理基元可以是一个像素，也可以是若干灰度上比较一致的像素点集合。纹理的表达可以是多层次的。纹理的结构分析方法认为复杂的纹理可由一些简单的纹理基元按某种规则重复排列组合而成。为了分析纹理结构，必须提取结构基元，并描述其特性和分布规则。第90页，共135页，编辑于2022年，星期五7.4.4 纹理的结构分析纹理结构的描述及排列 可以从像素或小块纹理一层一层地向上拼合。当然基元的排列可有不同的规则，第一级纹理排列为YXY，第二级排列为XYX等，其中X,Y代表基元或子纹理。第91页，共13

41、5页，编辑于2022年，星期五a)纹理基元b）由规则 生成的纹理模式c)由 和其它规则生成的二维纹理模式纹理结构分析图例 7.4.4 纹理的结构分析第92页，共135页，编辑于2022年，星期五频谱法借助于傅里叶频谱的特性来描述周期或近乎周期的2-D图像模式的方向性。常用的3个性质是：（1）傅里叶频谱中突起的峰值对应纹理模式的主方向；（2）这些峰在频域平面的位置对应模式的基本周期；（3）如果利用滤波把周期性成分除去，剩下的非周期性部分将可用统计方法描述。纹理的频谱法第93页，共135页，编辑于2022年，星期五小波变换法模型法：分形模型、随机场模型、同步自然回归模型。纹理的其他分析方法第94页

42、，共135页，编辑于2022年，星期五 基于图像灰度共生矩阵的特征提取算法基于图像灰度共生矩阵的特征提取算法1、将彩色图像转化为灰度图像；2、灰度级粗量化；由于求灰度共生矩阵的计算量较大，一般将灰度级粗量化，如256级量化为16级或8级。经过量化后的图像虽有失真，但 对纹理影响不大；3、求出四个方向上的共生矩阵；4、分别求出每个共生矩阵的特征。第95页，共135页，编辑于2022年，星期五 医学应用医学应用1、超声图像处理：肝脏、脾脏、2、X线：胸部如肺癌、肺结节等3、MRI：癫痫4、皮肤表面是人体生理健康状况的一面镜子，随着生活水平的提高，人们越来越关注人脸面部皮肤的美容及护理。皮肤粗糙度随

43、着年龄或外界的影响而变化，采用无创性定量描述皮肤纹理的分析方法，对皮肤生理和病理学研究或特殊类型化妆品的功效评价有重要意义。第96页，共135页，编辑于2022年，星期五超声图像的特征提取与分析超声图像的特征提取与分析T另外，近年来，随着生活水平的提高，饮食结构的变化及预防保健措施的相对滞后，脂肪肝的发病率呈上升趋势，且发病年龄越来越小。2000年美国加利福尼亚太平洋医学中心调查结果显示，成年人中脂肪肝已成为第三大疾病。这就越来越要求对脂肪肝的检查能够及时、准确、定量分析判断。虽然肝活检是目前诊断脂肪肝的金标准，但它是一种创伤性检查，故难以接受。因此，探索一种较为简便可靠的脂肪肝临床诊断方法，

44、对脂肪肝的研究、临床及预防工作具有积极的现实意义。第97页，共135页，编辑于2022年，星期五T一般情况下，脂肪肝可以分为轻度、中度和重度三个等级，它们的B超图像呈现各自不同的特征。多数文献公开研究报告中均注重区分正常肝与脂肪肝或是肝硬化、肝癌等，没有进一步量化脂肪肝的严重程度。超声影像检查是诊断脂肪肝的常用方法，但是由于超声图像与CT和MR等医学图像相比,质量较差、分辨率低等原因，诊断常以定性为主，受主观经验因素影响较大。超声图像的特征提取与分析超声图像的特征提取与分析第98页，共135页，编辑于2022年，星期五 超声图像的特征提取与分析超声图像的特征提取与分析 超声图像具有一定的纹理特

45、征，通过对图像进行纹理分析提取特征参数，为脂肪肝超声诊断提供量化指标。针对脂肪肝超声图像的特点，本文在图像预处理部分，对获得的图像进行降噪、增强，然后进行图像纹理特征提取，获取多方面的超声图像的参数指标，为医生提供诊断依据。第99页，共135页，编辑于2022年，星期五 超声图像的采集超声图像的采集 选用新西兰大白兔12 只（徐州医学院动物实验中心提供）,体质量1.52.5 kg,平均每只2.15 kg，雌雄不限，肝功能检查均正常。正常时采集一组超声图像，在普通饲料基础上加2%胆固醇、5%蔗糖、10%猪油喂养,每10天后各采集一组（此时脂肪肝程度越来越重），同时进行肝活检检查。超声图像采集（徐

46、州市医学科学研究所项目合作）采用上海麦迪逊SA-5000超声诊断仪，由于家兔皮层薄采用了7.5 MHz高频探头进行超声检查，并在相同机器参数下进行，有利于对比分析。第100页，共135页，编辑于2022年，星期五 特征提取的指标特征提取的指标 对采集的每幅超声图像选取6060的感兴趣区域。之所以选择6060的感兴趣区域，是因为家兔的肝脏较狭窄，太大则超出了肝脏组织，甚至可能包含肾脏组织，太小了则不能反映出组织的纹理结果，灰度直方图参数：灰度均值，扭曲度，灰度标准差，灰度熵；灰度共生矩阵参数：共生矩阵能量，熵，惯性动量，相关性与局部稳定性；行程长度参数：长行程因子，短行程因子，灰度分布不均匀因子

47、，行程长度不均匀性因子等。第101页，共135页，编辑于2022年，星期五第102页，共135页，编辑于2022年，星期五第103页，共135页，编辑于2022年，星期五 第104页，共135页，编辑于2022年，星期五表表1 1 直方图参数直方图参数-灰度均值数据灰度均值数据编号 正常 10d 20d 30d 40d 50d 60d001 36.010 36.367 45.911 46.324 48.297 34.896 48.240002 33.487 47.390 47.983 52.845 48.001 56.550 -003 38.594 39.739 31.470 49.659 5

48、5.153 58.929 71.932004 37.021 39.692 43.463 50.861 57.131 59.832 58.806005 39.124 39.345 40.471 57.165 43.041 43.173 48.506006 38.172 32.916 39.160 48.645 45.423 46.245 49.018007 46.052 48.043 50.575 -008 38.977 45.246 51.047 56.205 -009 47.866 -010 40.617 41.742 51.137 43.931 -011 43.449 -注：表中-表示在家

49、兔实验中已死亡 ;第105页，共135页，编辑于2022年，星期五共生矩阵能量参数结果共生矩阵能量参数结果编号 正常 10d 20d 30d 40d 50d 60d001 0.105 0.148 0.189 0.193 0.121 0.207 0.220002 0.140 0.075 0.060 0.121 0.212 0.238 -003 0.094 0.127 0.185 0.126 0.164 0.251 0.255004 0.074 0.130 0.183 0.182 0.182 0.237 0.265005 0.075 0.100 0.188 0.160 0.153 0.151 0.

50、153006 0.121 0.150 0.188 0.190 0.194 0.228 0.232共生矩阵能量度量的是图像灰度均匀性，图像呈现较粗的纹理，该值相应较大。从各组检查的数据看出，随着脂肪性程度加重，其值呈现上升趋势，说明该参数对于病情变化有所反映。第106页，共135页，编辑于2022年，星期五T首先创建一幅包含两个相互交迭的圆形对象的二进制图像。定义两个圆的圆心和半径。为了产生交迭，首先创建一幅包含两个相互交迭的圆形对象的二进制图像。定义两个圆的圆心和半径。为了产生交迭，半径必须大于两圆心距离的一半。半径必须大于两圆心距离的一半。TCenter1=-10;TCenter2=-cen