数字图像处理实验-形态学.doc

《数字图像处理实验-形态学.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西安邮电学院实验报告实验名称形态学图像处理课程名称数字图像处理A姓名方健成绩班级电子0802学号（01）日期2011-05-31地点3#531备注：仅供参考不予下载1.实验目的（1）了解并掌握膨胀、腐蚀及开运算、闭运算的基本原理；（2）编写程序使用开运算、闭运算处理图像，进一步理解开运算、闭运算的实质；（3）编写程序使用开运算、闭运算进行图像去噪处理，根据实验结果分析效果；（4）总结实验过程（实验报告，左侧装订）：方案、编程、调试、结果、分析、结论。2.实验环境（软件条件）WindowsXPMATLAB7.x3.实验方法对两幅受噪声干扰的数字图像curve_128.bmp（如图3.

2、1所示）和enoise.bmp（如图3.2所示）进行如下处理：图3.2实验图像enoise.bmp图3.1实验图像mili.bmp（1）对两幅图像进行腐蚀、膨胀处理，显示处理前、后图像：可以用不同尺度及形状的结构元素进行腐蚀、膨胀处理，分析结构元素对处理效果的影响；（2）分别对两幅图像设计相应的开、闭运算进行降低噪声的处理，显示处理前、后图像；（3）分析两幅图像为了降低噪声所设计的运算有无不同？如果有请分析为什么会有这种不同？4.实验分析A基本概念：A.膨胀已知二值图像A，如果Ab1,Ab2,…，Abn是由二值图像B={b1,b2,b3,…，bn}中像素值为1的点平移得到，则A由B

3、平移的并称为A被B膨胀。B.腐蚀腐蚀是膨胀的逆运算。二值图像A经二值图像B腐蚀后在p点仍为1的充分必要条件电子0802-01是：B平移到B后，B中的1像素也是A中的1像素。C.开运算用同一结构元腐蚀后在膨胀可去除比结构元小的所有区域像素点，而留下其余部分，这一顺序称为“开”运算。D.闭运算与开运算顺序相反的过程是先膨胀后再腐蚀，称为“关”运算或“闭”运算。B所用函数：构造结构元素：SE=strel(shape,parameters)；腐蚀运算：IM2=imerode(IM,SE)；膨胀运算：IM2=imdilate(IM,SE)；开运算：IM2=imopen(IM,SE)闭运算：I

4、M2=imclose(IM,SE)A.实验结论A.对mili.bmp进行腐蚀膨胀操作。1.使用半径为2像素的圆作为结构元素：2.使用2*2方形作为结构元素：B.对enoise.bmp进行腐蚀膨胀操作。1.使用半径为1像素的圆作为结构元素：2.使用边长为2像素的竖直线作为结构元素：分析：对“米粒图”经行腐蚀后，米粒边缘被腐蚀掉而变得细长，且灰白色噪声得到一定的消除。膨胀后，米粒边缘被扩张而变得饱满，但也将灰白色噪点放大了。“文字图”腐蚀后，黑色噪点消除，文字形态遭到破坏，膨胀后，黑色噪点张大，文字变得突出醒目。使用不同的结构元素对图像处理的结果是有一定差异的，从上图来看，使用面结构元

5、素还是线结构元素，都反映到噪点的形状上了。比如线结构元素使原本点状噪声变得细长。C.对mili.bmp经行开闭运算。A.对enoise.bmp经行开闭运算。分析：对于“米粒图”，开操作是先腐蚀再膨胀，因为腐蚀使白色噪点消除，膨胀又使米粒恢复原来的形态，因此平滑的效果好一些。先比开操作，闭操作并没有很好地消除噪声。对于二值“文字图”，由于目标文字和噪点都是黑色的，背景是白色的，这和典型的开闭所使用的图像是相反的，因此需要将二值原图经行反转，才能经行正常开闭操作，而得到的图像也要经行一次反转才是上图所示的。和“米粒图”类似，开操作很好地消除了噪声，但多次的开操作和一次的效果是一样的。开

6、和闭的线性组合使噪声进一步消除，目标进一步变得突出。附件电子0802-011.基本代码A.腐蚀膨胀A=imread('mili.bmp');subplot(231)imshow(A);title('原图');SE=strel('disk',2);A2=imerode(A,SE);subplot(232);imshow(A2);title('腐蚀');A3=imdilate(A,SE);subplot(233);imshow(A3);title('膨胀');B.开闭运算clc;clear;A=imread('enoise.bmp');A=im2bw(A);subplot(231);i

7、mshow(A);title('原图');A=abs(A-1);SE=strel('square',2);fo=imopen(A,SE);subplot(232);imshow(abs(fo-1));title('一次开运算');foc=imclose(A,SE);subplot(233);imshow(abs(foc-1));title('一次闭运算');foct=imclose(fo,SE);subplot(234);imshow(abs(foct-1));tit