基于MatlabGUI空域滤波增强的方案设计书.pdf

《基于MatlabGUI空域滤波增强的方案设计书.pdf》由会员分享，可在线阅读，更多相关《基于MatlabGUI空域滤波增强的方案设计书.pdf（53页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、基于 Matlab GUI 空域滤波增强的设计专业名称：电子信息工程班级：08 电子 2 班学生姓名：董浩文指导教师：章瑞平完成时间：2012 年 5 月目 录1 绪论 .11.1 课题研究背景 .11.2 图像去噪的研究现状 .21.3 本文主要工作 .3第 2 章 Matlab 简介 .42.1 Matlab 概况.42.1.1 Matlab 发展过程 .52.1.2 Matlab 的语言特点 .72.2 Matlab 图像处理常用函数 .7第 3 章 图像去噪算法 .103.1 图像噪声概述 .103.1.1 图像噪声的概念 .103.1.2 图像噪声的分类 . 113.2 邻域平均法图

2、像去噪 .123.2.1 模板操作和卷积运算 .123.2.2 邻域平均法原理 .153.3 中值滤波法图像去噪 .163.3.1 中值滤波原理 .173.3.2 中值滤波主要特性 .193.3.3 复合型中值滤波 .223.4 维纳滤波法图像去噪 .233.4.1 维纳滤波原理 .233.4.2 维纳滤波器特性 .25第 4 章 基于 Matlab GUI 的空域增强设计 .264.1 Matlab GUI 设计的介绍 .274.2 Matlab GUI 界面的操作 .274.2.1 启动 GUIDE.274.2.2 GUI 文件的介绍 .294.2.3 GUIDE 环境的设置 .304.3

3、 基于 Matlab GUI 空域增强软件平台的设计 .304.3.1 启动 GUIDE.304.3.2 加入菜单项编辑 .314.3.3 增加子菜单和布局 .324.4 各种方法的仿真 .344.4.1 邻域平均法的仿真 .344.4.2 中值滤波的仿真 .364.4.3 维纳滤波的仿真 .394.4.4 几种去噪方法的比较分析 .414.5 Matlab GUI 软件界面处理的结果 .43致 谢 .45参考文献.46摘要在许多情况下图像信息会受到各种各样噪声的影响，严重时会影响图像中的有用信息，所以对图像的噪声处理就显得十分重要。根据噪声的频谱分布的规律和统计特征以及图像的特点，出现了多种

4、多样的去噪方法。经典去噪方法有：空间域合成法、频域合成法和最优线性合成法等。与之相适应的出现了许多应用方法：如均值滤波器、中值滤波器、低通滤波器、维纳滤波器、最小失真法等。这些方法广泛地应用，促进数字信号处理的极大发展，并且能够显著提高图像质量。本文的主要工作就是研究四种常用去噪方法：邻域平均法、中值滤波法、维纳滤波法和模糊小波变换法的原理，并进行相应的仿真。本文首先介绍了本课题的研究背景以及国内外的研究现状，之后简单介绍了Matlab 软件的发展过程以及特点。第二，介绍了滤波器的工作原理以及分类。第三，在了解滤波器的基本工作原理的基础上，针对不同的数字滤波器所采用最合适的设计方法 ,根据基于

5、 matlab 空域滤波当中的均值 ,中值,维纳等滤波器来设计所需要的内容为了改善图像质量 .从图像中提取有效信息，必须对图像进行去噪预处理。第四，设计 Matlab GUI 的空域滤波软件 ,利用 Matlab GUI 操作界面设计了一个空域增强滤波器，该滤波器具有三种滤波方法，分别为均值滤波、中值滤波以及维纳滤波。最后根据不同的去噪算法（邻域平均法、中值滤波法、维纳滤波法），运用Matlab 软件编写代码，对一张含噪图片（含高斯噪声或椒盐噪声） 进行仿真去噪，并对结果分析讨论，比较几种方法的优缺点，并设计了相应的空域增强软件。通过 Matlab GUI 界面软件设计，可以方便地进行图像空域

6、增强处理，即使对Matlab 软件不熟练的人也可以较为方便和快速地进行操作。只要点击操作界面，载入图像，即可进行相应的图像处理。关键词:中值滤波；均值滤波；维纳滤波；小波变换法；Matlab GUIAbstractIn many casesimages information can be affected by various noises, seriouslyaffect the useful information of a image, so getting rid of noise of image is veryimportant.According to the noise sp

