阈值分割Word文档格式.docx

1、 关键词：图像分割 阈值选取 全局阈值 局部阈值 直方图 1 阈值分割 1.1 图像分割的概述 在对图像的研究和应用中，人们往往仅对图像中的某些部分感兴趣，这些部分称为目标或前景（其他部分称为背景），他们一般对应图像中特定的、具有独特性质的区域。为了辨识和分析目标，需要将他们分离提取出来，在此基础上才有可能对目标进一步利用。图像分割就是指把图像分成格局特性的区域并提取出感兴趣目标的技术和过程。这里特性可以是象素的灰度、颜色、纹理等，预先定义的目标可以对应单个区域，也可以对应多个区域。现有的图像分割算法有：阈值分割、边缘检测和区域提取法。本文着重研究基于阈值法的图像分割技术。 所谓图像分割是指根

2、据灰度、彩色、空间纹理、几何形状等特征把图像划分成若干个互不相交的区域，使得这些特征在同一区域内，表现出一致性或相似性，而在不同区域间表现出明显的不同。简单的讲，就是在一幅图像中，把目标从背景中分离出来，以便于进一步处理。图像分割是图像处理与计算机视觉领域低层次视觉中最为基础和重要的领域之一，它是对图像进行视觉分析和模式识别的基本前提。同时它也是一个经典难题，到目前为止既不存在一种通用的图像分割方法，也不存在一种判断是否分割成功的客观标准。 阈值法是一种传统的图像分割方法，因其实现简单、计算量小、性能较稳定而成为图像分割中最基本和应用最广泛的分割技术已被应用于很多的领域，例如，在红外技术应用中

3、，红外无损检测中红外热图像的分割，红外成像跟踪系统中目标的分割；在遥感应用中，合成孔径雷达图像中目标的分割等；在医学应用中，血液细胞图像的分割，磁共振图像的分割；在农业工程应用中，水果品质无损检测过程中水果图像与背景的分割。在工业生产中，机器视觉运用于产品质量检测等等。在这些应用中，分割是对图像进一步分析、识别的前提，分割的准确性将直接影响后续任务的有效性，其中阈值的选取是图像阈值分割方法中的关键技术。 1.2 阈值分割的基本原理 图像阈值化分割是一种最常用，同时也是最简单的图像分割方法，它特别适用于目标和背景占据不同灰度级范围的图像。它不仅可以极大的压缩数据量，而且也大大简化了分析和处理步骤

4、，因此在很多情况下，是进行图像分析、特征提取与模式识别之前的必要的图像预处理过程。图像阈值化的目的是要按照灰度级，对像素集合进行一个划分，得到的每个子集形成一个与现实景物相对应的区域，各个区域内部具有一致的属性，而相邻区域布局有这种一致属性。这样的划分可以通过从灰度级出发选取一个或多个阈值来实现。 阈值分割法是一种基于区域的图像分割技术，其基本原理是：通过设定不同的特征阈值，把图像像素点分为若干类常用的特征包括：直接来自原始图像的灰度或彩色特征；由原始灰度或彩色值变换得到的特征设原始图像为f（x，y），按照一定的准则在f（x，y）中找到特征值T，将图像分割为两个部分，分割后的图像为 : b0f

5、（x,y）0 2抖zz 和这些极小值点对应的灰度值就可用作分割阈值。 2.2 基于最优阈值的方法实现图像切割 有时目标和背景的灰度值有部分交错，用1个全局阈值并不能将它们绝对分开。这时常希望能减小误分割的概率，而选取最优阈值是一种常用的方法。设一幅图像仅包含2类主要的灰度值区域（目标和背景），它的直方图可看成灰度值概率密度函数p（z）的一个近似。这个密度函数实际上是目标和背景的2个单峰密度函数之和。如果已知密度函数的形式，那么就有可能选取1个最优阈值把图像分成2类区域而使误差最小。 设有这样1幅混有加性高斯噪声的图像，它的混合概率密度是 轾轾22（z-m）P（z-m）P犏犏1122p（z）=P

