function [contour_map,count_levels] = hist2_contour( counts, numContours )


tmp = linspace(0,100,numContours+2); 
prctile_ps = tmp(2:end-1);
count_levels = prctile(counts(:), prctile_ps );
first_good = min(find(count_levels>0));
tmp = linspace(0,100,numContours+2); 
prctile_ps = tmp(2:end-1);
count_levels = prctile(counts(counts > count_levels(first_good)), prctile_ps );

contour_map = zeros(size(counts));
for c = 1 : numContours
    level = count_levels(c);
    ZC = zero_crossings(counts-level);
    contour_map(find(ZC)) = c;
end