7、ectrum distribution and statistical characteristics andimage characteristics, appeared a variety of denoising method. Classical denoisingmethods are: spatial domain synthesis, frequency synthesis method and optimal linearsynthesis. To adapt to the emergence of many application methods: such as the m

8、eanfilter, median filter, low-pass filter, Wiener filter, minimum distortion method. Thesemethods are widely used, and promote the great developmentof digital signalprocessing, and can significantly improve the quality of images. The main work of thispaper is to study the four kinds of denoising met

9、hods: neighborhood average, medianfiltering, Wiener filtering method and fuzzy wavelet transform method, and thecorresponding simulation.This paper first introduces the research background of this topic and the presentresearch situation at home and abroad, then simply introduced the Matlab softwared

10、evelopment process and characteristics. Second, introduces the working principle andclassification of filter. In third, understanding filter based on the principle, accordingto different digital filter by using the appropriate design method, based on MATLABspatial filtering of the mean, median, Wien

11、er filter to design the content needed inorder to improve the image quality. Extracted from image information effectively,must the image denoising pretreatment. Fourth, design Matlab GUI spatial filteringsoftware, using Matlab GUI interface design of a spatial enhancement filter, the filterhas a thr

12、ee filtering methods, respectively mean filtering, median filtering and Wienerfiltering.Finally, according to the different denoising algorithm ( neighborhood average,median filtering, Wiener filtering method), using the software of Matlab code, to asheet containing noise images (including Gauss noi

13、se and salt and pepper noisedenoising) simulation,and the results analyzed,comparedthe advantagesanddisadvantages of several methods, and designs corresponding airspace enhancementsoftware.Through the Matlab GUI interface software design, can be easily carried outinspatial domain image enhancement p

14、rocessing, even for Matlab software unskilledperson can also be more convenient and rapid operation. Just click interface, loadimages, corresponding to image processing.分享到翻译结果 HYPERLINK http:/ l# 重试抱歉，系统响应超时，请稍后再试支持中英、中日在线互译支持网页翻译，在输入框输入网页地址即可提供一键清空、复制功能、支持双语对照查看，使您体验更加流畅Keywords:Median Filtering;M

15、ean Filtering;Wiener Filtering;WaveletTransform;Matlab GUI河北经贸大学毕业论文第1章 绪论1.1 课题研究背景21 世纪，人类已经进入了信息化时代，计算机在处理各种信息中发挥着重要作用。据统计，人类从自然界获取的信息中，视觉信息占75%85%。俗话说“百闻不如一见 ”，有些场景或事物，不管花费多少笔墨都难以表达清楚，然而，若用一幅图像描述，可以做到一目了然。可见，在当代高度信息化的社会中，图形和图像在信息传播中所起的作用越来越大，在图像处理领域，数字图像处理得到了飞速发展。图像是信息社会人们获取信息的重要来源之一。在通过图像传感器将现实

16、世界中的有用图像信号进行采集、量化、编码、传输、恢复的过程中，存在大量影响图像质量的因素。因此图像在进行使用之前，一般都要经过严格的预处理如去噪、量化、压缩编码等。噪声的污染直接影响着对图像边缘检测、特征提取、图像分割、模式识别等处理，使人们不得不从各种角度进行探索以提高图像的质量。所以采用适当的方法尽量消除噪声是图像处理中一个非常重要的预处理步骤。图像处理技术在 20 世纪首先应用于图像的远距离传送，而改善图像质量的应用开始于 1964 年美国喷气动力实验室用计算机对 “徘徊者七号 ”太空船发回的月球照片进行处理，并获得巨大成功。现在图像处理技术已深入到科学研究、军事技术、工农业生产、医学、

17、气象及天文学等领域。科学家利用人造卫星可以获得地球资河北经贸大学毕业论文源照片、气象情况；医生可以通过X 射线或 CT 对人体各部位的断层图像进行分析。但在许多情况下图像信息会受到各种各样噪声的影响，严重时会影响图像中的有用信息，所以对图像的噪声处理就显得十分重要。1.2 图像去噪的研究现状图像处理中，输入的是质量低的图像，输出的是改善质量后的图像。常用的图像处理方法有图像增强、复原、编码、压缩等。光学相干层析成像是近些年来发展较快的一种层析成像技术。因为其对生物组织无辐射损伤、具有微米级的分辨率、高探测灵敏度和越来越快的扫描速率等优点，在医学诊断病变组织方面，尤其是对生物组织活体检测具有诱人

