资源描述：
《数学形态学图像处理的基本运算实现及分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、数学形态学图像处理的基本运算实现及分析—、基本原理数学形态学是一种应用于图像处理和模式识别领域的新的方法。它的基本思想是用具有一定形态的结构元素去度量和提取图像中的对应形状以达到对图像进行分析和识别的目的。数学形态学的数学基础和所用语言是集合论。数学形态学的应用可以简化图像数据，保持它们基本的形状特性，并除去不相干的结构。另一方面，数学形态学的算法具有天然的并行实现的结构。1、基本运算数学形态学的基木运算有四个：膨胀、腐蚀、开启和关。如用A表示图像集合，B表示结构元素，形态学运算就是用B对A进行操作。A被B膨胀，记为A㊉B,㊉为膨胀算子，膨胀的定义为A㊉B={x
2、[(B)vnA]/0
3、}该式表明的膨胀过程是B首先做关于原点的映射，然后平移X。A被B的膨胀是B被所有x平移后与A至少有一个非零公共元素。A被B腐蚀，记为AOB,。为腐蚀算子，腐蚀的定义为AeB={x[(B)xnA]^0}也就是说，A被B的腐蚀的结果为所有使B被x平移后包含于A的点x的集合。换句话说，用B来腐蚀A得到的集合是B完全包括在A屮时B的原点位置的集合。膨胀和腐蚀并不互为逆运算，所以它们可以级连结合使用。例如，利用同一个结构元素B,先对图像腐蚀然后膨胀其结果，或先对图像膨胀然后痫蚀其结果,前一种运算称为开运算，后一种运算称为关运算。它们也是数学形态学中的重要运算。开启的运算符为o,A用B来开启写
4、作AoB,其定义为：AoB=(AQB)®B关的运算符为•，A用B来关写作A-B,其定义为：A•B=(A㊉3)03开和关两种运算都可以去除比结构元素小的特定图像细节，同时保证不产生全局的几何失真。开运算可以把比结构元素小的椒盐噪声滤除，切断细长搭接而起到分离作用。关运算可使比结构元素小的缺口或孔填补上，搭接短的间断而起到连通作用。2、实际应用近年来，数学形态学在图像处理方面得到了日益广泛的应用。下面主要就数学形态学在边缘检测、骨架提取等方面的应用做简要介绍。集合A的边界记为0(A),可以通过下述算法提取边缘：设B是一个合适的结构元素，首先令A被B腐蚀，然后求集合A和它的腐蚀的差。如下式
5、所示:0(A)*_(AG)B)下解释了边缘提取的过程。它表示了一个简单的二值图像，一个结构元素和用以上公式得出的结果。图(b)屮的结构元素是最常用的一种，但它决不是唯一的。如果采用一个5X5全“1”的结构元素，可得到一个二到三个像素宽的边缘。应注意的是，当集合B的原点处在集合的边界时，结构元素的一部分位于集合之外。这种条件下的通常的处理是约定集合边界外的值为0。(b)(d)利用形态学方法提取一个区域的骨格可以用腐蚀和开运算表示。也就是，A的骨骼记为S(A),骨骼化可以表示如下：S(A)=J5,(A)k=0和Ks(A)=U{(AQkB)-[AQkB)。B]}k=0其中B是结构元素，(A
6、ORB)表示对A连续腐蚀k次。下图的解释说明了以上讨论的概念。第一列显示了原始集合(顶部)和通过结构元素B两次腐蚀的图形。由于再多一次对A的腐蚀将产生空集，所以选取K=2o第二列显示了第一列通过B的开运算而得到的图形。以上结果可以通过以前讨论过的开运算拟合性质加以解释。第三列仅仅显示出第一列与第二列的差别。第四列包含两个部分骨骼及最后的结果(第四列的底部)。最后的骨骼不但比所要求的更粗，而且相比较更重要，它是不连续的。形态学给出了就特定图形侵蚀和空缺的描述。形态学给出了就特定图形侵蚀和空缺的描述。通常，骨骼必须最人限度的细化、相连、最小限度的腐蚀。第五列显示了S°(A)、&(A)㊉B
7、以及©(A)㊉2B)=S2(A)㊉B)㊉B。最后一列显示了图像A的重构。可知，A就是第五列中膨胀骨骼子集的“并”。厂4JJr二、源程序1、膨胀bw=imrcad(Jtext・png');se=strelCline，,11,90);%生成线形的结构元素bw2=imdilate(bw,se);%对图像进行膨胀subplot(121),imshow(bw),title(，原始图像');subplot(122),imshow(bw2),膨胀后的图像');2、腐蚀originalBW=imread(，circles.png');se=strelCdisk*,11);%生成圆盘形的结构元素ero
8、dedBW=imerode(originalBW,se);%对图像进行腐蚀subplot(121),imshow(origina1BW),title('原始图像');subplot(122),imshow(erodedBW);titlc(，腐蚀后的图像');3、开BW1=imread(，circbw.tif，);%读取图像subplot(121),imshow(BW1);title(，原始的图像');se二strel('rectangle',[4030])
