cvLabMATLAB/SAMPLE.m at main · ComputerVisionLaboratory/cvLabMATLAB · GitHub

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
clear
load('FACE_16x16x1000x10.mat')
whos
X = orzNormalize(orzReshape(DATA,1));
[nDim, nNum, nClass] = size(X);
X1 = X(:,1:100,:);
X2 = X(:,101:1000,:);
X2 = orzSlidingData(X2,25,25);
nData = size(X1,2);
[nDim,nNum2,nSet2,nClass] = size(X2);

nSubDim1 = 15;
nSubDim2 = 4;
nSigma = 1;
Lc = 3;
Label2 = orzLabel(nSet2,nClass);
LabelS = orzLabel(nNum2,nClass);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% [V1 D1] = orzBasisVector(X1,nSubDim1);
% SIM = orzSM(X2,V1);
% RSLT = OrzEval(SIM(:,:,end),orzLabel(size(X2,2)*size(X2,3),nClass));
% SM(1,:) = [RSLT.ER,RSLT.EER];

% SSS = reshape(SIM(:,:,end),nClass,nNum2,nSet2,nClass);
% RSLT = OrzEval(squeeze(mean(SSS,2)),Label2);
% SM(2,:) = [RSLT.ER,RSLT.EER];

% VVV=orzVoting(orzReshape(permute(SSS,[1,3,4,2]),2));
% RSLT = OrzEval(VVV,Label2);
% SM(3,:) = [RSLT.ER,RSLT.EER];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% [V1 D1] = orzBasisVector(X1,nSubDim1);
% [V2 D2] = orzBasisVector(X2,nSubDim2);
% SIM = orzCanonicalAngles(V1,V2);
% l = [1,2,3,4,5,6,7,8,9,10];
% RSLT = OrzEval(SIM(:,:,end,end),orzLabel(size(X2,3),nClass));
% MSM = [RSLT.ER,RSLT.EER];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[A1 D1 C1 K1] = orzKernelBasisVector(X1,nSubDim1,nSigma);
[A2 D2 C2 K2] = orzKernelBasisVector(X2,nSubDim2,nSigma);
SIM = orzKernelCanonicalAngles(X1,A1,X2,A2,nSigma);
RSLT = OrzEval(SIM(:,:,end,end),orzLabel(size(X2,3),nClass));
KMSM = [RSLT.ER,RSLT.EER];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% KC = OrzKCMSM(X1,nSubDim1,nSigma);
% DIFDIM = [KC.nOrthDim:-10:KC.nOrthDim-60];
% for iDifDim = 1:size(DIFDIM,2)
%     nDifDim = DIFDIM(iDifDim);
%     V1 = KC.TransformS(X1,nSubDim1,nDifDim);
%     V2 = KC.TransformS(X2,nSubDim2,nDifDim);
%     SIM = orzCanonicalAngles(V1,V2);
%     RSLT = OrzEval(SIM(:,:,end,end),orzLabel(size(X2,3),nClass));
%     KCMSM(iDifDim,:) = [RSLT.ER,RSLT.EER];
% end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
KO = CvtKOMSMSVD(X1,nSubDim1,nSigma,1,1e-5);
V1 = KO.TransformS(X1,nSubDim1);
V2 = KO.TransformS(X2,nSubDim2);
SIM = orzCanonicalAngles(V1,V2);

RSLT = OrzEval(SIM(:,:,end,end),orzLabel(size(X2,3),nClass));
KOMSM = [RSLT.ER,RSLT.EER];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
nK = 2;
% M = zeros(nDim, nK, nClass);
L = zeros(nClass, nData);
for I1=1:nClass
    [L(I1,:), ~] = kmeansPlusPlus(X1(:, :, I1),nK);
end
LC = histc(L',(1:nK))';
KKO = kmeansKOMSM(X1,Lc,nSubDim1,nSigma,1,1e-5);
V1 = KKO.TransformS(X1,nSubDim1);
V2 = KKO.TransformS(X2,nSubDim2);
SIM = orzCanonicalAngles(V1,V2);

RSLT = OrzEval(SIM(:,:,end,end),orzLabel(size(X2,3),nClass));
KKOMSM = [RSLT.ER,RSLT.EER];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% disp('SM')
% disp(SM);
% disp('MSM')
% disp(MSM);
disp('KMSM')
disp(KMSM);
% disp('KCMSM')
% disp(KCMSM);
disp('KOMSM')
disp(KOMSM);
% disp('KOMSM')
% disp(KKOMSM);