数学形态学与数字图像处理

数学形态学与数字图像处理
数学形态学与数字图像处理

数学形态学与数字图像处理

目录

一绪论 (1)

1.1 数学形态学简介 (1)

1.2 数学形态学与数字图像处理 (1)

1.3 本次课程设计的目的与要求 (2)

二数学形态学的基本运算 (3)

2.1 基本概念 (3)

2.1.1结构元素 (3)

2.1.2膨胀与腐蚀 (3)

2.2 二值形态学图像处理 (4)

2.2.1 膨胀 (4)

2.2.2 腐蚀 (6)

2.2.3 开运算 (7)

2.2.4 闭运算 (8)

2.3 灰度形态学图像处理 (9)

2.3.1 膨胀 (9)

2.3.2 腐蚀 (10)

2.3.3 开运算与闭运算 (11)

2.4 综述 (13)

三数学形态学滤波器去噪 (15)

3.1 概述 (15)

3.2噪声模型 (16)

3.2.1 高斯噪声 (16)

3.2.2 椒盐噪声 (16)

3.3形态学滤波器 (17)

3.4形态学图像去噪原理 (20)

3.5形态学图像去噪的应用 (20)

小结与体会 (21)

参考文献 (22)

附录 (23)

一绪论

1.1数学形态学简介

数学形态学作为一门新兴的图像处理与分析学科,1964年由法国的G.Mathern和J.Serra在积分几何的基础上首次创立。70年代初,采用数学形态学的学者们开拓了图像分析的一个新的领域。经过十多年的理论与实践探索,G.Mathern和J.Serra等人在研究中认识到,对图像先作开运算接着再作闭运算,可以产生一种幂等运算;采用递增尺寸的交变开闭序列作用于图像,可有效地消除图像的噪声,1982年他们正式提出了形态学滤波器的概念。90年代数学形态学有两个显著的发展趋势,第一个是致力于运动分析,包括编码与运动景物描述;第二个是算法与硬件结构的协调发展,用于处理数值函数的形态学算子的开发与设计。

目前国内许多有效的图像处理系统有的是基于数学形态学方法原理设计的,有的是把数学形态学算法纳入其基本软件,并以其运算速度作为系统性能的重要标志之一

1.2数学形态学与数字图像处理

数学形态学在图像处理中属于非线性滤波方法,现在数学形态学的方法已经发展成为图象处理技术的一个重要方面,并且被广泛的应用到图象处理的各个领域,利用数学形态学可以进行图像去噪、图象分割、增强、边缘检测、形态分析、图象压缩等各个方面。可以通过以下几个步骤来实现数学形态学算法对数字图像的处理:

步骤 1、提取图像的几何结构特征,也就是针对所要处理的图像找出相应的几何结构模式。

步骤 2、根据步骤 1 找出的几何结构模式选合适的结构元素,这里结构元素的选择标准择首先是要能最有效的展现该几何结构模式,其次该结构元素的形态还应该尽量的最简。

步骤 3、为了得到比原始图像更能显著突出物体特征信息的图像,用步骤 2 选取的结构元素对目标进行相应的数学形态学变换,如果能对结构元素给予合

适的变量,则还能够定量的表示出目标的几何结构模式。

步骤 4、通过上面的三个步骤,相对于我们的处理需求,目标图像会变得更加清晰、明了,并且更有利于我们提取出相应的图像信息。

现在,数学形态学处理图像已经发展成为一个专门的图像科学领域。该领域已经形成了一个理论概念、非线性滤波、设计算法以及应用系统相互连贯而有广阔的整体。与其他很多图像处理技术相比,数学形态学技术的理论框架完善、算法效率高、易于在专门硬件上使用并且适合处理很多与形状相关的问题。例如对于图像噪声去由于可以在去除噪声前有效的探究目标图像的几何结构模式,尽可能的解决去除噪声与保护图像边缘细节信息相冲突的基本矛盾。再如在提取图像边缘时,与其他算法相比,数学形态学方法提取的边缘更为连续,间断点也会少很多。所以很多学术机构及工业研究所在处理数字图形图像、计算机视觉、模式识别等很多问题时都会重点考虑数学形态学方法。

1.3本次课程设计的目的与要求

(1)通过形态学方面的知识处理各种图像。

(2)学会应用形态学知识处理加有高斯噪声与椒盐噪声的图像。

(3)理解不同的形态学运算在处理图像方面的应用。

(4)通过运用MATLAB软件实现仿真。

二数学形态学的基本运算

2.1基本概念

数学形态学是由一组形态学的代数运算子组成的,它的基本运算有4个:膨胀(或扩张)、腐蚀(或侵蚀)、开运算和闭运算,它们在二值图像和灰度图像中各有特点。基于这些基本运算还可推导和组合成各种数学形态学实用算法,用它们可以进行图像形状和结构的分析及处理,包括图像分割、特征抽取、边界检测、图像滤波、图像增强和恢复等。数学形态学方法利用一个称作结构元素的“探针”收集图像的信息,当探针在图像中不断移动时,便可考察图像各个部分之间的相互关系,从而了解图像的结构特征。数学形态学基于探测的思想,与人的FOA(Focus Of Attention)的视觉特点有类似之处。作为探针的结构元素,可直接携带知识(形态、大小、甚至加入灰度和色度信息)来探测、研究图像的结构特点。

2.1.1结构元素

所谓结构元素就是一定尺寸的背景图像,通过将输入图像与之进行各种形态学运算,实现对输入图像的形态学变换。结构元素没有固定的形态和大小,它是在设计形态变换算法的同时根据输入图像和所需信息的形状特征一并设计出来的,结构元素形状、大小及与之相关的处理算法选择得恰当与否,将直接影响对输入图像的处理结果。通常结构元素的形状有正方形、矩形、圆盘形、菱形、球形以及线形等。

2.1.2膨胀与腐蚀

膨胀在数学形态学中的作用是把图像周围的背景点合并到物体中。如果两个物体之间距离比较近,那么膨胀运算可能会使这两个物体连通在一起,所以膨胀对填补图像分割后物体中的空洞很有用。腐蚀在数学形态学运算中的作用是消除物体边界点,它可以把小于结构元素的物体去除,选取不同大小的结构元素可以去掉不同大小的物体。如果两个物体之间有细小的连通,当结构元素足够大时,通过腐蚀运算可以将两个物体分开。

2.2二值形态学图像处理

二值图像数学形态学的运算就是基于上述集合论的理论,进行击中与否变换(HMT),在定义了 HMT 及其基本运算膨胀(Dilation)和腐蚀(Erosion)后,再从积分几何和体视学移植一些概念和理论,根据图像分析的各种要求,构造出统一的、相同的或变化很小的结构元素进行各种形态变换。数学形态学中有两种最基本的操作即膨胀和腐蚀,其他的所有形态学操作都是基于这两个操作的组合或级联。如开运算和闭运算就都是膨胀和腐蚀的最基本组合。膨胀、腐蚀、开运算和闭运算构成了整个数学形态学变换的基础,下面分别对这四种基本形态学变换进行具体的分析。

二值图像

腐蚀膨胀

图2.1 膨胀与腐蚀示意图

2.2.1膨胀

膨胀操作是指一个集或对象目标从其原来的形状扩大的过程。该目标扩大的方式是由结构元素决定的。和待处理的对象相比较,结构元素的大小更小,一般用于膨胀的结构元素大小取到 3×3。膨胀的过程类似于卷积,结构元素在目标图像内从左到右、从上到下的移动,在每次移动的过程之中,都会寻找结构元素与目标对象之间重叠的像素,只要存在重叠的像素点,结构元素所在的中心位置点的像素值都会被标为 1。用集合论该过程可表示如下:A,B 为Z2中的集合,Φ为空集, A 被 B 的膨胀,记为 A ⊕ B,⊕为膨胀算子,膨胀的集合定义式为:

(2.1)

该式表示的膨胀过程首先是B做关于原点的映射,然后平移 x。A被B 膨

胀也就是被所有x平移后于 A至少有一个公共非零元素的集合。根据上述对膨胀过程的解释,公式(2.2.1)也可以被写作下面的形式:

(2.2)

和其他的形态学操作一样,公式中集合 B 在膨胀运算中一般被叫做结构元素。

