cvLabMATLAB/OrzKCMSM.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
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
classdef OrzKCMSM
    properties (SetAccess = public)
        nDim;
        nNum1;
        nClass;
        nSubDim1;
        nOrthDim;
        nSigma;

        nAlpha;
        nBeta;

        X1;
        A1;
        E1;
        C1;

        D;
        O;
        W;
    end% properties

    methods
        function OB = OrzKCMSM(X1, nSubDim1,nSigma, varargin)
            %function OB = OrzKOMSM(X1, nSubDim1,nSigma, varargin)
            % nDim:     次元
            % nNum1:    各クラスのデータ数
            % nClass:   クラス数
            % nSubDim1: 辞書部分空間の時限
            % nOrthDim: 空間の次元
            % nSigma:   ガウシアンカーネルパラメータ
            % nAlpha:   直交化行列の寄与率(第４パラメータ、デフォルトで1.0)
            %
            % X1:       辞書データ
            % A1:       X1のKPCA結合係数
            % E1:       固有値
            % C1:       寄与率
            %
            % D:        直交化カーネルグラム行列
            % O:        カーネル直交化行列
            % W:        Dの固有値
            %
            % [V2 E2 C2 A2 ] = TransformS(OB, X2,nSubDim2)
            % 部分空間に対する直交化変換
            % KPCA ⇒ 直交化変換 ⇒ グラムシュミット
            %
            % [V2 E2 C2 Y2 ] = TransformV(OB, X2, nSubDim2)
            % パターンベクトルに対する直交化変換し基底ベクトルを求める
            % 非線形直交化変換 ⇒ PCA
            %
            %function [Y2] = Transform(OB, X2)
            % パターンベクトルに対する直交化変換だけ
            % 非線形直交化変換


            if nargin == 3
                OB.nAlpha = 1;
            else
                OB.nAlpha = varargin{1};
            end

            OB.nSubDim1 = nSubDim1;
            OB.nSigma = nSigma;
            OB.X1 = X1;


            [OB.nDim , OB.nNum1, OB.nClass] = size(X1);
            OB.nOrthDim = OB.nClass*OB.nSubDim1;
            OB.A1 = zeros(OB.nNum1,OB.nSubDim1,OB.nClass);
            OB.E1 = zeros(OB.nSubDim1,OB.nClass);
            OB.C1 = zeros(1,OB.nClass);
            for I=1:OB.nClass
                [OB.A1(:,:,I) OB.E1(:,I) OB.C1(I)] = orzKPCA(OB.X1(:,:,I),OB.nSubDim1,OB.nSigma,'R');
                I
            end

            OB.D = zeros(OB.nSubDim1, OB.nSubDim1, OB.nClass,OB.nClass);
            for I1 = 1:OB.nClass
                for I2 = I1:OB.nClass
                    K = exp(-orzL2Distance(X1(:,:,I1),X1(:,:,I2))/nSigma);
                    if I1 == I2
                        OB.D(:,:,I1,I2) = eye(OB.nSubDim1,OB.nSubDim1);
                    else
                        OB.D(:,:,I1,I2) =  OB.A1(:,:,I1)'*K* OB.A1(:,:,I2);
                        OB.D(:,:,I2,I1) = OB.D(:,:,I1,I2)';
                    end
                end
            end
            OB.D = reshape(permute(OB.D,[1,3,2,4]), OB.nSubDim1 * OB.nClass, OB.nSubDim1 * OB.nClass);


            [B,OB.W] = eig(OB.D);
            [OB.W ind] = sort(diag(OB.W),'ascend');
            B = B(:,ind);
            OB.O = diag(1./sqrt(OB.W))*B';


        end
        function [V2 E2 C2 A2 ] = TransformS(OB, X2,nSubDim2,nSubDim3)
            % 部分空間に対する直交化変換
            % KPCA ⇒ 直交化変換 ⇒ グラムシュミット
            nSize = size(X2);
            X2=X2(:,:,:);
            [tmp, nNum2,nSet2] = size(X2);

            A2 = zeros(nNum2,nSubDim2,nSet2);
            E2 = zeros(nSubDim2,nSet2);
            C2 = zeros(1,nSet2);
            for I=1:nSet2
                [A2(:,:,I) E2(:,I) C2(I)] = orzKPCA(X2(:,:,I),nSubDim2,OB.nSigma,'R');
            end

            V2 = zeros(nSubDim3,nSubDim2,nSet2);
            for J = 1:nSet2
                a = zeros( OB.nSubDim1, nNum2,OB.nClass);
                for I = 1:OB.nClass
                    Z = exp(-orzL2Distance(OB.X1(:,:,I),X2(:,:,J))/OB.nSigma);
                    a(:,:,I) = OB.A1(:,:,I)' * Z;
                end
                a = permute(a,[1,3,2]);
                a = reshape(a,size(a,1)*size(a,2),size(a,3));
                V2(:,:,J) = cvtGramSchmidt(OB.O(1:nSubDim3,:)*a*A2(:,:,J));
            end
            V2 = reshape(V2,[nSubDim3,nSubDim2,nSize(3:end),1]);
        end


        function [Y2] = Transform(OB, X2)
            % パターンベクトルに対する直交化変換だけ
            % 非線形直交化変換
            nSize = size(X2);
            X2=X2(:,:,:);
            [tmp nNum2,nSet2] = size(X2);

            Y2 = zeros(OB.nOrthDim,nNum2,nSet2);
            for J=1:nSet2
                a = zeros( OB.nSubDim1,nNum2,OB.nClass);
                for I = 1:OB.nClass
                    Z = exp(-orzL2Distance(OB.X1(:,:,I),X2(:,:,J))/OB.nSigma);
                    a(:,:,I) = OB.A1(:,:,I)' * Z;
                end
                a = permute(a,[1,3,2]);
                a = reshape(a,OB.nSubDim1*OB.nClass,nNum2);
                Y2(:,:,J) = OB.O*a;
            end
            Y2 = reshape(Y2,[OB.nOrthDim,nNum2,nSize(3:end),1]);

        end

    end% methods
end% classdef