imshow函数matlab_大津阈值法matlab实现

clear all
I &＃61; imread(&＃39;1.jpg&＃39;);
%先把彩色图变灰度图%figure(1)
subplot(221)
imshow(I);
title(&＃39;原图&＃39;);
I&＃61;rgb2gray(I);
subplot(222)
imshow(I);
title(&＃39;灰度图&＃39;);
level &＃61; graythresh(I);
BW &＃61; imbinarize(I,level);
subplot(223)
%figure(2)
imshow(BW);
title(&＃39;matlab自带函数&＃39;);
IMAX &＃61; max(max(I));
IMIN &＃61; min(min(I));T&＃61;IMIN:IMAX;
ISIZE&＃61;size(I); %图像大小muxSize &＃61; ISIZE(1) * ISIZE(2);Tmp&＃61;0;%zeros(1,length(T));
mid&＃61;0;
for i&＃61;1:length(T) %从MIN到MAXTK &＃61;T(1,i);ifground&＃61;0;ibground&＃61;0;FgroundS&＃61;0;BgroundS&＃61;0;for j&＃61;1:ISIZE(1)for k&＃61;1:ISIZE(2)tmp &＃61; I(j,k);if(tmp>&＃61;TK)ifground&＃61;ifground&＃43;1;%像素个数FgroundS&＃61;FgroundS&＃43;double(tmp);%总的前景色灰度elseibground&＃61;ibground&＃43;1;Bground&＃61;BgroundS&＃43;double(tmp);endendendw0&＃61;ifground/muxSize;w1&＃61;ibground/muxSize;u0&＃61;FgroundS/ifground;u1&＃61;BgroundS/ibground;tmp&＃61;w0*w1*(u0-u1)*(u0-u1);if(Tmpendfor j&＃61;1:ISIZE(1)for k&＃61;1:ISIZE(2)if(I(j,k)>&＃61;mid)I(j,k)&＃61;255;elseI(j,k)&＃61;0;endend
end
subplot(224)
%figure(3)
imshow(I);
title(&＃39;编程实现大津阈值法&＃39;);