特征提取

clear all;clc;
root_path = 'I:\PathologicalSection\PulmonarySquamousCellCarcinoma\OpticalParameters\CrossValidation1\WithoutAugment\Data';
ImagePath=[root_path, '\train_image\'];
Images = dir([ImagePath, '*.tif']);
PictureNums = length(Images);
FeatureNums = 5;
gray_comatrix=zeros(PictureNums,FeatureNums);
label_features=zeros(PictureNums,FeatureNums*3+1);

for i=1:PictureNums
    i
    name = Images(i).name;
    I=imread([ImagePath,name]);
    
    for j=1:3
        B=I(:,:,j);
        
        H=graycomatrix(B,'GrayLimits',[],'NumLevels', FeatureNums,'Offset',[0 1;-1 1;0 -1;-1 -1]);%得到灰度共生矩阵
        stats=new_graycoprops(H,'all');
        
        Energy = [];
        Energy= [Energy,stats.Energy];
        Energy=mean(Energy);
        
        Contrast=[];
        Contrast=[Contrast,stats.Contrast];
        Contrast=mean(Contrast);
        
        Correlation=[];
        Correlation=[Correlation,stats.Correlation];
        Correlation=mean(Correlation);
        
        Dissimilarity=[];
        Dissimilarity=[Dissimilarity,stats.Dissimilarity];
        Dissimilarity=mean(Dissimilarity);
        
        Entropy=[];
        Entropy=[Entropy,stats.Entropy];
        Entropy=mean(Entropy);
        
        Homogeneity=[];
        Homogeneity=[Homogeneity,stats.Homogeneity];
        Homogeneity=mean(Homogeneity);
        
        Mean=[];
        Mean=[Mean,stats.Mean];
        Mean=mean(Mean);
        
        Variance=[];
        Variance=[Variance,stats.Variance];
        Variance=mean(Variance);
        

        label_features(i,(j-1)*5+2) = Energy;
        label_features(i,(j-1)*5+3) = Contrast;
        label_features(i,(j-1)*5+4) = Correlation;
        label_features(i,(j-1)*5+5) = Entropy;
        label_features(i,(j-1)*5+6) = Mean;
        
    end
    
    %类别标签
    if(name(1)=='C')
        label = 1;
    else
        label = 0;
    end
    label_features(i,1) = label;
    
end

dataset_train = label_features;
% save('mat文件名', '变量名')
save('dataset_train', 'dataset_train');

%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all;clc;

ImagePath=[root_path, '\val_image\'];
Images = dir([ImagePath, '*.tif']);
PictureNums = length(Images);
FeatureNums = 5;
gray_comatrix=zeros(PictureNums,FeatureNums);
label_features=zeros(PictureNums,FeatureNums*3+1);

for i=1:PictureNums
    i
    name = Images(i).name;
    I=imread([ImagePath,name]);
    
    for j=1:3
        B=I(:,:,j);
        
        H=graycomatrix(B,'GrayLimits',[],'NumLevels', FeatureNums,'Offset',[0 1;-1 1;0 -1;-1 -1]);%得到灰度共生矩阵
        stats=new_graycoprops(H,'all');
        
        Energy = [];
        Energy= [Energy,stats.Energy];
        Energy=mean(Energy);
        
        Contrast=[];
        Contrast=[Contrast,stats.Contrast];
        Contrast=mean(Contrast);
        
        Correlation=[];
        Correlation=[Correlation,stats.Correlation];
        Correlation=mean(Correlation);
        
        Dissimilarity=[];
        Dissimilarity=[Dissimilarity,stats.Dissimilarity];
        Dissimilarity=mean(Dissimilarity);
        
        Entropy=[];
        Entropy=[Entropy,stats.Entropy];
        Entropy=mean(Entropy);
        
        Homogeneity=[];
        Homogeneity=[Homogeneity,stats.Homogeneity];
        Homogeneity=mean(Homogeneity);
        
        Mean=[];
        Mean=[Mean,stats.Mean];
        Mean=mean(Mean);
        
        Variance=[];
        Variance=[Variance,stats.Variance];
        Variance=mean(Variance);
        

        label_features(i,(j-1)*5+2) = Energy;
        label_features(i,(j-1)*5+3) = Contrast;
        label_features(i,(j-1)*5+4) = Correlation;
        label_features(i,(j-1)*5+5) = Entropy;
        label_features(i,(j-1)*5+6) = Mean;
        
    end
    
    %类别标签
    if(name(1)=='C')
        label = 1;
    else
        label = 0;
    end
    label_features(i,1) = label;
    
end

dataset_val = label_features;
% save('mat文件名', '变量名')
save('dataset_val', 'dataset_val');


%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all;
clc;

ImagePath=[root_path, '\test_image\'];
Images = dir([ImagePath, '*.tif']);
PictureNums = length(Images);
FeatureNums = 5;
gray_comatrix=zeros(PictureNums,FeatureNums);
label_features=zeros(PictureNums,FeatureNums*3+1);

for i=1:PictureNums
    i
    name = Images(i).name;
    I=imread([ImagePath,name]);
    
    for j=1:3
        B=I(:,:,j);
        
        H=graycomatrix(B,'GrayLimits',[],'NumLevels', FeatureNums,'Offset',[0 1;-1 1;0 -1;-1 -1]);%得到灰度共生矩阵
        stats=new_graycoprops(H,'all');
        
        Energy = [];
        Energy= [Energy,stats.Energy];
        Energy=mean(Energy);
        
        Contrast=[];
        Contrast=[Contrast,stats.Contrast];
        Contrast=mean(Contrast);
        
        Correlation=[];
        Correlation=[Correlation,stats.Correlation];
        Correlation=mean(Correlation);
        
        Dissimilarity=[];
        Dissimilarity=[Dissimilarity,stats.Dissimilarity];
        Dissimilarity=mean(Dissimilarity);
        
        Entropy=[];
        Entropy=[Entropy,stats.Entropy];
        Entropy=mean(Entropy);
        
        Homogeneity=[];
        Homogeneity=[Homogeneity,stats.Homogeneity];
        Homogeneity=mean(Homogeneity);
        
        Mean=[];
        Mean=[Mean,stats.Mean];
        Mean=mean(Mean);
        
        Variance=[];
        Variance=[Variance,stats.Variance];
        Variance=mean(Variance);
        

        label_features(i,(j-1)*5+2) = Energy;
        label_features(i,(j-1)*5+3) = Contrast;
        label_features(i,(j-1)*5+4) = Correlation;
        label_features(i,(j-1)*5+5) = Entropy;
        label_features(i,(j-1)*5+6) = Mean;
        
    end
    
    %类别标签
    if(name(1)=='C')
        label = 1;
    else
        label = 0;
    end
    label_features(i,1) = label;
    
end

dataset_test = label_features;
% save('mat文件名', '变量名')
save('dataset_test', 'dataset_test');