function output=topover(input,sparsity);


input(logical(eye(size(input))))=nanmax(input(:))+1;
output=input;

path(path,'/home/jesse')
  for m=1:size(input,1)
    [B,i,j] = unique(input(m,:));
freq = histc(j,1:length(i));
B=B(end:-1:1);
freq=freq(end:-1:1);
spar=cumsum(freq)/sum(freq);
thresh=B(max(find((spar-sparsity)<=0)));
output(m,:)=real(input(m,:)>=thresh);
%if m/100==ceil(m/100)
%m
%end

end

  output=real(sparse(output & output'));