膨胀运算的实质是遍历待膨胀图像中的每个像素点,根据所选取的结构元素的值以及要处理像素点周围点的灰度值进行计算。比较局部范围内的像素点与结构元素中所对应点的灰度值之和。根据比较的结果,选取所计算的这些和中的最大值。所以经过膨胀,图像边缘的像素点灰度值会增加,图像边缘向外扩张,最终达到图像膨胀的视觉效果。不同的数学形态学文献对膨胀都有着不同的定义,公式(2.2.1)不是现在形态学文献中膨胀的唯一定义。然而,相比其它定义,这个定义存在一个明显的好处,即当把结构元素B 被当做卷积模板时,膨胀的概念会更加的形象化。因为虽然膨胀的本质是集合运算,而卷积本质上属于算术运算,但由于结构元素B 做相对于原点的“映射”后在集合 A(图像 A)上的平移是连续的,因此可以近似的将它滑过集合 A的整个过程近似看做卷积过程。

图2.2 膨胀操作

2.2.2腐蚀

腐蚀可以看做膨胀的逆运算或反过程。如果说膨胀是扩张了图像,那么腐蚀的作用则是使图像收缩。图像目标收缩的方式也是有结构元素决定的。和膨胀一样,腐蚀所选取的结构元素也要比目标要小,一般也是取 3×3 的大小。选取结构元素的尺寸较小的好处是可以减少腐蚀算法运行的时间。和膨胀相似,腐蚀操作也是将结构元素从左到右、从上到下在待处理图像中移动,以结构元素的中心点作为运算的中心,检验图像周围像素是否与结构元素完全重合。只要没有完全重叠,则该中心点像素就被标为 0。该过程同样可以用集合论的方法表示如下:

A ,

B 为Z2中的集合, A 被 B 腐蚀,记为 AΘB ,其定义公式为:

(2.3)

也就是说 A被B 的腐蚀的结果为所有使B 被 x平移后包含于 A的点 x的集合。和膨胀一样,腐蚀运算的实质也是遍历待腐蚀图像中的每个像素点,根据所选取的结构元素的值以及要处理像素点周围点的灰度值进行计算。比较局部范围内的像素点与结构元素中所对应点的灰度值之差。根据比较的结果,选取所计算的这些差中的最小值。经过腐蚀,图像边缘的像素点灰度值会降低,从而图像边缘会向内收缩,最终达到腐蚀的视觉效果。

膨胀和腐蚀的关系可以看做集合补和反转的对偶,可以用下面的公式表示:

(2.4)

图2.3 腐蚀操作

对一幅图像进行膨胀和腐蚀操作结果如下图:

原始图像阈值为0.8的图像腐蚀后图像1

膨胀后图像1腐蚀后图像2膨胀后图像2

图2.4

这两种运算具有对偶性,即一种运算对目标的操作相当于另一种运算对图像背景的操作。由图 2.2.4可知,膨胀操作后,图像相当于膨胀了一圈,而腐蚀操作后,图像则像被剥掉了一层。

2.2.3开运算

设 A是原始图像,B 是结构元素图像,则集合 A被结构元素B 做开运算,

记为:AoB,其公式为:

AoB=(AΘB)⊕B (2.5)

从公式(3.2.13)可以看出 A被结构元素B 做开运算就是 A被B 腐蚀后的结果再被 B 膨胀。同膨胀和腐蚀一样,我们也可以用用集合论的概念来定义开运算。A被结构元素B 做开运算就是B 在 A内的平移所得到的集合的并集,即:

(2.6)

开运算一般能平滑图像的轮廓,削弱狭窄的部分,去掉细的突出。根据开运算的这种作用,我们可以利用开运算来去除图像的噪声。

图2.5 开运算示意图

2.2.4闭运算

设 A是原始图像,B 是结构元素图像,则集合 A被结构元素B 做闭运算,记为:A?B,可用下面公式给出闭运算的定义:

A?B=(A⊕B)ΘB (2.7)

从公式(3.2.15)可以看出 A被结构元素B 做闭运算就是 A被B 膨胀后的结果再被 B 腐蚀。闭运算的效果基本上与开运算相反,它一般是能填充背景中狭窄以及凹陷的部分,消除小洞,还可以填充轮廓上空白的部分,使图像的轮廓得到一定的平滑。我们同样可以利用闭运算的这些作用来去除图像的噪声。

和膨胀和腐蚀一样,开运算和闭运算也是关于集合补和反转的对偶,可用下面公式表示这种关系:

(2.8)

对一幅二值图像进行开运算与闭运算的结果如下图:

原始图像阈值为0.7的二值图像开运算后图像1

闭运算后图像1开运算后图像2闭运算后图像2

图2.6

由图2.6可知,开运算能够去除孤立的小点、毛刺和小桥,而闭运算能够填平小孔,弥合小裂缝。

2.3灰度形态学图像处理

前面所提到的形态学方法都是基于二值图像的,下面把这些基本的形态学操作推广到灰度图像。这里把设为输入图像,为结构元素。

2.3.1膨胀

用b对函数f进行的灰度膨胀表示为 f⊕b,定义式为:

(2.9)其中,分别是 f 和b 的定义域,这里的 f 和b 表示的是函数而不再是二值形态学中所说的集合。

在灰度膨胀中,每个结构元素的位置上,在这一点的膨胀值是在跨度 b 的区间内 f 于b 之和的最大值。为通常对灰度图像进行膨胀处理会把图像中白色

的部分扩张,而且会带来双重的效果:(1)若使用的所有结构元素的值均为正,则输出图像会的亮度将会大于输入图像的亮度;(2) 膨胀过程中所用的结构元素的值和形状决定图像中暗的细节是部分全部减少了还是被消除。

2.3.2腐蚀

同样的道理,用b 对函数 f 进行的灰度腐蚀表示为 fΘb ,定义式为:

(2.10)和灰度膨胀一样,这里,分别是 f 和b 的定义域。在灰度腐蚀中,每个结构元素的位置上,在这一点的腐蚀值是在跨度 b 的区间内 f 于 b 之差的最小值。为通常对灰度图像进行腐蚀处理会把图像中黑色的部分扩张,而且会带来双重的效果:(1)若使用的所有结构元素的值均为正,则输出图像会的亮度将会小于输入图像的亮度;(2)在输入图像中如果结构元素的面积比亮的细节的面积大,则亮的效果将被削弱。环绕于亮细节周围的灰度值和结构元素自身的形态结构和幅值决定着最终消弱的程度。

灰度膨胀和腐蚀与二值膨胀和腐蚀一样也满足对偶关系。

对一幅灰度图像进行膨胀与腐蚀的结果如下图所示:

原始图像腐蚀后图像1膨胀后图像1

腐蚀后图像2膨胀后图像2

图2.7

由图2.7可知,灰值图像的腐蚀,如果结构元素的值都为正,则输出图像会比输入图像暗;如果输入图像中亮细节的尺寸比结构元素小,则其影响会被减弱,减弱的程度取决于这些亮细节周围的灰度值的结构元素的形状和幅值。灰度图像的膨胀运算,如果结构元素的值都为正,则输出图像会比输入图像亮,根据输入图像中暗细节的灰度值以及它们的形状相对于结构元素的关系,它们在膨胀中或被消减或被除掉。

2.3.3开运算与闭运算

用结构元素b 对图像 f 进行开操作即先用b 对 f 进行腐蚀操作再进行膨胀,记为fob,定义式为:

(2.11)

用结构元素b 对图像 f 进行闭操作即先用b 对 f 进行膨胀操作再进行腐蚀,记为f ? b,定义式为:

f?b = (f⊕b)Θb (2.12)

用几何的方法可以解释灰度形态学的开、闭运算。用b 对 f 进行开操作的原理可以解释为:顺着曲面的下侧面推动小球上下滚动,使球体能够在曲面的整个下侧面来回移动。这时球体表面的所有接触到曲面部分的最高点就组成了开操作的曲面。相对应的,闭操作就是在曲面的上侧面滚动,以便球体能在曲面的整个上侧面来回移动。当球体的任何部分接触到曲面的最高点就构成了闭操作的曲面。如下图2.8就详细的解释了这一过程:

图2.8

灰度形态学的开、闭运算和二值形态学一样都满足对偶关系。并且灰度开、

闭运算的作用效果也基本上于二值形态学中的一样。其中灰度开运算一般能平

滑图像的轮廓,消减图像中狭窄的部分,抹平突出的细节。灰度闭运算的效果

基本上与灰度开运算相反,它一般是能填充背景中狭窄以及凹陷的部分,消除

小洞,还可以填充轮廓上空白的部分,使图像的轮廓得到一定的平滑。我们同

样可以利用灰度开、闭运算的这些作用来去除图像中不希望看到的掺杂细节(如

噪声)

