北航数图实验报告四图像腐蚀和膨胀.docx

北 京 航 空 航 天 大 学 数字图像处理实验报告 实验四：图像分割处理 学 院 专业方向 班 级 学 号 学生姓名 指导教师 实验四 图像分割处理实验 1．实验目的 （1）了解图像分割的基本原理，并利用图像分割算法进行图像分割处理； （2）掌握数学形态学的基本运算。 2．实验内容 （1）利用类间方差阈值算法实现图像的分割处理； （2）利用形态学处理进行处理结果修正。 3. 实验要求 （1）实验用图: （2）对输入图像进行平滑处理，以减小噪声对分割处理的影响； （3）利用类间方差阈值算法对滤波处理后图像进行分割处理，获取分割图像； （4）利用数学形态学中的腐蚀和膨胀运算处理，剔除分割处理结果中的一些细小的残余误分割点，在进行腐蚀和膨胀运算时可采用半径为r的圆形结构元素，注意比较选取不同r值时的处理结果。 四、实验代码 function STshiyan4_OpeningFcn(hObject, eventdata, handles, varargin) I=imread(4.bmp); I=rgb2gray(I); subplot(3,3,1); imshow(I); title(原图像); IM=medfilt2(I);%中值滤波 subplot(3,3,2); imshow(IM); title(中值滤波后图像); function pushbutton1_Callback(hObject, eventdata, handles) I1=imread(4.bmp); I1=rgb2gray(I1); T=Otsu(I1); IM1=medfilt2(I1); s=size(IM1); for m=1:s(1) for n=1:s(2) if IM1(m,n)=T IM1(m,n)=255; else IM1(m,n)=0; end end end subplot(3,3,3); imshow(IM1); title(最大类间方差法); function pushbutton2_Callback(hObject, eventdata, handles) I = imread(4.bmp); %--------------r=1腐蚀-------------------------- se1 = strel(disk,1,8); I1 = imerode(I,se1); subplot(3,3,4); imshow(I1); title(r=1的腐蚀运算); %--------------r=2腐蚀-------------------------- se2 = strel(disk,2,8); I2 = imerode(I,se2); subplot(3,3,5); imshow(I2); title(r=2的腐蚀运算); %--------------r=3腐蚀-------------------------- se3 = strel(disk,3,8); I3 = imerode(I,se3); subplot(3,3,6); imshow(I3); title(r=3的腐蚀运算); function pushbutton3_Callback(hObject, eventdata, handles) I = imread(4.bmp); %--------------r=1膨胀-------------------------- se1 = strel(disk,1,8); I1 = imdilate(I,se1); subplot(3,3,7); imshow(I1); title(r=1的膨胀运算); %--------------r=2膨胀-------------------------- se2 = strel(disk,2,8); I2 = imdilate(I,se2); subplot(3,3,8); imshow(I2); title(r=2的膨胀运算); %--------------r=3膨胀-------------------------- se3 = strel(disk,3,8); I3 = imdilate(I,se3); subplot(3,3,9); imshow(I3); title(r=3的膨胀运算); 五、实验结果 程序运行后，进入到GUI页面，选择相应的按钮就可以显示相应的变换图像结果，如下: