图像阈值分割及去噪的实现毕业论文目录摘要 (1)Abstract (2)目录 (3)引言 (4)第一章图像噪音 (5)第二章图像缩放和灰度变换处理 (6)2.1图像缩放处理方法 (6)2.2图像灰度变换处理 (6)第三章图像阈值分割 (8)3.1 图像分割技术概要 (8)3.2图像阈值分割原理 (8)3.3图像阈值分割方法 (9)第四章图像去噪 (12)4.1 滤波原理 (12)4.2滤波实现方法 (12)第五章仿真实验结果和讨论 (16)5.1图像二值化算法对比 (16)5.2图像去噪效果对比 (17)结论 (21)参考文献 (22)致谢语 (23)引言数字图像处理是从 20 世纪 60 年代以来随着计算机技术和 VLSI 的发展而产生、发展和不断成熟起来的一个新兴技术领域，它在理论上和实际应用上都取得了巨大的成就，并引起各方面人士的广泛重视[1]。

首先，视觉是人类最重要的感知手段，图像又是视觉的基础。

因此数字图像成为心理学、生物医学、计算机科学等诸多方面的学者研究视觉感知的有效工具。

其次，数字图像处理在军事、遥感、工业图像处理等大型应用中也有不断增长的需求。

为适用特殊的场合和获得较好的视觉效果，常常需要一种有效的方法来对图像进行处理。

数字图像处理技术从广义上可看作是各种图像加工技术的总称。

它包括利用计算机和其他电子设备完成的一系列工作，如图像分割、图像变换、图像去噪等。

本文主要是在整合各种优秀的阈值分割和滤波算法的基础上，实现对图像进行分割和去噪，达到处理和读取图像的目的。

在MATLAB仿真的基础上，比对各种分割和去噪方法的优缺点。

第一章图像噪声大量的实验研究发现，由摄像机拍摄得到的图像受离散的脉冲、椒盐噪声和零均值的高斯噪声的影响较严重。

噪声给图像处理带来很多困难，对图像分割、特征提取、图像识别等具有直接影响。

因此，实时采集的图像需进行滤波处理。

消除图像中的噪声成份叫做图像的平滑化或滤波操作。

按信号与噪声的关系，噪声主要分为加性噪声和乘性噪声[2]：1)加性噪声，即与信号存在与否无关，是独立于信号之外的噪声，叠加方式为[3]：()()()= + (1.1)x t s t n t其中()表示噪声信号。

比如说放大器噪声。

n tx t表示输出信号，()s t表示源信号，()2)乘积性噪声又称为卷积噪声，只有在信号出现在信道中才表现出来，它不会主动对信号形成干扰。

叠加方式为[4]=+(1.2)x t s t n t()()[1()]乘积性噪声可以通过同态变换成为加性噪声。

比如说胶片颗粒噪声。

本文讨论通过对图像施加高斯噪声、椒盐噪声和乘性噪声，来模拟实际应用中图像损毁的部分情况。

其中，高斯噪声和椒盐噪声属于加性噪声。

第二章图像缩放和灰度变换处理2.1 图像缩放处理方法在计算机图像处理中，图像缩放（image scaling）是指对数字图像的大小进行调整的过程。

图像放大是从一幅低分辨率图像获得其高分辨率版本的一种图像处理技术[5]。

各种插值技术是实现图像放大的最常用的方法，其实质是对源图像建立连续的数学模型，然后按缩放要求进行重采样得到目的图像，如最近邻插值、双线性插值、双三次插值等。

最近邻插值是最简便的插值，在这种算法中，每一个插值输出像素的值就是在输入图像中与其最临近的采样点的值，运算量非常小。

当图像中包含像素之间灰度级变化的细微结构时，最近邻插值法会在图像中产生人工的痕迹。

具体原理为如下[5]：最近邻点插值的一个隐含的假设条件是任一网格点 p(x,y) 的属性值都使用距它最近的位置点的属性值，用每一个网格节点的最邻点值作为待定节点值。

当数据已经是均匀间隔分布，要先将数据转换为SURFER网格文件，可以应用最近邻点插值法；或者在一个文件中,数据紧密完整，只有少数点没有取值，可用最近邻点插值法来填充无值的数据点。

最近邻点插值网格化法没有选项,它是均质且无变化的,对均匀间隔的数据进行插值很有用,同时,它对填充无值数据的区域很有效2.2图像灰度变换处理图像的灰度变换(Gray-Scale Transformation GST) 处理是图像增强处理技术中一种非常基础、直接的空间域图像处理方法，也是图像数字化软件和图像显示软件的一个重要组成部分[6]。

灰度变换是指根据某种目标条件按一定变换关系逐点改变原图像中每一个像素灰度值的方法。

目的是为了改善画质，使图像的显示效果更加清晰。

灰度变换有时又被称为图像的对比度增强或对比度拉伸。

例如为了显示出图像的细节部分或提高图像的清晰度，需要将图像整个围的灰度级或其中某一段(,)a b 灰度级扩展或压缩到(,)a b ''，这些都要求采用灰度变换方法。

从图像输入装置得到的图像数据，以浓淡表示，各个像素与某一灰度值相对应。

设原图像像素的灰度值(,)D f x y =，处理后图像像素的灰度值(,)D g x y '=。

