diedaifa.zip_迭代法

clear;close;clc; I=imread('cameraman.tif'); I=double(I); T=(min(I(:))+max(I(:)))/2; %给定初始阈值 done=false; while ~done for m=1:256 for n=1:256 if I(m,n)<=T R1(m,n)=I(m,n); else R2(m,n)=I(m,n); end end end Tnew=(mean(R1(:))+mean(R2(:)))/2; %确定新的阈值 done=abs(Tnew-T)<0.1; %判断 T=Tnew; end for m=1:256 for n=1:256 if I(m,n)<=T J(m,n)=0; %J是二值化后的图像，小于阈值为0，大于阈值为1 else J(m,n)=1; end end end figure('Name','迭代法二值化','NumberTitle','off') subplot(121) imshow(I,[]),title('原图') subplot(122) imshow(J,[]),title('二值化后的图像')