%圖像OTSU閾值分割的程序設(shè)計
clear;??%清除系統(tǒng)變量
close?all;?%關(guān)閉所有打開的窗口
[name?path]=uigetfile（‘*.jpg‘）;?%出現(xiàn)打開文件界面?
filename=[path?name];???????????%獲取文件的路徑與文件名
I=rgb2gray（imread（filename））;?%??讀圖并轉(zhuǎn)化為黑白圖片
%以下為閾值計算程序
Ni=imhist（I）;%計算直方圖數(shù)組
N=sum（Ni）;???%總像素點(diǎn)個數(shù)
delamax=0;???%類間方差最大值
threshold=0;?%閾值
for?k=2:255
??u=dot（[0:255]Ni/N）;?%圖像的總平均灰度級
??w0=sum（Ni（1:k）/N）;??%C0類像素所占面積的比例
??w1=1-w0;????????????%C1類像素所占面積的比例
??if?w0==0|w0==1??????%當(dāng)w0為1或0時提前結(jié)束本次循環(huán)
?????continue?
??end
??u0=dot（[0:k-1]Ni（1:k）/N）/w0;???%C0類像素的平均灰度
??u1=dot（[k:255]Ni（k+1:256）/N）/w1;?%C1類像素的平均灰度
????dela（k）=w0*（u-u0）^2+w1*（u-u1）^2;??%類間方差公式
??%求出