matlab練習程式（最大類間方差二值化映像）

　　這個過去也實現過，不過現在回頭看，當時寫的還真是糟糕，所以現在決定用matlab重寫了。而且當時的方法現在來看還真是不怎麼樣呢。

　　這裡用的方法是《特徵提取與影像處理》3.3.4裡的。

　　下面是代碼：

function main    img=imread('lena.jpg');    imshow(img);    img=double(img);    [m n]=size(img);        Hist=zeros(1,256);    for i=1:m        for j=1:n            Hist(img(i,j)+1)=Hist(img(i,j)+1)+1;    %求長條圖        end    end    p=Hist/(m*n);       %長條圖機率分布    uT=sum((1:256).*p(1:256));  %映像亮度均值，其實比真正的均值要大1，所以後面減了1    sigma_2=zeros(1,256);    for k=1:256        sigma_2(k)=(uT*w(k,p)-u(k,p))^2/(w(k,p)*(1-w(k,p)));    %類間方差     end    [tmp index]=max(sigma_2);       %求最大類間方差的索引    index=index-1;                  %這裡索引是1-256，實際映像灰階是0-255，所以減1        imgn=img>index;    figure;    imshow(imgn);        function re=w(k,p)      %長條圖前k個亮度級的0階累積矩        re=sum(p(1:k));       end    function re=u(k,p)      %長條圖前k個亮度級的1階累積矩        re=sum((1:k).*p(1:k));    end        end

效果：

原圖

二值化後