则灰度增强[5]可表示为：(,)[((,)] g x y T f x y = (2.2.1)或 () D T D ' = (2.2.2)要求D 和D ′都在图像的灰度围之。

函数T 称为灰度变换函数，它描述了输入灰度值和输出灰度值之间的转换关系。

根据(2.2.1)式可以将灰度变换分为线性变换和非线性变换。

本文采取rgb2gray 函数来实现灰度变换。

其原理如下：以R 、G 、B 为轴建立空间直角坐标系，则RGB 图的每个象素的颜色可以用该三维空间的一个点来表示，而Gray 图的每个象素的颜色可以用直线R=G=B 上的一个点来表示。

于是RGB 图转Gray 图的本质，就是寻找一个三维空间到一维空间的映射[6]，最直接的方法是射影（即过RGB 空间的一个点向直线R=G=B 做垂线），matlab 中也是如此实现的。

并且有0.29900 * 0.58700 * 0.11400 * Gray R G B =++(2.2.3)第三章图像阈值分割3.1 图像分割技术概要图像分割是一种重要的图像技术，在理论研究和实际应用中都得到了人们的广泛重视。

图像分割的方法和种类有很多，值得指出的是，没有惟一的标准的分割方法。

许多不同种类的图像或景物都可作为待分割的图像数据，不同类型的图像，已经有相对应的分割方法对其分割，同时，某些分割方法也只是适用于某些特殊类型的图像分割。

可以说，图像分割结果的好坏直接影响对图像的理解。

在图像研究中，图像分割的方法主要可以分成两大类[1]。

一类是边界方法，这种方法的假设是图像分割结果的某个子区域在原来图像中一定会有边缘存在；一类是区域方法，这种方法的假设是图像分割结果的某个子区域一定会有相同的性质，而不同区域的像素则没有共同的性质，阈值分割方法属于区域分割。

不可否认，这两类方法都有优缺点，限于时间和篇幅，本文将采用阈值法来分割图像。

从某种角度讲，图像分割的目的之一就是目标提取―将人类感兴趣的区域提取出来。

从这一点说，区域分割是最直接的方法，利用这种分割方法可以直接得到感兴趣的目标区域。

本文介绍其中的一种方法-----阈值分割方法。

阈值分割这种方法，作为一种图像预处理方式，图像分割并不考虑图象降质，只将图象中感兴趣的特征有选择地突出，衰减其不需要的特征，故预处理后的输出图象并不需要去逼近原图象。

3.2 图像阈值分割原理阈值分割的基本想法是确定一个阈值，然后把每个像素点的像素值和阈值相比较，根据比较的结果把该像素划分成两类——前景或者背景。

一般阈值分割可以分成以下3步[7]。

（1）确定阈值T。

（2）将阈值T和像素值比较。

（3）把像素归类。

上面的3个步骤，关键的是第一步，如果能确定一个合适的阈值，就可以对图像进行正确、方便的分割。

即：设定某一阈值T，用T将图像的数据分成两大部分：大于T的像素群和小于T 的像素群。

这是研究灰度变换最特殊的方法，称为图像的二值化（binarization）。

阈值处理的操作过程是先由用户指定或通过算法生成一个阈值，如果图像中某像素的灰度值小于该阈值，则将该像素的灰度值设置为0或255，否则灰度值设置为255或0。

阈值方法分全局阈值和局部阈值两种，如果分割过程中对图像上每个像素所使用的阈值都相等，则为全局阈值方法；如果每个像素所使用的阈值可能不同，则为局部阈值方法。

前者较简单，但抑制噪声和不均匀光照的能力较差[8]。

后者则是根据区域灰度分布特性，自适应调节阈值，效果一般好于前者。

可以认为局部阈值是全局阈值的一般情况，全局阈值则是最简单的图像分割方法。

本文中介绍了三种经典的图像区域阈值分割方法。

3.3 图像阈值分割方法一、Otsu算法Otsu是一种全局阈值算法，是在最小二乘法原理基础上推导出来的。

这种方法的基本思想是将直方图在某一阈值处分割成两组，一组对应于背景文字部分，当被分成的两组的组方差最小，组间方差最大时，决定阈值。

这种基于两组间最佳分类而决定阈值的方法也成为最大类间方差阈值分割法。

设给定图像具有L级灰度值，对1 < i < L 中的每个i 将1，2，……，L分成两组，计算落在组1的像素个数，平均灰度，方差；落在组2的像素个数，平均灰度，方差。

设给定图像的灰度级L ，灰度值为i 的像素数 i n ，总的像素数为 N ，每一个灰度值出现的概率为 i P ，则有i 0 = Li N n =∑(3.3.2) 0 Li i i N n P ==∑(3.3.3)对1 < i < L 中的每个 i 将1，2，… …，L 分成两组，以 (1)t t L <<为阈值，将像素分为两组，组1为灰度值小于 t 的像素集合。

分别计算落在组1和组2的像素数，平均灰度，方差。

设组1出现的概率为 1q ，平均灰度为 1 μ，方差为 21σ 。

组2像素灰度值的均值为2μ，方差为 22σ 。

2W σ 是组间方差，2B σ 是组间方差，则有1i 0 t iq P ==∑ (3.3.4)2221 W N q σσ= (3.3.5) 2221212()B q q N μμσ=-(3.3.6)对于一幅给定的图像可以证明22 W B C σσ+=(3.3.7)其中，C 为常数。

因此只需求出2max()B σ ，则 2W σ 自然达到最小[9]。

二、Niblack 算法 Niblack 二值化算法是局部阈值法，其基本思想是根据点( x , y)邻域 N ×N 块所有像素的灰度值来决定一个阈值。

若点( x , y)处的灰度值大于该阈值时被判为白色,否则为黑色。

点( x , y)处的阈值这样来计算[10]:(,) (,)(,) T x y m x y k s x y =+•(3.3.8)其中 m( x , y)和 s ( x , y)依次是点( x , y)处的局部邻域的样本平均值和标准差。