1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
| figure(1);subplot(2,1,1);
I=imread('Image_ombre.png');
imshow(I);
I1 = im2bw(I, graythresh(I)) ;
M=size(I,1);N=size(I,2);
hist=zeros(1,256);
for i=1:M
for j=1:N
hist(I(i,j)+1)=hist(I(i,j)+1)+1;
end;
end;
hist= double(hist);
for i=1:256
Prb(i)=hist(i)/(M*N);
end;
omega(1)=Prb(1);
myu(1)=0;
for i=2:256
omega(i)=omega(i-1)+Prb(i);
myu(i)=myu(i-1)+ i*Prb(i);
end;
Thr=0;
max_sigma=0;
for i=1:255
if(omega(i)~= 0 && omega (i)~=1)
sigma(i)=power(myu(255)*omega(i)-myu(i),2)/(omega(i)*(1- omega(i)));
else sigma(i)=0;
end;
if(sigma(i)>max_sigma)
max_sigma=sigma(i);
Thr=i;
end;
end;
for i=1:M
for j=1:N
if (I(i,j)> Thr)
I(i,j,:)=0;
else I(i,j,:)=255;
end;
end;
end;
figure(1);subplot(2,1,2);imshow(I); |
Partager