图像分割matlab阈值自动计算,求matlab图像直方图阈值分割算法（包含预处理步骤）...

在此提供两种方法,一,最大阈值分割(大律法)：

I=imread('test.jpg');

subplot(1,3,1);imshow(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

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('阈值分割图');