基于灰度共生矩阵的图像分割

基于灰度共生矩阵的图像分割－Matlab算法

clear;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   n=1 3*3窗口
%   n=2 5*5窗口
%   n=3 7*7窗口
%  Contrast 对比度; Correlation 相关性; Energy 角二阶矩;  Homogeneity 熵
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = 3;
% numlevel为灰度级别
numlevel = 8;
% 步长
D = 1;
% 读取灰度图像
gray = imread（‘aa.bmp‘）;
gsize = size（gray）;
graycopy = zeros（gsize（1）+n*2gsize（2）+n*2）;
gcsize = size（graycopy）;
graycopy（1+n:gcsize（1）-n1+n:gcsize（2）-n） = gray;
% 填充行
for index=1:n
    graycopy（index:） = graycopy（n+1:）;
    graycopy（gcsize（1）-（index-1）:） = graycopy（gcsize（1）-n:）;
end
% 填充列
for index = 1:n
    graycopy（:index） = graycopy（:n+1）;
    graycopy（:gcsize（2）-（index-1）） = graycopy（:gcsize（2）-n）;
end

%- 定义特征值矩阵 -%
Contrast = zeros（gsize（1）gsize（2））;
Correlation = zeros（gsize（1）gsize（2））;
Energy = zeros（gsize（1）gsize（2））;
Homogeneity = zeros（gsize（1）gsize（2））;

for row = 1:gsize（1）
    for col = 1:gsize（2）
        copyRow = row + n;
        copyCol = col + n;
        newMatrix = graycopy（copyRow-n:copyRow+ncopyCol-n:copyCol+n）;
        GLCM = graycomatrix（newMatrix‘GrayLimits‘[min（newMatrix（:）） max（newMatrix（:））]...
            ‘NumLevels‘numlevel‘Offset‘[0 D;-D D;-D0;-D -D]‘Symmetric‘true）;
        stats = graycoprops（GLCM）;
        Contrast（rowcol） = mean（stats.Contrast）;
        Correlation（rowcol） = mean（stats.Correlation）;
        Energy（rowcol） = mean（stats.Energy）;
        Homogeneity（rowcol） = mean（stats.Homogeneity）;
    end
end
% save（‘7.mat‘‘Contrast‘‘Correlation‘‘Energy‘‘Homogeneity‘）;
% % 将图像转换为uint8图，并进行直方图均衡化
% Contrast_uint = uint8（Contrast）;
% Contrast_uint = histeq（Contrast_uint）;
% Correlation_uint = uint8（Correlation）;
% Correlation_uint = histeq（Correlation_uint）;
% Energy_uint = uint8（Energy）;
% Energy_uint = histeq（Energy_uint）;
% Homogeneity_uint = uint8（Homogeneity）;
% Homogeneity_uint = histeq（Homogeneity_uint）;
% % 保存图像
% imwrite（Contrast_uint‘Contrast.tif‘‘tiff‘）;
% imwrite（Correlation_uint‘Correlation.tif‘‘tiff‘）;
% imwrite（Energy_uint‘Energy.tif‘‘tiff‘）;
% imwrite（Homogeneity_uint‘Homogeneity.tif‘‘tiff‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2174  2011-01-03 15:14  f.m

----------- ---------  ---------- -----  ----

                 2174                    1