大津阈值分割法.docx

clc
clear
y=imread('mat4.ico');
a=rgb2gray(y);
L=256;
Ps=zeros(L,1);
nk=zeros(L,1);
[row,col]=size(y);
n=row*col;
fori=1:row
forj=1:col
num=double(y(i,j))+1;
nk(num)=nk(num)+1;end
end
fori=1:L
Ps(i)=nk(i)/n;
end
%figure(4),imhist(hist);
%p=imhist(hist);
figure
imshow(a)
count=imhist(a);
[m,n]=size(a);
A=m*n;
L=256;
count=count/A;
fori=1:L
ifcount(i)~=0
st=i-1;
break;
end
end
fori=L:-1:1
ifcount(i)~=0
nd=i-1;
break;
end
end
f=count(st+1:nd+1);
p=st;q=nd-st;
u=0;
fori=1:q
u=u+f(i)*(p+i-1);
ua(i)=u;
end;

fori=1:q
w(i)=sum(f(1:i));
end;

d=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);
th=tp+p;
fori=1:m
forj=1:n
ifa(i,j)>th
a(i,j)=255;
else
a(i,j)=0;
end
end
end
figure
imshow(a);

处理前照片

处理后照片：