多阈值算法分割算法

%-----------------------------基于最大内间方差思想的多阈值图像分割算法-------------------------
tic
clc
clear all;
close all;
img = imread（‘lenna.bmp‘）;
I_gray=rgb2gray（img）;%将三原色彩色图像转化为灰度图像
I_double=double（I_gray）;    %转化为双精度，因为大多数的函数和操作都是基于double的
[widlen]=size（I_gray）;     %wid为行数，len为列数b
colorlevel=256;     %灰度级
hist=zeros（colorlevel1）;   %直方图，256×1的0矩阵
threshold_num=3;%代表有三个阈值
threshould_vector=zeros（threshold_num1）;%初始化一个向量数组用于存储阈值
W=zeros（threshold_num+11）;%用于存储每个部分的w的值
miu=zeros（threshold_num+11）;%用于存储阈值分割出来的图像各处的平均值
best_threshold=zeros（1threshold_num）;%用于保存最好的阈值
% file_ceache = fopen（‘test.txt‘‘w‘）;%用于存储命令行输出的数据
%计算直方图，统计灰度值的个数
for i=1:wid
    for j=1:len
        m=I_gray（ij）+1;    %因为灰度为0-255所以+1
        hist（m）=hist（m）+1;
    end
end
%直方图归一化
hist=hist/（wid*len）;
%miuT为总的平均灰度，hist[m]代表像素值为m的点个数
miuT=0;
for m=1:colorlevel
    miuT=miuT+（m-1）*hist（m）;
end
max_fangcha=0;  %初始方差大小
%---------------------程序的开始-----------------------------------
for t1=1:colorlevel-2
    for t2=t1+1:colorlevel-1
        for t3=t2+1:colorlevel
            miu=zeros（threshold_num+11）;
            W=zeros（threshold_num+11）;
            %计算第一个目标所占的比例
            for m=1:t1
                W（1）=W（1）+hist（m）;
            end
            %计算第二个目标所占的比例
            for m=t1+1:t2
                W（2）=W（2）+hist（m）;
            end
            %计算第三个目标所占的比例
            for m=t2+1:t3
                W（3）=W（3）+hist（m）;
            end
            %计算第四个目标所占的比例
            W（4）=1-W（1）-W（2）-W（3）;
            %             fprintf（file_ceache‘%d---%d---%d目标W的值：%f%f%f%f\n‘t1t2t3W（1）W（2）W（3）W（4））;
            %计算目标的平均灰度
            for m=1:colorlevel
                if（m<=t1）
                    miu（1）=miu（1）+（m-1）*hist（m）;
                else if（m>t1&&m<=t2）
                        miu（2）=miu（2）+（m-1）*hist（m）;
                    else if（m>t2&&m<=t3）
                            miu（3）=miu（3）+（m-1）*hist（m）;
                        else
                            miu（4）=miu（4）+（m-1）*hist（m）;
                        end
                    end
                end
            end
            %    f