18、的应用前景。由于噪声干扰的存在，这些生理信号可能失真。甚至面目全非，这给图像处理带来了难度。因此，必须对含噪的图像进行处理，从而改善图像质量，最大程度地显现图像本来的特点。在所处理的图像中，相邻像素的灰度之间大多具有很高的相关性， 也就是说，一幅图像中大多数像素的灰度差别不大。因为这种灰度相关性的存在，一般图像的能量主要集中在低频区域中，只有图像的细节部分的能量才处于高频区域中。因为在图像的数字化和传输中常有噪声出现，而这部分干扰信息主要集中在高频区域内，所以消除噪声的一般方法是衰减高频分量或称低通滤波，但与之同时带来的负面影响是图像的细节也有一定的衰减，从视觉效果上来看图像比处理前模糊。一个

19、较好的去噪方法应该是既能消去噪声对图像的影响又不使图像细节变模糊。为了改善图像质量，从图像中提取有效信息，必须对图像进行去噪预处理。河北经贸大学毕业论文根据噪声的频谱分布的规律和统计特征以及图像的特点，出现了多种多样的去噪方法。经典去噪方法有：空间域合成法、频域合成法和最优线性合成法等。与之相适应的出现了许多应用方法：如均值滤波器、中值滤波器、低通滤波器、维纳滤波器、最小失真法等。这些方法广泛应用，促进数字信号处理的极大发展，显著提高了图像处理的质量。近年来，小波变换去除噪声的方法得到广泛的应。与传统的去噪方法相比，它利用非线性域值，在时间域和频率域同时具有良好的局部化性质，而且时窗和频窗的宽

20、度可以调节。对高频成分采用逐渐精细的时域或空域取样步长，从而可以聚焦到对象的任意细节。因此可以提高散斑高散射特性的噪声对比度，很好的消除散斑噪。小波变换去除噪声的方法在不断地发展，去噪方法很多，如非线性小波变换阈值法去噪、小波变换模极大值去噪及基于小波变换域的尺度相关性去噪法等。1.3 本文主要工作图像在获取和传输过程中，往往受到噪声的干扰，而降噪的目的是尽可能保持原始信号主要特征的同时，除去信号中的噪声。目前的图像去噪方法可以将图像的高频成分滤除，虽然能够达到降低噪声的效果，但同时破坏了图像细节。边缘特性是图像最为有用的细节信息，本文对邻域平均法、中值滤波法、维纳滤波法及模糊小波变换法的图像

21、去噪算法进行了研究分析和讨论。第 1 章 绪论：主要介绍本课题的研究背景以及国内外的研究现状。第 2 章 Matlab 简介：介绍软件 Matlab 的概况、发展历程、语言特点、图像处理应用时的常用函数及其用法。河北经贸大学毕业论文第 3 章 图像去噪算法：简要说明了图像噪声的概念及分类， 详细阐述了邻域平均法、中值滤波法、维纳滤波法及模糊小波变换法的去噪原理及特点。第 4 章 基于 Matlab GUI 的空域增强设计：根据邻域平均法、中值滤波法、维纳滤波法的原理，运用 Matlab 仿真软件编写代码，对一张含噪图片（含高斯噪声或椒盐噪声）进行仿真去噪，并对结果分析讨论，比较几种方法的优缺点

22、。最后设计了一个 Matlab GUI 空域增强软件，方便他人的操作。第2章 Matlab 简介2.1 Matlab 概况Matlab(Matrix Laboratory) 为美国 Mathworks 公司 1983 年首次推出的一套高性能的数值分析和计算软件，其功能不断扩充，版本不断升级，目前的最新版本河北经贸大学毕业论文为 7.1 版。Matlab 将矩阵运算、数值分析、图形处理、编程技术结合在一起，为用户提供了一个强有力的科学及工程问题的分析计算和程序设计工具，它还提供了专业水平的符号计算、文字处理、可视化建模仿真和实时控制等功能，是具有全部语言功能和特征的新一代软件开发平台。Matla