6、p（z）+Pp（z）=exp-+exp- 犏犏 1122222s2s2ps2ps臌臌1212mmss其中和分别是背景和目标区域的平均灰度值，和分别是关于均值1212PP的均方差，和分别是背景和目标区域灰度值的行验概率。根据概率定义有12PP+=1，所以混合概率密度中有5个未知的参数。如果能求出这些参数就可以12确定混合概率密度。 mm假设，需定义1个阈值T使得灰度值小于T的像素分割为背景而使得12灰度值大于T的像素分割为目标。这时错误地将1个目标像素划分为背景的概率和1个背景像素错误地划分为目标的概率分别是 TE（T）=p（z）dz12- E（T）=p（z）dz21T总的误差概率是 E（T）=

7、PE（T）+PE（T）2112为求得使该误差最小的阈值可将E（T）对T求导并令导数为零，这样得到： Pp（T）=Pp（T）1122将这个结果用于高斯密度，可得到解一元二次方程的根判别式的系数 22A=s-s12（）22 B=2ms-ms1221nsp2222222=-+21Csmsm2ssXln（） 122112ispi=112该二次在一般情况下有2个解。如果2个区域的方差相等，刚只有1个最优阈值 2m+mps122T=+ln（） optimal2m-mp121 2.3 基于迭代的方法实现图像切割 迭代法是基于逼近的思想，其步骤如下： （1）求出图象的最大灰度值和最小灰度值，分别记为Zmax和

8、Zmin，令初始阈值T0=（Zmax+Zmin）/2； （2）根据阈值TK将图象分割为前景和背景，分别求出两者的平均灰度值ZO和ZB； （3）求出新阈值TK+1=（ZO+ZB）/2； （4）若TK=TK+1，则所得即为阈值；否则转2，迭代计算。 迭代所得的阈值分割的图象效果良好。基于迭代的阈值能区分出图像的前景和背景的主要区域所在，但在图像的细微处还没有很好的区分度。但令人惊讶的是，对某些特定图象，微小数据的变化却会引起分割效果的巨大改变，两者的数据只是稍有变化，但分割效果却反差极大 经试验比较，对于直方图双峰明显，谷底较深的图像，迭代方法可以较快地获得满意结果。但是对于直方图双峰不明显，或图

9、像目标和背景比例差异悬殊，迭代法所选取的阈值不如最大类间方差法。 2.4最大类间方差的方法实现图像切割 由Otsu于1978年提出的最大类间方差法以其计算简单、稳定有效，一直广为使用。从模式识别的角度看，最佳阈值应当产生最佳的目标类与背景类的分离性能，此性能我们用类别方差来表征，为此引入类内方差 、类间方差 和总体方差。最大类间方差法计算简单、稳定有效，一直广为使用，是一种受到普遍欢迎的阈值选取方法。其基本思路是将直方图在某一阈值处分割成两组，当被分成的两组的方差为最大时，得到阈值。因为方差是灰度分布均匀性的一种量度，方差值越大，说明构成图像的两部分差别越大，当部分目标错分为背景或部分背景错分

10、为目标都会导致两部分差别变小，因此使类间方差最大的分割意味着错分概率最小。 图像的灰度级范围是0，1，2，L-1，设灰度级i的象素点个数为mi，l-1图像的象素点的总数为，则灰度级i的出现概率定义为pM=mii0L-1mi在Ostu方法中，阈值t把图像的象素分为C0=（0，1，t）pi=,pi=1M0和C1=（t+1，t+2，L-1）两类（分别代表门标与背景）。 CC类出现概率及均值分别为:0和1tw= p=w（t） 0ii=0L-1w=P=1-w（t）1ii=t+1tipm（t）im= 0ww（t）i=00L-1ipm（t）-m（t）iTm=11-w（t）wi=t+11tL-1邋其中:m（t

11、）=ip,m（t）ip iTii=0i=t+1 CC类的方差:0和12t（i-m）p20is（t）=0wi=002L-1（i-m）p21is（t）=1wi=t+11222类间方差为:s（t）=ws+wsm0011222类内方差为:s（t）=w（m-m）+w（m-m）n00T11T222总体方差为:s（t）=s+sTmn2sm引入关于t的等价判决准则函数:h（t）=2sn*最优阈值t通过等价判决准则的最大值得到:*t=argmaxh（t） t=0,L-1即阈值T将图像分成目标，背景两部分，使得两类总方差取得最大值的t，即为最佳分割阈值 2.5基于最大熵的方法实现图像切割 八十年代以来，许多学者将