对一幅加有椒盐噪声的灰度图像进行开运算与闭运算的结果如下:

原始图

像加入椒盐噪

声开运算后图像

1

闭运算后图像

1开运算后图像

2闭运算后图像2

图2.9 开运算操作消除与结构元素相比尺寸较小的亮细节,而保持图像整体灰度

值和大的亮区域基本不受影响。闭运算操作消除与结构元素相比尺寸较小的暗

细节,而保持图像整体灰度值和大的暗区域基本不受影响。所以,开-闭运算可

用于图像的去噪,特别是对于加有椒盐噪声的图像的去噪效果更为明显。详细

过程见下一章。

对一幅加有高斯噪声的图像进行开运算与闭运算的结果如下:

原始图

像加入高斯噪

声开运算后图像

1

闭运算后图像

1开运算后图像

2闭运算后图像2

图3.0 单独运用开运算或闭运算对一幅加有高斯噪声的图像进行去噪效果并不理

想。

2.4综述

本章首先对数学形态学的大致发展历程及基本原理进行了概括性的介绍。

说明了数学形态学在图像处理中的应用领域、处理步骤及其在图像处理中的独

有优势。数学形态学是用数学集合论中的语言来描述的,也即集合论是数学形

态学的基础,所以本章对集合论的主要内容做了介绍。然后基于集合论的基本

原理,本章重点讨论了形态学中最基本的部分——二值形态学。用集合论的基

本理论给出了二值形态学的四种基本运算——膨胀、腐蚀、开运算和闭运算的

公式;说明了这四种基本二值形态学变换在处理图像时的运算原理;通过实验

(主要是去噪实验)给出了四种形态学变换对图像的处理效果并对其处理噪声

的效果进行了讨论。分别对膨胀和腐蚀、开运算和闭运算的关系以集合论公式

的形式给予了介绍。然后本章把二值形态学的基本理论推广到灰度图像,以交叉对比的形式把灰度形态学中膨胀、腐蚀、开运算及闭运算的基本集合论公式、原理等进行了介绍。并同样通过实验说明了膨胀和腐蚀对图像的处理效果以及开运算和闭运算在去除图像噪声时所表现出的特点。通过本章的介绍,可以明确数学形态学的基本理论;对数学形态学在图像处理,特别是去除噪声时的优势有了进一步的说明。

本章是下一章关于形态学去噪的理论基础,特别是其中的开-闭运算,在图像去噪方面有着十分重要的应用。

三数学形态学滤波器去噪

3.1概述

通过图像去噪(也可称之为图像清洁)来对图像进行增强是图像处理中的一个基本问题。于这个问题,研究者们已经进行了长达多年的研究。大多数的图像处理书籍中都会有涉及到图像去噪的章节。从本质上讲,图像去噪属于图像恢复问题中的一类,目的是将一个受到噪声侵蚀的图像版本恢复到其原先完善的版本。这个问题看似很简单,但是由于我们常常无法对恢复后图像的实际结构作出合理的分析,而会变得异常复杂。而且不同的去噪方法往往会根据图像的不同类型以及去噪的不同目标来分析推断恢复图像的具体结构特征。因此,在图像去噪领域没有一种适用于所有去噪问题的图像去噪技术。在大多数情况下,要针对具体情况制定专门的去噪算法。在图像处理中,数学形态学是一种对图像进行分析的数学工具,其分析的基础是图像的形态。运用数学形态学处理图像时,要根据不同的图像类别单独的设计出整套的变换算法,当然这些算法也都是基于前文提到的膨胀、腐蚀、开运算、闭运算等基本变换。通过这些基本变换而设计出的算法应该能有效地描述待处理图像各个元素与部分之间的关系。通过数学形态学的应用,应该最终达到在简化图像数据、取出图像不相干结构的同时,保留住图像基本形状特征并且更加适合人眼视觉特征的目的。

图3.1 形态学处理图像的基本过程

3.2噪声模型

数字图像的噪声主要来源于图像的获取(数字化过程)和传输过程。图像

传感器的工作情况受各种因素的影响,如图像获取中的环境条件和传感元器件

自身的质量。例如,使用CCD 摄像机获取图像,关照程度和传感器温度是生成

图像中产生大量噪声的主要因素。图像在传输过程中主要由于所用的传输信道

的干扰受到噪声污染。比如通过无线网络传输的图像可能因为光或其其它大气

因素的感染被污染。

3.2.1高斯噪声

在现实中这种噪声比较普遍,所以我们对其考虑甚多。事实上,这种易处

理性比较方便,考虑时这种模型经常居于临界情况下。

高斯随即变量Z 的PDF 由下式给出: p(z)=2

