cosmo clusterize hdr

function clusters = cosmo_clusterize(sample, nbrhood_mat)