【图像分割】基于最大类间法otsu实现图像分割附matlab代码

原创

matlab科研助手 2022-02-23 18:52:13 博主文章分类：图像处理 ©著作权

©著作权归作者所有：来自51CTO博客作者matlab科研助手的原创作品，请联系作者获取转载授权，否则将追究法律责任

1 简介

在现代战争中,快速而有效的对攻击目标进行识别和跟踪对获取战争主动权起着很重要的作用,而要达到这个目的,就需要从目标图像中准确的分割出目标.在图像分割中,阈值的选取至关重要.最大类间方差法是一种常用而有效的图像分割算法,并已在许多实时场合中采用.为满足高速场合的要求, 本文采用 Altera 公司的 Cyclone II 系列的 FPGA 实现类间方差的计算.实验结果表明,本设计能够实时稳定的对目标分割提取,分割效果良好.

2 完整代码

clear
close all
clc
I=imread('rice.png');    
[m,n]=size(I);   
N=m*n;   
L=256;   
for i=1:L  
    count(i)=length(find(I==(i-1)));  
    f(i)=count(i)/(N);  %每个灰度对应的概率，i=1,对应灰度值为0（i-1）
end  
for i=1:L   
    if count(i)~=0   
        st=i-1;   %开始的灰度值
        break;   
    end   
end   
for i=L:-1:1   
    if count(i)~=0   
        nd=i-1;   %结束的灰度值
        break;   
    end   
end   
p=st;   q=nd-st+1;   
u=0;   
for i=1:q   
    u=u+f(p+i)*(p+i-1);  %u是像素的平均值    
    ua(i)=u;           %ua（i）是前i+p个像素的平均灰度值  (前p个无取值) 
end;   
for i=1:q   
    w(i)=sum(f(1+p:i+p));  %w（i）是前i个像素的累加概率,对应公式中P0 
end;   
w=w+eps;  
   %对照sigmaB的公式写出目标函数。实际是遍历所有值
d=(w./(1-w)).*(ua./w-u).^2;
[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点   
th=tp+p;  
figure;imshow(im2bw(I,th/255),[]);  title('最大类间方差');
%% matlab自带函数
figure;imshow(im2bw(I,graythresh(I)),[]);  title('matlab自带');