%calculate features
clear
cohort = '';
path = '';

seg_num = [];
CT_phase = 'A';%A/V/p

total = dir([path '\' cohort '\' CT_phase]);
for i=3:1:length(total)
    load([path '\' cohort '\' CT_phase '\' total(i).name]);
    
    feature=main_get_feature_all(data.ori_slice,data.seg_slice,data.ReconstructionDiameter,data.dls);
    data.feature=feature;
    
    save([path '\' cohort '\' CT_phase '\' total(i).name], 'data');
end
