f=imread（‘45.tif‘）;
subplot（221）;imshow（f）;title（‘原圖像‘）;%顯示原圖像

%全局閾值分割
count=0;
T=mean2（f）;
done=false;
while?~done
???count=count+1;
???g=f>T;
???Tnext=0.5*（mean（f（g））+mean（f（~g）））;
???done=abs（T-Tnext）<0.5;
???T=Tnext;
end
disp（[‘count=‘num2str（count）]）?%打印輸出count的值
disp（[‘T=‘num2str（T）]）??%打印輸出T的值
g=im2bw（fT/255）;
subplot（222）;imhist（f）;title（‘f直方圖‘）;
subplot（223）;imshow（g）;title（‘閾值后圖像‘）;

%otsu閾值分割
[T?SM]?=?graythresh（f）;
g1?=?im2bw（f?T）;??
subplot（224）;?imshow（g1）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????455566??2005-09-17?21:06??數字圖像處理之全局閾值分割與Otsu閾值分割MATLAB實現\45.tif
?????文件????????1368??2019-05-25?19:03??數字圖像處理之全局閾值分割與Otsu閾值分割MATLAB實現\T.m
?????文件?????????549??2019-05-25?19:11??數字圖像處理之全局閾值分割與Otsu閾值分割MATLAB實現\otsu.m
?????目錄???????????0??2019-05-26?18:48??數字圖像處理之全局閾值分割與Otsu閾值分割MATLAB實現\