23、b 已发展成为适合众多学科，多种工作平台、功能强大的大型软件。在欧美等国家的高校， Matlab 已成为线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具。成为攻读学位的本科、硕士、博士生必须掌握的基本技能。 在设计研究单位和工业开发部门，Matlab 被广泛的应用于研究和解决各种具体问题。在中国， Matlab 也已日益受到重视，短时间内就将盛行起来， 因为无论哪个学科或工程领域都可以从Matlab 中找到合适的功能。当今的信息化社会，图像是人类赖以获取信息的最重要的来源之一。随着计算机技术的迅猛发展，图像技术与计算机技术不断融合，产生了一系列

24、图像处理软件， 如 VC、 Matlab， 这些软件的广泛应用为图像技术的发展提供了强大的支持。Matlab 已成为国际公认的最优秀的科技应用软件之一，具有编程简单、数据可视化功能强、可操作性强等特点，而且配有功能强大、专业函数丰富的图像处理工具箱，是进行图像处理方面工作必备的软件工具。2.1.1Matlab 发展过程Matlab 名字由 MATrix 和 LABoratory 两词的前三个字母组合而成。那是20世纪七十年代后期的事：时任美国新墨西哥大学计算机科学系主任的Cleve Moler教授出于减轻学生编程负担的动机，为学生设计了一组调用 LINPACK 和河北经贸大学毕业论文EISPA

25、CK 库程序的 “通俗易用 ”的接口，此即用FORTRAN 编写的萌芽状态的Matlab。经几年的校际流传，在 Little 的推动下，由 Little、Moler、Steve Bangert 合作，于 1984 年成立了 MathWorks 公司，并把 MATLAB 正式推向市场。从这时起，Matlab 的内核采用 C 语言编写，而且除原有的数值计算能力外，还新增了数据图视功能。Matlab 仅短短几年，就以其良好的开放性和运行的可靠性，使原先控制领域里的封闭式软件包。在时间进入20 世纪九十年代的时候， Matlab 已经成为国际控制界公认的标准计算软件。MathWorks 公司于 199

26、3 年推出 Matlab4.0 版本，从此告别 DOS 版。4.X 版在继承和发展其原有的数值计算和图形可视能力的同时， 出现了以下几个重要变化：（1）推出了 SIMULINK。这是一个交互式操作的动态系统建模、仿真、分析集成环境。它的出现使人们有可能考虑许多以前不得不做简化假设的非线性因素、随机因素，从而大大提高了人们对非线性、随机动态系统的认知能力。（2）开发了与外部进行直接数据交换的组件，打通了Matlab 进行实时数据分析、处理和硬件开发的道路。（3）推出了符号计算工具包。 1993 年 MathWorks 公司从加拿大滑铁卢大学购得 Maple 的使用权，以 Maple 为“引擎”开

27、发了 Symbolic Math Toolbox 1.0 。（4）构作了 Notebook 。MathWorks 公司瞄准应用范围最广的Word，运用DDE 和 OLE，实现了 MATLAB 与 Word 的无缝连接，从而为专业科技工作者创造了融科学计算、图形可视、文字处理于一体的高水准环境。河北经贸大学毕业论文1997 年仲春，Matlab5.0 版问世，现今的 Matlab 拥有更丰富的数据类型和结构、更友善的面向对象、更加快速精良的图形可视、更广博的数学和数据分析资源、更多的应用开发工具。诚然，到 1999 年底，Mathematica 也已经升到 4.0 版，它特别加强了以前欠缺的大规

28、模数据处理能力。在国际学术界， Matlab 已经被确认为准确、可靠的科学计算标准软件。在许多国际一流学术刊物上，（尤其是信息科学刊物），都可以看到 Matlab 的应用。在设计研究单位和工业部门， Matlab 被认作进行高效研究、开发的首选软件工具。2.1.2Matlab 的语言特点Matlab 语言有如下特点：（1）编程效率高（2）用户使用方便（3）扩充能力强（4）语句简单，内涵丰富（5）高效方便的矩阵和数组运算（6）方便的绘图功能由于 Matlab 的上述特点， 该软件已经广泛地应用于各行各业， 例如数值分析、图像处理、模拟仿真等等。2.2 Matlab 图像处理常用函数Matlab7

