用matlab实现自适应图像阈值分割最大类间方差法

用matlab实现自适应图像阈值分割（最大类间方差法） %本程序是利用最大类间方差算法求解自适应阈值，对图像进行分割

%本程序是利用最大类间方差算法求解自适应阈值，对图像进行分割
clear;
warning off;

web -browser http://www.ilovematlab.cn/thread-10282-1-1.html
SE = strel（‘diamond‘4）;
BW1 = imread（‘img1.bmp‘）;
BW2 = imerode（BW1SE）;
BW3 = imdilate（BW2SE）;
BW4 = BW1-BW3;

%rgb转灰度
if isrgb（BW4）==1
    I_gray=rgb2gray（BW4）;
else
    I_gray=BW4;
end
figureimshow（I_gray）;
I_double=double（I_gray）;%转化为双精度
[widlen]=size（I_gray）;
colorlevel=256;  %灰度级
hist=zeros（colorlevel1）;%直方图
%threshold=128; %初始阈值

%计算直方图
for i=1:wid
    for j=1:len
        m=I_gray（ij）+1;
        hist（m）=hist（m）+1;
    end
end
hist=hist/（wid*len）;%直方图归一化
miuT=0;
for m=1:colorlevel
    miuT=miuT+（m-1）*hist（m）;
end
xigmaB2=0;

for mindex=1:colorlevel
    threshold=mindex-1;
    omega1=0;
    omega2=0;
    for m=1:threshold-1
          omega1=omega1+hist（m）;
    end
    omega2=1-omega1;
    miu1=0;
    miu2=0;

    for m=1:colorlevel
        if m           miu1=miu1+（m-1）*hist（m）;
        else
           miu2=miu2+（m-1）*hist（m）;
        end
    end
    miu1=miu1/omega1;
    miu2=miu2/omega2;
    xigmaB21=omega1*（miu1-miuT）^2+omega2*（miu2-miuT）^2;
    xigma（mindex）=xigmaB21;
    if xigmaB21>xigmaB2
        finalT=threshold;
        xigmaB2=xigmaB21;
    end
end
fT=finalT/255 %阈值归一化
T=graythresh（I_gray）%matlab函数求阈值
for i=1:wid
    for j=1:len
        if I_double（ij）>finalT
            bin（ij）=1;
        else
            bin（ij）=0;
        end
    end
end
figureimshow（bin）;
figureplot（1:colorlevelxigma）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1609  2008-05-23 18:02  tuxiangfenge.m

     文件        122  2008-05-14 18:20  Matlab中文论坛--助努力的人完成毕业设计.url

     文件       3434  2008-05-15 08:26  使用帮助：新手必看.htm

----------- ---------  ---------- -----  ----

                 5165                    3