灰度-梯度共生矩阵数字特征的研究

1、PAGE PAGE 9灰度-梯度共生矩阵数字特征的研究0引言灰度共生矩阵概念最早由Haralick于1973年提出，在1992年Ohanian P. P.通过实验证明了基于灰度共生矩阵的统计特征的有效性。洪继光于1984年在灰度共生矩阵的基础上提出了灰度-梯度共生矩阵1，并应用该模型上的15个特征对五类白血球样本进行了分类识别，其实验结果表明，对于像细胞核边界不清晰的图像，该方法较灰度共生矩阵的分类结果好。灰度-梯度共生矩阵模型集中反映了图像中两种最基本的信息，即像素的灰度和梯度（或边缘）的相互关系。各像素的灰度是构成一幅图像的基础，而梯度则是构成图像边缘轮廓的要素。灰度-梯度空间可以很清晰地

2、表现图像内像素灰度与梯度的分布规律，同时也体现了各像素与其邻域像素的空间关系，对图像的纹理能很好地描绘。1灰度-梯度共生矩阵模型灰度-梯度共生矩阵纹理分析方法是利用图像的灰度和梯度的综合信息提取纹理特征2。灰度-梯度共生矩阵的元素定义为在归一的灰度图像及其归一的梯度图像中具有灰度值和梯度值的像素数，即在集合中元素的个数。其中，。对灰度-梯度共生矩阵进行归一化处理，使其各元素之和为1。如式（1）所示：（1）而，所以上式可以表示为式（3-2）：（2）该灰度-梯度共生矩阵的原点在左上角，向右梯度值增加，向下灰度值增加。对于粗纹理的图像，在灰度轴附近集中分布，仅少量边界点远离灰度轴分布。对于细纹理图像

3、，离开灰度轴而沿梯度轴散开分布3。灰度-梯度共生矩阵常用的数字特征有：（1）小梯度优势（2）大梯度优势（3）灰度分布的不均匀性（4）梯度分布的不均匀性（5）能量（6）灰度平均（7）梯度平均（8）灰度均方差（9）梯度均方差（10）相关（11）灰度熵（12）梯度熵（13）混合熵（14）惯性（15）逆差矩因梯度表达了图像灰度的最大变化率，图像各点梯度的大小则表达了该点处等灰度线的密集程度，所以大梯度优势和小梯度优势可从一定程度上反映图像的灰度变化剧烈程度。当图像的灰度变化平缓时，小梯度像素数多，小梯度优势大。而当图像的灰度变化剧烈时，大梯度像素数多，大梯度优势大。2 仿真实验及分析因为灰度-梯度共生

4、矩阵是定义于多灰度级图像上的，故采用灰度化预处理方法将半调图像“下采样”变换成灰度图像。具体做法如图1所示，先将半调图像按分块，并依光栅扫描这4个像素，将其分别记为a、b、c、d，再按式（3-3）便可将的二值块映射成一个十进制。这样，原来仅有两级灰度值的半调图像就缩小成原图像1/4大小的16级灰度图像。（3）图1 灰度化处理示意图灰度-梯度共生矩阵是在归一化灰度图像及其归一化梯度图像上构造的，将灰度图像和梯度图像均归一化为16级，而应用灰度化预处理后半调图像就变成16级的灰度图像，因此在这里不需再进行灰度归一化。而由半调图像获得对应梯度图像，则是利用sobel算子，该算子包含两组3x3的矩阵，

5、分别为横向及纵向，将之与半调图像作平面卷积，即可分别得出横向及纵向的亮度差分近似值。如果以A代表原始图像，Gx及Gy分别代表经横向及纵向边缘检测的图像，其公式如下:那么，最终求得的梯度图像可用下面的公式算出：（4）对于利用sobel算子求得梯度图像可通过式（5）求得归一化梯度图像。然后，便可通过1中的数字特征计算公式求得特征值。（5）其中，。最后，使用公式（2）便得到灰度-梯度共生矩阵。实验选用一幅lena_hFloyd 256256半调图作为实验对象。见图2。程序使用C+语言实现算法。最终计算出的15个数字特征值如下：小梯度优势 0.0435289大梯度优势 71.8826灰度分布的不均匀性

6、 6.64258梯度分布的不均匀性 6.29035能量 0.0229605灰度平均 5.97595 图2 实验用半调图梯度平均 7.74176灰度均方差 6.21491梯度均方差 3.45656相关 5.3755灰度熵 0.250095梯度熵 0.373588混合熵 1.86179惯性 1.1107e+019逆差矩 0.0931894该结果各项数据经与使用Matlab处理的结果完全一致。3 总结本次SSRT计划项目，我熟悉了BMP文件的文件物理与逻辑结构，明白了什么是灰度图、梯度图、半调图，学习了运用编程语言实现灰度-梯度共生矩阵数字特征计算的重要算法与核心思想。信控学院的孔月萍老师在此期间给与了我很多必要的支持，使得我得以克服困难，完成此次SSRT实验计划，在此我衷心地表达我的感谢。参考文献1 洪继光, 灰度一梯度共生矩阵纹理分析方法. 自动化学报,Jan 1984. Vol. 10 No.12 夏德深,金盛,王健, 基于分数维与灰度梯度共生矩