[ci,cp,out]=thresh(I,rmin,rmax)