29、.0 为用户在进行图像处理中提供了一些常用的函数：（1）文件的读入与显示河北经贸大学毕业论文函数 imread 作用是读入文件，其调用格式如下：X,map=imread(filename,fmt) ， 其中，filename 为需要读入图像的文件名。 fmt为图像格式。函数 imshow 作用是显示文件，其语法格式如下：imshow(BW)imshow(X,map)（2）计算二维卷积函数 conv2 格式：C=conv2(A,B) 作用是算矩阵 A 和 B 的卷积。（3）噪声及其噪声的Matlab 实现函数 imnoise 格式：J=imnoise(I,type)J=imnoise(I,typ

30、e,parameter)返回对图像 I 添加典型噪声后的有噪图像 J， 参数type 和 parameter 用于确定噪声的类型和相应的参数。（4）二维离散小波变换函数 dwt2 格式：cA,cH,cV,cD=dwt2(X,wname)是使用指定的小波基函数“wname” 对二维信号 X 进行二维离散小波变换， cA，cH，cV，cD 分别为近似分量、水平细节分量、垂直细节分量和对角细节分量。cA,cH,cV,cD=dwt2(X,Lo_D,Hi_D)是使用指定的分解低通和高通滤波器Lo_D 和 Hi_D 分解信号 X。河北经贸大学毕业论文（5）二维信号的多层小波分解函数 wavedec2 格式

31、：C,S=wavedec2(X,N,wname)使用小波基函数“wname” 对二维信号 X 进行N 层分解。C,S=wavedec2(X,N,Lo_D,Hi_D)使用指定的分解低通和高通滤波器Lo_D和 Hi_D 分解信号 X。（6）二维离散小波反变换函数 idwt2 格式：X=idwt2(cA,cH,cV,cD,wname)由信号小波分解的近似信号cA 和细节信号cH、cH、cV、cD 经小波反变换重构原信号 X。X=idwt2(cA,cH,cV,cD,Lo_R,Hi_R)使用指定的重构低通和高通滤波器Lo_R和 Hi_R 重构原信号 X。X=idwt2(cA,cH,cV,cD,wname

32、,S)返回中心附近的S 个数据点。（7）二维信号的多层小波重构函数 waverec2 格式：X=waverec2(C,S,wname)由多层二维小波分解的结果 C、 S 重构原始信号 X。X=waverec2(C,S,Lo_R,Hi_R)使用重构低通和高通滤波器 Lo_R 和 Hi_R 重构原信号。河北经贸大学毕业论文第3章 图像去噪算法3.1 图像噪声概述3.1.1图像噪声的概念噪声可以理解为 “妨碍人们感觉器官对所接收的信源信息理解的因素”。例如一幅黑白图片，其平面亮度分布假定为f (x, y)， 那么对其接收起干扰作用的亮度分布R(x, y)即可称为图像噪声。但是，噪声在理论上可以定义为

33、“不可预测，只能用概率统计方法来认识的随机误差 ”。 因此将图像噪声看成是多维随机过程是合适的，因而描述噪声的方法完全可以借用随机过程的描述，即用其概率分布函数和概率密度分布函数。但在很多情况下，这样的描述方法是很复杂的，甚至是不可能的。而实际应用往往也不必要。通常是用其数字特征，即均值方差，相关函数等。因为这些数字特征都可以从某些方面反映出噪声的特征。目前大多数数字图像系统中，输入图像都是采用先冻结再扫描方式将多维图河北经贸大学毕业论文像变成一维电信号，再对其进行处理、存储、传输等加工变换。最后往往还要在组成多维图像信号，而图像噪声也将同样受到这样的分解和合成。在这些过程中电气系统和外界影响

34、将使得图像噪声的精确分析变得十分复杂。另一方面图像只是传输视觉信息的媒介，对图像信息的认识理解是由人的视觉系统所决定的。不同的图像噪声，人的感觉程度是不同的，这就是所谓人的噪声视觉特性课题。图像噪声在数字图像处理技术中的重要性越来越明显，如高放大倍数航片的判读，X 射线图像系统中的噪声去除等已经成为不可缺少的技术步骤。3.1.2图像噪声的分类图像噪声按其产生的原因可以分为：外部噪声， 即指系统外部干扰以电磁波或经电源串进系统内部而引起的噪声。如电气设备，天体放电现象等引起的噪声。内部噪声：一般又可分为以下四种：（1） 由光和电的基本性质所引起的噪声。 如电流的产生是由电子或空穴粒子的集合，定向

