clear;
warning?off;
I=imread（‘D:\mytest\true1.jpg‘）;
I_gray=rgb2gray（I）;%rgb轉(zhuǎn)灰度
subplot（121）;
imshow（I_gray）title（‘灰度圖像‘）;
Y3=medfilt2（I_gray[3?3]）;
I_double=double（Y3）;%轉(zhuǎn)化為雙精度
[widlen]=size（I_gray）;%wid為行數(shù)，len為列數(shù)
colorlevel=256;??%灰度級(jí)
hist=zeros（colorlevel1）;%直方圖?256*1的0矩陣
%threshold=128;?%初始閾值

%計(jì)算直方圖
for?i=1:wid
????for?j=1:len
????????m=I_gray（ij）+1;
????????hist（m）=hist（m）+1;
????end
end
hist=hist/（wid*len）;%直方圖歸一化
%miuT為總的平均灰度，hist[m]代表像素值為m的點(diǎn)個(gè)數(shù)
miuT=0;
for?m=1:colorlevel
????miuT=miuT+（m-1）*hist（m）;
end
xigmaB2=0;%用于保存每次計(jì)算的方差，與下次計(jì)算的方差比較大小

for?mindex=1:colorlevel
????threshold=mindex-1;
????omega1=0;?%前景點(diǎn)所占比例
????omega2=0;?%背景點(diǎn)所占比例
????for?m=1:threshold-1
??????????omega1=omega1+hist（m）;%計(jì)算前景比例
????end
????omega2=1-omega1;%計(jì)算背景比