matlab otsu算法代码,OTSU算法matlab程序代码

function th=thresh_md(a);

x=imread(a);

a=rgb2gray(x);

imshow(a)

count=imhist(a);

[m,n]=size(a);

N=m*n;

L=256;

count=count/N;

for i=1:L

if count(i)~=0

st=i-1;

break;

end

end

for i=L:-1:1

if count(i)~=0

nd=i-1;

break;

end

end

f=count(st+1:nd+1); %fê??????ò?è3???μ????ê

p=st; q=nd-st;

u=0;

for i=1:q

u=u+f(i)*(p+i-1); %uê?????μ????ù?μ

ua(i)=u; %ua￡¨i￡?ê??°i??????μ????ù?ò?è?μend;

for i=1:q

w(i)=sum(f(1:i)); %w￡¨i￡?ê??°i??????μ?à??ó???êend;

d=(u*w-ua).^2./(w.*(1-w));

[y,tp]=max(d); %?éò?è?3?êy×éμ?×?′ó?μ?°è?×?′ó?μμ?μ?th=tp+p;

for i=1:m

for j=1:n

if a(i,j)>th

a(i,j)=0;

else

a(i,j)=255;