function?newx=diedai（x）
ZMax=max（max（x））;
%取出最大灰度值
ZMin=min（min（x））;
%取出最小灰度值
TK=（ZMax+ZMin）/2;
bCal=1;
[MN]=size（x）;
%圖像大小
????iForeground=0;
????iBackground=0;
????%定義前景和背景數(shù)
????ForegroundSum=0;
????BackgroundSum=0;
????%定義前景和背景灰度總和
while（bCal）???
????for?i=1:M
?????????for?j=1:N
???????????????tmp=x（ij）;
???????????????if（tmp>=TK）
????????????????????iForeground=iForeground+1;
????????????????????ForegroundSum=Foregro