35、运动所形成。因这些粒子运动的随机性而形成的散粒噪声；导体中自由电子的无规则热运动所形成的热噪声；根据光的粒子性，图像是由光量子所传输，而光量子密度随时间和空间变化所形成的光量子噪声等。（2）电器的机械运动产生的噪声。如各种接头因抖动引起电流变化所产生的噪声；磁头、磁带等抖动或一起的抖动等。（3）器材材料本身引起的噪声。如正片和负片的表面颗粒性和磁带磁盘表面缺陷所产生的噪声。随着材料科学的发展，这些噪声有望不断减少，但在目前来讲，还是不可避免的。河北经贸大学毕业论文（4）系统内部设备电路所引起的噪声。如电源引入的交流噪声；偏转系统和箝位电路所引起的噪声等。图像噪声从统计理论观点可以分为平稳和非平

36、稳噪声两种。在实际应用中，不去追究严格的数学定义，这两种噪声可以理解为：其统计特性不随时间变化的噪声称其为平稳噪声。其统计特性随时间变化而变化的称其为非平稳噪声。3.2邻域平均法图像去噪3.2.1模板操作和卷积运算模板操作是数字图像处理中常用的一种运算方式， 图像的平滑、锐化、细化、边缘检测等都要用到模板操作。例如，有一种常见的平滑算法是将原图中的一个像素的灰度值和它周围邻近 8 个像素的灰度值相加，然后将求得的平均值作为新图像中该像素的灰度值。可用如下方法来表示该操作：1111（3-1）1 1* 19111上式有点类似矩阵，通常称之为模板（ Template） ，带星号的数据表示该元素为中心

37、元素，即这个元素是将要处理的元素。如果模板为：1* 1 1111 1（3-2）911 1该操作的含义是：将原图中的一个像素的灰度值和它右下相邻的8 个像素值相加，然后将求得的平均值作为新图像中该像素的灰度值。河北经贸大学毕业论文图 3-1 卷积处理过程模板操作实现了一种邻域运算，即某个像素点的结果不仅和本像素灰度有关，而且和其邻域点的值有关。模板运算的数学含义是卷积（或互相关）运算。卷积是一种用途很广的算法，可用卷积来完成各种变换，图3-1 说明了卷积的处理过程。卷积运算中的卷积核就是模板运算中的模板，卷积就是做加权求和的过程。邻域中的每个像素（假定邻域为33 大小，卷积核大小与邻域相同），分

38、别和卷积核中的每一个元素相乘，乘积求和所得的结果即为中心像素的新值。卷积核中的元素称作加权系数 （亦称为卷积系数），卷积核中的系数大小及排列顺序， 决定了对图像进行区处理的类型。改变卷积核中的加权系数，会影响到总和的数值与符号，从而影响到所求像素的新值。河北经贸大学毕业论文在模板或卷积的加权运算中，还存在一些具体问题需要解决：首先是图像边界问题，当在图像上移动模板（卷积核）至图像边界时，在原图像中找不到与卷积核中的加权系数相对应的 9 个像素，即卷积核悬挂在图像缓冲区的边界上，这种现象在图像的上下左右四个边界上均会出现。例如，当模板为1111（3-3）1 1*19111设原图像为1111122

39、222（3-4）3333344444经过模板操作后的图像为222（3-5）333“”表示无法进行模板操作的像素点。解决这个问题可以采用两种简单的方法：一种方法是忽略图像边界的数据，另一种方法是在图像四周复制原图像边界像素值，从而使卷积核悬挂在图像四周河北经贸大学毕业论文时可以进行正常的计算。实际应用中，多采用第一种方法。3.2.2邻域平均法原理邻域平均法是一种利用 Box 模板对图像进行模板操作 （卷积操作）的图像平滑方法， 所谓 Box 模板是指模板中所有系数都取相同值的模板， 常用的 33 和 55模板如下：1 111 11 111 1111111 1*11 1 1* 1 1（3-6）25

