Otsus 图像分割

图像分割，双峰法.显示直方图，根据直方图确定阈值，并根据次阈值进行图像的风格

I=imread（‘cell.bmp‘）;
subplot（131）imshow（I）;
title（‘原始图像‘）
%下面使用MATLAB函数计算阈值
level=graythresh（I）;
%大津法计算全局图像I的阈值
BW=im2bw（Ilevel）;
%阈值分割
subplot（132）imshow（BW）
title（‘graythresh计算阈值‘）
disp（strcat（‘graythresh计算灰度阈值:‘num2str（uint8（level*255））））

%Matlab程序实现简化计算阈值
iMax=max（max（I））;
iMin=min（min（I））;
%计算最大和最小值

T = zeros（26）;
for j =1:length（iMin）
    T（1j） = double（iMin（::j））;
    T（1j+3） = double（iMax（::j））;
end
T 
iSize=size（I）;
muxSize=iSize（1）*iSize（2）;
for i=1:length（T）
      %从最小灰度值到最大值分别计算方差
      TK=T（1i）;
      iForeground=0;
      iBackground=0;
      %定义前景和背景数
      ForegroundSum=0;
      BackgroundSum=0;
      %定义前景和背景灰度总和
      for j=1:iSize（1）
          for k=1:iSize（2）
              tmpData=I（jk）;
              if（tmpData>=TK）

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

     文件       1842  2010-05-27 19:28  Otsus.m

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

                 1842                    1