文档

graycomatrix

从图像创建灰度共现矩阵

语法

灰度矩阵(I)
glcms = graycomatrix(I,Name,Value,…)
[glcms,SI] = graycomatrix(___

描述

glcms= graycomatrix (从图像中创建灰度共现矩阵(GLCM).灰度共现矩阵的另一个名称是a灰度空间依赖矩阵.同样,在没有连字符的文献中,单词cooccurrence经常被使用。

graycomatrix通过计算像素的灰度(灰度强度)值的频率来创建GLCM在具有值的像素附近水平发生j.属性指定其他像素空间关系“补偿”参数)。每个元素(我,我)glcm指定像素与值相匹配的次数发生在具有值的像素的水平附近j

例子

glcms= graycomatrix(我名称,值,……)返回一个或多个灰度共现矩阵,具体取决于可选名称/值对的值。参数名可以缩写,大小写无关紧要。

例子

glcms如果] = graycomatrix(___返回缩放后的图像,如果,用于计算灰度共现矩阵。中的值如果之间1而且NumLevels

例子

全部折叠

将灰度图像读入工作区。

I = imread(“circuit.tif”);imshow(我)

计算灰度图像的灰度共生矩阵(GLCM)。默认情况下,graycomatrix根据像素的水平接近度计算GLCM:[0 1]。这是同一行中相邻的像素。这个例子指定了一个不同的偏移量:同一列上相隔两行。

glcm =灰度矩阵(I,“抵消”[2 0])
glcm =8×814205 2107 126 00 00 0⋯2242 14052 3555 400 00 00 191 3579 7341 1505 37 00 00 0 683 1446 7184 1368 00 00 7 116 1502 10256 1124 00 00 00 2 1153 1435 00 00 00 00 00 00 00 00 00 00 00 00 0

创建一个简单的3 × 6示例数组。

I = [1 1 5 6 8 8;2 3 5 7 0 2;0 2 3 5 6 7]
我=3×61 1 5 6 8 8 2 3 5 7 0 2 0 2 3 5 6 7

计算灰度共生矩阵(GLCM)并返回用于计算的缩放图像。方法指定空括号GrayLimits参数,示例使用输入图像中的最小和最大灰度值作为限制。

[glcm,SI] = graycomatrix(I,“NumLevels”9“GrayLimits”[])
glcm =9×90 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
如果=3×62 2 6 7 9 9 3 4 6 8 1 3 1 3 4 6 7 8

将灰度图像读入工作区。

I = imread(“cell.tif”);imshow(我)

定义四个偏移量。

偏移量= [0 1;-1 0;-1 1];

计算glcm,并返回缩放后的图像。显示缩放后的图像,将数据值额外重新缩放到范围[0,1]。

[glcms,SI] = graycomatrix(I,“抵消”、补偿);imshow(重新调节(SI))

注意该函数如何返回一个包含四个glcm的数组。

名称大小字节类别属性I 159x191 30369 uint8 SI 159x191 242952 double glcms 8x8x4 2048 double offset 4x2 64 double

将灰度图像读入工作区。

I = imread(“circuit.tif”);imshow(我)

方法计算GLCM对称的选项,返回缩放后的图像。设置时创建的GLCM对称的真正的在对角线上是对称的,等价于Haralick(1973)描述的GLCM。

[glcm,SI] = graycomatrix(I,“抵消”(2 0),“对称”,真正的);glcm
glcm =8×828410 4349 317 0 0 0 0 0⋯4349 28104 7134 1083 7 0 0 0 0 0 0 0 1083 317 7134 14682 2951 153 2951 14368 2870 2 0 0 0 7 153 2870 20512 2277 0 0 0 0 0 2 2277 2870 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

显示缩放后的图像,将数据值额外重新缩放到范围[0,1]。

imshow(重新调节(SI))

输入参数

全部折叠

输入图像,指定为二维、实数、非稀疏、数字或逻辑数组。

例子:

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

名称-值对参数

的可选逗号分隔对名称,值参数。的名字参数名称和价值对应的值。的名字必须出现在单引号内(' ').您可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家

例子:

全部折叠

Range使用缩放输入图像到灰度级别,指定为两元素向量(低高).如果N灰度等级的数量(见参数“NumLevels”)用于缩放,范围(低高)分为N等宽的箱子和值在一个箱子被映射到一个单一的灰度。灰度值小于或等于缩放到1。灰度值大于或等于被缩放到“NumLevels”.如果“GrayLimits”设置为[]graycomatrix中的最小和最大灰度值的限制,我(:)[min()马克斯(我(:))),例如,[0 1]对于类double和(-32768 32767)为类int16

例子:

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

灰度级别数,以整数形式指定。例如,如果NumLevels是8,graycomatrix缩放中的值它们是1到8之间的整数。灰度的数目决定了灰度共生矩阵的大小(glcm).

例子:

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

感兴趣的像素与其邻居之间的距离,指定为ap-by-2整数数组。数组中的每一行都是两个元素的向量,[row_offset, col_offset],指定关系,或者抵消,表示一对像素。row_offset感兴趣的像素与其相邻像素之间的行数。col_offset感兴趣的像素与其相邻像素之间的列数。由于偏移量通常表示为一个角度,下表列出了指定公共角度的偏移量值,给定像素距离D。

抵消

0

[0 D]

45

[- D D]

90 (0 - d)
135 [- d - d]

下图说明了数组:偏移量= [0 1;1 1;1 0;1 (1)

例子:

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

考虑值的顺序,指定为布尔值真正的.例如,当“对称”设置为真正的graycomatrix在计算值1与值2相邻的次数时,同时计算1,2和2,1对。当“对称”设置为graycomatrix只计算1,2或2,1,取决于的值“抵消”

例子:

数据类型:逻辑

输出参数

全部折叠

灰度共现矩阵(或矩阵),返回为NumLevels——- - - - - -NumLevels——- - - - - -P类数组,在那里P偏移量在吗抵消

用于计算GLCM的缩放图像,返回为a矩阵的大小与输入图像相同。

算法

graycomatrix从图像的缩放版本计算GLCM。默认情况下,如果是二值图像,graycomatrix将图像缩放为两个灰度级。如果是一个强度图像,graycomatrix将图像缩放到8个灰度级别。可以指定灰度的个数graycomatrix控件用于缩放图像“NumLevels”参数,以及graycomatrix属性缩放值“GrayLimits”参数。

如下图所示graycomatrix计算4 × 5图像的GLCM中的几个值.GLCM中的元素(1,1)包含该值1因为在图像中只有一个实例,其中两个水平相邻的像素具有值1而且1.元素(1、2)在GLCM中包含该值2因为在图像中有两个实例,其中两个水平相邻的像素具有值1而且2graycomatrix继续该处理,以填充GLCM中的所有值。

graycomatrix忽略像素对,如果其中一个像素包含,替换为正用值NumLevels,并替换为负用值1graycomatrix忽略边界像素,如果对应的邻居像素落在图像边界之外。

GLCM创建于何时“对称”设置为真正的在对角线上是对称的,等价于Haralick(1973)描述的GLCM。由以下语法生成的GLCM“对称”设置为真正的

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

等于下列表述产生的两个glcm的和,在哪里“对称”设置为

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

参考文献

[1] Haralick, r.m., K. Shanmugan,和I. Dinstein,“图像分类的纹理特征”,IEEE系统、人与控制论学报,卷SMC-3, 1973,第610-621页。

[2]哈拉利克,r.m.和L.G.夏皮罗。计算机与机器人视觉:第1卷,Addison-Wesley, 1992,第459页。

另请参阅

R2006a之前介绍

这个话题有用吗?