基于直方图双峰特性的图像分割Matlab代码

假设某图像的灰度直方图具有二峰性（ f(Ta)=Ha; f(Tb)=Hb ），表明这个图像较亮的区域和较暗的区域可以较好地分离。取二峰间的谷点为阈值Th，可以得到好的二值处理的效果。

clc 
clear 
I=imread（‘ct.bmp‘）;
[m n]=size（I）;
a=141;b=173;
imhist（I）;%显示图像直方图
I=double（I）;
h=zeros（1256）;%定义直方图向量h
%计算每个灰度级的像素个数，将其存入h向量
for i=1:m 
    for j=1:n
        h（I（ij）+1）=h（I（ij）+1）+1;
    end
end
%寻找直方图双峰之间谷底的灰度级Th
%hmin=m*n;for k=a:b  if h（k+1）h1=h（a:b）; 
Th=a+fin