MATLAB图像处理.rar_86X_图像处理_机器人视觉_机器视觉

大津法 I=imread('coins.png'); subplot(131),imshow(I); title('原始图像') level=graythresh(I); BW=im2bw(I,level); subplot(132),imshow(BW) title('graythresh计算阈值') disp(strcat('graythresh计算灰度阈值：',num2str(uint8(level*255)))) iMax=max(max(I)); iMin=min(min(I)); T=double(iMin:iMax); iSize=size(I); muxSize=iSize(1)*iSize(2); for i=1:length(T) TK=T(1,i); iForeground=0; iBzckground=0; ForegroundSum=0; BzckgroundSum=0; for j=1:iSize(1) for k=1:iSize(2) tmpData=I(j,k); if(tmpData>=TK) iForeground=iForeground+1; ForegroundSum=ForegroundSum+double(tmpData); else iBackground=iBackground+1; BackgroundSum=BackgroundSum+double(tmpData); end end end w0=iForeground/muxSize; w1=iBackground/muxSize; u0=ForegroundSum/iForeground; u1=BackgroundSum/iBackground; T(2,i)=w0*w1*(u0-u1)*(u0-u1); end oMax=max(T(2,:)); idx=find(T(2,:)>=oMax); T=uint8(T(1,idx)); disp(strcat('简化大津法计算灰度阈值：',num2str(T))) BW=im2bw(I,double(T)/255); subplot(133),imshow(BW) title('简化大津法计算灰度阈值')