22

)(21σμπσ--z (3.1)

其中z 表示灰度值,μ表示z 的平均值或期望值,σ表示z 的标准差。标

准差的平方σ2称为z 的方差。高斯函数的曲线如图3-1(a )所示。当z 服从

式(3-1)的分布时候,其值有70%落在[(μ-σ),(μ+σ)]内,且有95%落在

[(μ-2σ),( μ+2σ)]范围内。

3.2.2椒盐噪声

椒盐噪声的PDF 可由下式给出:

P(z)=?????==其他0b z P a z P b a (3.2)

如果b>a ,灰度值b 在图像中将显示为一个亮点,相反,a 的值将显示为一

个暗点。若a p 或b p 为零,则脉冲噪声称为单极脉冲。如果a p 和b p 均不可能为

零,尤其是它们近似相等时,脉冲噪声值将类似于随机分布在图像上的胡椒和

盐粉微粒。由于这个原因,双极脉冲声也称为椒盐噪声。同时,它们有时也称为散粒和尖峰噪声。

3.3形态学滤波器

前文提到了四种数学形态学的基本变换:膨胀、腐蚀、开运算以及闭运算。其中开运算和闭运算经常被用来去除图像的噪声。而这两种运算最典型的特征是它们都具有幂等性,因此作为滤波器,开运算和闭运算能够完全的完成滤波,没有必要进行重复操作。而很多其他滤波器不具有这一性质。例如中值滤波器会产生震荡,在有限次迭代之后会进一步的改变图像并不能保证得到一个稳定的结果。

开运算可以有效地滤除图像中的正噪声,即滤除目标中的噪声部分,特别是一些小的结构。而闭运算可以被用来滤除负噪声,既是能填充背景中狭窄以及凹陷的部分,尤其是小的孔洞。通过对开运算和闭运算的组合运用,可以得到四种新的形态学滤波器:开运算后闭运算、闭运算后开运算、开运算后闭运算再开运算、闭运算后开运算再闭运算。这四种新的形态学变换具有幂等性,因此在进行其他的组合并不会得到新的形态学滤波。这四种变换可以作为滤波器来去除正噪声以及负噪声,例如可以使用它们来替代中值滤波法来滤除图像中的斑点噪声。而且在去除图像噪声方面,上面的四种新的形态学组合运算一般会比单纯使用形态学开运算和闭运算会取得更好的去噪峰值信噪比以及视觉效果,即具有更佳的去噪效果。

形态滤波器是用一个结构元素B对初始图像串联地使用开、闭操作。这样图像中比结构元素小的游离的噪声将被滤除。若初始图像为A,结构元素为B,则形态滤波器可以这样来构成:

OC(A,B)=C(O(A,B),B)或CO(A,B)=O(C(A,B),B) (3.3) 形态滤波器的详尽描述如下:

(((AΘB)⊕B)⊕BΘB或(((A⊕B)ΘB)ΘB)⊕B (3.4) 如果结构元素包括原点(0,0),则腐蚀和膨胀满足以下性质:

性质1: AΘB?A?A⊕B (3.5) 这一性质表明,在B包括原点的前提下,腐蚀后的结果只会使A的点数减

少或者不变,而膨胀则使A的点数增加或者不变。利用前一点,可以通过设计

适当的结构元素B,使得腐蚀后得以消除A中的微小颗粒,即噪声点。利用后一

点,又可以对腐蚀结果再用B进行膨胀,以恢复有用信息(细节部分)。

性质2:对开运算和闭运算,恒有

O(A,B)?A?C(A,B) (3.6)

即开运算使原图形缩小而闭运算使原图形增大。

根据上面的讨论以及开闭运算的性质不难证明形态开一闭(OC)和形态闭一

开(CO)滤波器具有如下一些重要性质:

(l)平移不变性:

OC(A+x,B)=OC(A,B)+x

CO(A+x,B)=CO(A,B)+x (3.7)

(2)递增性:

如果1A是2A的子集,则:

OC(A1,B)?OC(A2,B)

CO(A1,B)?CO(A2,B) (3.8)

(3)幂等性:

OC(CO(A,B))=CO(OC(A,B))

CO(OC(A,B))=OC(CO(A,B)) (3.9) (4)对偶性:

(OC(A,B))c=CO(A,B)

(CO(A,B))c=OC(A,B) (3.10) 形态滤波器的输出不仅取决于变换的形式,而且取决于结构元素的尺寸和形

状,一般只有与结构元素的尺寸和形状相匹配的基元才能被保留。

如我们在同一个图像中都加入 6%椒盐噪声,分别使用开运算、闭运算、

闭运算后开运算、开运算后闭运算四种形态学变换进行去噪处理,如下图:

数学形态学的基本运算

第二章数学形态学的基本运算 2.1二值腐蚀和膨胀 二值图象是指那些灰度只取两个可能值的图象,这两个灰度值通常取为0和1。习惯上认为取值1的点对应于景物中的点,取值为0的点构成背景。这类图象的集合表示是直接的。考虑所有1值点的集合(即物体)X,则X与图象是一一对应的。我们感兴趣的也恰恰是X集合的性质。 如何对集合X进行分析呢?数学形态学认为,所谓分析,即是对集合进行变换以突出所需要的信息。其采用的是主观“探针”与客观物体相互作用的方法。“探针”也是一个集合,它由我们根据分析的目的来确定。术语上,这个“探针”称为结构元素。选取的结构元素大小及形状不同都会影响图象处理的结果。剩下的问题就是如何选取适当的结构元素以及如何利用结构元素对物体集合进行变换。为此,数学形态学定义了两个最基本的运算,称为腐蚀和膨胀即1。 2.1 .1二值腐蚀运算 腐蚀是表示用某种“探针”(即某种形状的基元或结构元素)对一个图象进行探测,以便找出图象内部可以放下该基元的区域。它是一种消除边界点,使边界向内部收缩的过程。可以用来消除小且无意义的物体。腐蚀的实现同样是基于填充结构元素的概念。利用结构元素填充的过程,取决于一个基本的欧氏空间概念—平移。我们用记号A二表示一个集合A沿矢量x平移了一段距离。即: 集合A被B腐蚀,表示为AΘB,其定义为: 其中A称为输入图象,B称为结构元素。AΘB由将B平移x仍包含在A内的所有点x组成。如果将B看作模板,那么,AΘB则由在将模板平移的过程中,所有可以填入A内部的模板的原点组成。根据原点与结构元素的位置关系,腐蚀后的图象大概可以分为两类: (1)如果原点在结构元素的内部,则腐蚀后的图象为输入图象的子集,如图2.1所示。 (2)如果原点在结构元素的外部,那么,腐蚀后的图象则可能不在输入图象的内部,如图2.2所示。 图2.1腐蚀类似于收缩

数字图像处理实验五

数字图像处理 实验 实验五:图像增强-空域滤波 学院:信息工程学院 姓名: 学号: 专业及班级: 指导教师:

一、 实验目的 进一步了解MatLab 软件/语言,学会使用MatLab 对图像作滤波处理,使学生有机会掌握滤波算法,体会滤波效果。 了解几种不同滤波方式的使用和使用的场合,培养处理实际图像的能力,并为课堂教学提供配套的实践机会。 二、 实验内容 (1)学生应当完成对于给定图像+噪声,使用平均滤波器、中值滤波器对不同强度的高斯噪声和椒盐噪声,进行滤波处理;能够正确地评价处理的结果;能够从理论上作出合理的解释。 (2)利用MATLAB 软件实现空域滤波的程序: I=imread('electric.tif'); J = imnoise(I,'gauss',0.02); %添加高斯噪声 J = imnoise(I,'salt & pepper',0.02); %添加椒盐噪声 ave1=fspecial('average',3); %产生3×3的均值模版 ave2=fspecial('average',5); %产生5×5的均值模版 K = filter2(ave1,J)/255; %均值滤波3×3 L = filter2(ave2,J)/255; %均值滤波5×5 M = medfilt2(J,[3 3]); %中值滤波3×3模板 N = medfilt2(J,[4 4]); %中值滤波4×4模板 imshow(I); figure,imshow(J); figure,imshow(K); figure,imshow(L); figure,imshow(M); figure,imshow(N); 三、实验具体实现 a) 调入并显示原始图像Sample2-1.jpg 。 b) 利用imnoise 命令在图像Sample2-1.jpg 上加入高斯(gaussian) 噪声 c)利用预定义函数fspecial 命令产生平均(average)滤波器 111191111---????--????---? ? d )分别采用3x3和5x5的模板,分别用平均滤波器以及中值滤波器,对加入噪声的图像进行处理并观察不同噪声水平下,上述滤波器处理的结果; e )选择不同大小的模板,对加入某一固定噪声水平噪声的图像进行处理,观察上述滤波器处理的结果。 f )利用imnoise 命令在图像Sample2-1.jp g 上加入椒盐噪声(salt & pepper)

数字图像处理教学大纲(2014新版)

数字图像处理 课程编码:3073009223 课程名称:数字图像处理 总学分: 2 总学时:32 (讲课28,实验4) 课程英文名称:Digital Image Processing 先修课程:概率论与数理统计、线性代数、C++程序设计 适用专业:自动化专业等 一、课程性质、地位和任务 数字图像处理课程是自动化专业的专业选修课。本课程着重于培养学生解决智能化检测与控制中应用问题的初步能力,为在计算机视觉、模式识别等领域从事研究与开发打下坚实的理论基础。主要任务是学习数字图像处理的基本概念、基本原理、实现方法和实用技术,并能应用这些基本方法开发数字图像处理系统,为学习图像处理新方法奠定理论基础。 二、教学目标及要求 1.了解图像处理的概念及图像处理系统组成。 2.掌握数字图像处理中的灰度变换和空间滤波的各种方法。 3.了解图像变换,主要是离散和快速傅里叶变换等的原理及性质。 4.理解图像复原与重建技术中空间域和频域滤波的各种方法。 5. 理解解彩色图像的基础概念、模型和处理方法。 6. 了解形态学图像处理技术。 7. 了解图像分割的基本概念和方法。 三、教学内容及安排 第一章:绪论(2学时) 教学目标:了解数字图像处理的基本概念,发展历史,应用领域和研究内容。通过大量的实例讲解数字图像处理的应用领域;了解数字图像处理的基本步骤;了解图像处理系统的组成。 重点难点:数字图像处理基本步骤和图像处理系统的各组成部分构成。 1.1 什么是数字图像处理 1.2 数字图像处理的起源

1.3.1 伽马射线成像 1.3.2 X射线成像 1.3.3 紫外波段成像 1.3.4 可见光及红外波段成像 1.3.5 微波波段成像 1.3.6 无线电波成像 1.3.7 使用其他成像方式的例子 1.4 数字图像处理的基本步骤 1.5 图像处理系统的组成 第二章:数字图像基础(4学时) 教学目标:了解视觉感知要素;了解几种常用的图像获取方法;掌握图像的数字化过程及其图像分辨率之间的关系;掌握像素间的联系的概念;了解数字图像处理中的常用数学工具。 重点难点:要求重点掌握图像数字化过程及图像中像素的联系。 2.1 视觉感知要素(1学时) 2.1.1 人眼的构造 2.1.2 眼镜中图像的形成 2.1.3 亮度适应和辨别 2.2 光和电磁波谱 2.3 图像感知和获取(1学时) 2.3.1 用单个传感器获取图像 2.3.2 用条带传感器获取图像 2.3.3 用传感器阵列获取图像 2.3.4 简单的图像形成模型 2.4 图像取样和量化(1学时) 2.4.1 取样和量化的基本概念 2.4.2 数字图像表示 2.4.3 空间和灰度级分辨率 2.4.4 图像内插 2.5 像素间的一些基本关系(1学时) 2.5.1 相邻像素 2.5.2 临接性、连通性、区域和边界 2.5.3 距离度量 2.6 数字图像处理中所用数学工具的介绍 2.6.1 阵列与矩阵操作

基于数学形态学的图像噪声处理.

