基于OSTU的阈值分割算法研究

function H_total=ostu（threshold）
% fxy=fstatus（xy）;
% 适应值函数 xy为多维数组
% 对约束化问题，这里用惩罚函数的外点法，
% 惩罚函数法的基本思想是借助惩罚函数把约束化问题转化为列约束化问题．
% 外点法　对违反约束的点在目标函数中加入相应的惩罚，而对可行点不惩罚．

% 求:max f（xy）

global G_min G_max
fvarsize=size（threshold）;
H_total=zeros（fvarsize（1）1）;
global p

w1=0;%定义参数
w2=0;
u1=0;
u2=0;
u=0;
m=1;

%f1=csvread（‘d:\f1.txt‘）;
%p=f1;
%for a=1:254
%s=f1（a）+s;    %计算总的像素数
%end
%for a=1:254
 %   p（a）=f1（a）/s;
 %end           %计算各灰度值的概率
thr