毕业设计（论文）-基于Matlab的图像预处理算法实现.doc

《毕业设计（论文）-基于Matlab的图像预处理算法实现.doc》由会员分享，可在线阅读，更多相关《毕业设计（论文）-基于Matlab的图像预处理算法实现.doc（32页珍藏版）》请在三一文库上搜索。

1、I 基于 Matlab 的图像预处理算法实现 （吉首大学张家界学院，湖南 吉首 416000） 摘 要 在图像处理的过程中，由于获取图像的工具或手段的影响，使获取图像无法完 全体现原始图像的全部信息。因此，以改善图像数据、抑制不需要的变形或者增强 某些对于后续处理来说比较重要的图像特征为目的的图像预处理在图像处理的过程 中就显得非常重要。 常用的图像预处理手段包括图像大小的修改、图像浓度的扩展，图像由彩色差 到灰度的转变、图像二值化、图像的锐化处理及图像平滑处理等方法。本文在简单 介绍数字图像处理技术的基础上，着重讨论了灰度直方图的构建、直方图均衡、图 像平滑、锐化和噪声处理的原理及实现方法，

2、并用 matlab 开发工具实现了上述图 像预处理算法并给出了各种算法的处理结果。 关键词：图像处理；灰度直方图；图像平滑；噪声处理；图像锐化；matlab Image preprocessing algorithm Based on the Matlab He Jinze (Zhangjiajie College of Jishou University,Jishou,Hunan 416000) Abstract In the process of image processing, due to the affection of image tools, the acquisition i

3、mage cant reflect the original image information completely. Therefore, to improve image data, inhibit deformation and enhance image characteristics is very important in image processing. Common image preprocessing method include the image size modification, the image thickness expansion, changing c

4、olor image to gray, image gray-scale binary, image sharpening and the image smooth processing. Based on the brief introduction of digital image processing technique, this paper discussed the construction of gray histogram, equilibrium of histogram, the image smoothing， sharpening and noise treatment

5、s principle and method. The realizations of all those image preprocessing algorithm which is II realized with matlab has been discussed in this paper, and all those result has been given after the Implementations. Key word:Image processing; grayscale histogram; image smoothing; noise processing; ima

6、ge sharpening ;matlab III 目 录 第一章 绪 论 1 1.1 何谓数字图像处理.1 1.2 数字图像处理的特点及其应用 1 1.2.1 数字图像处理的特点.1 1.2.2 图像预处理的内容 .2 1.2.3 数字图像处理的应用.3 1.3 MATLAB 4 1.3.1 matlab 简述 4 1.3.2 matlab 处理图像的特点 5 第二章 数字图像处理的灰度直方图 .6 2.1 灰度的定义.6 2.2 直方图定义.6 2.2.1 直方图的典型用途 .6 2.2.2 灰度直方图的计算 .7 2.2.3 图像直方图实现代码 .7 2.3 直方图均衡.8 2.3.1

7、直方图均衡原理.8 2.3.2 直方图均衡的实现 .8 第三章 图像平滑与图像锐化 12 3.1 图像的平滑12 3.1.1 领域平均法基础理论 12 3.1.2 算法实现 13 3.2 图像锐化15 3.2.1 图像锐化的目的和意义 15 3.2.2 图像锐化算法 15 3.2.3 图像锐化的实现代码 16 第四章 图像噪声与噪声的处理 18 4.1 噪声的概念18 4.2 图像噪声对图像的影响18 4.3 噪声来源 .18 IV 4.4 噪声图像模型及噪声特性 .19 4.4.1 含噪模型19 4.4.2 噪声特性20 4.5 图像二值化 .20 4.5.1 理论基础 20 4.5.2 图

8、像二值化的实现代码 20 4. 6 二值图像的去噪.21 4.6.1 理论基础 21 4.6.2 二值图像去噪的实现代码 22 第五章 结论 24 参考文献 .25 图像预处理系统设计 绪论 1 第一章 绪 论 1.1 何谓数字图像处理 数字图像处理（Digital Image Processing） ，就是利用数字计算机或则其 他数字硬件，对从图像信息转换而得到的电信号进行某些数学运算，以提高图 像的实用性。例如从卫星图片中提取目标物的特征参数，三维立体断层图像的 重建等。总的来说，数字图像处理包括点运算、几何处理、图像增强、图像复 原、图像形态学处理、图像编码、图像重建、模式识别等。目前数

