PATHS_GROUPS = {};
tmp = [];
len = length(PATHS);
for i=1:len
element = PATHS(i);
if element ~= ';'
tmp = [tmp element];
elseif tmp(end) == 'n'
tmp = [ tmp '\'];
PATHS_GROUPS = [ PATHS_GROUPS; tmp ];
tmp = [];
else
tmp = [];
end
end
num_targets = length(PATHS_GROUPS);
UV_mus_musp_g = cell(num_targets, 4);
for i=1:num_targets
i
target_path=cell2mat( PATHS_GROUPS(i) );
root=[];
len_target_path = length(target_path);
for j=length('H:\Breast_')+1 : len_target_path-length('\target_region\')
ele = target_path(j);
if ele == '\'
ele = '_';
elseif ele == '-'
ele = '_';
end
root=[root ele];
end
load( [target_path, 'mus_map2.mat'] )
load( [target_path, 'musp_map2.mat'] )
load( [target_path, 'g_map2.mat'] )
UV = imread( [target_path, 'uv.ppm'] );
%% UV归一化
a=double(UV(100:2400,100:2400,3));
originalMinValue = double(min(min(a)));
originalMaxValue = double(max(max(a)));
originalRange = originalMaxValue - originalMinValue;
desiredMin = 0;
desiredMax = 1;
desiredRange = desiredMax - desiredMin;
UV_01 = desiredRange * (double(a) - originalMinValue) / originalRange + desiredMin;
%% 相乘
UV_musmap2 = mus_map2 .* UV_01;
UV_muspmap2 = musp_map2 .* UV_01;
UV_gmap2 = g_map2 .* UV_01;
%% 选择是否保存二维图
save2D = 0
if save2D == 1
mus2D_name=[root, '_', 'UVmus'];
musp2D_name=[root, '_', 'UVmusp'];
g2D_name=[root, '_', 'UVg'];
% eval([mus2D_name,'=UV_musmap2;'])
% eval([musp2D_name,'=UV_muspmap2;'])
% eval([g2D_name,'=UV_gmap2;'])
save(['H:\UV_Params2D\', mus2D_name], 'UV_musmap2')
save(['H:\UV_Params2D\', musp2D_name], 'UV_muspmap2')
save(['H:\UV_Params2D\', g2D_name], 'UV_gmap2')
end
%% 平均值
meanUV_musmap2 = mean(mean(UV_musmap2));
meanUV_muspmap2 = mean(mean(UV_muspmap2));
meanUV_gmap2 = mean(mean(UV_gmap2));
UV_mus_musp_g(i,1) = cellstr(root);
UV_mus_musp_g(i,2) = num2cell(meanUV_musmap2);
UV_mus_musp_g(i,3) = num2cell(meanUV_muspmap2);
UV_mus_musp_g(i,4) = num2cell(meanUV_gmap2);
end
xlswrite('H:\UV_multi_Params\UV_mus_musp_g.xlsx', UV_mus_musp_g);