40、91 111 11111 111 1中间的黑点表示以该像素为中心元素，即该像素是要进行处理的像素。Box模板对当前像素及其相邻的像素点都一视同仁，统一进行平均处理，这样就可以滤去图像中的噪声。例如，用33Box 模板对一幅数字图像处理结果，如图3-2所示（图中计算结果按四舍五入进行了调整，对边界像素不进行处理） 。图 3-2 33Box 模板平滑处理示意图最简单的平滑滤波是将原图中一个像素的灰度值和它周围邻近像素的灰度值相加，然后将求得的平均值作为新图中该像素的灰度值。它采用模板计算的思河北经贸大学毕业论文想，模板操作实现了一种邻域运算， 即某个像素点的结果不仅与本像素灰度有关，而且与其邻域点

41、的像素值有关。设f (i, j)为给定的含有噪声的图像，经过邻域平均处理后的图像为g(x, y)，则邻域平均法也可以用数学公式表达：1（3-7）g ( x , y ) 式中：x， y=0，1，2，N-1；S 是以(x, y)为中心的邻域的集合， M 是S 内的点数。邻域平均法的思想是通过一点和邻域内的像素点求平均来去除突变的像素点，从而滤掉一定的噪声，其主要在实际应用中，也可以根据不同的需要选择M( i , j ) sf (i, j )使用不同的模板尺寸，如33、55、77、99 等。邻域平均处理方法是以图像模糊为代价来减小噪声的，且模板尺寸越大，噪声减小的效果越显著。如果f (i, j)是噪

42、声点，其邻近像素灰度与之相差很大，采用邻域平均法就是用邻近像素的平均值来代替它，这样能明显消弱噪声点，使邻域中灰度接近均匀，起到平滑灰度的作用。因此，邻域平均法具有良好的噪声平滑效果， 是最简单的一种平滑方法。3.3 中值滤波法图像去噪中值滤波是一种非线性信号处理方法， 与其对应的中值滤波器也就是一种非线性滤波器。中值滤波器于 1971 提出并应用在一维信号时间序列分析中， 后来被二维图像信号处理技术所引用。它在一定条件下可以克服线性滤波器（如邻域平滑滤波等）所带来的图像细节模糊，而且对滤除脉冲干扰及图像扫描噪声最为有效。在实际运算过程中并不需要图像的统计特性，这也带来不少方便。但是对一河北经

43、贸大学毕业论文些细节多，特别是点、线、尖顶细节多的图像不宜采用中值滤波。由于中值滤波是一种非线性运算，对随机输入信号的严格数学分析比较复杂，下面采用直观的方法简要介绍中值滤波的原理。3.3.1中值滤波原理中值滤波就是用一个奇数点的移动窗口， 将窗口中心点的值用窗口内个点的中值代替。假设窗口内有五点，其值为80、90、200、110、120，那么此窗口内各点的中值即为 110。设有一个一维序列f1, f2, fn，取窗口长度（点数）为 m（m 为奇数） ，对其进行中值滤波，就是从输入序列中相继抽出 m 个数fiv, fi1, fi, fi1, fiv（其中fi为窗口的中心点值，v m 1/ 2）

44、 ，再将这 m 个点按其数值大小排序，取其序号为中心点的那个数作为滤波输出。用数学公式表示为： vii vii（3-8）yMedf, , f , fm 12i N , v 例如，有一序列 0，0，3，4，7，则 Med0，0，3，4，7=3。此列若用平滑滤波，窗口也是取 5，那么平滑滤波输出为 (0+3+4+0+7)/5=2.8 。图 3-3 使用内含 5 个像素的窗口对离散阶跃函数、斜坡函数、脉冲函数以及三角函数进行中值滤波和平均值滤波的示例。左边一列为原波形，中间为平均值滤波结果，右边为中值滤波结果。可以看出，中值滤波器不影响阶跃函数和斜坡函数。周期小于 m/2（窗口之半）的脉冲受到抑制，

45、另外三角函数的顶部变平。河北经贸大学毕业论文图 3-3 中值滤波和平均值滤波比较a)阶跃；b)斜坡；c)单脉冲；d)双脉冲； e)三脉冲；f)三角波河北经贸大学毕业论文二维中值滤波可由下式表示：yij Med fij（3-9）式中：A 为窗口；fij为二维数据序列。二维中值滤波的窗口形状和尺寸对滤波效果影响较大，不同的图像内容和不同的应用要求，往往采用不同的窗口形状和尺寸。常用的二维中值滤波窗口有线状、方形、圆形、十字形以及圆环形等。窗口尺寸一般先用33，再取 55 逐渐增大，直到滤波效果满意为止。就一般经验来讲，对于有缓变的较长轮廓线物体的图像，采用方形或圆形窗口为宜。 对于包含有尖顶物体的

