THESE/codes/seuilhisto.m

   1 function [ seuil ] = seuilhisto( Iin, Tini, dT)
   2 %UNTITLED3 Summary of this function goes here
   3 %   Detailed explanation goes here
   4 [nb, gl] = imhist(Iin) ;
   5 Ti = Tini ;
   6
   7 while abs((moyhisto(nb,gl,0,Ti-1)+moyhisto(nb,gl,Ti,255))/2-Ti) > dT,
   8     Ti = round((moyhisto(nb,gl,0,Ti-1)+moyhisto(nb,gl,Ti,255))/2)
   9
  10 end
  11 seuil = Ti ;