graycomatrix 计算(图像)灰度共生矩阵(CLCM)

matlab函数: graycomatrix()

功           能：创建灰度共生矩阵

Gray-level co-occurrence matrix from an image

图像的灰度共生矩阵

灰度共生矩阵是像素距离和角度的矩阵函数，它通过计算图像中一定距离和一定方向的两点灰度之间的相关性，来反映图像在方向、间隔、变化幅度及快慢上的综合信息。

使用方法：
glcm = graycomatrix(I)
glcms = graycomatrix(I,param1,val1,param2,val2,...)
[glcms,SI] = graycomatrix(...)

描述：
glcms = graycomatrix(I) 产生图像I的灰度共生矩阵GLCM。它是通过计算两灰度值 i,j 在图像 I 中水平相邻的次数而得到的 (你也可以通过调整' Offsets' 参数来指定其它的像素空间关系），GLCM中的每一个元素（i,j）代表灰度 i 与灰度 j 在图像 I 中水平相邻的次数。

graycomatrix()先将图像 I 归一化到指定的灰度级，再计算GLCM；这是因为动态地求取图像的GLCM区间代价过高。如果I是一个二值图像，那么灰度共生矩阵就将图像转换到二值灰度级（黑和白）。如果I是一个灰度图像， 那将转换到8灰度级(默认)。灰度的级数决定了GLCM的大小尺寸，假设灰度级为L，则GLCM的尺寸是L x L。你可以通过设定参数“NumLevels”来指定灰度级数目，还可以通过设置“GrayLimits"参数来设置灰度共生矩阵的转换方式。

下图在一个4x5的图像I中显示了如何求解灰度共生矩阵，以（1，1）点为例，在图像 I 中水平相邻的像素对的灰度值都为1的情况只出现了1次，所以GLCM(1,1)的值是1。，同理，在图像 I 中水平相邻的像素对的灰度值分别为 1和2 的情况出现了2次，所以GLCM(1,2)的值是2。 graycomatrix迭代以上过程，就可以计算出GLCM的所有位置(L^2)的取值。

  

glcms = graycomatrix(I,param1,val1,param2,val2,...) 返回一个或多个灰度灰度共生矩阵，根据指定的参数对的值。参数可以简写，并且对大小写不敏感。

参数
下面按照字母的顺序列写了参数：
     'GrayLimits'  是两个元素的向量[low,high]，指明了图像 I 中的灰度值如何线性归一化到灰度级别。低于或等于low的灰度值置成1，大于或等于high的灰度值置成NumLevels。如果其设为[]，灰度共生矩阵将使用图像I的最小和最大灰度值分别作为GrayLimits的low和high，即[min(I(:) , max(I(:)))]。

     'NumLevels'    一个整数，指定灰度级的数目。例如，如果NumLevels为8，意思就是将图像I的灰度映射到1到8之间，它也决定了灰度共生矩阵的大小。默认值是8。

     'Offset'   一个p*2的整数矩阵，指定了感兴趣像素对之间的距离和方向。矩阵中的每一行是一个两元素的向量，[row_offset , col_offset]，它指定了一对像素之间的关系，或者说是位移。row_offset是感兴趣像素对间隔的行的数目；col_offset是感兴趣像素对间隔的列的数目。offset通常表示一个角度，下面列写的offset的值指定了常见角度。D代表是当前像素与邻居的距离。

Angle        Offset
  0              [0 D]
 45             [-D D]
 90             [-D 0]
135            [-D -D]
     下图说明了数组：offset = [0 1; -1 1; -1 0; -1 -1]

 

     'Symmetric'  一个布尔型数(逻辑型)，指定创建GLCM时像素对中的两像素的顺序是否考虑。例如，当 'Symmetric' 是true时，graycomatrix计算1连接2的次数时，(1,2)和(2,1)这两种数对都计算在内。当'Symmetric'是false时，graycomatrix只是计算(1,2)或(2,1).

[glcm,SI] = graycomatrix(....) 返回归一化(灰度级的)图像,SI,它被用来计算灰度共生矩阵(GLCM)，SI图像的取值范围是[1,NumLevels]。

支持类型

      I可以是数字型或逻辑型，但必须是二维的，实数的，非稀疏的矩阵。SI是一个double型矩阵，它和I的尺寸相同。glcms是一个‘NumLevels’ x ‘NumLevels’ x P的double型矩阵，P是offsets的数目(即‘Offset’参数值的列数)。
说明：

      灰度共生矩阵(GLCM)的另一个名字是灰度空间相关矩阵(gray-level spatial dependence matrix)。另一方面，co-occurrence在文献中使用时经常不带连字符，即cooccurrence。

     如果像素对中的一个像素值为NaN，graycomatrix忽略该像素对。

     graycomatrix用NumLevels值替代positive Inf，用1代替negative Inf。

      如果边界像素的邻居落在图像边界的外边，graycomatrix忽略该边界像素。

     当'Symmetric'设置成'true'时，GLCM 是关于对角线对称的，就是Haralick (1973)描述的GLCM。下面句法（1）使用'Symmetric'为'true'时创建了GLCM等于句法（2）和句法（3）使用'Symmetric'为‘false’时产生的GLCM的和。

 graycomatrix(I, 'offset', [0 1], 'Symmetric', true)    (1)

 g