matlab图片固定阙值分割

close all;
clc;
clear;
 
% 1.计算均值、方差
Thre=128;                 %手动设置阈值
I=imread('填充后平滑处理.bmp');   %载入真彩色图像
%I=rgb2gray(I);            %转换为灰度图像
figure;imhist(I);
a=find(I=Thre);
e1=mean(I(a));
e2=mean(I(b));
std1 = std(im2double(I(a)) * 255, 0);
std2 = std(im2double(I(b)) * 255, 0);
 
% 2.绘制两个高斯分布
x=0:1:300;
%y1=gaussmf(x,[std1 e1]);
y1=normpdf(x,e1,std1);
figure,
plot(x,y1)
hold on
%y2=gaussmf(x,[std2 e2]);
y2=normpdf(x,e2,std2);
plot(x,y2)
xlabel('gaussmf, P1=[e1 v1],P2=[e2 v2]')
legend('灰度小于128的灰度值正态分布','灰度大于128的灰度值正态分布')
title('两个高斯分布')

% 3.固定阈值分割方法
I1=imbinarize(I,Thre/256);
figure,imshow(I1),title('固定阈值分割方法（取经验值为128）');
imwrite(I1,'阈值二值化.bmp');

运行：