数学形态学去噪要点.docx
《数学形态学去噪要点.docx》由会员分享,可在线阅读,更多相关《数学形态学去噪要点.docx(15页珍藏版)》请在冰豆网上搜索。
数学形态学去噪要点
一绪论.......................................................
1
数学形态学简介..............................................
1
数学形态学与数字图像处理....................................
1
本次课程设计的目的与要求....................................
2
二数学形态学的根本运算........................................
3
根本概念....................................................
3
结构元素..................................................
3
膨胀与腐蚀................................................
3
二值形态学图像处理..........................................
4
膨胀......................................................
4
腐蚀......................................................
6
开运算....................................................
7
闭运算....................................................
8
灰度形态学图像处理..........................................
9
膨胀......................................................
9
腐蚀.....................................................
10
开运算与闭运算...........................................
11
综述.......................................................
13
三数学形态学滤波器去噪.......................................
15
概述.......................................................
15
噪声模型...................................................
16
高斯噪声.................................................
16
椒盐噪声.................................................
16
形态学滤波器...............................................
17
形态学图像去噪原理.........................................
20
形态学图像去噪的应用.......................................
20
小结与体会....................................................
21
参考文献.......................................................
22
附录
..........................................................
23
一绪论数学形态学简介数学形态学作为一门新兴的图像处理与分析学科,1964年由法国的和在积分几何的根底上首次创立。
70年代初,采用数学形态学的学者们开拓了图像分析的一个新的领域。
经过十多年的理论与实践探索,和等人在研究中认识到,对图像先作开运算接着再作闭运算,可以产生一种幂等运算;采用递增尺寸的交变开闭序列作用于图像,可有效地消除图像的噪声,1982年他们正式提出了形态学滤波器的概念。
90年代数学形态学有两个显著的开展趋势,第一个是致力于运动分析,包括编码与运动景物描述;第二个是算法与硬件结构的协调开展,用于处理数值函数的形态学算子的开发与设计。
目前国内许多有效的图像处理系统有的是基于数学形态学方法原理设计的,有的是把数学形态学算法纳入其根本软件,并以其运算速度作为系统性能的重要标志之一数学形态学与数字图像处理数学形态学在图像处理中属于非线性滤波方法,现在数学形态学的方法已经开展成为图象处理技术的一个重要方面,并且被广泛的应用到图象处理的各个领域,利用数学形态学可以进行图像去噪、图象分割、增强、边缘检测、形态分析、图象压缩等各个方面。
可以通过以下几个步骤来实现数学形态学算法对数字图像的处理:
步骤1、提取图像的几何结构特征,也就是针对所要处理的图像找出相应的几何结构模式。
步骤2、根据步骤1找出的几何结构模式选适宜的结构元素,这里结构元素的选择标准择首先是要能最有效的展现该几何结构模式,其次该结构元素的形态还应该尽量的最简。
步骤3、为了得到比原始图像更能显著突出物体特征信息的图像,用步骤2选取的结构元素对目标进行相应的数学形态学变换,如果能对结构元素给予适宜
1
的变量,那么还能够定量的表示出目标的几何结构模式。
步骤4、通过上面的三个步骤,相对于我们的处理需求,目标图像会变得更加清晰、明了,并且更有利于我们提取出相应的图像信息。
现在,数学形态学处理图像已经开展成为一个专门的图像科学领域。
该领域已经形成了一个理论概念、非线性滤波、设计算法以及应用系统相互连贯而有广阔的整体。
与其他很多图像处理技术相比,数学形态学技术的理论框架完善、算法效率高、易于在专门硬件上使用并且适合处理很多与形状相关的问题。
例如对于图像噪声去由于可以在去除噪声前有效的探究目标图像的几何结构模式,尽可能的解决去除噪声与保护图像边缘细节信息相冲突的根本矛盾。
再如在提取图像边缘时,与其他算法相比,数学形态学方法提取的边缘更为连续,间断点也会少很多。
所以很多学术机构及工业研究所在处理数字图形图像、计算机视觉、模式识别等很多问题时都会重点考虑数学形态学方法。
本次课程设计的目的与要求〔1〕通过形态学方面的知识处理各种图像。
〔2〕学会应用形态学知识处理加有高斯噪声与椒盐噪声的图像。
〔3〕理解不同的形态学运算在处理图像方面的应用。
〔4〕通过运用MATLAB软件实现仿真。
2
二数学形态学的根本运算根本概念数学形态学是由一组形态学的代数运算子组成的,它的根本运算有4个:
膨胀〔或扩张〕、腐蚀〔或侵蚀〕、开运算和闭运算,它们在二值图像和灰度图像中各有特点。
基于这些根本运算还可推导和组合成各种数学形态学实用算法,用它们可以进行图像形状和结构的分析及处理,包括图像分割、特征抽取、边界检测、图像滤波、图像增强和恢复等。
数学形态学方法利用一个称作结构元素的“探针〞收集图像的信息,当探针在图像中不断移动时,便可考察图像各个部分之间的相互关系,从而了解图像的结构特征。
数学形态学基于探测的思想,与人的FOA(FocusOfAttention)的视觉特点有类似之处。
作为探针的结构元素,可直接携带知识〔形态、大小、甚至参加灰度和色度信息〕来探测、研究图像的结构特点。
结构元素所谓结构元素就是一定尺寸的背景图像,通过将输入图像与之进行各种形态学运算,实现对输入图像的形态学变换。
结构元素没有固定的形态和大小,它是在设计形态变换算法的同时根据输入图像和所需信息的形状特征一并设计出来的,结构元素形状、大小及与之相关的处理算法选择得恰当与否,将直接影响对输入图像的处理结果。
通常结构元素的形状有正方形、矩形、圆盘形、菱形、球形以及线形等。
膨胀与腐蚀膨胀在数学形态学中的作用是把图像周围的背景点合并到物体中。
如果两个物体之间距离比拟近,那么膨胀运算可能会使这两个物体连通在一起,所以膨胀对填补图像分割后物体中的空洞很有用。
腐蚀在数学形态学运算中的作用是消除物体边界点,它可以把小于结构元素的物体去除,选取不同大小的结构元素可以去掉不同大小的物体。
如果两个物体之间有细小的连通,当结构元素足够大时,通过腐蚀运算可以将两个物体分开。
3
二值形态学图像处理二值图像数学形态学的运算就是基于上述集合论的理论,进行击中与否变换HMT〕,在定义了HMT及其根本运算膨胀〔Dilation〕和腐蚀(Erosion)后,再从积分几何和体视学移植一些概念和理论,根据图像分析的各种要求,构造出统一的、相同的或变化很小的结构元素进行各种形态变换。
数学形态学中有两种最根本的操作即膨胀和腐蚀,其他的所有形态学操作都是基于这两个操作的组合或级联。
如开运算和闭运算就都是膨胀和腐蚀的最根本组合。
膨胀、腐蚀、开运算和闭运算构成了整个数学形态学变换的根底,下面分别对这四种根本形态学变换进行具体的分析。
二值图像腐蚀膨胀图膨胀与腐蚀示意图
膨胀膨胀操作是指一个集或对象目标从其原来的形状扩大的过程。
该目标扩大的方式是由结构元素决定的。
和待处理的对象相比拟,结构元素的大小更小,一般用于膨胀的结构元素大小取到3×3。
膨胀的过程类似于卷积,结构元素在目标图像内从左到右、从上到下的移动,在每次移动的过程之中,都会寻找结构元素与目标对象之间重叠的像素,只要存在重叠的像素点,结构元素所在的中心位置点的像素值都会被标为1。
用集合论该过程可表示如下:
A,B为Z2中的集合,Φ为空集,A被B的膨胀,记为A⊕B,⊕为膨胀算子,膨胀的集合定义式为:
〔〕该式表示的膨胀过程首先是B做关于原点的映射,然后平移x。
A被B膨胀
4
也就是被所有x平移后于A至少有一个公共非零元素的集合。
根据上述对膨胀过程的解释,公式(2.2.1)也可以被写作下面的形式:
〔〕和其他的形态学操作一样,公式中集合B在膨胀运算中一般被叫做结构元素。
膨胀运算的实质是遍历待膨胀图像中的每个像素点,根据所选取的结构元素的值以及要处理像素点周围点的灰度值进行计算。
比拟局部范围内的像素点与结构元素中所对应点的灰度值之和。
根据比拟的结果,选取所计算的这些和中的最大值。
所以经过膨胀,图像边缘的像素点灰度值会增加,图像边缘向外扩张,最终到达图像膨胀的视觉效果。
不同的数学形态学文献对膨胀都有着不同的定义,公式(2.2.1)不是现在形态学文献中膨胀的唯一定义。
然而,相比其它定义,这个定义存在一个明显的好处,即当把结构元素B被当做卷积模板时,膨胀的概念会更加的形象化。
因为虽然膨胀的本质是集合运算,而卷积本质上属于算术运算,但由于结构元素B做相对于原点的“映射〞后在集合A(图像A)上的平移是连续的,因此可以近似的将它滑过集合A的整个过程近似看做卷积过程。
图膨胀操作
5
腐蚀腐蚀可以看做膨胀的逆运算或反过程。
如果说膨胀是扩张了图像,那么腐蚀的作用那么是使图像收缩。
图像目标收缩的方式也是有结构元素决定的。
和膨胀一样,腐蚀所选取的结构元素也要比目标要小,一般也是取3×3的大小。
选取结构元素的尺寸较小的好处是可以减少腐蚀算法运行的时间。
和膨胀相似,腐蚀操作也是将结构元素从左到右、从上到下在待处理图像中移动,以结构元素的中心点作为运算的中心,检验图像周围像素是否与结构元素完全重合。
只要没有完全重叠,那么该中心点像素就被标为0。
该过程同样可以用集合论的方法表示如下:
A,B为Z2中的集合,A被B腐蚀,记为AΘB,其定义公式为:
〔〕也就是说A被B的腐蚀的结果为所有使B被x平移后包含于A的点x的集合。
和膨胀一样,腐蚀运算的实质也是遍历待腐蚀图像中的每个像素点,根据所选取的结构元素的值以及要处理像素点周围点的灰度值进行计算。
比拟局部范围内的像素点与结构元素中所对应点的灰度值之差。
根据比拟的结果,选取所计算的这些差中的最小值。
经过腐蚀,图像边缘的像素点灰度值会降低,从而图像边缘会向内收缩,最终到达腐蚀的视觉效果。
膨胀和腐蚀的关系可以看做集合补和反转的对偶,可以用下面的公式表示:
〔〕图腐蚀操作
6
对一幅图像进行膨胀和腐蚀操作结果如下列图:
原始图像阈值为的图像腐蚀后图像1
膨胀后图像1腐蚀后图像2膨胀后图像2图这两种运算具有对偶性,即一种运算对目标的操作相当于另一种运算对图像背景的操作。
由图可知,膨胀操作后,图像相当于膨胀了一圈,而腐蚀操作后,图像那么像被剥掉了一层。
开运算
设A是原始图像,B是结构元素图像,那么集合A被结构元素B做开运算,记为:
AoB,其公式为:
AoB=(AΘB)⊕B(2.5)从公式(3.2.13)可以看出A被结构元素B做开运算就是A被B腐蚀后的结果再被B膨胀。
同膨胀和腐蚀一样,我们也可以用用集合论的概念来定义开运算。
A被结构元素B做开运算就是B在A内的平移所得到的集合的并集,即:
〔〕开运算一般能平滑图像的轮廓,削弱狭窄的局部,去掉细的突出。
根据开运算的这种作用,我们可以利用开运算来去除图像的噪声。
7
图开运算示意图
闭运算
设A是原始图像,B是结构元素图像,那么集合A被结构元素B做闭运算,记为:
A?
B,可用下面公式给出闭运算的定义:
A?
B=(A⊕B)ΘB〔〕从公式(3.2.15)可以看出A被结构元素B做闭运算就是A被B膨胀后的结果再被B腐蚀。
闭运算的效果根本上与开运算相反,它一般是能填充背景中狭窄以及凹陷的局部,消除小洞,还可以填充轮廓上空白的局部,使图像的轮廓得到一定的平滑。
我们同样可以利用闭运算的这些作用来去除图像的噪声。
和膨胀和腐蚀一样,开运算和闭运算也是关于集合补和反转的对偶,可用下面公式表示这种关系:
〔〕
对一幅二值图像进行开运算与闭运算的结果如下列图:
8
原始图像阈值为的二值图像开运算后图像1
闭运算后图像1开运算后图像2闭运算后图像2
图由图可知,开运算能够去除孤立的小点、毛刺和小桥,而闭运算能够填平小孔,弥合小裂缝。
灰度形态学图像处理前面所提到的形态学方法都是基于二值图像的,下面把这些根本的形态学操
作推广到灰度图像。
这里把设为输入图像,为结构元素。
膨胀
用b对函数f进行的灰度膨胀表示为f⊕b,定义式为:
〔〕其中,分别是f和b的定义域,这里的f和b表示的是函数而
不再是二值形态学中所说的集合。
在灰度膨胀中,每个结构元素的位置上,在这一点的膨胀值是在跨度b的区间内f于b之和的最大值。
为通常对灰度图像进行膨胀处理会把图像中白色的局部扩张,而且会带来双重的效果:
(1)假设使用的所有结构元素的值均为正,
9
那么输出图像会的亮度将会大于输入图像的亮度;
(2)膨胀过程中所用的结构元素的值和形状决定图像中暗的细节是局部全部减少了还是被消除。
腐蚀
同样的道理,用b对函数f进行的灰度腐蚀表示为fΘb,定义式为:
〔〕
和灰度膨胀一样,这里,分别是f和b的定义域。
在灰度腐蚀中,每个结构元素的位置上,在这一点的腐蚀值是在跨度b的区间内f于b之差的最小值。
为通常对灰度图像进行腐蚀处理会把图像中黑色的局部扩张,而且会带来双重的效果:
(1)假设使用的所有结构元素的值均为正,那么输出图像会的亮度将会小于输入图像的亮度;〔2〕在输入图像中如果结构元素的面积比亮的细节的面积大,那么亮的效果将被削弱。
环绕于亮细节周围的灰度值和结构元素自身的形态结构和幅值决定着最终消弱的程度。
灰度膨胀和腐蚀与二值膨胀和腐蚀一样也满足对偶关系。
对一幅灰度图像进行膨胀与腐蚀的结果如下列图所示:
原始图像腐蚀后图像1膨胀后图像1
腐蚀后图像2膨胀后图像2
图
10
由图可知,灰值图像的腐蚀,如果结构元素的值都为正,那么输出图像会比输入图像暗;如果输入图像中亮细节的尺寸比结构元素小,那么其影响会被减弱,减弱的程度取决于这些亮细节周围的灰度值的结构元素的形状和幅值。
灰度图像的膨胀运算,如果结构元素的值都为正,那么输出图像会比输入图像亮,根据输入图像中暗细节的灰度值以及它们的形状相对于结构元素的关系,它们在膨胀中或被消减或被除掉。
开运算与闭运算
用结构元素b对图像f进行开操作即先用b对f进行腐蚀操作再进行膨胀,记为fob,定义式为:
〔〕用结构元素b对图像f进行闭操作即先用b对f进行膨胀操作再进行腐蚀,记为f?
b,定义式为:
f?
b=(f⊕b)Θb〔〕用几何的方法可以解释灰度形态学的开、闭运算。
用b对f进行开操作的原理可以解释为:
顺着曲面的下侧面推动小球上下滚动,使球体能够在曲面的整个下侧面来回移动。
这时球体外表的所有接触到曲面局部的最高点就组成了开操作的曲面。
相对应的,闭操作就是在曲面的上侧面滚动,以便球体能在曲面的整个上侧面来回移动。
当球体的任何局部接触到曲面的最高点就构成了闭操作的曲面。
如下列图就详细的解释了这一过程:
11
图灰度形态学的开、闭运算和二值形态学一样都满足对偶关系。
并且灰度开、闭运算的作用效果也根本上于二值形态学中的一样。
其中灰度开运算一般能平滑图像的轮廓,消减图像中狭窄的局部,抹平突出的细节。
灰度闭运算的效果根本上与灰度开运算相反,它一般是能填充背景中狭窄以及凹陷的局部,消除小洞,还可以填充轮廓上空白的局部,使图像的轮廓得到一定的平滑。
我们同样可以利用灰度开、闭运算的这些作用来去除图像中不希望看到的掺杂细节〔如噪声〕对一幅加有椒盐噪声的灰度图像进行开运算与闭运算的结果如下:
原始图像参加椒盐噪声开运算后图像1
闭运算后图像1开运算后图像2闭运算后图像2图开运算操作消除与结构元素相比尺寸较小的亮细节,而保持图像整体灰度值和大的亮区域根本不受影响。
闭运算操作消除与结构元素相比尺寸较小的暗细
12
节,而保持图像整体灰度值和大的暗区域根本不受影响。
所以,开-闭运算可用于图像的去噪,特别是对于加有椒盐噪声的图像的去噪效果更为明显。
详细过程见下一章。
对一幅加有高斯噪声的图像进行开运算与闭运算的结果如下:
原始图像参加高斯噪声开运算后图像1
闭运算后图像1开运算后图像2闭运算后图像2图单独运用开运算或闭运算对一幅加有高斯噪声的图像进行去噪效果并不理想。
综述本章首先对数学形态学的大致开展历程及根本原理进行了概括性的介绍。
说明了数学形态学在图像处理中的应用领域、处理步骤及其在图像处理中的独有优势。
数学形态学是用数学集合论中的语言来描述的,也即集合论是数学形态学的根底,所以本章对集合论的主要内容做了介绍。
然后基于集合论的根本原理,本章重点讨论了形态学中最根本的局部——二值形态学。
用集合论的根本理论给出了二值形态学的四种根本运算——膨胀、腐蚀、开运算和闭运算的公式;说明了这四种根本二值形态学变换在处理图像时的运算原理;通过实验〔主要是去噪实验〕给出了四种形态学变换对图像的处理效果并对其处理噪声的效果进行了讨论。
分别对膨胀和腐蚀、开运算和闭运算的关系以集合论公式的形式给予了介绍。
然后本章把二值形态学的根本理论推广到灰度图像,以交叉比照的形式把灰度形
13
态学中膨胀、腐蚀、开运算及闭运算的根本集合论公式、原理等进行了介绍。
并同样通过实验说明了膨胀和腐蚀对图像的处理效果以及开运算和闭运算在去除图像噪声时所表现出的特点。
通过本章的介绍,可以明确数学形态学的根本理论;对数学形态学在图像处理,特别是去除噪声时的优势有了进一步的说明。
本章是下一章关于形态学去噪的理论根底,特别是其中的开-闭运算,在图像去噪方面有着十分重要的应用。
14
三数学形态学滤波器去噪概述通过图像去噪〔也可称之为图像清洁〕来对图像进行增强是图像处理中的一个根本问题。
于这个问题,研究者们已经进行了长达多年的研究。
大多数的图像处理书籍中都会有涉及到图像去噪的章节。
从本质上讲,图像去噪属于图像恢复问题中的一类,目的是将一个受到噪声侵蚀的图像版本恢复到其原先完善的版本。
这个问题看似很简单,但是由于我们常常无法对恢复后图像的实际结构作出合理的分析,而会变得异常复杂。
而且不同的去噪方法往往会根据图像的不同类型以及去噪的不同目标来分析推断恢复图像的具体结构特征。
因此,在图像去噪领域没有一种适用于所有去噪问题的图像去噪技术。
在大多数情况下,要针对具体情况制定专门的去噪算法。
在图像处理中,数学形态学是一种对图像进行分析的数学工具,其分析的根底是图像的形态。
运用数学形态学处理图像时,要根据不同的图像类别单独的设计出整套的变换算法,当然这些算法也都是基于前文提到的膨胀、腐蚀、开运算、闭运算等根本变换。
通过这些根本变换而设计出的算法应该能有效地描述待处理图像各个元素与局部之间的关系。
通过数学形态学的应用,应该最终到达在简化图像数据、取出图像不相干结构的同时,保存住图像根本形状特征并且更加适合人眼视觉特征的目的。
图形态学处理图像的根本过程
15
噪声模型数字图像的噪声主要来源于图像的获取〔数字化过程〕和传输过程。
图像传感器的工作情况受各种因素的影响,如图像获取中的环境条件和传感元器件自身的质量。
例如,使用CCD摄像机获取图像,照顾程度和传感器温度是生成图像中产生大量噪声的主要因素。
图像在传输过程中主要由于所用的传输信道的干扰受到噪声污染。
比方通过无线网络传输的图像可能因为光或其其它大气因素的感染被污染。
高斯噪声在现实中这种噪声比拟普遍,所以我们对其考虑甚多。
事实上,这种易处理性比拟方便,考虑时这种模型经常居于临界情况下。
高斯随即变量Z的PDF由下式给出:
2
1
(z)2
2
p(z)=
2
〔〕
其中z表示灰度值,μ表示z的平均值或期望值,σ
表示z的标准差。
标
准差的平方σ2称为z的方差。
高斯函数的曲线如图
3-1〔a〕所示。
当z服从
式〔3-1〕的分布时候,其值有
70%落在[(μ-σ),(
μ+σ)]内,且有95%落在
[(μ-2σ),(
μ+2σ)]范围内。
椒盐噪声
椒盐噪声的PDF可由下式给出:
Pa
z
a
Pb
z
b
P(z)=
0
其他
(3.2)
如果b>a,灰度值b在图像中将显示为一个亮点,相反,a的值将显示为一
个暗点。
假设pa或pb为零,那么脉冲噪声称为单极脉冲。
如果
pa和pb均不可能为
零,尤其是它们近似相等时,脉冲噪声值将类似于随机分布在图像上的胡椒和盐粉微粒。
由于这个原因,双极脉冲声也称为椒盐噪声。
同时,它们有时也称为散粒和尖峰噪声。
16
形态学滤波器前文提到了四种数学形态学的根本变换:
膨胀、腐蚀、开运算以及闭运算。
其中开运算和闭运算经常被用来去除图像的噪声。
而这两种运算最典型的特征是它们都具有幂等性,因此作为滤波器,开运算和闭运算能够完全的完成滤波,没有必要进行重复操作。
而很多其他滤波器不具有这一性质。
例如中值滤波器会产生震荡,在有限次迭代之后会进一步的改变图像并不能保证得到一个稳定的结果。
开运算可以有效地滤除图像中的正噪声,即滤除目标中的噪声局部,特别是一些小的结构。
而闭运算可以被用来滤除负噪声,既是能填充背景中狭窄以及凹陷的局部,尤其是小的