Remote.Sensing/testNmfASCL1.m at master · yuehniu/Remote.Sensing · 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
%% hyperNmfASCL1_2 test with more than 3 endmembers
clear all
expTimes = 1;
dataSize = 900;
noiseLevel = 0.0;
bandNum = 3;
endNum = 3;
HTrue = zeros(endNum, bandNum);
HI = zeros(endNum, bandNum);
HASCL1 = zeros(endNum, bandNum);
WTrue = zeros(dataSize, endNum);
WI = zeros(dataSize, endNum);
WASCL1 = zeros(dataSize, endNum);
% for i = 1:exp_times
    % generate true W, H, V
    HTrue = abs( randn( endNum, bandNum ) );
    [V, W_true] = create4(dataSize, HTrue);
%     max_V = max(V);
%     for i = 1:bandNum
%        V(:,i) = V(:,i) / (max_V(i));
%     end
    figure;
    scatter(V(:,1), V(:,2));
    xlabel('band 1');
    ylabel('band 2');
%     xlim([0,1])
%     ylim([0,1])

    %% find initial H using n_findr
    HInitIndx = nFindr(V, endNum);
    HI = V(HInitIndx, :);

    alpha = 1;
    tol = 0.1;
    maxIter = 5000;
    [WI, EI] = nmfAbundance(V, endNum, HI,...
                        alpha, tol, maxIter);
    VASCL1 = WI * HI;

    figure;
    scatter(VASCL1(:,1), VASCL1(:,2), 'c', 'full'); hold on
    scatter(V(:,1), V(:,2), 'r')
    scatter(HI(:, 1), HI(:, 2), 'full');

    % factorize V using nmf_MDC_simple method

%% ASCL1_2 test
    random_ = 1;
    if random_
        HI = abs(randn(size(HI)));
        WI = abs(randn(size(WI)));
        WI = WI ./ ( repmat( sum(WI,2),1, endNum ) );
        [ WASCL1, HASCL1, SRcL1, errRcL1, objRcL1] = ...
            hyperNmfASCL1(...
                V', HI', WI',...
                0.001,... % tolObj
                20000,... % maxIter
                10 ... %fDelta
                );
        WMdcRandomInit = WASCL1;
        HMdcRandomInit = HASCL1;
    else
        [ WASCL1, HASCL1, SRcL1, errRcL1, objRcL1] = ...
            hyperNmfASCL1(...
                V', HI', WI',...
                0.001,... % tolObj
                20000,... % maxIter
                5 ... %fDelta
                );
        WMdcWellInit = WASCL1;
        HMdcWellInit = HASCL1;
    end
    %% visualize rest
    VASCL1 = HASCL1 * WASCL1;
    figure;
    bandIndx1 = 1;
    bandIndx2 = 2;
    scatter(V(:,bandIndx1), V(:,bandIndx2), 'c' ); hold on ;
    scatter(HTrue(:, bandIndx1), HTrue(:, bandIndx2), 'filled', 'r');
    scatter(HI(:, bandIndx1), HI(:, bandIndx2), 'filled', 'b');
    scatter( VASCL1(bandIndx1,:), VASCL1(bandIndx2,:), 5, 'k' );
    scatter( HASCL1(bandIndx1, :), HASCL1(bandIndx2,:) , 'filled','k')
    plot(SRcL1(1, :,bandIndx1), SRcL1(1, :,bandIndx2), 'r-.', 'MarkerSize', 5);
    plot(SRcL1(2, :,bandIndx1), SRcL1(2, :,bandIndx2), 'g-.', 'MarkerSize', 5);
    plot(SRcL1(3, :,bandIndx1), SRcL1(3, :,bandIndx2), 'b-.',  'MarkerSize', 5);

    figure;
    hold on;
    plot(errRcL1, 'r');
    plot(objRcL1, 'k')
    plot(objRcL1-errRcL1, 'm')