x=imread ('C:\Users\mengbin\Desktop\171.jpg')
%figure;
%imshow(x);
[m,n]=size(x);
N=m*n;
num=zeros(1,256);
p=zeros(1,256);
for i=1:m
for j=1:n
num(x(i,j)+1)=num(x(i,j)+1)+1;
end
end
for i=0:255;
p(i+1)=num(i+1)/N;屠波
end
totalmean=0;
for i=0:255;
totalmean=totalmean+i*p(i+1);
poco2007
end羌族锅庄
maxvar=0;
for k=0:255
kk=k+1;
zerosth=sum(p(1:kk));
firsth=0;
for h=0:k
firsth=firsth+h*p(h+1);
end
var=totalmean*zerosth-firsth;
var=var*var;
var=var/(zerosth*zerosth+ 0.1);
var=sqrt(var);
maxvar=var;
point=k;
end
end
threshold=point;
for i=1:m
for j=1:n
if x(i,j)>二等残废point
bin(i,j)=1;
else
bin(i,j)=0;
end
end
end
李芸昕imshow(bin);