灰度共生矩阵提取图像特征matlab代码

% 灰度共生矩阵提取图像特征代码

% 读入图像
img = imread('lena.jpg');
% 转换为灰度图像
gray = rgb2gray(img);

% 定义灰度共生矩阵的方向和距离
directions = [0, 1; -1, 1; -1, 0; -1, -1];
distances = [1, 2, 3, 4];

% 初始化灰度共生矩阵
glcm = zeros(256, 256, length(distances), length(directions));

% 计算灰度共生矩阵
for d = 1:length(distances)
for dir = 1:length(directions)
glcm(:, :, d, dir) = graycomatrix(gray, 'NumLevels', 256, 'Offset', distances(d)*directions(dir, :));
glcm(:, :, d, dir) = glcm(:, :, d, dir) / sum(glcm(:, :, d, dir), 'all');
end
end

% 计算灰度共生矩阵的统计特征
stats = graycoprops(glcm, {'Contrast', 'Correlation', 'Energy', 'Homogeneity'});

% 显示灰度共生矩阵的统计特征
disp(stats);