基于数学形态学的图像噪声处理 摘要 本文首先介绍了数学形态学的发展简史及其现状,紧接着详细的阐述了数学形态学在图像处理和分析中的理论基础。并从二值数学形态 学出发着重研究了数学形态学的膨胀、腐蚀、开运算、闭运算等各种 运算和性质,然后根据已有的运算,接着引入了形态滤波器设计、形态学图像处理的实用算法。由于在图像的获取中存在各种可能的噪声,比 如高斯噪声、瑞利噪声、伽马噪声、指数噪声、均匀噪声以及椒盐等 噪声,由于这些噪声的普遍存在,因此,利用数学形态学的腐蚀、膨胀、开启、闭合设计出了一种比较理想的(闭和开)形态学滤波器,并且用MATLAB语言编写程序,反复的使用这种开闭、闭开来处理图像中存在的噪声,其效果比较满意。 关键词:数学形态学图像处理腐蚀膨胀滤波Studies on Mathematical Morphology for Image Processing ABSTRACT In this paper ,we first introduced the brief history and development of mathematical morphology some general theory of mathematical morphology analysis and many experiment results are https://www.360docs.net/doc/e02644189.html,ter ,from the aspect of morphology of dual value, special emphasis on various operations and properties including dilation, erosion,open operation and close operation etc.In addition, morphology analysis method of the dual value image is also discussed and the practical and improved operations of the morphological image processing such as electric filter design, marginal pattern testing are introduced. As the image of the acquisition in the range of possible noise, such as Gaussian noise, Rayleigh noise, Gamma noise, Uniform noise Salt and Pepper noise and so on. As the prevalence of such noise, so using mathematical morphology of erosion,dilation, opening, closing designed a more ideal (open and closed morphological filter, And repeated to use opening and closing, closing and opening handle image processing in the noise. It is satisfied with its results.And the simulation results is more satisfactory after the use of MATLAB language programming. Keyword:mathematical morphology image processing erosion dilation

数字图像处理实验5 冈萨雷斯

一、 实验目的: 实验五 彩色图像处理 1、使用 RGB 、Indexed 和 Gray 图像间转换函数 2、掌握彩色图像平滑与锐化的算法 3、彩色图像的分割 二、 实验内容 1、使用 RGB 、Indexed 和 Gray 图像间转换函数 对图像“Fig0630(01)(strawberries_fullcolor).tif ”使用 rgb2ind 分别产生 8 色抖 动和非抖动图像;使用 rgb2gray 实现图像转换,并使用函数 dither 产生其抖 动形式的图像。分别显示如下: 2、彩色图像平滑与锐化 对图像“ F ig0604(a)(iris).tif ” 在 RGB 空 间 实 现 彩色图像的平滑 ( w =ones (25)./(25*25))与锐化(w= [1 1 1 1 1; 1 1 1 1 1; 1 1 -24 1 1 ; 1 1 1 1 1; 1 1 1 1 1];)(函数 imfilter );输出结果分别如下:

在 HSI 空间实现彩色图像的平滑(w=ones (25)./(25*25)),(函数 imfilter, rgb2hsi, cat ),观察仅平滑亮度分量(intensity )和平滑全部三个分量结果的差 别。 3、彩色图像分割 用 colorseg 函 数 实 现 基 于 肤 色 的 人 脸 分 割 。 对 图 像 “'Fig0636(woman_baby_original).tif ”在 HSI 空间,用函数 colorseg 分别用参 数“euclidean ”和“mahalanobis ”实现人脸分割。结果类似如下:

数字图像处理实验报告5

数字图像处理与分析 实 验 报 告 学院: 班级: 姓名: 学号:

实验五(1)图像压缩 一、实验目的 1. 理解图像压缩的基本定义和常见方法 2. 掌握在MATLAB中进行图像压缩的方法 3. 掌握利用DCT进行图像压缩的方法 4. 进一步熟悉了解MATLAB语言的应用,为进行综合性图像处理实验打下基础 二、实验内容 %利用离散余弦变换进行JPEG图像压缩 I=imread('1.bmp'); %读入原图像; I=im2double(I); %将原图像转为双精度数据类型; T=dctmtx(8); %产生二维DCT变换矩阵 B=blkproc(I,[8 8],'P1*x*P2',T,T'); %计算二维DCT,矩阵T及其转置T’是DCT 函数P1*x*P2的参数 Mask=[ 1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; %二值掩膜,用来压缩DCT系数,只留下数中 左上角的10个 B2=blkproc(B,[8 8],'P1.*x',Mask); %只保留DCT变换的10个系数 I2= blkproc(B2,[8,8],'P1*x*P2',T',T); %逆DCT,重构图像 Subplot(1,2,1); Imshow(I);title('原图像'); %显示原图像 Subplot(1,2,2); Imshow(I2);title('压缩图像'); %显示压缩后的图像。虽然舍弃了85%的DCT 系数,但图像仍然清晰(有一些质量损失)得到图像:

数字图像处理应用论文数字图像处理技术论文

数字图像处理应用论文数字图像处理技术论文 关于数字图像处理及其应用的研究 摘要:首先对数字图像处理的关键技术以及相应的处理设备进行详细的探讨,然后对数字图像处理的应用领域以及发展趋势进行详尽论述。 关键词:数字图像处理:关键技术;应用领域 0 引言 人类通过眼、耳、鼻、舌、身接受信息,感知世界。约有75%的信息是通过视觉系统获取的。数字图象处理是用数字计算机处理所获取视觉信息的技术,上世纪20年代Bartlane电缆图片传输系统(纽约和伦敦之间海底电缆)传输一幅图片所需的时间由一周多减少到小于3个小时;上世纪50年代,计算机的发展,数字图像处理才真正地引起人们的巨大兴趣;1964年,数字图像处理有效地应用于美国喷气推进实验室(J.P.L)对“徘徊者七号”太空船发回的大批月球照片的处理;但是直到上世纪六十年代末至七十年代扔,由于离散数学理论的创立和完善,使之形成了比较完整的理论体系,成为一门新兴的学科。数字图像处理的两个主要任务:如何利用计算机来改进图像的品质以便于人类视觉分析;对图像数据进行存储、传输和表示,便于计算机自动化处理。图像处理的范畴是一个受争论的话题,因此也产生了其他的领域比如图像分析和计算机视觉等等。

1 数字图像处理主要技术概述 不论图像处理是基于什么样的目的,一般都需要通过利用计算机图像处理对输入的图像数据进行相关的处理,如加工以及输出,所以关于数字图像处理的研究,其主要内容可以分为以下几个过程。图像获取:这个过程基本上就是把模拟图像通过转换转变为计算机真正可以接受的数字图像,同时,将数字图像显示并且体现出来(例如彩色打印)。数据压缩和转换技术:通过数据压缩和数据转换技术的研究,减少数据载体空间,节省运算时间,实现不同星系遥感数据应用的一体化。图像分割:虽然国内外学者已提出很多种图像分割算法,但由于背景的多变性和复杂性,至今为止还没有一种能适用于各种背景的图像分割算法。当前提出的小波分析、模糊集、分形等新的智能信息处理方法有可能找到新的图像分割方法。图像校正:在理想情况下,卫星图像上的像素值只依赖于进入传感器的辐射强度;而辐射强度又只与太阳照射到地面的辐射强度和地物的辐射特性(反射率和发射率)有关,使图像上灰度值的差异直接反映了地物目标光谱辐射特性的差异,从而区分地物目标。图像复原,以图像退化的数学模型为基础,来改善图像质量表达与描述,图像分割后,输出分割标记或目标特征参数;特征提取:计算描述目标的特征,如目标的几何形状特征、统计特征、矩特征、纹理特征等。图像增强:显示图像中被模糊的细节。或是突出图像中感兴趣的特征。图像识别:统计模式识别、模糊模式识别、人工神经网络等。

数字图像处理复习重点整理