46、图像， 用十字形窗口，而窗口大小则以不超过图像，用十字形窗口，而窗口大小则以不超过图像中的最小有效物体的尺寸为宜。如果图像中点、线、尖角细节较多，则不宜采用中值滤波。3.3.2中值滤波主要特性（1）对某些输入信号中值滤波的不变性对某些特定的输入信号，如在窗口内单调增加或减少的序列，中值滤波输出信 号 仍 保 持 输 入 信 号 不 变 ， 即 ：fin fi fin或fin fi fin，则yifi。河北经贸大学毕业论文a)b)图 3-4 中值滤波不变性示例a)原始图像； b)中值滤波输出一维中值滤波这种不变性可以从图3-3 中 a)和 b)上看出来。二维中值滤波的不变性如图 3-4 所示。它

47、不但与输入信号有关，而且还与窗口形状有关。一般与窗口对顶角连线垂直的边缘线保持不变性。利用这个特点，可以使中值滤波既能去除图像中的噪声，又能保持图像中一些物体的边缘。对于一些周期性的数据序列，中值滤波对此序列保持不变性。例如，下列一维周期性的数序列fi,1,1,1,1,1,1,1,1,河北经贸大学毕业论文若设窗口长度为 9，则中值滤波对此序列保持不变性。对于二维周期序列不变性，如周期网状结构图案，分析起来就更复杂了，可以通过试验改变窗口形状和尺寸来获取。（2）中值滤波去噪声性能2med对于零均值正态分布的噪声输入，中值滤波输出的噪声方差近似为i214mf2(m)m 122（3-10）2med2

48、i式中：为输入噪声功率（方差），m为中值滤波窗口长度（点数），m为输入噪声均值，f (m)为输入噪声密度函数。20而均值滤波的输出噪声方差为0212im（3-11）比较两公式，可以看出，中值滤波的输出与输入噪声的密度分布有关。对随机噪声的抑制能力，中值滤波比均值滤波要差一些。但对脉冲干扰，特别是脉冲宽度小于 m/2、相距较远的窄脉冲干扰，中值滤波的效果较好。（3）中值滤波的频谱特性设 G 为输入信号频谱， F 为输出信号频谱，定义中值滤波的频率响应特性为H （3-12）试验表明，中值滤波频谱特性起伏不大，其均值比较平坦。可以认为信号经GF河北经贸大学毕业论文中值滤波后，频谱基本不变。这一特点对

49、设计和使用中值滤波器很有意义。3.3.3复合型中值滤波对一些内容复杂的图像，可以使用复合型中值滤波。如中值滤波线性组合、高阶中值滤波组合、加权中值滤波以及迭代中值滤波等。（1）中值滤波的线性组合是将几种窗口尺寸大小和形状不同的中值滤波器复合使用，只要个窗口都与中心对称，滤波输出可保持几个方向上的边缘跳变，而且跳变幅度可调节。其线性组合方程如下：k（3-13）AkYijmaxMed fij式中：Ak为窗口。（2）高阶中值滤波组合如下式所示：NYijakMed fijk 1Ak（3-14）式中：ak为不同中值滤波的系数，Ak为窗口。这种中值滤波可以使输入图像中任意方向的线条保持不变。例如，可选择图

50、3-5 中的 4 种线性窗口A1 A4，用上式组合的中值滤波，可以使输入图像中各种方向线条保持不变，而且又有一定的噪声平滑性能。河北经贸大学毕业论文图 3-5 几种线性窗口（3）其他类型的中值滤波：为了在一定的条件下尽可能去除噪声，又有效保持图像细节，可以对中值滤波器参数进行修正，如加权中值滤波，也就是对输入窗口进行加权。也可以是对中值滤波器的使用方法进行变化， 保证滤波的效果，还可以和其他滤波器联合使用。3.4维纳滤波法图像去噪3.4.1维纳滤波原理从连续的（或离散的）输入数据中滤除噪声和干扰以提取有用信息的过程称为滤波，这是信号处理中经常采用的主要方法之一，具有十分重要的应用价值，而相应的