提取图像特征，包括灰度梯度，灰度距等相关特征

 

%% 
while 1==1
    choice=menu('Texture Analysis',...
                'Select Directory',...
                'Select Measure File',...
                'Co-comatrix',...
                'Exit');
    if choice == 1
        disp('Read from Directory');
        inputdirectory_name = uigetdir
        allfiles = dir(inputdirectory_name);
    end
%%
    if choice == 2
        disp('Write to a measure files');
        [FileName,PathName] = uiputfile('*.mes');
        outputfilename=[PathName,FileName];
    end
%%
 if choice == 3
     %Setup measurement file
        file = fopen(outputfilename,'w+');
        fprintf(file,'Texture Measture File \t\n');
        offset=[0 1;-1 1;-1 0;-1 -1];
        fprintf(file,'the order of the statistics is 0° 45° 90° 135°\n');
        count=1;
        for fileindex=1:length(allfiles)
          tmpfile=allfiles(fileindex);
          if tmpfile.isdir
              continue;
          end         
          %Find image name
          tempfile = strcat(inputdirectory_name,'\');
          filename = strcat(tempfile,tmpfile.name);
          Image = imread(filename);
          % Calculate the co-comatrix
          
                 
          glcms=graycomatrix(Image,'Offset',offset);
          % Calculate the chateratics of co-matrix
          stats=graycoprops(glcms,{'contrast','Correlation','Energy','Homogeneity'});
          %put the statistics into file.
          fprintf(file,'%s\n',filename);
          fprintf(file,'%f\t%f\t%f\t%f\t\n',stats.Contrast,stats.Correlation,stats.Energy,stats.Homogeneity);
          disp(count);disp('is OK');count=count+1;
          
        end
        disp('all is OK!');
      fclose(file);
 end
   if choice == 4
        save texturestructure;
        clear all;
        clc;
        close all;
        return;
    end
end

D00007