MATLAB求灰度共生矩阵特征的代码

MATLAB中求灰度共生矩阵特征的代码如下：

% 读取灰度图像
I = imread('image.png');
% 将图像转化为灰度图像
I = rgb2gray(I);
% 灰度级数
grayLevels = 256;

% 定义灰度共生矩阵
glcm = zeros(grayLevels, grayLevels);

% 定义方向和距离
offsets = [0 1; -1 1; -1 0; -1 -1];
distances = [1];

% 计算灰度共生矩阵
for i = 1:size(offsets, 1)
    glcm = glcm + graycomatrix(I, 'Offset', offsets(i, :), 'NumLevels', grayLevels, 'GrayLimits', [], 'Symmetric', true, 'Distance', distances);
end

% 将灰度共生矩阵归一化
glcm = glcm ./ sum(glcm(:));

% 定义灰度共生矩阵特征
contrast = 0;
dissimilarity = 0;
homogeneity = 0;
energy = 0;
entropy = 0;

% 计算灰度共生矩阵特征
for i = 1:grayLevels
    for j = 1:grayLevels
        contrast = contrast + (i - j)^2 * glcm(i, j);
        dissimilarity = dissimilarity + abs(i - j) * glcm(i, j);
        homogeneity = homogeneity + glcm(i, j) / (1 + abs(i - j));
        energy = energy + glcm(i, j)^2;
        if glcm(i, j) ~= 0
            entropy = entropy - glcm(i, j) * log2(glcm(i, j));
        end
    end
end

在这个代码中，我们首先读入一张彩色图像，然后将其转化为灰度图像。接着我们定义了灰度共生矩阵的大小和方向，然后使用 `graycomatrix` 函数计算灰度共生矩阵，并将其归一化处理。最后，我们定义了灰度共生矩阵的特征，包括对比度(contrast)、不相似度(dissimilarity)、同质性(homogeneity)、能量(energy)和熵(entropy)，并计算这些特征。