12、Shannon信息熵的概念应用于图像阈值化，其基本思想都是利用图像的灰度分布密度函数定义图像的信息熵，根据假设的不同或视角的不同提出不同的熵准则，最后通过优化该准则得到阈值。Pun16通过使后验熵的上限最大来确定阈值。此方法又称为KSW熵方法。 2.6 其它方法实现阈值分割 对于基于点的全局阈值选取方法，除上述主要几种之外还许多，但大多都是以上述基本方法为基础，做出的改进方法或者对算法的优化，如使用递推方法以降低算法复杂性。例如一种使目标和背景差距最大的阈值求取方法，类似于最大类间方差阈值法。是它的一种简化算法。又如1984年Dunn等人23提出了均匀化误差阈值选取方法，这种方法实质上是要使将

13、背景点误分为目标点的概率等于将目标点误分为背景点的概率。类似于最小误差阈值法。近年来有一些新的研究手段被引入到阈值选取中。比如人工智能，描述了如何用人工智能的方法，寻找直方图的谷底点，作为全局阈值分割。其它如神经网络，数学形态学，小波分析与变换等等。 总的来说，基于点的全局阈值算法，与其它几大类方法相比，算法时间复杂度较低，易于实现，适合应用于在线实时图像处理系统。基于区域的全局阈值选取方法对一幅图像而言，不同的区域，比如说目标区域或背景区域，同一区域内的象素，在位置和灰度级上同时具有较强的一致性和相关性。而在上述基于点的全局阈值选取方法中，有一个共同的弊病，那就是它们实际上只考虑了直方图提供

14、的灰度级信息，而忽略了图像的空间位置细节，其结果就是它们对于最佳阈值并不是反映在直方图的谷点的情况会束手无策，不幸我们通常遇到的很多图像恰恰是这种情况。另一方面，完全不同的两幅图片却可以有相同的直方图，所以即使对于峰谷明显的情况，这些方法也不能保证你得到合理的阈值。于是，人们又提出了很多基于空间信息的阈值化方法。可以说，局部区域的全局阈值选取方法，是基于点的方法，再加上考虑点领域内象素相关性质组合而成，所以某些方法常称为“二维xxx方法”。由于考虑了象素领域的相关性质，因此对噪声有一定抑止作用。 3 图像切割的实现 3.1 图像切割实现代码 clc clear all; %极小值图像切割% I

15、=imread（C:Documents and SettingsAdministrator桌面课程设计题目 2.bmp）; figure（1）; subplot（121）,imhist（I）; % 观察灰度直方图， 灰度80处有谷，确定阈值T=80 title（直方图 subplot（1,2,2）,imshow（I）; title（原图） I1=im2bw（I,80/255）; % im2bw函数需要将灰度值转换到0,1范围内 figure（2） subplot（3,2,1）,imshow（I1）;极小值点阈值切割 %迭代法图像切割% I=imread（ x,y=size（I）; a=imhi

16、st（I）; I=double（I）; max=1; for i=2:x if a（max）a（i） min=i; end end z0=max z1=min T=（z0+z1）/2; TT=0; S0=0; n0=0; S1=0; n1=0; allow=0.5; d=abs（T-TT）; count=0; while（d=allow） count=count+1; for i=1:x for j=1:y if （I（i,j）=T） S0=S0+I（i,j）; n0=n0+1; end if （I（i,j）=T） Seg（i,j）=1; end end end subplot（322）,im

17、show（Seg）;迭代阈值分割1 SI=1-Seg; subplot（323）,imshow（SI）;迭代阈值分割2 %最优法切割% I=imread（ m,n=size（I）; dt = 1; d = 0.3; I9 = I; I9=double（I9）; while （abs（dt-d）0.001） background = 0; object = 0; bm = 0;%背景象素数 om = 0;%前景象素数 for i=1:m for j=1:n if （I9（i,j）d） background = background +I9（i,j）; bm = bm+1; else object = object + I9（i,j）; om = om+1; end end end dt = d; d = （background/bm+object/om）/2; %假设先验概率相同 end tmax3=d %tmax3 =79.9102 for i = 1:m for j = 1:n if I9（i,j）tp） y1（i,j）=x1（i,j）; else y1（i,j）=0; end end end %上面一段代码实现分割 subplot（325）,imshow（y1）;最大方差法分割的图像 %最大熵法图像切割% I=imread（ X=double（I）;