局部阈值分割法matlab，求matlab图像直⽅图阈值分割算法
（包含预处理步骤）
在此提供两种⽅法,⼀，最⼤阈值分割(⼤律法)：
I=imread('test.jpg');
subplot(1,3,1);(I);title('原图');
I1=rgb2gray(I);
subplot(1,3,2);imhist(I1);title('');
level=graythresh(I);
g=im2bw(I,level);%最佳阈值level
subplot(1,3,3);imshow(g);title('阈值分割');
%在Command窗⼝输出灰度阈值
disp(strcat('graythresh 计算灰度阈值：',num2str(uint8(level*255))))
第⼆种：迭代求阈值
I=imread('C:test.jpg');
Zmax=max(max(I));
Zmin=min(min(I));
TK=(Zmax+Zmin)/2;%初始阈值
flag=1;
[m,n]=size(I);
while(flag)
fg=0;
bg=0;
fgsum=0;
bgsum=0;
for i=1:m
for j=1:n
tmp=I(i,j);
if(tmp>=TK)
fg=fg+1;
fgsum=fgsum+double(tmp);
else
bg=bg+1;
bgsum=bgsum+double(tmp);
end
end
end
u1=fgsum/fg;
u2=bgsum/bg;
TKTmp=uint8((u1+u2)/2);
if(TKTmp==TK)
flag=0;
else
TK=TKTmp;
end
matlab直方图
end
disp(strcat('迭代后的阈值：',num2str(TK))); newI=im2bw(I,double(TK)/255);
subplot(1,2,1);imshow(I);title('原图')
subplot(1,2,2);imshow(newI);tltle('阈值分割图');

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。