%call thjis after go to proces the result

r=res;  %copy to play with

thresh=800;
r(find(r<thresh))=0;
r(find(r>thresh))=1;

r=wither(r);

out(:,:,2) = cf_rescale(gpic, min(gpic(:)), max(gpic(:)), 0, 1);
out(:,:,1) = cf_rescale(r, min(r(:)), max(r(:)), 0, 1);

imagesc(out)