graycomatrix函数生成灰度共生矩阵

灰度共生矩阵的接触知识，http://www.cnblogs.com/tiandsp/archive/2013/03/23/2977673.html

初学者再看如此多的公式，可能看也看不懂。那就简单的办法就是直接看图（matlab help文件中的）：

刚开始看的人，肯定也摸不着头脑。

我这里解释下两边矩阵的含义：

（1）左边的矩阵，你可以想成某图像的像素（当然不是rgb的，可以认为是灰度图）

（2）右边的矩阵（定义为GLCM）：是生成的灰度共生矩阵。关于维数：这里是8*8的，但实际情况中，你可以是2*2,4*4，16*16等等（在matlab中，如果你是二值图形，维数默认为2，如果是灰度图，默认为8）。

    现在再讲讲GLCM里面数字的含义，如GLCM（1,1）就是从左往右水平的像素是分别都是1,1的次数。从左边的矩阵可以看出，只有一次，同理GLCM（1,2）=2，因为有从左往右的像素分别是1,2的情况有两种，再比如GLCM（2,2）=0，因为没有从左往右像素分别是2,2的情况。（这里的从左往右是相邻的）。

   我相信在这里，初学者都明白了，GLCM的如何生成的吧。

（1）下面讲讲关于matlab自带的graycomartix函数吧（原理其实很简单，你可以自己写个函数）：

我上面讲的例子，其实可以只用一个函数就可以完成。

gray = [1 1 5 6 8
        2 3 5 7 1
        4 5 7 1 2
        8 5 1 2 5];
GLCM = graycomatrix(gray,'GrayLimits',[]);

（注意，必须用'GrayLimtis',[]）,不然你得到的GLCM是很奇葩的

GLSM =
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0    16
原因是，GrayLimits参数的默认值是改种类型的范围，如double的范围[0 1]   int16的范围[-32768 32767]

（2）另外，关于graycomatrix函数，你可以写成

GLCM = graycomatrix(gray,'G',[]);

因为不论参数不论是‘G’还是‘Graycomatrix’还是‘g’都可以被识别的。

（3）最后不要混淆，GLCM的矩阵大小和你原始的gray是不一样的。（GLCM的大小只和参数‘graylimits’有关！）

[GLCM, SI] = graycomatrix(gray,'G',[]);  这里的得到的矩阵SI才是和gray矩阵一样的。

SI的解释：对矩阵gray的像素进行分类，如果你设的参数“NumLvels”是8，那么灰度共生矩阵就将gray线性分类成8类，对应的类的数值（1-8）就是SI矩阵。

for example：

gray = [ 1 1 5 6 8 8; 2 3 5 7 0 2; 0 2 3 5 6 7];
[glcm,SI] = graycomatrix(gray,'NumLevels',9,'G',[])

得到的结果：

glcm =


     0     0     2     0     0     0     0     0     0
     0     1     0     0     0     1     0     0     0
     0     0     0     2     0     0     0     0     0
     0     0     0     0     0     2     0     0     0
     0     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     2     1     0
     0     0     0     0     0     0     0     1     1
     1     0     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0     1

SI =


     2     2     6     7     9     9
     3     4     6     8     1     3
     1     3     4     6     7     8

gray就是3*6的矩阵，SI就是gray对应的分类，因为参数“NumLevels”设为9，所以SI的数值就是1-9，大小和gray一样。而glcm和和前面讲的一样！

关于Graycomatrix函数的其他参数：
NumLevels:(numeric）数值的默认值是8 ， (binary)默认是是2，我们上面的例子就是默认值为8，所以产生的GLSM是8*8的矩阵。
OffSet：表明像素对比的顺序。     [0 D]:表示从左往右像素距离为D（如果D=1，就是相邻）
symmetric:默认值为false，比如从左往右就是从左往右，而当该参数设定为true是，从左往右和从右往左都被算上的。
更多的参数表示可以看help文件吧。这里贴个截图：