MATLAB图像分割代码.pdf

[matlab图像处理]阈值分割

%迭代式阈值分割?otsu阈值分割?二值化

closeall;%关闭所有窗口

clear;%清除变量的状态数据

clc;%清除命令行

%迭代式阈值分割?

zmax=max(max(I));%取出最大灰度值

zmin=min(min(I));%取出最小灰度值

tk=(zmax+zmin)/2;

bcal=1;

[m,n]=size(I);

while(bcal)

???%定义前景和背景数

???iforeground=0;

???ibackground=0;

???%定义前景和背景灰度总和

???foregroundsum=0;

???backgroundsum=0;

???fori=1:m

???????forj=1:n

???????????tmp=I(i,j);

???????????if(tmp=tk)

???????????????%前景灰度值

???????????????iforeground=iforeground+1;

???????????????foregroundsum=foregroundsum+double(tmp);

???????????else

???????????????ibackground=ibackground+1;

???????????????backgroundsum=backgroundsum+double(tmp);

???????????end

???????end

???end

???%计算前景和背景的平均值

???z1=foregroundsum/iforeground;

???z2=foregroundsum/ibackground;

???tktmp=uint8((z1+z2)/2);

???if(tktmp==tk)

???????bcal=0;

???else

???????tk=tktmp;

???end

???%当阈值不再变化时,说明迭代结束

end

disp(strcat(迭代的阈值:,num2str(tk)));%在commandwindow里显示出:

迭代的阈值:阈值

newI=im2bw(I,double(tk)/255);%函数im2bw使用阈值（threshold）变换法把

灰度图像（grayscaleimage）

????%转换成二值图像。所谓二值图像，一般意义上是指只有纯黑（0）、纯白

（255）两种颜色的图像。

????%语法

???%BW=im2bw(I,level)

???%BW=im2bw(X,map,level)

???%BW=im2bw(RGB,level)

???%其中level就是设置阈值的。level取值范围[0,1]。

subplot(2,2,2);

imshow(newI);

title(2rice的迭代法分割效果图);

%otsu阈值分割

bw=graythresh(I);

disp(strcat(otsu阈值分割的阈值:,num2str(bw*255)));%在commandwindow

里显示出:迭代的阈值:阈值

newII=im2bw(I,bw);

subplot(2,2,3);

imshow(newII);

title(3rice的otsu阈值分割);

%二值化阈值为135

[width,height,bmsize]=size(I);

fori=1:width

???forj=1:height

???????ifI(i,j)135

???????????I(i,j)=255;

???????else?

???????????I(i,j)=0;

???????end

???end

end???

subplot(2,2,4);

imshow(I);

title(4rice的二值阈值分割);