9、字图像处理 的应用越来越广泛，已经渗透到工业、医疗保健、航空航天、军事等各个领域， 在国民经济中发挥越来越大的作用。 1.2 数字图像处理的特点及其应用 在计算机出现之前，模拟图像处理占主导地位。随着计算机的发展，数字 图像处理发展速度越来越快。尽管目前一般采用顺序处理的计算机，对大数据 量的图像处理速度不如光学方法快，但是其处理的精度高，实现多种功能的、 高度复杂的运算求解非常灵活方便。在其短短的历史中，它却成功的应用于几 乎所有与成像有关的领域，并正发挥相当重要的作用。 1.2.1 数字图像处理的特点 同模拟图像处理相比，数字图像处理有很多优点。主要表现在： 1. 精度高 不管是对 4bi

10、t 还是 8bit 和其他比特图像的处理，对计算机程序来说几乎 是一样的。即使处理图像变大，只需改变数组的参数，而处理方法不变。所以 从原理上不管处理多高精度的图像都是可能的。而在模拟图像处理中，要想使 精度提高一个数量级，就必须对处理装置进行大幅度改进。 2. 再现性好 不管是什么图像，它们均用数组或集合表示。将它们输入到计算机内，用 计算机容易处理的方式表示。在传送和复制图像时，只在计算机内部进行处理， 这样数据就不会丢失或遭破坏，保持了完好的再现性。而在模拟图像处理中， 就会因为各种干扰及设备故障而无法保持图像的再现性。 3. 通用性、灵活性高 图像预处理系统设计 绪论 2 不管是可视图

11、像还是 X 线照片、红外线热成像、超声波图像等不可见光成 像，尽管这些图像成像体系中的设备规模和精度各不相同，但当把图像信号直 接进行 A/D 变换，或记录成照片再数字化，对于计算机来说都能用二维数组表 示，不管什么样的图像都可以用同样的方法进行处理，这就是计算机处理的通 用性。另外，对处理程序自由加以改变，可进行各种各样的处理。如上下滚动、 漫游、拼图、合成、变换、放大、缩小和各种逻辑运算等，所以灵活性很高。 1.2.2 图像预处理的内容 在图像分析中，对输入图像进行特征抽取、分割和匹配前所进行的处理。 图像预处理的主要目的是消除图像中无关的信息，恢复有用的真实信息，增强 有关信息的可检测性

12、和最大限度地简化数据，从而改进特征抽取、图像分割、 匹配和识别的可靠性。预处理过程一般有数字化、几何变换、归一化、平滑、 复原和增强等步骤。 1.数字化 一幅原始照片的灰度值是空间变量（位置的连续值）的连续函数。在 MN 点阵上对照片灰度采样并加以量化（归为 2b 个灰度等级之一），可以得到计算 机能够处理的数字图像。为了使数字图像能重建原来的图像,对 M、N 和 b 值的 大小就有一定的要求。在接收装置的空间和灰度分辨能力范围内,M、N 和 b 的 数值越大,重建图像的质量就越好。当取样周期等于或小于原始图像中最小细节 周期的一半时，重建图像的频谱等于原始图像的频谱，因此重建图像与原始图 像

13、可以完全相同。由于 M、N 和 b 三者的乘积决定一幅图像在计算机中的存储 量，因此在存储量一定的条件下需要根据图像的不同性质选择合适的 M、N 和 b 值，以获取最好的处理效果。 2.几何变换 用于改正图像采集系统的系统误差和仪器位置的随机误差所进行的变换。 对于卫星图像的系统误差，如地球自转、扫描镜速度和地图投影等因素所造成 的畸变,可以用模型表示,并通过几何变换来消除。随机误差如飞行器姿态和高 度变化引起的误差，难以用模型表示出来，所以一般是在系统误差被纠正后， 通过把被观测的图和已知正确几何位置的图相比较，用图中一定数量的地面控 制点解双变量多项式函数组而达到变换的目的。 3.归一化

