channel_distances.m

% kymograph channels -- check CONSTANTS numbers may vary
% ERK_CHANNEL = 1;
% AKT_CHANNEL = 2;
% VELOCITY_CHANNEL = 4; 
function res = channel_distances(flist,classList,channel1,channel2,clipLimits,tClip)
res=table();

if ~exist('tClip','var')
    tClip = [];
end

if ~exist('clipLimits','var') || isempty(clipLimits)
    clipLimits = SSF.getClipLimits(flist,[channel1,channel2],0,tClip);
end

parfor ff = 1:length(flist)    
    class = classList{ff};
    
    nr = table();
    
    strDB = fullfile(flist(ff).folder,flist(ff).name);
    im1 = SSF.loadImage(strDB,channel1,0,tClip);  
    im2= SSF.loadImage(strDB,channel2,0,tClip);   
    
    im1 = SSF.quantize8(im1,clipLimits{ff,1});
    im2 = SSF.quantize8(im2,clipLimits{ff,2});
    
    d1 = SSF.ncd_ssf_volume(im1,im2);
    d2 = SSF.ncd_ssf_volume(im2,im1);
    d = min(d1,d2);
    
    nr.idx = ff;
    nr.class = {class};
    nr.d = d;
    res(ff,:) = nr;
end
res.Properties.VariableNames = {'idx','class','d'};
res.class=cellfun(@upper,res.class,'UniformOutput',false);
boxplot(res.d,res.class,'notch','on');
ylabel('NCD( Akt_{SSF}, Velocity_{SSF} )','FontSize',12)
xlabel('experimental condition','FontSize',12)
set(gcf,'color','w')
4;

% 
% classes = unique(res.class);
% mu=[];stderr=[];
% for i = 1:length(classes)
%     idx = find(strcmp(res.class,classes(i)));
%     mu(i) = mean(res.d(idx));
%     stderr(i) = std(res.d(idx))/sqrt(length(idx));
% end