算法

1  OTSU阈值分割方法

num_bins = 256;
  counts = imhist(I,num_bins);%列向量
  
  % Variables names are chosen to be similar to the formulas in
  % the Otsu paper.
  p = counts / sum(counts);
  omega = cumsum(p);
  mu = cumsum(p .* (1:num_bins)');
  mu_t = mu(end);
  
  sigma_b_squared = (mu_t * omega - mu).^2 ./ (omega .* (1 - omega));

上面的代码是以灰度图像（I）来说明的，omega对应于w(k),而mu对应于u(k)，mu_t对应于u。sigma_b_squared对应于方差（在这里是一个列向量），然后求 sigma_b_squared的最大值即可