《数字图像处理》复习 第一章绪论 数字图像处理技术的基本内容:图像变换、图像增强、图象恢复、图像压缩编码、图像分割、图像特征提取(图像获取、表示与描述)、彩色图像处理和多光谱及高光谱图像处理、形态学图像处理 第二章数字图像处理基础 2-1 电磁波谱与可见光 1.电磁波射波的成像方法及其应用领域: 无线电波(1m-10km)可以产生磁共振成像,在医学诊断中可以产生病人身体的横截面图像☆微波(1mm-1m)用于雷达成像,在军事和电子侦察领域十分重要 红外线(700nm-1mm)具有全天候的特点,不受天气和白天晚上的影响,在遥感、军事情报侦察和精确制导中广泛应用 可见光(400nm-700nm)最便于人理解和应用最广泛的成像方式,卫星遥感、航空摄影、天气观测和预报等国民经济领域 ☆紫外线(10nm-400nm)具有显微镜方法成像等多种成像方式,在印刷技术、工业检测、激光、生物学图像及天文观测 X射线(1nm-10nm)应用于获取病人胸部图像和血管造影照片等医学诊断、电路板缺陷检测等工业应用和天文学星系成像等 伽马射线(0.001nm-1nm)主要应用于天文观测 2-2 人眼的亮度视觉特征 2.亮度分辨力——韦伯比△I/I(I—光强△I—光照增量),韦伯比小意味着亮度值发生较小变化就能被人眼分辨出来,也就是说较小的韦伯比代表了较好的亮度分辨力 2-3 图像的表示 3. 黑白图像:是指图像的每个像素只能是黑或白,没有中间的过渡,一般又称为二值图像 (黑白图像一定是二值图像,二值图像不一定是黑白图像) 灰度图像:是指图像中每个像素的信息是一个量化了的灰度级的值,没有彩色信息。 彩色图像:彩色图像一般是指每个像素的信息由R、G、B三原色构成的图像,其中的R、B、G是由不同的灰度级来描述的。 4.灰度级L、位深度k L=2^k 5.储存一幅M×N的数字图像所需的比特 b=M×N×k 例如,对于一幅600×800的256灰度级图像,就需要480KB的储存空间(1KB=1024Byte 1Byte=8bit) 2-4 空间分辨率和灰度级分辨率 6.空间分辨率是图像中可分辨的最小细节,主要由采样间隔值决定,反映了数字化后图像的实际分辨率。一种常用的空间分辨率的定义是单位距离内可分辨的最少黑白线对数目(单位是每毫米线对数),比如每毫米80线对。对于一个同样大小的景物来说,对其进行采样的空间分辨率越高,采样间隔就越小,图片的质量就越高。 7.灰度级分辨率是指在灰度级别中可分辨的最小变化,通常把灰度级级数L称为图像的灰度级分辨率(灰度级通常是2的整数次幂) 8.在图像空间分辨率不变的情况下,采样数越少,图像越小。同时也证实了,在景物大小不变的情况下,图像阵列M×N越小,图像的尺寸就越小; 随着空间分辨率的降低,图像大小尺寸不变,图像中的细节信息在逐渐损失,棋盘格似的粗颗粒像素点变得越来越明显。由此也说明,图像的空间分辨率越低,图像的视觉效果越差;随着灰度分辨率的降低,图像的细节信息在逐渐损失,伪轮廓信息在逐渐增加。由于伪轮

形态学图像处理小结

一.形态学基础知识理解 形态学图像处理基本的运算包括:二值腐蚀和膨胀、二值开闭运算、骨架抽取、极限 腐蚀、击中击不中变换、形态学梯度、Top-hat变换、颗粒分析、流域变换、灰值腐蚀和膨胀、灰值开闭运算、灰值形态学梯度等。 1.膨胀与腐蚀 最基本的形态学操作有二种:膨胀与腐蚀(Dilation与Erosion)。膨胀是在二值图像中“加长”和“变粗”的操作。这种方式和变粗的程度由一个结构元素组成的集合来控制。 腐蚀是“收缩”或“细化”二值图像中的对象。同样,收缩的方式和程度由一个结构元素 控制。腐蚀和膨胀是对白色部分(高亮部分)而言的,不是黑色部分。膨胀就是图像中的 高亮部分进行膨胀,“领域扩张”,效果图拥有比原图更大的高亮区域。腐蚀就是原图中 的高亮部分被腐蚀,“领域被蚕食”,效果图拥有比原图更小的高亮区域。 常用的三种膨胀与腐蚀的组合:开运算、闭运算、击中或击不中变换。(1)开运算 和闭运算: A被B的形态学开运算是A被B腐蚀后再用B来膨胀腐蚀结果。其几何解释为:B在A内完全匹配的平移的并集。形态学开运算完全删除了不能包含结构元素的对象区域,平滑了对象的轮廓,断开了狭窄的连接,去掉了细小的突出部分;(2)闭运算: A被B的形态学闭运算是先膨胀再腐蚀的结果,其几何解释为:所有不与A重叠的B的平移的并集。形态学闭运算会平滑对象的轮廓,与开运算不同的是,闭运算一般会将狭窄的缺口连接起 来形成细长的弯口,并填充比结构元素小的洞。(3)击中击不中变换: 击中与击不中变换先对目标图像进行目标结构元素的腐蚀操作;后对目标图像的对偶进行背景结构元素的腐 蚀操作;最后取两次结果的交集。 2.重构 重构是一种涉及到两幅图像和一个结构元素的形态学变换。一幅图像,即标记(marker),是变换的开始点。另一幅图像是掩模(mask),用来约束变换过程。结构 元素用于定义连接性。 3.灰度图像形态学

基于数学形态学的图像边缘检测方法研究文献综述

文献综述 课题:基于数学形态学的图像边缘检测方法研究 边缘检测是图像分割的核心容,而图像分割是由图像处理到图像分析的关键步骤,在图像工程中占据重要的位置,对图象的特征测量有重要的影响。图像分割及基于分割的目标表达、特征提取和参数测量等将原始图像转化为更抽象更紧凑的形式,使得更高层的图像分析和理解成为可能。从而边缘检测在图像工程中占有重要的地位和作用。因此对边缘检测的研究一直是图像技术研究中热点,人们对其的关注和研究也是日益深入。 首先,边缘在边界检测、图像分割、模式识别、机器视觉等中有很重要的作用。边缘是边界检测的重要基础,也是外形检测的基础。同时,边缘也广泛存在于物体与背景之间、物体与物体之间,基元与基元之间,是图像分割所依赖的重要特征。其次,边缘检测对于物体的识别也是很重要的。第一,人眼通过追踪未知物体的轮廓而扫视一个未知的物体。第二,如果我们能成功地得到图像的边缘,那么图像分析就会大大简化,图像识别就会容易得多。第三,很多图像并没有具体的物体,对这些图像的理解取决于它们的纹理性质,而提取这些纹理性质与边缘检测有极其密切的关系。 理想的边缘检测是能够正确解决边缘的有无、真假、和定向定位。长期以来,人们一直关心这一问题的研究,除了常用的局

部算子及以后在此基础上发展起来的种种改进方法外,又提出了许多新的技术,其中,比较经典的边缘检测算子有 Roberts cross算子、Sobel算子、Laplacian算子、Canny算子等,近年来又有学者提出了广义模糊算子,形态学边缘算子等。这些边缘检测的方法各有其特点,但同时也都存在着各自的局限性和不足之处。 本次研究正是在已有的算法基础上初步进行改进特别是形 态学边缘算子,以期找到一个更加简单而又实用的算子,相信能对图像处理中的边缘检测方法研究以及应用有一定的参考价值。 一、课题背景和研究意义: 伴随着计算机技术的高速发展,数字图像处理成为了一门新兴学科,并且在生活中的各个领域得以广泛应用。图像边缘检测技术则是数字图像处理和计算机视觉等领域最重要的技术之一。在实际图像处理中,图像边缘作为图像的一种基本特征,经常被用到较高层次的图像处理中去。边缘检测技术是图像测量、图像分割、图像压缩以及模式识别等图像处理技术的基础,是数字图像处理重要的研究课题之一。 边缘检测是图像理解、分析和识别领域中的一个基础又重要的课题, 边缘是图像中重要的特征之一,是计算机视觉、模式识别等研究领域的重要基础。图像的大部分主要信息都存在于图像的边缘中,主要表现为图像局部特征的不连续性,是图像中灰度变化比较强烈的地方,也即通常所说的信号发生奇异变化的地

数字图像处理实验五复习课程

数字图像处理实验五

数字图像处理 实验 实验五:图像增强-空域滤波 学院:信息工程学院 姓名: 学号: 专业及班级: 指导教师:

一、实验目的 进一步了解MatLab软件/语言,学会使用MatLab对图像作滤波处理,使学生有机会掌握滤波算法,体会滤波效果。 了解几种不同滤波方式的使用和使用的场合,培养处理实际图像的能 力,并为课堂教学提供配套的实践机会。 二、实验内容 (1)学生应当完成对于给定图像+噪声,使用平均滤波器、中值滤波器对不同强度的高斯噪声和椒盐噪声,进行滤波处理;能够正确地评价处理的结果;能够从理论上作出合理的解释。 (2)利用MATLAB软件实现空域滤波的程序: I=imread('electric.tif'); J = imnoise(I,'gauss',0.02); %添加高斯噪声 J = imnoise(I,'salt & pepper',0.02); %添加椒盐噪声 ave1=fspecial('average',3); %产生3×3的均值模版 ave2=fspecial('average',5); %产生5×5的均值模版 K = filter2(ave1,J)/255; %均值滤波3×3 L = filter2(ave2,J)/255; %均值滤波5×5 M = medfilt2(J,[3 3]); %中值滤波3×3模板 N = medfilt2(J,[4 4]); %中值滤波4×4模板 imshow(I); figure,imshow(J); figure,imshow(K); figure,imshow(L); figure,imshow(M); figure,imshow(N); 三、实验具体实现

