A=imread('image.jpg'); %读入图像
B=rgb2gray(A);
B=double(B);
hist(B) %求图像的灰度直方图
[m,n]=size(B);
for i=1:m %根据直方图进行阈值分割
for j=1:n
if B(i,j)>70&B(i,j)<130%阈值
B(i,j)=1;
else
B(i,j)=0;
end
end
end
subplot(121),imshow(A) %显示分割结果
subplot(122),imshow(B)