资源描述：
《实验四 图像分割》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验四图像分割一实验目的1．掌握图像分割的目的。2．掌握基于阈值选取的图像分割方法（灰度阈值分割、直方图阈值、最大熵阈值等方法）3．掌握基于区域的图像分割方法。4．掌握基于边缘检测的图像分割方法。5．了解运动图像目标分割方法。二实验原理三实验内容1、测试微小结构：%图像分割测试图像中的微小结构I=imread('cell.tif');%读入原始图像到I变量subplot(2,2,1),imshow(I),title('原始图像');Ic=imcomplement(I);%使用imcomplement函数对图像求反色BW=im2bw(Ic,gray
2、thresh(Ic));%使用im2bw函数对图像进行阈值分割subplot(2,2,2),imshow(BW),title('阈值截取分割后图像');se=strel('disk',6);%创建形态学结构元素，选择一个半径为6个像素的圆盘%形结构元素BWc=imclose(BW,se);%图像形态学关闭运算BWco=imopen(BWc,se);%图像形态学开启运算subplot(2,2,3),imshow(BWco),title('对小图像进行删除后图像');mask=BW&BWco;%对两幅图进行逻辑与操作subplot(2,2,4),im
3、show(mask),title('检测结果图像');2、双峰法%直方图双峰法阈值分割图像程序clearI=imread('cameraman.bmp');%读入灰度图像并显示imshow(I);figure;imhist(I);%显示灰度图像直方图Inew=im2bw(I,140/255);%图像二值化，根据140/255确定的阈值，划分目标与背景figure;imshow(Inew);3、最大熵%基于最大信息熵算法程序clear;closeall;I=imread('goldhill.bmp');%输入原图像figure,imshow(I);
4、%显示原始彩色图像iflength(size(I))==3%如果是彩色图像转换为灰度图像I=rgb2gray(I);%将RGB图像转换为灰度图像end[X,Y]=size(I);V_max=max(max(I));V_min=min(min(I));T0=(V_max+V_min)/2;%初始分割阈值h=imhist(I);%计算图像直方图figure,plot(h);grayp=imhist(I)/numel(I);%求图像像素概率I=double(I);H0=-sum(grayp(find(grayp(1:end)>0)).*log(gray
5、p(find(grayp(1:end)>0))));cout=100;%是指迭代次数为100次while(cout>0)Tmax=0;%初始化graypd=0;graypd=0;Hd=0;Hb=0;T1=T0;A1=0;A2=0;B1=0;B2=0;fori=1:X%计算灰度平均值forj=1:Yif(I(i,j)<=T1)A1=A1+1;B1=B1+I(i,j);elseA2=A2+1;B2=B2+I(i,j);endendendM1=B1/A1;M2=B2/A2;T2=(M1+M2)/2;TT=round(T2);graypd=sum(gra
6、yp(1:TT));%计算分割区域G1的概率和ifgraypd==0graypd=eps;endgraypb=1-graypd;ifgraypb==0graypb=eps;endHd=-sum((grayp(find(grayp(1:TT)>0))/graypd).*log((grayp(find(grayp(1:TT)>0))/graypd)));%计算分割后区域G1的信息熵Hb=-sum(grayp(TT+(find(grayp(TT+1:end)>0)))/graypb.*log(grayp(TT+(find(grayp(TT+1:end)
7、>0)))/graypb));%计算分割后区域G2的信息熵H1=Hd+Hb;cout=cout-1;if(abs(H0-H1)<0.0001)
8、(cout==0)Tmax=T2;break;elseT0=T2;H0=H1;endendTmaxcoutfori=1:X%根据所求阈值Tmax转换图像forj=1:Yif(I(i,j)<=Tmax)I(i,j)=0;elseI(i,j)=1;endendendfigure,imshow(I);%输出图像分割处理后的结果运行结果：Tmax=128.9812cout=96Tmax=137.6849cout=
9、924、全局阈值%MATLAB对图像基于全局阈值的分割程序%functiontest2(a)clearall;a=150;%awie可设