数字图像处理技术及其应用_李红俊

·620· 计算机测量与控制.2002.10(9)  Computer Measurement &Control 设计与应用 收稿日期:2001-12-04。作者简介:李红俊(1974-),男,山西省平遥县人,硕士研 究生,主要从事机械电子方向的研究。 文章编号:1671-4598(2002)09-0620-03 中图分类号:T P391.41 文献标识码:B 数字图像处理技术及其应用 李红俊,韩冀皖 (太原理工大学机械工程学院,山西太原 030024) 摘要:介绍了数字图像处理的基本概念、基本原理,对其中一些算法进行了详细的说明,对不同算法进行了比较。同时,在对现有图像处理方法进行应用的同时,对滤波做了一些新的尝试。最后,将像素细分算法应用于实际生产中, 获得了较好的效果。 关键词:数字图像处理;边缘检测;滤波;像素细分算法 Digital Image Processing and Its Application LI Hong -jun ,HAN Ji -w an (Taiy uan University of T echnolo gy ,T aiyuan 030024,China ) Abstract :T he basic co ncepts and basic principals of digital imag e processing are introduced .Some arithmetics and compari -so n between different arithme tics are expounded .New methods of sieve are adopted when existing image processing methods is being applied .A t last ,the arithmetic of subpixel is applied into practice and obtains effect preferably . Key words :digital image processing ;edge detecting ;sieve ;arithmetic of subpixel 1 序言 图像处理技术基本可以分成两大类:模拟图像处 理(Analog Image Processing )和数字图像处理(Dig -ital Image Processing )。数字图像处理,通俗地讲就是利用计算机对图像进行处理。因此也称之为计算机图像处理(Computer Image Processing )。其优点是处理精度高,处理内容丰富,可进行复杂的非线性处理,有灵活的变通能力,一般来说只要改变软件就可以改变处理内容。存在的问题主要在于处理速度,特别是进行复杂的处理更是如此。数字图像处理概括地说主要包括如下几项内容:几何处理(Geometrical Pro -cessing )、算术处理(Arithmetic Processing )、图像增强(Image Enhancement )、图像复原(Image Restora -tion )、图像重建(Image Reconstruction )、图像编码(Image Encoding )、图像识别(Im age Recognition )、图像理解(Image Understanding )。图像处理技术的发展涉及越来越多的基础理论知识,雄厚的数理基础及相关的边缘学科知识对图像处理科学的发展有越来越大的影响。总之,图像处理科学是一项涉及多学科的综合性科学。 2 边缘检测 所谓边缘应是物体的轮廓或物体不同表面之间的交界在图像中的反映。它的形成是由于物体的材料不同或表面的朝向不同,引起在图像中的边缘处存在明暗、色彩、纹理的变化。因此反过来在图像中检查不 同灰度、色彩等特性区域的交界处就可得到边缘。边缘轮廓是人类识别物体形状的重要因素,也是图像处理中重要的处理对象。 图1 边缘和灰度值模型示意图 如上所述,边缘常常发生在灰度突然变化的部 位,如图1(a )所示,两边为不同的灰度级g 1、g 2,则x 0处为边缘。但实际上由于物体表面交界处灰度常常缓慢变化,在图像中表现为边缘是有一定宽度的,如图1(b )所示,而且由于物体表面的曲折变化加上噪声干扰,边缘时常显得模糊不清,这给边缘的检测带来一定的困难。另外,有的物体本身为条状的区域,例如河流、道路或物体表面的裂缝,它们的边缘表现为狭长的平行线(1~2个像元宽度),如图1(c )所示,而且两边灰度相同或相近,因此检查的方法也有所不同。 边缘检测主要采用各种算法来发现、强化图像中那些可能存在边缘的像素点。边缘检测算子可分为微分(梯度)法、模板匹配法和区域拟合法3种基本方法。对于边缘检测影响较大的是图像中的噪声、退化、模糊等因素,这些都需要特殊的算法来解决。 3 滤波 当图像输入到计算机的时候,由于输入转换器件(如光敏器件、A /D 转换器等性质的差别)及周围环 DOI :10.16526/j .cn ki .11-4762/tp .2002.09.022

数字图像处理实验 数学形态学其应用

