-
Notifications
You must be signed in to change notification settings - Fork 0
/
runMetricsComparisonCLBP_v2.m
106 lines (94 loc) · 2.8 KB
/
runMetricsComparisonCLBP_v2.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
% /* --------------------------------------------------------------------------------------
% * File: runMetricsComparison.m
% * Date: 01/06/2015
% * Author: David Pastor Escuredo, [email protected]
% * Version: 0.2
% * License: BSD
% * --------------------------------------------------------------------------------------
% Copyright (c) 2013-2017, David Pastor Escuredo
% with Biomedical Image Technology, UPM (BIT-UPM)
% with BioEmergences, CNRS
% with LifeD lab
% All rights reserved.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%% SPATIO-TEMPORAL CLUSTERS %%%%%%%%%%%%%%%%%%%%%%%
addpath('colormaps')
cs=1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%% Calculate distance to LBPs of the reference %%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%D=zeros(size(clus,2), size(clus,2));
%Df=zeros(size(clus,2), size(clus,2));
D=zeros(size(clus,2)-5, 5);
Df=zeros(size(clus,2)-5, 5);
%WE COMPARE AT LEAST 3 TIME STEPS
%We skip t=1 as the lagrangian is all 1
for t=1:size(clus,2)
% for t2=1:size(lbpRef,2)
%min of distance to ref lbps
lbpt=lbpRef(:,t);
clust=clus(:,t);
for j=1:size(clus,1)
clustj=clust(j,:);
for b=1:size(lbpt,1)
dclustj(b)=pdist2(clustj,lbpt(b,:),'euclidean');
end
dclust(j)=min(dclustj);
end
size(clust)
dref=prctile(dclust,25);
dref2=prctile(dclust,75);
drefmin=min(dclust);
drefmax=max(dclust);
drefm=mean(dclust);
D(t,:)=[drefm dref dref2 drefmin drefmax];
end
dlmwrite([desc 'CLBP_score' '-' tagDesc '_v2.csv'],D)
% maxsD=[0 0 0 0 0.0002 0.006 0.006 0.007 0]
% max(D(:))
% max(Df(:))
%
% h=figure
% colormap('bone')
% %Df=Df./maxsD(DescriptorIndex);
%
% if normType==1
%
% D=D./maxsD(DescriptorIndex);
% D=1-D;
% tagn='1'
% else
% 'hols'
% end
%
% imagesc(D, [0 1])
% axis off
% screen_size = get(0, 'ScreenSize');
% set(h, 'Position', [0 0 screen_size(3) screen_size(4) ] );
% if savePlots
% saveas(h,[StatsPath dataset tagfolder filesep 'past_score' tagn '-' datref '-' tagDesc '-' selCtag tag_modes '-modes' num2str(mxCluster) '.png'],'png')
% end
%
% h2=figure
% colormap('pink')
% %Df=max(Df(:))-Df;
%
% if normType==1
% Df=Df./maxsD(DescriptorIndex);
% Df=1-Df;
% tagn=''
% else
% 'hola'
% end
% imagesc(Df, [0 1])
% axis off
% screen_size = get(0, 'ScreenSize');
% set(h2, 'Position', [0 0 screen_size(3) screen_size(4) ] );
% if savePlots
% saveas(h2,[StatsPath dataset tagfolder filesep 'future_score' tagn '-' datref '-' tagDesc '-' selCtag tag_modes '-modes' num2str(mxCluster) '.png'],'png')
% end
%
% max(D(:))
% max(Df(:))
%
% pause