I=imread（‘cell.bmp‘）;
subplot（131）imshow（I）;
title（‘原始圖像‘）
%下面使用MATLAB函數(shù)計(jì)算閾值
level=graythresh（I）;
%大津法計(jì)算全局圖像I的閾值
BW=im2bw（Ilevel）;
%閾值分割
subplot（132）imshow（BW）
title（‘graythresh計(jì)算閾值‘）
disp（strcat（‘graythresh計(jì)算灰度閾值:‘num2str（uint8（level*255））））

%Matlab程序?qū)崿F(xiàn)簡(jiǎn)化計(jì)算閾值
iMax=max（max（I））;
iMin=min（min（I））;
%計(jì)算最大和最小值

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）
??????%從最小灰度值到最大值分別計(jì)算方差
??????TK=T（1i）;
??????iForeground=0;
??????iBackground=0;
??????%定義前景和背景數(shù)
??????ForegroundSum=0;
??????BackgroundSum=0;
??????%定義前景和背景灰度總和
??????for?j=1:iSize（1）
??????????for?k=1:iSize（2）
??????????????tmpData=I（jk）;
??????????????if（tmpData>=TK）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

?????????????????1842????????????????????1