利用聚类技术实现纹理图像分割

MATLAB实现的利用聚类技术实现纹理图像分割，纹理特征提取采用的是灰度共生矩阵，聚类采用的是k-means方法。

% Get the texture feature vector from the image given.
%
% Usage:  texture = get_texture_features（ image_name ）
%
% Parameters:
%   image_name    the name of the image file like ‘Test_Img_1.jpg‘
%
% Output:
%   texture         the texture vector of each pixel whose size is N*P where N is the total number of the pixels column by column 
%
% This code was written by Liu Ruoyu. Contact: 12112062@bjtu.edu.cn
% Last change: December 11th 2012
%

function texture = get_texture_features（ image_name ）

        window_size = 35;
        dir_image = [‘./data/‘ image_name];
        I = imread（ dir_image ）;
        [rsize csize] = size（I）;
        texture = [];
        
        padding_size = fix（ window_size/2 ）;
        I_padding = padarray（ I [ padding_size padding_size] ‘replicate‘ ‘both‘ ）;    % pad the boundary of the image use the value of the edge pixel
        
        for c = 1:csize
            for r = 1:rsize
                    % get the texture column by column because the fuction ‘reshape‘ will put the elements into one column first
                
                    slidWin = I_padding（ r : r+window_size-1  c : c+window_size-1 ）;
                    glcm = graycomatrix（ slidWin ‘NumLevels‘ 16 ‘Offset‘ [0 1] ）;
                    stats = graycoprops（ glcm ‘Contrast Correlation Energy ‘ ）;
                    texture_v = [ stats.Contrast stats.Correlation stats.Energy ];
                    texture = [ texture ; texture_v ];
                    
            end
        end
        

end

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

     文件      49708  2012-11-28 15:07  Texture_mosica\data\Texture_mosaic_1.jpg

     文件      43587  2012-11-28 15:08  Texture_mosica\data\Texture_mosaic_2.jpg

     文件      41155  2012-11-28 15:11  Texture_mosica\data\Texture_mosaic_3.jpg

     文件      47552  2012-11-28 15:13  Texture_mosica\data\Texture_mosaic_4.jpg

     文件       1579  2012-12-11 19:34  Texture_mosica\demo\get_texture_features.m

     文件       1169  2012-12-11 19:38  Texture_mosica\demo\segment_based_on_grouping.m

     文件        581  2012-12-11 19:38  Texture_mosica\demo\segment_based_on_grouping_main.m

     文件      23406  2012-12-11 15:28  Texture_mosica\result\result_1_sqEuclidean.jpg

     文件      20917  2012-12-11 15:29  Texture_mosica\result\result_2_sqEuclidean.jpg

     文件      21231  2012-12-11 15:30  Texture_mosica\result\result_3_sqEuclidean.jpg

     文件      22610  2012-12-11 18:42  Texture_mosica\result\result_3_sqEuclidean_5_parts.jpg

     文件      23495  2012-12-11 15:33  Texture_mosica\result\result_4_sqEuclidean.jpg

     目录          0  2012-12-09 16:28  Texture_mosica\data

     目录          0  2012-12-11 19:38  Texture_mosica\demo

     目录          0  2012-12-11 19:26  Texture_mosica\result

     目录          0  2012-12-11 19:26  Texture_mosica

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

               296990                    16