Otsu算法实现（matlab）

function threshold = Otsu(normalizedHist)
%normalizedHist为归一化直方图
%threshold为自动分割的阈值
wk = 0;
uk = 0;
L = length(normalizedHist);
uT = sum((0:L-1).*normalizedHist');
maximum = 0;

for k=1:L
    wk = wk + normalizedHist(k); %到阈值k处的累加直方图分量
    uk = uk + (k-1) * normalizedHist(k); %到阈值k处的灰度均值
    eta = (wk*uT - uk)^2 / (wk*(1-wk)); %类间方差/总体方差，目标是使它最大
    if ( eta >= maximum )
        threshold = k;
        maximum = eta;
    end
end

end

测试代码

img = imread('图片路径');
[counts x] = imhist(img);
[m n] = size(img);
normalizedHist = counts./(m*n);
level = Otsu(normalizedHist);
output = img;
output(output<level) = 1;
output(output>=level) = 0;
imagesc(output);