Mercurial > hg > smallbox
view DL/two-step DL/dico_color.m @ 236:5f4e47b78f2b ver_2.0_alpha1
Added example.
author | luisf <luis.figueira@eecs.qmul.ac.uk> |
---|---|
date | Thu, 19 Apr 2012 17:59:08 +0100 |
parents | 485747bf39e0 |
children |
line wrap: on
line source
function colors = dico_color(dico, mu) % DICO_COLOR cluster a dictionary in pairs of high correlation atoms. % Called by dico_decorr. % % Parameters: % -dico: the dictionary % -mu: the correlation threshold % % Result: % -colors: a vector of indices. Two atoms with the same color have a % correlation greater than mu numAtoms = length(dico); colors = zeros(numAtoms, 1); % compute the correlations G = abs(dico'*dico); G = G-eye(size(G)); % iterate on the correlations higher than mu c = 1; maxCorr = max(max(G)); while maxCorr > mu % find the highest correlated pair x = find(max(G)==maxCorr, 1); y = find(G(x,:)==maxCorr, 1); % color them colors(x) = c; colors(y) = c; c = c+1; % make sure these atoms never get selected again G(x,:) = 0; G(:,x) = 0; G(y,:) = 0; G(:,y) = 0; % find the next correlation maxCorr = max(max(G)); end % complete the coloring with singletons index = find(colors==0); colors(index) = c:c+length(index)-1; end