• 大小: 5KB
    文件类型: .m
    金币: 1
    下载: 0 次
    发布日期: 2021-06-16
  • 语言: Matlab
  • 标签: 多阈值  

资源简介

代码用matlab编写的,三个阈值分割的示例,本算法用的是最大内间方差算法改造而成的

资源截图

代码片段和文件信息

%-----------------------------基于最大内间方差思想的多阈值图像分割算法-------------------------
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

评论

共有 条评论