14、使图像的某些特征在给定变换下具有不变性质的一种图像标准形式。图像 的某些性质，例如物体的面积和周长，本来对于坐标旋转来说就具有不变的性 质。在一般情况下，某些因素或变换对图像一些性质的影响可通过归一化处理 得到消除或减弱，从而可以被选作测量图像的依据。例如对于光照不可控的遥 图像预处理系统设计 绪论 3 感图片，灰度直方图的归一化对于图像分析是十分必要的。灰度归一化、几何 归一化和变换归一化是获取图像不变性质的三种归一化方法。 4.图像平滑 消除图像中随机噪声的技术。对平滑技术的基本要求是在消去噪声的同时 不使图像轮廓或线条变得模糊不清。常用的平滑方法有中值法、局部求平均法 和 k 近邻平均法

15、。局部区域大小可以是固定的，也可以是逐点随灰度值大小变 化的。此外，有时应用空间频率域带通滤波方法。 5.复原 校正各种原因所造成的图像退化，使重建或估计得到的图像尽可能逼近于 理想无退化的像场。在实际应用中常常发生图像退化现象。例如大气流的扰动， 光学系统的像差，相机和物体的相对运动都会使遥感图像发生退化。基本的复 原技术是把获取的退化图像 g(x，y)看成是退化函数 h(x，y)和理想图像 f(x，y)的卷积。它们的傅里叶变换存在关系 G(u，vH(u，v)F(u，v)。根据 退化机理确定退化函数后，就可从此关系式求出 F(u，v)，再用傅里叶反变换 求出 f(x，y)。通常把 称为反向滤

16、波器。实际应用时，由于 H(u，v)随离开 uv 平面原点的距离增加而迅速下降，为了避免高频范围内噪声的强化，当 u2+v2 大于某一界限值 W 娿时，使 M(u，v)等于 1。W0 的选择应使 H(u，v)在 u2+v2W 范围内不会出现零点。图像复原的代数方法是以最小二乘法最佳准则 为基础。寻求一估值,使优度准则函数值最小。这种方法比较简单，可推导出最 小二乘法维纳滤波器。当不存在噪声时，维纳滤波器成为理想的反向滤波器。 6.增强 对图像中的信息有选择地加强和抑制，以改善图像的视觉效果，或将图像 转变为更适合于机器处理的形式，以便于数据抽取或识别。例如一个图像增强 系统可以通过高通滤波器来

17、突出图像的轮廓线，从而使机器能够测量轮廓线的 形状和周长。图像增强技术有多种方法，反差展宽、对数变换、密度分层和直 方图均衡等都可用于改变图像灰调和突出细节。实际应用时往往要用不同的方 法，反复进行试验才能达到满意的效果。 1.2.3 数字图像处理的应用 计算机图像处理和计算机、多媒体、智能机器人、专家系统等技术的发展 紧密相关。近年来计算机识别、理解图像的技术发展很快，也就是图像处理的 目的除了直接供人观看（如医学图像是为医生观看作诊断）外，还进一步发展 了与计算机视觉有关的应用，如邮件自动分检，车辆自动驾驶等。下面仅罗列 了一些典型应用实例，而实际应用更广。 1. 在生物医学中的应用 图像

18、预处理系统设计 绪论 4 主要包括显微图像处理；DNA 显示分析；红、白血球分析计数；虫卵及组 织切片的分析；癌细胞的识别；染色体分析等等。 2. 遥感航天中的应用 军事侦察、定位、导航、指挥等应用；多光谱卫星图像分析；地形、地图、 国土普查；地质、矿藏勘探；天文、太空星体的探测及分析等。 3. 工业应用 CAD 和 CAM 技术用于模具、零件制造、服装、印染业；零件、产品无损检 测，焊缝及内部缺陷检查；交通管制、机场监控；火车车皮识别等。 4. 军事公安领域中的应用 巡航导弹地形识别；指纹自动识别；警戒系统及自动火炮控制；反伪装侦 察；手迹、人像、印章的鉴定识别；过期档案文字的复原；集装箱的

19、不开箱检 查等。 5. 其他应用 图像的远距离通信；多媒体计算机系统及应用；电视电话；服装试穿显示； 理发发型预测显示；电视会议；办公自动化、现场视频管理等。 1.3 MATLAB 1.3.1 matlab 简述 MATLAB 是矩阵实验室（Matrix Laboratory）的简称，是美国 MathWorks 公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计 算的高级技术计算语言和交互式环境，主要包括 MATLAB 和 Simulink 两大部分。 MATLAB 是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交 互式程序设计的高科技计算环境。它将数值

20、分析、矩阵计算、科学数据可视化 以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗 环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提 供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言 （如 C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应 用软件中在数值计算方面首屈一指。MATLAB 可以进行矩阵运算、绘制函数和数 据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程 计算、控制设计、信号处理与通讯、图像处理

21、、信号检测、金融建模设计与分 析等领域。 图像预处理系统设计 绪论 5 MATLAB 的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形 式十分相似，故用 MATLAB 来解算问题要比用 C，FORTRAN 等语言完成相同的事 情简捷得多，并且 MATLAB 也吸收了像 Maple 等软件的优点,使 MATLAB 成为一个 强大的数学软件。在新的版本中也加入了对 C，FORTRAN，C+ ，JAVA 的支持。 可以直接调用,用户也可以将自己编写的实用程序导入到 MATLAB 函数库中方便 自己以后调用，此外许多的 MATLAB 爱好者都编写了一些经典的程序，用户可以 直接进行下载就可以

22、用。 1.3.2 matlab 处理图像的特点 第一，MATLAB 的基本数据单位是矩阵，它的指令表达式与数学,工程中常 用的形式十分相似,故用 MATLAB 来解算问题要比用 C,FORTRAN 等语言完成相同 的事情简捷得多。 第二，友好的工作平台和编程环境。MATLAB 由一系列工具组成。这些工具 方便用户使用 MATLAB 的函数和文件，其中许多工具采用的是图形用户界面。包 括 MATLAB 桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于 用户浏览帮助、工作空间、文件的浏览器。 第三，简单易用的程序语言。Matlab 一个高级的距阵/阵列语言，它包含 控制语句、函数、数据

23、结构、输入和输出和面向对象编程特点。用户可以在命 令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用 程序(M 文件)后再一起运行。 第四，强大的科学计算机数据处理能力。MATLAB 是一个包含大量计算算法 的集合。其拥有 600 多个工程中要用到的数学运算函数，可以方便的实现用户 所需的各种计算功能。 第五，出色的图形处理功能，MATLAB 自产生之日起就具有方便的数据可视 化功能，以将向量和距阵用图形表现出来，并且可以对图形进行标注和打印。 高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用 于科学计算和工程绘图。 第六，应用广泛的模块集合工具箱，MAT

24、LAB 对许多专门的领域都开发了功 能强大的模块集和工具箱。一般来说，他们都是由特定领域的专家开发的，用 户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。 文章主要分为四部分，其中，文章的开头简单介绍了数字图像应用的发展 过程，还对本文应用的 matlab 开发工具的发展进行简单讲述。文章的第二部分 主要针对图像预处理系统中怎么做到灰度转换及灰度直方图的构建。第三部分 讲的是图像的预处理技术，其中包括图像平滑、锐化等。第四部分讲的是噪声 处理，噪声特性以及 matlab 去除噪声的方法。 图像预处理系统设计 数字图像处理的灰度直方图 6 第二章 数字图像处理的灰度直方图 2

25、.1 灰度的定义 灰度使用黑色调表示物体。 每个灰度对象都具有从 0%（白色）到 图 2.1 灰度条 100%（黑色）的亮度值。 使用黑白或 灰度扫描仪生成的图像通常以灰度显 示。使用灰度还可将彩色图稿转换为高质量黑白图稿。 灰度级变换(点运算)的定义 对于输入图象 f(x，y)，灰度级变换 T 将产生一个输出图像 g(x，y)，且 g(x，y)的每一个像素值都是由 f(x，y)的对应输入像素点的值决定的，g(x，y) =T(f(x，y)。 对于原图象 f(x，y)和灰度值变换函数 T(f(x，y)，由于灰度值总是有 限个(如：O2 5 5)，非几何变换可定义为：R=T(r)。 2.2 直方图

26、定义 直方图：图像直方图是图像处理中一种十分重要的图像分析工具，它描述 了一幅图像的灰度级内容，任何一幅图像的直方图都包含了丰富的信息，它主 要用在图象分割，图像灰度变换等处理过程中。 从数学上来说图像直方图是图像各灰度值统计特性与图像灰度值的函数， 它统计一幅图像中各个灰度级出现的次数或概率；从图形上来说，它是一个二 维图，横坐标表示图像中各个像素点的灰度级，纵坐标为各个灰度级上图像各 个像素点出现的次数或概率。直方图的纵坐标都对应着该灰度级在图像中出现 的概率。 2.2.1 直方图的典型用途 图像预处理系统设计 数字图像处理的灰度直方图 7 对图像进行数字化时，利用直方图可以检查输入图像的

27、灰度值在可利用的 灰度范围内分配得是否适当。 在医学方面，为了改善 X 射线操作人员的工作条件，可采用低辐射 X 射线 曝光，但这样获得的 X 光片灰度级集中在暗区，导致某些图像细节无法看清， 判读困难。通过直方图修正使灰度级分布在人眼合适的亮度区域，便可使 X 片 中的细节清晰可见。 可以根据直方图确定二值化的阈值；当物体部分的灰度值比其它部分的灰 度值大时，可以用直方图求出物体的面积(实际上是象素数=灰度大于和等于 q 的象素的总和)；当物体部分的灰度值比其它部分的灰度值大时，可以用直方图 求出物体的面积(实际上是象素数=灰度大于和等于 q 的象素的总和)；利用色彩 直方图可以进行基于颜色

28、的图象分割。 2.2.2 灰度直方图的计算 若图象具有 L 级灰度(通常 L=256，即 8 位灰度级)，则大小为 m(n 的灰度 图象 f(x，Y)的灰度直方图 Hk，k=0L-1，可按如下步骤计算获得： 1) 初始化数据，清空表中所以内容； 2) 对图像中的各个像素中的灰度级进行统计，并把其相各个灰度级中所对 应的像素点统计下来； 3) 用横坐标表示图像中各个像素点的灰度级，纵坐标表示各个灰度级上图 像各个像素点出现的次数或概率。 2.2.3 图像直方图实现代码 A,map=imread(f:234.bmp); imshow(A,map), image=double(A); for i=1

29、:256, a(i)=0; end, for i=1:256, for j=1:256, b=image(i,j)+1; a(b)=a(b)+1; end, end, for k=1:256, p(k)=a(k)/(2562); 图像预处理系统设计 数字图像处理的灰度直方图 8 end, figure,stem(p); 用上述代码生成图 2.2 的直方图，结果如图 2.3 所示。 050100150200250300 0 0.005 0.01 0.015 0.02 0.025 图 2.2 原图 图 2.3 直方图 2.32.3 直方图均衡直方图均衡 中文名称：直方图均衡 英文名称：histog

30、ram equalization 定 义：使原直方图变换为具有均匀密度分布的直方图，然后按该直方图调 整原图像的一种图像处理技术。 2.3.1 直方图均衡原理 直方图均衡化是灰度变换的一个重要应用，广泛应用在图像增强处理中， 它是以累计分布函数变换为基础的直方图修正法，可以产生一幅灰度级分布具 有均匀概率密度的图像，扩展了像素的取值动态范围。若像素点的原灰度为 R，变换后的灰度为 S，需要注意的是 R、S 是归一化后的灰度值，其灰度变换 函数 T(R)为： （2- 1 , 010)()( 0 klR n n RpRTS j k oj j j k j r ； 1） 式中， 是第 j 级灰度值的概

31、率，是图像中 j 级灰度的像素总数，)( j RPr j n 是图像中灰度级的总数目，是图象中像素的总数。对变换后的 S 值取最靠近ln 的一个灰度级的值，建立灰度级变换表，将原图像变换为直方图均衡的图像。 下面是实现图像直方图均衡化函数的源代码和效果图： 2.3.2 直方图均衡的实现 (1) 对给定的待处理图像统计其直方图，求出. NnrP kkr / 图像预处理系统设计 数字图像处理的灰度直方图 9 (2) 根据统计出的直方图采用累积分布函数作变换, 求变换后的新灰度； （2-)()( 0 kj k j rk RpRTS 2） (3) 用新灰度代替旧灰度，求出 Ps(s)，这一步是近似过程

32、，应根据处理 目的尽量做到合理，同时把灰度值相等或近似地合并到一起。 A,map=imread(f:456.bmp); imshow(A,map), title(原图); image=double(A); for i=1:256, n(i)=0; end, for i=1:256, for j=1:256, s=image(i,j)+1; n(s)=n(s)+1; end, end, for k=1:256, p(k)=n(k)/(2562); end, figure,bar(p,r); title(直方图); for k=1:256, q(k)=0; end, for k=1:256, fo

33、r j=1:k, q(k)=q(k)+p(j); end, end, figure,bar(q,y); title(累积直方图); 图像预处理系统设计 数字图像处理的灰度直方图 10 N=256; for k=1:256, o(k)=round(N-1)*q(k)+0.5); end, for i=1:256, for j=1:256, rimage(i,j)=(image(i,j); end end figure,imshow(uint8(rimage); title(均衡化后的图); for i=1:256, n(i)=0; end, for i=1:256, for j=1:256, s

34、=rimage(i,j)+1; n(s)=n(s)+1; end, end, for k=1:256, p(k)=n(k)/(2562); end, figure,bar(p,b); title(均衡化的直方图); 图像预处理系统设计 数字图像处理的灰度直方图 11 050100150200250300 0 0.005 0.01 0.015 0.02 0.025 图 2.4 原图 图 2.5 原始图象的直方图 050100150200250300 0 0.005 0.01 0.015 0.02 0.025 图 2.6 直方图均衡化后的效果图 图 2.7 均衡化后的直方图 从上述效果图可以看出，

35、经过直方图均衡化处理后，图像变的清晰了，从 直方图来看，处理后的 LENA 的图像直方图分布更均匀了，在每个灰度级上图像 都有像素点。但是直方图均衡化存在着两个缺点： （1）变换后图像的灰度级减少，某些细节消失； （2）某些图像，如直方图有高峰，经处理后对比度不自然的过分增强。 图像预处理系统设计 图像平滑与图像锐化 12 第三章 图像平滑与图像锐化 3.1 图像的平滑 图像的平滑方法是一种实用的图像处理技术，能减弱或消除图像中的高频 率分量，但不影响低频率分量。因为高频率分量主要对应图像中的区域边缘等 灰度值具有较大较快变化的部分，平滑滤波将这些分量滤去可减少局部灰度起 伏，使图像变得比较平

36、滑。实际应用中，平滑滤波还可用于消除噪声，或者在 提取较大目标前去除过小的细节或将目标内的小间断连接起来。它的主要目的 是消除图像采集过程中的图像噪声，在空间域中主要利用邻域平均法、中值滤 波法和选择式掩模平滑法等来减少噪声；在频率域内，由于噪声主要存在于频 谱的高频段，因此可以利用各种形式的低通滤波器来减少噪声。本节主要介绍 常用的空间域图像平滑方法。 3.1.1 领域平均法基础理论 最简单的平滑滤波是将原图中一个像素的灰度值和它周围邻近 8 个像素的 灰度值相加，然后将求得的平均值（除以 9）作为新图中该像素的灰度值。它 采用模板计算的思想，模板操作实现了一种邻域运算，即某个像素点的结果不

37、 仅与本像素灰度有关，而且与其邻域点的像素值有关。模板运算在数学中的描 述就是卷积运算，这里不再赘述。邻域平均法也可以用数学公式表达： 设为给定的含有噪声的图像，经过邻域平均处理后的图像为，),(jif),(jig 则: , （3-Mji N jif jig),( , ),( ),( 1） 所取邻域中各邻近像素的坐标，是邻域中包含的邻近像素的个数。邻MN 域平均法的模板为：，中间的黑点表示以该像素为中心元素，即该 111 111 111 9 1 像素是要进行处理的像素。在实际应用中，也可以根据不同的需要选择使用不 同的模板尺寸，如 33、55、77、99 等。 邻域平均处理方法是以图像模糊为代

38、价来减小噪声的，且模板尺寸越大， 噪声减小的效果越显著。如果是噪声点，其邻近像素灰度与之相差很大，),(jif 采用邻域平均法就是用邻近像素的平均值来代替它，这样能明显消弱噪声点， 使邻域中灰度接近均匀，起到平滑灰度的作用。因此，邻域平均法具有良好的 图像预处理系统设计 数字图像处理的灰度直方图 13 噪声平滑效果，是最简单的一种平滑方法。 图像预处理系统设计 图像平滑与图像锐化 14 3.1.2 算法实现 程序中采用的是 33 的模板，即把当前图像和它周围 8 个像素的灰),(jif 度值相加，然后将求得的平均值（除以 9）作为该点的像素值。 A,map=imread(f:234.bmp);

39、 figure(1); imshow(A,map), title(原图); image=double(A); for i=1:256, a(i)=0; end, for i=2:255, for j=2:255, b=image(i-1,j-1)+image(i-1,j)+image(i-1,j+1) +image(i,j-1)+image(i,j)+image(i,j+1) +image(i+1,j-1)+image(i+1,j)+image(i+1,j+1); b=b/9; image(i,j)=b; end, end, figure(2); imshow(image,map), titl

40、e(平滑); 使用上述代码对含有高斯噪声的图像 3.1 分别利用邻域平均法的不同尺寸 模板进行平滑后的实验结果，图 3.2、图 3.3 显示的是分别使用了 33、55 模板平滑后的图像。从实验结果可以看出，当所用平滑模板尺寸增大时，对噪 声的消除效果也有所增强，但同时会带来图像的模糊，边缘细节逐步减少，且 运算量增大。在实际应用中，可以根据不同的应用场合选择合适的模板大小。 图像预处理系统设计 图像平滑与图像锐化 15 二 二 图 3.1 含有噪声的原始图像 二 二 3*3 图 3.2 33 邻域平均法的平滑图像 二 二 5*5 图 3.3 55 邻域平均法的平滑图像 图像预处理系统设计 图像

41、平滑与图像锐化 16 3.2 图像锐化 3.2.13.2.1 图像锐化的目的和意义 图像锐化是一种补偿轮廓、突出边缘信息以使图像更为清晰的处理方法.锐 化的目标实质上是要增强原始图像的高频成分.常规的锐化算法对整幅图像进行 高频增强,结果呈现明显噪声.为此,在对锐化原理进行深入研究的基础上,提出 了先用边缘检测算法检出边缘,然后根据检出的边缘对图像进行高频增强的方法. 实验结果表明,该方法有效地解决了图像锐化后的噪声问题 锐化的目的在于使图像中对象轮廓上的像素灰度大的更大，小的更小，但 对轮廓外的像素不起作用。由于这一原因，图像的锐化对孤立点或对孤立线条 的 边缘增强作用十分明显，但在一定程度

42、上也会对噪声信号产生增强作用。 图像锐化处理的主要目的是突出图像中的细节或者增强被模糊了的细节。 模糊可能是由于错误操作，或者是由于图像获取方法的固有影响所导致的。例 如，当图像的分辨率有限时，所获得的像素值不是一点的亮度，而是周围景物 亮度的平均值。 图像平滑往往使图像中的边界、轮廓变的模糊，为了减少这类不利效果的 影响，这就需要利用图像鋭化技术，使图像的边缘变的清晰。图像銳化处理的 目的是为了使图像的边缘、轮廓线以及图像的细节变的清晰，经过平滑的图像 变得模糊的根本原因是因为图像受到了平均或积分运算，因此可以对其进行逆 运算（如微分运算）就可以使图像变的清晰。从频率域来考虑，图像模糊的实

43、质是因为其高频分量被衰减，因此可以用高通滤波器来使图像清晰。 3.2.2 图像锐化算法 为了要把图像中间任何方向伸展的的边缘和轮廓线变得清晰，我们希望对 图像的某种运算是各向同性的。可以证明偏导平方和的运算是各向同性的，既: （3- 2 2 2 2 y f x f y f x f 2） 式中是图像旋转前的坐标，是图像旋转后的坐标。梯度运算就是在yx,yx , 这个式子的基础上开方得到的。图像（x，y）点的梯度值： （3-fkfg 2 3） 为了突出物体的边缘，常常采用梯度值的改进算法，将图像各个点的梯度 图像预处理系统设计 图像平滑与图像锐化 17 值与某一阈值作比较，如果大于阈值，该像素点的

44、灰度用梯度值表示，否则用 一个固定的灰度值表示。我们在对图像增强的过程中，采用的是一种简单的高 频滤波增强方法： （3-1, 1, 2 1 2 2 yxfyxfyxfyxf y f x f yxG 4） 式中 f，g 分别为锐化前后的图像，是与扩散效应有关的系数。表kf 2 示对图像 f 进行二次微分的拉普拉斯算子。这表明不模糊的图像可以由模糊的 图像减去乘上系数的模糊图像拉普拉斯算子来得到。可以用下面的模板 2 H=1，4，1，4，-20，4，1，4，1来近似。在具体实现时，上述模板 H 中的各个系数可以改变，这个系数的选择也很重要，太大了会使图像的轮廓k 过冲，太小了则图像锐化不明显。实验

45、表明，选取 2-8 之间往往可以达到比k 较满意的效果。下面给出等于 4 的情况下的实现代码和效果图。k 3.2.33.2.3 图像锐化的实现代码 A,map=imread(f:234.bmp); figure(1); imshow(A,map), title(原图); image=double(A); for i=2:255, for j=2:255, b=image(i,j)-2*1*image(i-1,j-1)+4*image(i-1,j)+1*image(i- 1,j+1)+4*image(i,j-1)-20*image(i,j)+4*image(i,j+1)+1*image(i+1,

46、j- 1)+4*image(i+1,j)+1*image(i+1,j+1); image(i,j)=b; end, end, figure(2); imshow(image,map), title(锐化); 图像预处理系统设计 图像平滑与图像锐化 18 二 二 图 3.5 原图 图 3.6 锐化图 从上述效果图可以看出，经过拉普拉斯锐化处理后，图像变的色彩浓密了， 从拉普拉斯锐化图来看，处理后的 LENA 的图像像素点对比度更明显了，图像的 边缘和轮廓都变的清晰了，每个像素点都被体现了出来。 图像预处理系统设计 图像噪声与噪声的处理 19 第四章 图像噪声与噪声的处理 4.1 噪声的概念 目前

47、大多数数字图像系统中，输入图像都是采用先冻结再扫描方式将多维 图像变成一维电信号，再对其进行处理、存储、传输等加工变换。最后往往还 要在组成多维图像信号，而图像噪声也将同样受到这样的分解和合成。在这些 过程中电气系统和外界影响将使得图像噪声的精确分析变得十分复杂。另一方 面图像只是传输视觉信息的媒介，对图像信息的认识理解是由人的视觉系统所 决定的。不同的图像噪声，人的感觉程度是不同的，这就是所谓人的噪声视觉 特性课题。 图像噪声在数字图像处理技术中的重要性越来越明显，如高放大倍数航片的判 读，射线图像系统中的噪声去除等已经成为不可缺少的技术步骤。 4.2 图像噪声对图像的影响 人类获取外界信息有视觉、听觉、触觉、味觉等多种方法，但绝大部分 (约 80%