图像纹理分析灰度共生矩阵matlab源码

function  finalmat = batchglcm（newdirfilenamenmlabelDtxtname）
%newdir为需要处理的图像所在的文件夹路径，filename是图像名字，为字符串，n是图像总数，m为特征向量的维数
if ~exist（newdir‘dir‘） 
mkdir（newdir）; 
end 
paths=[pwd‘\‘newdir‘\‘];
labelmat = zeros（1n）;
labelmat（1:）=label;
idmat = 1:1:n;
featuremat = zeros（mn）;
finalmat = zeros（m+2n）;
for k=1:n; 
    img = imread（[pathsfilenamenum2str（k）‘.jpg‘]）;
    glcms1 = graycomatrix（img‘numlevels‘64‘offset‘[0 D;-D D;-D 0;-D -D]）;
    stats = graycoprops（glcms1{‘contrast‘‘correlation‘‘energy‘‘homogeneity‘}）;
    ga1 = glcms1（::1）;
    ga2 = glcms1（::2）;
    ga3 = glcms1（::3）;
    ga4 = glcms1（::4）;
    energya1 = 0;energya2 = 0;energya3 = 0;energya4 = 0;
    for i = 1:64
        for j = 1:64
            energya1 = energya1 + sum（ga1（ij