function?T=ostu1（I）
%I?=?im2uint8（I（:））;
counts=imhist（I）;
rank=256;
p=counts/sum（counts）;
omega=cumsum（p）;
mu?=?cumsum（p?.*?（1:rank）‘）;
mu_t?=?mu（end）;
%sigma_2=zeros（1256）;
%for?i=1:256
????%if（1-omega（i）==0）
????%????continu