实验四:数学形态学及其应用 1.实验目的 1.了解二值形态学的基本运算 2.掌握基本形态学运算的实现 3.了解形态操作的应用 2.实验基本原理 腐蚀和膨胀是数学形态学最基本的变换,数学形态学的应用几乎覆盖了图像处理的所有领域,给出利用数学形态学对二值图像处理的一些运算。 膨胀就是把连接成分的边界扩大一层的处理。而收缩则是把连接成分的边界点去掉从而缩小一层的处理。 二值形态学 I(x,y), T(i,j)为 0/1图像Θ 腐蚀:[]),(&),(),)((),(0 ,j i T j y i x I AND y x T I y x E m j i ++=Θ== 膨胀:[]),(&),(),)((),(0,j i T j y i x I OR y x T I y x D m j i ++=⊕== 灰度形态学 T(i,j)可取10以外的值 腐蚀: []),(),(min ),)((),(1 ,0j i T j y i x I y x T I y x E m j i -++=Θ=-≤≤ 膨胀: []),(),(max ),)((),(1 ,0j i T j y i x I y x T I y x D m j i +++=⊕=-≤≤ 1.腐蚀Erosion: {}x B x B X x ?=Θ: 1B 删两边 2B 删右上 图5-1 剥去一层(皮) 2.膨胀Dilation: {}X B x B X x ↑⊕:= 1B 补两边

2B 补左下 图5-2 添上一层(漆) 3.开运算open :B B X ⊕Θ=)(X B 4.闭close :∨ Θ⊕=B B X X B )( 5.HMT(Hit-Miss Transform:击中——击不中变换) 条件严格的模板匹配 ),(21T T T =模板由两部分组成。1T :物体,2T :背景。 {} C x x i X T X T X T X ??=?21, 图5-3 击不中变换示意图 性质: (1)φ=2T 时,1T X T X Θ=? (2))()()(21T X T X T X C Θ?Θ=? C T X T X )()(21Θ?Θ= )/()(21T X T X ΘΘ= 6.细化/粗化 (1)细化(Thin ) C T X X T X XoT )(/??=?= 去掉满足匹配条件的点。 图5-4 细化示意图 系统细化{}n B oB XoB T Xo ))(((21=, i B 是1 -i B 旋转的结果(90?,180?,270?) 共8种情况 适于细化的结构元素 1 1110 00d d I = d d d L 1011 00= (2)粗化(Thick ) X 21 1 1 2 3 T ? XoT X ? X X ?T X ΘT T ⊕

数学形态学去噪

目录 一绪论 (1) 1.1 数学形态学简介 (1) 1.2 数学形态学与数字图像处理 (1) 1.3 本次课程设计的目的与要求 (2) 二数学形态学的基本运算 (3) 2.1 基本概念 (3) 2.1.1结构元素 (3) 2.1.2膨胀与腐蚀 (3) 2.2 二值形态学图像处理 (4) 2.2.1 膨胀 (4) 2.2.2 腐蚀 (6) 2.2.3 开运算 (7) 2.2.4 闭运算 (8) 2.3 灰度形态学图像处理 (9) 2.3.1 膨胀 (9) 2.3.2 腐蚀 (10) 2.3.3 开运算与闭运算 (11) 2.4 综述 (13) 三数学形态学滤波器去噪 (15) 3.1 概述 (15)

3.2噪声模型 (16) 3.2.1 高斯噪声 (16) 3.2.2 椒盐噪声 (16) 3.3形态学滤波器 (17) 3.4形态学图像去噪原理 (20) 3.5形态学图像去噪的应用 (20) 小结与体会 (21) 参考文献 (22) 附录 (23) 一绪论 1.1数学形态学简介 数学形态学作为一门新兴的图像处理与分析学科,1964年由法国的G.Mathern和J.Serra在积分几何的基础上首次创立。70年代初,采用数学形态学的学者们开拓了图像分析的一个新的领域。经过十多年的理论与实践探索,G.Mathern和J.Serra等人在研究中认识到,对图像先作开运算接着再作闭运算,可以产生一种幂等运算;采用递增尺寸的交变开闭序列作用于图像,可有效地消除图像的噪声,1982年他们正式提出了形态学滤波器的概念。90年代数学形态学有两个显著的发展趋势,第一个是致力于运动分析,包括编码与运动景物描述;第二个是算法与硬件结构的协调发展,用于处理数值函数的形态学算子的开发与设计。

数字图像处理实验5 图像压缩

实验5 图像压缩 一.实验目的: 1.掌握图像压缩的原理——编码冗余,压缩比C R的计算等。 2.了解并掌握霍夫曼编码的原理、实现步骤。 3.掌握JPEG标准——通用的图像压缩/解压缩编码标准。 二.实验内容: 1.利用已给出的MATLAB自编函数库matlab_function文件夹,实现压缩比的计算。 2.对信号源符进行霍夫曼编码,以消除信源的冗余数据。 3.练习JPEG标准的压缩/解压缩技术。 三.实验原理: 1.图像压缩比C R的计算 函数imratio(f1, f2),计算图像压缩比C R,该函数来自MATLAB自编函数库matlab_function文件夹,语法如下: imratio(imread(‘filename’), ‘filename.jpg’) //第二个参数‘filename.jpg’仅是文件名,实际上是一个结构,内含压缩 //后的各种压缩信息,并不代表图像本身 >>f = imread(‘E:\医学图像处理实验讲义\实验五\car_lady.jpg’) >>imfinfo E:\医学图像处理实验讲义\实验五\car_lady.jpg //查看图像文件的详细信息 >>imwrite(f, ‘car_lady25.jpg’, ‘quality’, 25) //将压缩后的图像存到MATLAB默认路径中 >>imfinfo car_lady25.jpg //可依据图像信息计算出压缩率 >>f25 = imread(‘car_lady25.jpg’) >>Cr = imratio (f25, ‘car_lady25.jpg’) 2.霍夫曼编码 符号概率 a1 0.1875 a2 0.5 a3 0.125 a4 0.1875 函数huffman(p)进行霍夫曼编码,语法: huffman(p) //p为向量符号 >>p = [0.1875 0.5 0.125 0.1875] >>c = huffman(p)

基于数学形态学的图像分割方法研究

基于数学形态学的图像分割方法研究 专业:电子信息科学与技术 班级:2005级1班 姓名:杨晓琦

引言 3 1 图像分割基本理论7 1.1 图像分割的概念7 1.2 传统的图像分割方法9 1.3 特殊理论工具的图像分割方法 12 1.4 图像分割的评价 13 2 数学形态学基本理论16 2.1 形态学的概念 16 2.2 结构元素的选取 16 2.3 二值形态学理论 18 2.4 灰值形态学理论 20 2.5 形态学重建 21 2.6 形态学边缘检测 22 3 Matlab在图像分割处理中的应用 24 3.1 Matlab简介.24 3.2 Matlab在图像处理方面的应用.24 3.3 基于Matlab的图像分割.26 4 车牌图像分割的相关理论研究28 4.1 车牌定位算法简介 28 4.2 车牌的字符图像分割 37 5 基于数学形态学车牌图像分割42 5.1 形态学车牌定位 42 5.2 形态学字符图像分割 53 5.3 本章小结 61 结论62 致谢63 参考文献64附录1 源程序清单 68

本文运用形态学方法对车牌定位算法和车牌字符分割算法进行了系统的研究。这两种算法的研究为车牌识别做了先期准备,是智能交通中非常重要的组成部分。在车牌定位算法部分,提出了一种基于二值面积形态学的车牌定位算法。首先将车牌的灰度图像二值化,然后逐步缩小车牌候选区的面积,计算车牌图像中连通区域的面积,并根据车牌图像的实际情况确定面积阈值,并用形态学的方法对车牌图像进行处理运算,以实现车牌的精确定位。仿真实验结果表明此算法定位精度高,而且能适应复杂天气环境,能达到满意的定位效果。在字符分割算法部分,对投影分割算法进行了改进,将其与数学形态学分割算法进行了结合。首先将车牌图像二值化,然后用数学形态学分割方法结合水平与垂直投影分割方法,确定车牌字符宽和高并校正车牌实际位置,通过两次投影,校正车牌角度,去除车牌边框,确定车牌上下边界及中心点,最后分割提取车牌字符。实验结果表明该算法能有效的保持车牌字符边缘,获得较好的分割效果。 关键词: 图像分割;数学形态学;车牌定位;字符分割; Matlab

实验五 图像复原

信息工程学院实验报告 课程名称:数字图像处理Array 实验项目名称:实验五图像复原实验时间: 班级:姓名:学号: 一、实验目的 1.了解图像退化/复原处理的模型; 2. 掌握图像复原的原理及实现方法; 3. 通过本实验掌握利用MATLAB编程实现图像的恢复。 4. 掌握matlab代码的调试方法,熟悉常见代码错误及改正方法。 二、实验步骤及结果分析 MATLAB图像处理工具箱包含四个图像复原函数,请参照教材第126页例6.8编程实现图像复原。 1.用点扩散(PSF)函数创建运动模糊图像,修改参数改变模糊程度。 a) 无噪声运动模糊图像 b) 有噪声运动模糊图像 程序代码: I=imread('cameraman.tif'); %读取图像 subplot(1,3,1); imshow(I,[]);%显示图像 title('原始图像'); PSF=fspecial('motion',25,11); %运动模糊函数,运动位移是25像素,角度是11 Blurred=imfilter(I,PSF,'conv','circular'); %对图像运动模糊处理 subplot(1,3,2); imshow(Blurred,[]);title('无噪声运动模糊图像'); %显示无噪声运动模糊图像 Noise=0.05*randn(size(I)); %正态分布的随机噪声 BlurredNoisy=imadd(Blurred,im2uint8(Noise));%对退化后的图像附加噪声 subplot(1,3,3); imshow(BlurredNoisy,[]);title('有噪声运动模糊图像'); %显示运动模糊且加噪声后图像 执行结果:

图形学实验5

图形学实验5

《3D游戏图形学》实验报告书 (实验五)

姓名: 学号: 班级: 浙江理工大学 二0一二年十二月

实验五纹理映射实验 一、实验目的和要求 掌握纹理映射的基本原理,利用VC++ OpenGL实现纹理映射技术。 二、实验原理 纹理映射是真实感图形制作的一个重要部分,运用纹理映射可以方面地制作真实感图形,而不必花更多的时间去考虑物体的表面纹理。如一张木制桌子其表面的木纹是不规范的,看上去又是那么自然,如果在图形制作中不用纹理映射,那么只是这张桌面纹理的设计,就要花费很大精力,而且设计结果也未必能像现实中那么自然。如果运用纹理映射就非常方便,可以用扫描仪将这样的一张桌子扫成一个位图。然后的具体的操作中,只需把桌面形状用多边形画出来,把桌面纹理贴上去就可以了。 另外,纹理映射能够在多边形进行变换时仍保证纹理的图案与多边形保持一致性。例如,以透视投影方式观察墙面时,远端的砖会变小,而近处的砖就会大一些。 此外,纹理映射也可以用于其他方面。例如,使用一大片植被的图像映射到一些连续的多边

形上,以模拟地貌,或者以大理石、木纹等自然物质的图像作为纹理映射到相应的多边形上,作为物体的真实表面。 在OpenGL中提供了一系列完整的纹理操作 函数,用户可以用它们构造理想的物体表面,可以对光照物体进行处理,使其映射出所处环境的景象,可以用不同方式应用到曲面上,而且可以随几何物体的几何属性变换而变化,从而使制作的三维场景和三维物体更真实更自然。 在OpenGL中要实现纹理映射,需要经历创建纹理、指定纹理应用方式、启用纹理映射、使用纹理坐标和几何坐标绘制场景几个过程。 用于指定一维、二维和三维纹理的函数分别为: Void glTexImage1D(GLenum target, Glint level, Glint components, GLsizei width, Glint border, GLenum format, GLenum type, const GLvoid *texels); Void glTexImage2D(GLenum target, Glint level, Glint components, GLsizei width, GLsizei height, GLint border, GLenum format, GLenum type, const GLvoid *texels);

相关文档
最新文档