1

我有一组 16 位的图像数据集，我想从中创建 GLCM 矩阵以提取 GLCM 特征。

但是，生成的矩阵显示一个值（如下图所示），我想知道为什么。

16glcm

我尝试使用相同的图像但转换为 8 位，结果 GLCM 显示几个值。

8glcm

注意：我使用了以下 Matlab 函数：

glcm_matrix = graycomatrix(image.tif);

这是 16 位图像的裁剪示例：

rescaled_crop

注意：计算中使用的图像可以从这里下载。原始图像的对比度非常低，看起来很暗。上面显示的图像的对比度被拉伸，仅用于可视化目的。

编辑：

我用了

glcm_matrix = graycomatrix(image.tif, 'GrayLimits', []);

它给了我以下结果：

16glcm固定

4

2 回答 2

3

这是一个分箱/缩放问题。

让我们看看里面：

edit graycomatrix

在这种情况下，我们对“NumLevels”和“GrayLimits”这两个选项感兴趣

%   'NumLevels'      An integer specifying the number of gray levels to use
%                    when scaling the grayscale values in I. For example,
%                    if 'NumLevels' is 8, GRAYCOMATRIX scales the values in
%                    I so they are integers between 1 and 8.  The number of
%                    gray levels determines the size of the gray-level
%                    co-occurrence matrix (GLCM).
%
%                    'NumLevels' must be an integer. 'NumLevels' must be 2
%                    if I is logical.
%  
%                    Default: 8 for numeric
%                             2 for logical
%
%   'GrayLimits'     A two-element vector, [LOW HIGH], that specifies how 
%                    the values in I are scaled into gray levels. If N is
%                    the number of gray levels (see parameter 'NumLevels')
%                    to use for scaling, the range [LOW HIGH] is divided
%                    into N equal width bins and values in a bin get mapped
%                    to a single gray level. Grayscale values less than or
%                    equal to LOW are scaled to 1. Grayscale values greater
%                    than or equal to HIGH are scaled to NumLevels. If
%                    'GrayLimits' is set to [], GRAYCOMATRIX uses the
%                    minimum and maximum grayscale values in I as limits,
%                    [min(I(:)) max(I(:))].

因此，换句话说，该函数将您的数据分箱到 8x8 箱中，并假设缩放范围是完整的 uint16 范围 (0-65535)。但是，您提供的示例图像最小为 305，最大为 769，使其落入第一个 bin（0-8192 左右）。当我调用A = graycomatrix(I)它时，它给了我以下矩阵：

A =

    6600           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

然而，当A = graycomatrix(I,'GrayLimits', [])被调用时，缩放范围被视为 min(I) - max(I)，并且函数按预期工作：

A =
       4           2           1           0           0           0           0           0
       1           1           2           2           0           0           0           0
       2           2           4           7           1           0           0           0
       0           1           7         142          72           1           0           0
       0           0           0          65        1711         252           0           0
       0           0           0           0         230        3055         178           0
       0           0           0           0           0         178         654           8
       0           0           0           0           0           0           8           9

在您的原始示例中，单个值很可能位于 8x8 矩阵的中间，因为您的原始图像是 int16 而不是 uint16，因此 graycomatrix 是对称的，以考虑到负值的可能性。

当然，您也可以缩放原始图像以适合其数据类型。例如，如果您期望异常值等，百分位缩放可能是一个好主意。

于 2019-02-07T18:42:56.993 回答
2

我只想以@Tapio 的出色答案为基础。

graycomatrix当您在函数调用中使用名称/值对时产生的 GLCMGrayLimits', []看起来不错。但是，这种方法可能对您的应用程序无效。如果以这种方式计算一组图像的 GLCM，则对应于两个不同图像的两个不同 GLCM 的相同元素可能具有不同的含义。实际上，由于每个图像的强度重新缩放不同，GLCM 的组件实际上是在对从一个图像到另一个图像的不同共现进行编码。

为避免这种情况，您可以首先计算整个图像数据集的最小和最大强度（例如minImgs和maxImgs），然后使用这些值以完全相同的方式重新调整构成数据集的所有图像的强度：

glcm_matrix = graycomatrix(image_tif, 'GrayLimits', [minImgs maxImgs]);
于 2019-02-09T19:29:26.973 回答