clear;close all;
% I=imread('分割0.jpg');
% I=imread('分割1.jpg');
% I=imread('分割2.jpg');
I=imread('分割3.jpg');
figure
imshow(I);
I=rgb2gray(I);
% figure
% imshow(I);
% % 自编程序
Ni=imhist(I); %计算直方图数组
N=sum(Ni); %总像素点个数
delamax=0; %类间方差最大值
threshold=0; %阈值
for k=2:255
u=dot((0:255),Ni/N); %图像的总平均灰度级
w0=sum(Ni(1:k)/N); %C0类像素所占面积的比例
w1=1-w0; %C1类像素所占面积的比例
if w0==0||w0==1 %当w0为1或0时提前结束本次循环
continue
end
u0=dot((0:k-1),Ni(1:k)/N)/w0; %C0类像素的平均灰度
u1=dot((k:255),Ni(k+1:256)/N)/w1; %C1类像素的平均灰度
dela(k)=w0*(u-u0)^2+w1*(u-u1)^2; %类间方差公式
%求出类间方差的最大值，最大时的那个值对应的k值存入delamax
if dela(k)>delamax
delamax=dela(k);
threshold=k-1;
end
end
BW1=im2bw(I,threshold/255); %阈值分割
% figure
% imshow(BW1),title('自编程序运行结果') %显示图片
disp('自编程序求的的阈值：') %显示提示信息
disp(threshold) %显示结果
% %调用系统函数OTSU阈值分割
level=graythresh(I); %调用系统函数计算阈值
BW2=im2bw(I,level); %阈值分割
% figure
% imshow(BW2),title('调用库函数运行结果')
level=round(level*255); %四舍五入计算调用系统函数的阈值
disp('调用库函数求得的阈值：')
disp(level)
%%手工设定阈值
[width,height]=size(I); %获取图片宽高
th=input('请输入阈值:');
BW3=BW1;
for i=1:width
for j=1:height
if(I(i,j)<th) %灰度小于阈值时则为黑色
BW3(i,j)=0;
else % 灰度大于等于阈值时则为白色
BW3(i,j)=1;
end
end
end
% figure
% imshow(BW3),title('全局阈值')
% % 显示
figure
subplot(2,2,1);imshow(I);title('灰度图片');
subplot(2,2,2);imshow(BW1);title('自编程序运行结果');
subplot(2,2,3);imshow(BW2),title('调用库函数运行结果');
subplot(2,2,4);imshow(BW3),title('全局阈值');
