《迭代阈值法》由会员分享，可在线阅读，更多相关《迭代阈值法（16页珍藏版）》请在金锄头文库上搜索。
1、摘要数字图像处理的目的之一是图像识别 , 而图像分割是图像识别工作的基础。 图像分割是指把图像分解成具有特性的区域并提取出感兴趣目标的技术和过程， 是计算机视觉领域的一个重要而且基本的问题，分割结果的好坏将直接影响到视 觉系统的性能。因此从原理，应用和应用效果的评估上深入研究图像分割技术具 有十分重要的意义。本课题主要介绍了图像分割的基本知识。图像分割的算法有阈值分割法，边 缘检测法，区域分割等，本设计重点介绍了基于最小点阈值方法，基于最优阈值 分割方法，基于迭代图像分割方法，最大类间方差法（OTSU）的图像分割法的 原理和他们的MATLAB的实现代码与运行结果。关键词：图像分割； MATLA
2、B； 阈值分割；目录1 课程设计目的 32 课程设计要求 33 相关知识 33.1 图像分割的概述33.2 阈值分割的基本原理43.3 阈值分割方法的分类 53.3.1 基于点的全局阈值方法 63.3.2 基于区域的全局阈值方法 63.3.3 局部阈值法和多阈值法64 程设计分析 64.1 基于迭代的方法实现图像切割64.2 最大类间方差的方法实现图像切割75 程序设计 75.1 程序简单介绍75.2 程序代码86 结果与分析 11结束语13参考文献14迭代阈值法1 课程设计目的本设计的课题任务是掌握图像阈值分割算法研究，实现对图像的分割。了解图像 分割的应用及基本方法，理解阈值化图像分割原理
3、，理解三类典型的阈值化分割 算法，并利用之进行图像分割，给出实验结果并做出分析。2 课程设计要求 查阅相关资料； 理解基于各像素值的阈值分割算法，基于区域性质的阈值分割算法 , 基于坐 标位置的阈值分割算法1；软件编程实现利用基于各像素值的阈值分割算法 进行图像分割，要求完成如下内容：包括极小值点阈值、最优阈值、迭代阈 值，基于最大方差的阈值，基于最大熵的阈值等方法，利用之实现图像分割， 这里的图像可以针对核磁共振图像（MRI） 2用MATLAB实现，并观察各算法之间的区别。3 相关知识3.1 图像分割的概述在对图像的研究和应用中，人们往往仅对图像中的某些部分感兴趣，这些 部分称为目标或前景（
4、其他部分称为背景），他们一般对应图像中特定的、具有独 特性质的区域。为了辨识和分析目标，需要将他们分离提取出来，在此基础上才 有可能对目标进一步利用。图像分割就是指把图像分成格局特性的区域并提取出 感兴趣目标的技术和过程。这里特性可以是象素的灰度、颜色、纹理等，预先定 义的目标可以对应单个区域，也可以对应多个区域。现有的图像分割算法有： 阈值分割、边缘检测和区域提取法。本文着重研究基于阈值法的图像分割技术。所谓图像分割是指根据灰度、彩色、空间纹理、几何形状等特征把图像划分 成若干个互不相交的区域,使得这些特征在同一区域内，表现出一致性或相似性， 而在不同区域间表现出明显的不同。简单的讲，就是在
5、一幅图像中，把目标从背 景中分离出来，以便于进一步处理。图像分割是图像处理与计算机视觉领域低层 次视觉中最为基础和重要的领域之一，它是对图像进行视觉分析和模式识别的基 本前提。同时它也是一个经典难题，到目前为止既不存在一种通用的图像分割方 法，也不存在一种判断是否分割成功的客观标准。阈值法是一种传统的图像分割方法，因其实现简单、计算量小、性能较稳定 而成为图像分割中最基本和应用最广泛的分割技术。已被应用于很多的领域， 例如，在红外技术应用中，红外无损检测中红外热图像的分割，红外成像跟踪系 统中目标的分割；在遥感应用中，合成孔径雷达图像中目标的分割等；在医学应 用中，血液细胞图像的分割，磁共振图
6、像的分害；在农业工程应用中，水果品 质无损检测过程中水果图像与背景的分割。在工业生产中，机器视觉运用于产品 质量检测等等。在这些应用中，分割是对图像进一步分析、识别的前提，分割的 准确性将直接影响后续任务的有效性，其中阈值的选取是图像阈值分割方法中的 关键技术。3.2 阈值分割的基本原理图像阈值化分割是一种最常用，同时也是最简单的图像分割方法，它特别适 用于目标和背景占据不同灰度级范围的图像。它不仅可以极大的压缩数据量，而 且也大大简化了分析和处理步骤，因此在很多情况下，是进行图像分析、特征提 取与模式识别之前的必要的图像预处理过程。图像阈值化的目的是要按照灰度级, 对像素集合进行一个划分，得
7、到的每个子集形成一个与现实景物相对应的区域， 各个区域内部具有一致的属性，而相邻区域布局有这种一致属性。这样的划分 可以通过从灰度级出发选取一个或多个阈值来实现。阈值分割法是一种基于区域的图像分割技术，其基本原理是：通过设定不同 的特征阈值，把图像像素点分为若干类.常用的特征包括：直接来自原始图像的 灰度或彩色特征；由原始灰度或彩色值变换得到的特征.设原始图像为f(x, y)， 按照一定的准则在f(x，y)中找到特征值T,将图像分割为两个部分，分割后的 图像为：g(x, y)=乃o, f(x, y) t若取b0=0 (黑)，久=1 (白)，即为我们通常所说的图像二值化。一般意义下，阈值运算可以
8、看作是对图像中某点的灰度、该点的某种局部特性以及该点在图像中的位置的一种函数，这种阈值函数可记作T = T(x, y, n(x, y), f(x, y)(322)式中f(x，y)是点(x，y)的灰度值；n(x，y)是点(x，y)的局部邻域特性.根据 对T的不同约束，可以得到3种不同类型的阈值，即(1) 点相关的全局阈值T = T(f(x,y)(3.2.3)只与点的灰度值有关(2) 区域相关的全局阈值T = T(n(x, y),f(x,y)(3.2.4)与点的灰度值和该点的局部邻域特征有关(3) 局部阈值或动态阈值T = T(x,yn(x, y),f(x, y)(3.2.5)它与点的位置、该点的
9、灰度值和该点邻域特征有关因此本文分三大类对阈值选取技术进行综述：1) 基于点的全局阈值方法；2) 基于区域的全局阈值方法3) 局部阈值方法和多阈值方法3.3 阈值分割方法的分类全局阈值法指利用全局信息对整幅图像求出最优分割阈值，可以是单阈值， 也可以是多阈值；局部阈值法是把原始的幅图像分为几个小的子图像，再对每个 子图像应用全局阈值法分别求出最优分割值。其中全局阈值法又可分为基于点的 阈值法和基于区域的阈值法。阈值分割法的结果很大程度上依赖于阈值的选择， 因此该方法的关键是如何选择合适的阈值。由于局部阈值法中仍要用到全局阈值法，因此本文主要对全局阈值法中基于点的阈值法和基于区域的阈值法分别进行
10、 了研究。根据阈值法的原理可以将阈值选取技术分为三大类3.3.1 基于点的全局阈值方法基于点的全局阈值算法与其他几大类方法相比，算法时间复杂度较低，易于 实现，适合应用于在线实时图像处理系统。3.3.2 基于区域的全局阈值方法对一幅图像而言，不同的区域，比如说目标区域或背景区域，同一区域内的 象素，在位置和灰度级上同时具有较强的一致性和相关性。3.3.3 局部阈值法和多阈值法局部阈值（动态阈值）当图像中有如下一些情况：有阴影，照度不均匀，各 处的对比度不同，突发噪声，背景灰度变化等，如果只用一个固定的全局阈值对 整幅图像进行分割，则由于不能兼顾图像各处的情况而使分割效果受到影响。有 一种解决办
11、法就是用与象索位置相关的一组阈值（即阈值使坐标的函数）来对图 像各部分分别进行分割。这种与坐标相关的阈值也叫动态阈值，此方法也叫变 化阈值法，或自适应阈值法。这类算法的时间复杂性和空间复杂性比较大，但是 抗噪能力强，对一些用全局阈值不易分割的图像有较好的效果。多阈值法很显然，如果图像中含有占据不同灰度级区域的几个目标，则需要 使用多个阈值才能将他们分开。其实多域值分割，可以看作单阈值分割的推广。4程序设计分析4.1 基于迭代的方法实现图像切割迭代法是基于逼近的思想，其步骤如下：（1）求出图象的最大灰度值和最小灰度值，分别记为z护“和zm.，令初始阈值ma天M/cT0 =（Z叱+ Z伽）72（4
12、丄1）（2）根据阈值兀将图象分割为前景和背景，分别求出两者的平均灰度值Z。和Z；(3) 求出新阈值:Tki =(Z Z严22(4丄2)(4) 若Tk=Tk1，则所得即为阈值；否则转到第二步，迭代计算。 迭代所得的阈值分割的图象效果良好。基于迭代的阈值能区分出图像的前景和背 景的主要区域所在，但在图像的细微处还没有很好的区分度。 但令人惊讶的是，对某些特定图象，微小数据的变化却会引起分割效果的巨大改 变，两者的数据只是稍有变化，但分割效果却反差极大经试验比较，对于直方图双峰明显，谷底较深的图像，迭代方法可以较快地 获得满意结果。但是对于直方图双峰不明显，或图像目标和背景比例差异悬殊， 迭代法所选
13、取的阈值不如最大类间方差法。4.2 最大类间方差的方法实现图像切割由Otsu于1978年提出的最大类间方差法以其计算简单、稳定有效，一直广 为使用。从模式识别的角度看，最佳阈值应当产生最佳的目标类与背景类的分 离性能，此性能我们用类别方差来表征，为此引入类内方差、类间方差和总体 方差。最大类间方差法计算简单、稳定有效，一直广为使用，是一种受到普遍欢 迎的阈值选取方法。其基本思路是将直方图在某一阈值处分割成两组，当被分 成的两组的方差为最大时，得到阈值。因为方差是灰度分布均匀性的一种量度， 方差值越大，说明构成图像的两部分差别越大，当部分目标错分为背景或部分背 景错分为目标都会导致两部分差别变小
14、，因此使类间方差最大的分割意味着错分 概率最小。5程序设计5.1 程序简单介绍程序除了有题目中要求的自编程序实现最大类间方差法外，还扩展了迭代法 阈值分割，功能较多，为了使程序有良好的人机界面，主程序中提供了简单的菜 单界面。同时为了增加程序的可读性，分模块编写，然后在主函数中调用。5.2 程序代码function main %主程序clear;close all; %清除变量，关闭所有窗口name path=uigetfile(C:Documents and SettingsAll UsersDocumentsMy Pictures示例图片Winter.jpg);%获取图片路径及文件名=im
15、read(path name);% 读图I=rgb2gray(I);%彩色转黑白while 1%循环执行disp(0 原图 1 OTSU 2 全局阈值 3 迭代法 4 退出);%提示信息select=input(请选择：)；选择相应功能switch selectcase 0yuantu(I); %输入 0 显示原图case 1%输入 1 OTSU 法fun1(I);case 2%输入 2 全局阈值法fun2(I);case 3fun3(I) %输入 3 迭代法case 4return; %输入 4 退出otherwise%其他值不处理endend% OTSUfunction fun1(I) %阈值计算程序Ni=imhist(I);%计算直方图数组N=sum(Ni);%总像素点个数delamax=0;%类间方差最大值threshold=0; %阈值for k=2:255 u=dot(0:255,Ni/N); %图像的总平均灰度级 w0=sum(Ni(1:k)/N); %C0 类像素所占面积的比例 w1=1-w0;%C1 类像素所占面积的
