getcubes.m

Information Fusion of Flash and Non-Flash Images (Project)

function tcubes = getCubes(tdata,tol);
%learns multiple categories
[N D] = size(tdata);
fdata = [tdata 1-tdata];
%start out w/ one box
tcubes = fdata(1,:);
for i = 2:N
    
    %find # of nodes
    numNodes = size(tcubes,1);
    pt = ones(numNodes,1)*fdata(i,:);
    %find node w/ min distance
    tmp1 = abs(min(tcubes, pt)-tcubes);
    conmap = sum(tmp1,2);
    [minval, ind] = min(conmap);
    
    if (minval == 0)    %CASE 1: falls inside box: do nothing!
        %disp('in box');
    elseif (minval <= tol) %CASE 2: falls close enough to box: stretch box
        %disp('stretch box');
        tcubes(ind,:) = min(tcubes(ind,:), fdata(i,:));
        
    else %CASE 3: make new box!
        %disp('new box');
        tcubes = [tcubes; fdata(i,:)];
    end
end