图像二值化算法设计研究和实现

专业资料专业资料WORDWORD完美格式 下载可编辑图像二值化算法研究与实现摘要：图像二值化是图像预处理中的一项重要技术，在模式识别、光学字符识别、医学成像等方面都有重要应用。论文介绍了图像及数字图像处理技术的一些概念和相关知识；对VC++软件的发展和软件在图像处理中的应用做了简要介绍；还介绍了图像二值化算法以及利用VC+软件工具进行算法的实现。论文重点实现了图像分割技术中常用灰度图像二值化算法，如Otsu算法、Bernsen算法，并对这些算法运行的实验结果进行分析与比较。关键词：图像处理；二值化；VC++.引言图像与数字图像图像就是用各种观测系统观测客观世界获得的且可以直接或间接作用与人眼而产生视觉的实体。视觉是人类从大自然中获取信息的最主要的手段。 拒统计，在人类获取的信息中，视觉信息约占60%听觉信息约占20%其他方式加起来才约占20%由此可见，视觉信息对人类非常重要。同时，图像又是人类获取视觉信息的主要途径，是人类能体验的最重要、最丰富、信息量最大的信息源。通常，客观事物在空间上都是三维的 （3D）的,但是从客观景物获得的图像却是属于二维（2D）平面的。数字图像：数字图像是将连续的模拟图像经过离散化处理后得到的计算机能够辨识的点阵图像。在严格意义上讲，数字图像是经过等距离矩形网格采样，对幅度进行等间隔量化的二维函数。因此，数字图像实际上就是被量化的二维采样数组。数字图像处理技术内容与发展现状数字图像处理就是采用一定的算法对数字图像进行处理，以获得人眼视觉或者某种接受系统所需要的图像处理过程。图像处理的基础是数字，主要任务是进行各种算法设计和算法实现。图像处理技术的发展大致经历了初创期、 发展期、普及期和实用化期4个阶段。初创期开始与20世纪60年代，当时的图像采用像素型光栅进行少秒显示，大多采用中、大型机对其处理。在这一时期，由于图像存储成本高、处理设备昂贵，其应用面很窄。进入20世纪70年代的发展期，开始大量采用中、小型机进行处理，图像处理也逐渐改用光栅扫描方式，特别是CT和卫星遥感图像的出现，对图像处理技术的发展起到了很好的推动作用。到了20世纪80年代，图像处理技术进入普及期，此时的微机已经能够担当起图形图像处理的任务。超大规模集成电路（VeryLargeScaleIntegration,VLSI）的出现更使处理速度大大提高，设备造价也进一步降低，极大地促进了图形图像系统的普及和应用。20世纪90年代是图像处理技术的实用化时期，图像处理的信息量巨大，对处理速度的要求极高。图像二值化原理及意义图像二值化是指用灰度变换来研究灰度图像的一种常用方法，即设定某一阈值将灰度图像的像素分成大于阈值的像素群和小于阈值的像素群两部分。例如输人灰度图像函数为f(x,y),输出二值图像函数为g(x,y),则0f(x,y)：：Thresholdg(x,y) (1)255f(x,y).Threshold阈值(threshold)是把目标和背景区分开的标尺，选取适当的阈值就是既要尽可能保存图像信息，又要尽可能减少背景和噪声的干扰，这是选择阈值的原则。图像的二值化处理就是将图像上的点的灰度置为 0或255,也就是使整个图像呈现出明显的黑白效果。即将256个亮度等级的灰度图像通过适当的阈值选取而获得仍然可以反映图像整体和局部特征的二值化图像。在数字图像处理中，二值图像占有非常重要的地位，特别是在实用的图像处理中，以二值图像处理实现而构成的系统是很多的， 要进行二值图像的处理与分析，首先要把灰度图像二值化，得到二值化图像，这样子有利于再对图像做进一步处理时， 图像的集合性质只与像素的值为0或255的点的位置有关，不再涉及像素的多级值，使处理变得简单，而且数据的处理和压缩量小。二值图像在图像分析中应用非常广泛，二值图像就是指只有两个灰度级的图像，二值图像具有存储空间小，处理速度快，可以方便地对图像进行布尔逻辑运算等特点。更重要的是，在二值图像的基础上，还可以进一步对图像处理，获得该图像的一些几何特征或者其他更多特征。.软件工具一VC++软件概述MicrosoftVisualC++HMicrosoft公司推出的开发Win32环境程序，面向对象的可视化集成编程系统。它不但具有程序框架自动生成、灵活方便的类管理、代码编写和界面设计集成交互操作、可开发多种程序等优点，而且通过简单的设置就可使其生成的程序框架支持数据库接口、 OLE2WinSock网络、控制界面。它以拥有“语法高亮”，IntelliSense (自动编译功能)以及高级除错功能而著称。比如，它允许用户进行远程调试，单步执行等。还有允许用户在调试期间重新编译被修改的代码，而不必重新启动正在调试的程序。其编译及建置系统以预编译头文件、最小重建功能及累加连结著称。这些特征明显缩短程式编辑、编译及连结的时间花费，在大型软件计划上尤其显著。VC++SC+钠言、MFC的封装类、IDE构成。MFC(microsoftfoundationclasses,微软基础类库)是微软为Windows程序员提供的一个面向对象的 Windows编程接口，它大大简化了 Windows编程工作。?MFC®供了一个标准化的结构?提供了大量的代码，指导用户编程时实现某些技术和功能。VC在单个工具中，能完成如下功能，因此被称为 IDE(IntegretedDevelopmentEnviroment,集成开发环境)：?生成不同编写代码的起始应用程序?通过许多不同方法查看项目文件?编辑源文件和包含文件?建立应用程序的可视界面?编译及链接?运行及调试应用程序VC++开发环境VisualC++的集成开发环境VisualStudio 提供了大量的实用工具以支持可视化编程特性，它们包括： 项目工作区、ClassWizard、AppWizard、WizardBar等。卜面对它们作逐一介绍。 VisualC++的用户界面如图1所示向导栏菜单条工具条I的豆件拿iXVH3工唆MtiT*9D||CAb«viDlg6二||抑Idas*eheIieeI>币IZJJ♦SJ正中支Vt---「1•，”・皿P1TTFTF¥$ClMAWlza^dya%a■尸nt0。ui尸I:山mlFtancticn卜尸〔口¥/《\hFj<_U9FtTiiALt 亡曲上务p^lJE»Ct»d!luirEiidlu«id[i■口。斯=ripgGfCDQtdC量cvi3ngb，pCX>;^/JiHFj4_U!lftTiJftLd三IW， 等/iIIe®ie*t£lfli&SiC*・YCAbddlDIg«・：CMMgFe■£：甘叱CTesiA^p«","Z；CTediODC："FCTtsWlu收[»-_।GlcbalaFTFh想宣"tbK/Di"IR^l titinqnFJ4_lti;C4£ADOUt[>lg、//HaatimqoIwaniiltr-G//J>hFJ4_HSGDEDLhRE^HESSHG£_l1flF()tEAboutbl^:LtAbOiUtblyf>[C'DidilagqCAEi(hJtDi.gLiIiDil>^1小m4HFJ4_DA1ft_[hlT(eADOUtl>lg>

二IMT由二1HITCalSB■二项目工作区 调试区图1VisualC++的用户界面2.2.1项目工作区VisualC++使用项目工作区来组织项目、元素以及项目信息在屏幕上出现的方式。在桌面上，工作区以窗口方式组织项目、文件和项目设置。每个项目视图都有一个相应的文件夹，包含了关于该项目的各种元素。展开该文件夹可以显示该视图方式下工作区的详细信息。项目视图区包含三个选项卡，单击各选项卡，可以在其间切换。这三个选项卡从左到右依次是：ClassView、ResourceView、FileView。（1）ClassView（类视图）显示项目中定义的C+谈，展开文件夹显示工程中所定义的所有类，展开类可查看类的数据成员和成员函数以及全局变量、函数和类型定义。 ClassView显示所有已定义的类以及这些类中的数据成员、成员变量。⑵ResourceView（资源视图）显示项目中所包含的资源文件。展开文件夹可显示所有的资源类型，如图3。_孑ItrBlfBSOuf+L_JAcjcclcr^lor=i-*-JDIulQglinnAJBOUTnox.+I |Icuri+_iiMeou十L_JStringTabic+LJToolibnr+I1Ve-rsinn图2类视图图3资源视图⑶FileView（文件视图）显示所创建的工程。展开文件夹可以查看工程中所包含的文件，如图42.2.2向导栏I男■:Clas…]妁 MfIIeVI…]"图4文件视图jcExdlglessDoc▼]jcExdlglessDoc▼][AlldassmembersjjJlCExdlglessDoc图5向导栏WizardBar会自动跟踪用户程序的上下文一一比如，当文本编辑器中的光标从一个函数移动到另一个函数时，Wizard的显示会自动更新。WizardBar工具条包含了三个相关的下拉列表框： 类（Class）、过滤器（Filter）和成员（MemberWizardBar最右边是一个ActionControl下拉选项，单击ActionControl 的向下箭头符号会弹出一个菜单，用于执行跳到函数定义、增加消息处理函数等操作。.图像二值化的算法及实现论文研究对象论文主要研究BM哈式的灰度图像文件。BMP（BitmapPicture）文件格式是Windows系统交换图形、图像数据的一种标准格式。BMPS像的数据由四个部分组成，如表1所示。表1BMP图像文件结构文件部分属性说明BITMAPFILEHEADER（位图文件头）bfType文件类型，必须是0x424D,即字符串“BMbfSize指定文件大小，包括这14个字节bfReservered1r保留字，不用考虑bfReservered2保留字，不用考虑bfOffBits从文件头到实际位图数据的偏移字节数BITMAPINFOHEADE（位图信息头）biSizeP该结构长度，为40biWidth图像的宽度，单位是像素biHeight图像的高度，单位是像素Rbiplanesr位平面数，必须是1,不用考虑biBitCount指定颜色位数，1为二值，4为16色，8为256色，16,24、32为真彩色biCompression指定是否压缩，有效的值为bi_RGBBI_RLE8,BI_RLE4,BI_BITFIELDSbiSizeImage实际的位图数据占用的字节数biXPelsPerMeter目标设备水平分辨率，单位是每米的biYPelsPerMeter目标设备垂直分辨率，单位是每米的biClrUsed实际使用的颜色数，若该值为0,则使用颜色数为2的biBitCount次力种biClrImportant图像中重要的颜色数，若该值为0,则所有的颜色都是重要的Palette（调色板）rgbBlue该颜色的蓝色分量rgbGreen该颜色的绿色分量rgbRed该颜色的红色分量rgbReservedr保留字ImageData（位图数据）按像素按行优先顺序排序，每一行的字节数必须是4的整倍数第一部分为位图文件头BITMAPFILEHEADE定是个结构体，其定义如下:typedefstructtagBITMAPFILEHEADER{WORD bfType;DWORDbfSize；WORD bfReserved1 ;WORD bfReserved2 ;DWORDbfOffBits;}BITMAPFILEHEADER这个结构的长度是固定的，为14个字节（WOR为无符号16位二进制整数，DWORD无符号32位二进制整数）。第二部分为位图信息头BITMAPINFOHEADER是一个结构，其定义如下：typedefstructtagBITMAPINFOHEADER{DWORDbiSize；LONG biWidth;LONG biHeight;WORD biPlanes；WORDbiBitCount;DWORDbiCompressionDWORDbiSizeImage；LONG biXPelsPerMeter；LONG biYPelsPerMeter；DWORDbiClrUsed;DWORDbiClrImportant;}BITMAPINFOHEADER这个结构的长度是固定的，为40个字节（LONG为32位二进制整数）。其中，biCompression的有效值为BI_RGBBI_RLE&BI_RLE4BI_BITFIELDS,这者B是一些Windows定义好的常量。由于RLE4ffiRLE8的正缩格式用工勺不多，一般仅讨论biCompression的有效值为BI_RGB即不压缩的情况。第三部分为调色板(Palette),当然，这里是对那些需要调色板的位图文件而言的。真彩色图像是不需要调色板的，BITMAPINFOHEAD球直接是位图数据。调色板实际上是一个数组，共有biClrUsed个元素(如果该值为零，则有2的biBitCount次方个元素)。数组中每个元素的类型是一个RGBQUAD构，占4个字节，其定义如下：typedefstructtagRGBQUAD{BYTErgbBlue;BYTErgbGreen；BYTErgbRed;BYTErgbReserved;}RGBQUAD第四部分就是实际的图像数据。对于用到调色板的位图，图像数据就是该像素颜色在调色板中的索引值，对于真彩色图像，图像数据就是实际的 RGB值。下面就2色、16色、256色和真彩色位图分别介绍。对于2色位图，用1位就可以表示该像素的颜色(一般0表示黑色，1表示白色)，所以一个字节可以表示8个像素。对于16色位图，用4位就可以表示一个像素的颜色，所以一个字节可以表示2个像素。对于256色位图，一个字节刚好可以表示1个像素。下面有两点值得注意：(1)每一行的字节数必须是4的整倍数，如果不是，则需要补齐。BMPC件的数据存放是从下到上，从左到右的，也就是说，从文件中最先读到的是图像最下面的一行的左边的第一个像素，然后是左边的第二个像素，接下来是倒数第二行左边第一个像素，左边第二个像素。依次类推，最后得到的是最上面的最右边的一个像素。灰度变换黑白图片的黑白变换叫灰度变换，彩色图片的色彩变换也叫灰度变换。一幅彩色图像的象素矩阵中每个象素由RGB3种颜色按一定的比例混合形成一种颜色来表示， 比如黑色使RGB(0,0,0),纯红色是RGB(0,1,0)…。在处理相片时，有时可能因为环境光源太暗，使RG电的值偏小，就会使图形太暗看不清，如果环境光源太光，又使图像泛白，通过灰度变换，就可以使RGBS调和到合适的程度，使相片变得漂亮。灰度变换原理：首先提取一幅图像的一个象素，在BM昭式中一个象素由8位红色亮度值，8位绿色亮度值和8位蓝色亮度值组成，只要按一定的变换函数去变换。灰度变换公式： GRAY=0.11*R+0.59*G+0.3*B (2)图像二值化常用算法简单阈值通过直接对图像的灰度直方图曲线进行分析，判断极小值所在的灰度级，再对所有极小值点进行相应的分类合并来最终确定阈值的位置，避免出现过分割的现象。在预处理阶段本文采用了基于各向异性扩散的平滑方法对图像进行滤波，以消除噪声，同时还可以使图像的灰度直方图曲线更为连续，从而加速算法的执行效率。简单阈值是指用灰度变换来研究灰度图像的一种常用方法， 即设定某一阈值将灰度图像的像素分成大于阈值的像素群和小于阈值的像素群两部分。例如输人灰度图像函数为，输出二值图像函数为，则Threshold=127 (3)阈值(threshold)是把目标和背景区分开的标尺，选取适当的阈值就是既要尽可能保存图像信息，又要尽可能减少背景和噪声的干扰，这是选择阈值的原则。

图6 图6 简单阈值效果图像由于简单阈值法过于简单，丢失部分信息。随着时代的发展人们对于图像的要求也越来越高，这种简单的方法已经不能够满足人们的要求。Bernsen法局部阈值法把灰度阈值选取为随像素位置变化而变化的函数，它是一种动态选择阈值的自适应方法。设图像在像素点(i,j)处的灰度值为f(i,j),考虑以像素点(i,j)为中心的(20+1)父(2◎+1)窗口，(2：+1表示窗口的边长)，则Bernsen算法可以描述如下：计算图像中各个像素点(i,j)的阈值T(i,j)T(i,j)=0.5(maxf(im,jn)minf(im,jn))；.--,：：m_- -m-■一四一. 一四：；b(i,j)=广对图像中各像素点(i,j)用b(i,j)b(i,j)=广f(i,j)；T(i,j)f(i,j)-T(i,j)用I存储灰度图像的值，设I为NxM,把I边界扩展成(N+2WM+2)extend矩阵首先读取原图像I的大小为NmM。由于I中的元素不是每个都是在3M3窗口的中心，所以需要对灰度图像I进行扩展。首先创建一个(N+2)m(M+2)的矩阵extend,把矩阵I中的像素extend(i-1,j,1)=I(i,j),而第一行和最后一行，第一列和最后一列的填充依据是以它靠近的行或列为对称轴进行填充。遍历从extend(2,2)到extend(N+1,M+1)的像素，并取以当前像素为中心的3父3窗口的最大像素max和最小像素min,依据公式t=0.5父(max+min)求出阈值t。把灰度图像矩阵I赋值于另一矩阵B,以免改变当前得到的灰度图像矩阵。遍历该矩阵B,对当前灰度值与t比较，如果大于赋予1,判为目标像素类，否则赋予0,作为背景像素类。显示得到的二值图像BoBernsen方法流程图如图7所示。专业资料专业资料图7Bernsen方法流程图 Bernsen方法结果分析专业资料专业资料WORDWORD完美格式 下载可编辑(d) (e) (f)图8原图(a)lenna原始图像(256x256);(b)barbara原始图像(256x256);(c)peppers原始图像(256X256);(d)lenna原始图像(512x512)；(e)barbara原始图像(512x512)；(f)peppers原始图像(512X512)通过调研参考文献可知，对图8这六幅图像都采用Bernsen算法进行图像二值化处理。其结果如图所示。lenna(256X256)、lenna(512X512)原始图像通过Bernsen方法二值化后的所得到的二值图像如9所示。barbara(256x256)>barbara(512x512)原始图像通过Bernsen方法二值化后的所得到的二疝图像如图10所示。peppers(256X256)、peppers(512x512)原始图像通过Bernsen方法二值化后的所得到的二值图像如图 11所图9图9Bernsen方法二值化(a)二值lenna图像(256乂256);lenna图像(b)二值lenna图像(512X512)barbarabarbara图像二值barbara图像(512X512)(a) (b)图10Bernsen方法二值化(a)二值barbara图像(256乂256); (b)(a) (b)图11Bernsen方法二值化peppers图像(a)二值peppers图像(256乂256); (b)二值peppers图像(512乂512)实验结果分析如表112120表2Bernsen方法处理不同像素的不同图像性能指标图9a图9b图10a图10b图11a图11b阈值(T)时间(S)/s4.20328.75002.29708.59403.18708.8440嫡值(H)0.91310.85330.93390.92380.87150.8922因为Bernsen方法不存在预取阈值，得到的都是局部阈值，随像素的变化而变化，所以没有固定的阈值。图像越大用Bernsen算法对其进行二值化处理的时间越长。3.3.3Otsu法Otsu在1979年提出的最大类间方差法(有时也称之为大津法)，该方法的基本思想是：设阈值将图像分割成两组，一组灰度对应目标，另一组灰度对应背景，则这两组灰度值的类内方差最小，两组的类间方差最大。对图像 Image,记t为目标与背景的分割阈值，目标像素数占图像比例为剑，平均灰度为巳；背景像素数占图像比例为以，平均灰度为明。图像的总平均灰度为：N=00。)%(。+管。)也(t)。从最小灰度值到最大灰度值遍历t,当t使得值g=%(5-N)2+与(丹-N)2最大时t即为分割的最佳阈值。对大津法可作如下理解：该式实际上就是类间方差值，阈值 t分割出的目标和背景两部分构成了整幅图像，而目标取值为,概率为80,背景取值匕，概率为孙，总均值为N,根据方差的定义即得该式。因方差是灰度分布均匀性的一种度量，方差值越大说明构成图像的两部分差别越大，当部分目标点错分为背景或部分背景点错分为目标点都会导致两部分差别变小，因此使类间方差最大的分割意味着错分概率最小。记f(i,j)为NMM图像(i,j)点处的灰度值，灰度级为 N,不妨假设f(i,j)取值[0,m-1]0记p(k)为灰度值为k的频率,则有：TOC\o"1-5"\h\zp(k) P (6)MNf(i,j)小假设用灰度值t为阈值分割出的目标与背景分别为：｛f(i,j)Wt｝和｛f(i,j)〉t｝,于是目标部分比例：■0(t)-p(i) (7)0mM目标部分点数：N0(t)=MN'、p(i) (8)0父』背景部分比例：TOC\o"1-5"\h\z」(t)="P(i) (9)t::i理二背景部分点数：Ni(t)=MN'、P(i) (10)t：：i•型1目标均值：」0(t)='ip⑴/o(t) (11)0出炉背景均值：L(t)= 、ip(i)/i(t) (12)t：-iCm__1总均值： ,一口=o(t)」o(t) '1(t)」1(t) (13)大津法指出求图像最佳阈值g的公式为：g=Arg段盟1[so(t)(均⑴一"了+叫⑴(邑小)一邑)2] (14)该式右边括号内实际上就是主问方差值，阈值 g分割出的目标和背景两部分构成了整幅图像，而目标值七(t),概率为叫(t),背景取值也⑴，概率为,(t),总均值为4根据方差的定义即得该式。因方差是灰度分布均匀性的一种度量，方差值越大，说明构成图像的两部分差别越大，当部分目标错分为背景或部分背景错分为目标都会导致两部分差别变小，因此使类间方差最大的分割意味着错分概率最小，这便是大津法的真正含义。其实，换一种思想也能很好的理解大津方法。 ％(t)和为(t),可以分别代表目标和背景的中心灰度，口则代表整幅图像的中心灰度，要使目标和背景得到最好的分割，当然希望分割出的目标尽量远离图像中心，即 (N0(t)-N)2或|也。)-打尽量大，背景也尽量远离中心，即(七(。-与2或|方。)-阿尽量大，由于希望两者都大，于是有：(1)两者之加权和最大：g=ArgM*「o(t)(」o(t)-」)2 ，1(t)(L(t)-」)2] (15)o'、•【'、m-1(2)两者之积最大：g=ArgMax[(」0(t)-」)2(L(t)-」)2] (16)o父im-1注意到有■=S°(t)艮o(t)+%(t)t(t),且艮o(t)«艮4艮1。),因此有：So。)(乜⑴-均2+叫。)(艮1(t)一邑)2=(艮o(t)—艮)2(地⑴一均2。可见是二者等价的。Otsu方法流程图如图7所示。

图12Otsu方法流程图Otsu方法实验结果分析对图8这六幅图像都采用Otsu算法进行图像二值化处理。其结果如图所示。lenna(256X256)、lenna(512X512)原始图像通过Otsu方法二值化后的所得到的二值图像如图 13所示。barbara(256x256)>barbara(512x512)原始图像通过Otsu方法二值化后的所得到的二值图像如图14所示。peppers(256X256)、peppers(512X512)原始图像通过Otsu方法二值化后的所得到的二值图像如图 15所示。(a) (b)图13Otsu方法二值化lenna图像(a)二值lenna图像(256x256);(b)(a)二值lenna图像(256x256);(b)(a)图14Otsu二值lenna图像(512X512)(b)方法二值化barbara图像(a)二值barbara图像(256x256); (b)二值barbara图像(512x512)(a) (b)图15Otsu方法二值化barbara图像(a)二值peppers图像(256乂256); (b)二值peppers图像(512乂512)实验结果分析如表3。表3Otsu方法处理不同像素的不同图像性能指标图13a图13b图14a图14b图15a图15b阈值(T)10912911411896106时间(S)/s1.09403.87501.18703.67201.04703.7810嫡值(H)0.99931.00000.99970.99440.93560.9998图(a)为256X256的图像，图(b)为512X512的图像。以图13为例，lenna512X512的图像时间为3.8750,阈值为129;lenna256X256的时间为1.0940,阈值为109。可知，512X512图像的执行时间要比256X256图像的时间要长，说明图像越大，用Otsu方法对其进行二值化处理所需时间越长；大图像的阈值也比小图像要大。嫡值是衡量图像有序化程度的一个度量，嫡值越低图像越有序，图13b即lenna(512x512)在进行比较的图像中是最混乱的。0tsu 方法和Bernsen方法实验结果比较比较其它参考文献可知，在表中用一些变量保存了有用的实验数据， T表示Otsu方法、Bernsen方法的运行时间，单位为秒。H表示由Otsu方法、Bernsen方法处理得到的二值图像的嫡值，单位为比特/符号。嫡值的定义为：MH(x)八p(x。10g21/pxi (17)i1

图像的嫡值反映了整幅图像的效果。信息嫡是信息论中用于度量信息量的一个概念。一幅图像越是有序，信息嫡就越低；反之，一幅图像越是混乱，信息嫡就越高。所以，信息嫡也可以说是衡量图像有序化程度的一个度量。程序的运行时间代表算法运行效率，也是算法的一项性能指数。两种性能比较如表4所示【⑵表4两种方法性能比较二值化算法Otsu方法Bernsen方法性能指标阈值(T)时间(S)/s嫡值(H)时间(S)/s嫡值(H)lenna256X2561091.09400.99934.20320.9131512X5121293.87501.00008.75000.8533barbara256X2561141.18700.99972.29700.9339512X5121183.67200.99448.59400.9238peppers256X256961.04700.93563.18700.8715512X5121063.78100.99988.84400.8922Bernsen方法由于不存在预取阈值，得到的都是局部阈值，随像素的变化而变化，没有固定的阈值。程序运行时间T:T(Bernsen)>T(Otsu),可见Bernsen算法要寻找局部极大、极小值，因此速度较慢；TOtsu(256x256)<TOtsu(512X512),TBernsen(256乂256)<TBernsen(512X512),可见图像越大，二值化程序运行所需时间就越长。二值图像的嫡值：H(Otsu)>H(Bernsen),可以见Otsu方法对于光照不均、噪声干扰大的图像，其二值化效果较差，其反映了整个图像的整体灰度分布情况。3.3.4迭代法迭代法是另外一种全局二值化方法，它要求图像分割阈值的算法是基于逼近的思想，首先选择一个近似阈值作为估计值的初始值，然后进行分割，产生子图像，并根据子图像的特性来选取新的阈值，再利用新的阈值分割图像，经过几次循环，使错误分割的图像像素点降到最少。这样做的效果好于用初始阈值直接分割图像的效果。具体算法步骤如下所(1)求出图像中的最小灰度值和最大灰度值，分别记为Zmin和Zmax,则阈值初值T。为：Zmin•ZmaxToZmin•ZmaxTo二 2(2)根据阈值Tk将图像分割成目标和背景两部分，求出两部分的平均灰度值'z i, j Ni,jZ二也^ 乙0 _' N i,jZi，j 五、zi,jNi,jzi,jTkZG二 '、Ni,jzi，jTkZo和Zg：(18)(19)(20)为z(i,j)式中z(i,j)为图像上(i,j)点的灰度值；N(i,j)为(i,j)点的权重系数，一般为z(i,j)图16迭代法算法实现图图16迭代法算法实现图17迭代法效果图像的个数。(3)利用下式求出新阈值：T*豆三 (21)(4)如果Tk=Tk+i,则结束，否贝UK=K+1,车$向(2)。经实验测试发现，对于直方图双峰明显、谷底较深的图像，迭代法可以较快地获得满意结果；但是对于直方图双峰不明显或图像目标和背景比例差异悬殊， 迭代法所选取的阈值不如OTSUt。迭代法算法实现如图16示。(钎算图0废度宜方图)一1 一二~出明*的.大次瀛、他和•小灰度值.分别记班卜4K和Ms,令初始网回''\^7a=CZii*Jt+/in*)I二/^*值相♦♦，分、(为前景和背景，分别求出)3者的平均塔值4和/ 拿出金懒值~kn*l-<z»+z()/2对某些特定图像，某些微小的变化却会引起分割效果的巨大变化，两者数据只是稍有变化，但分割效果却反差极大，其中原因还有待进一步研究。4.总结Bernsen法虽然能够根据局部灰度特性来自适应地选取阀值,有较大的灵活性，但是仍然存在缺点和问题：实现速度慢、有伪影现象、有笔画断裂现象。简单阈值法能减少背景和噪声的干扰迭代所得的阀值分割的图像效果不好， 丢失部分图像信息，所以我主要采用Otsu法和迭代法。Otsu法是经典的非参数、无监督自适应阈值方法，是一种直方图技术。Otsu法实现简单，对于具有明显双峰直方图的图像效果明显，但对于低对比度和光照不均匀的图像效果不佳，抗噪声能力差，因而应用范围受到极大限制。对于直方图双峰明显、谷底较深的图像，迭代法可以较快地获得满意结果；但是对于直方图双峰不明显或图像目标和背景比例差异悬殊，迭代法所选取的阈值不如OTSUt。几种常用的图像二值化算法各有优缺点，由于时间的问题本次设计没有逐一研究。LiuJieping,YuYinglin.Aflexiblemethodforimagenoiseremoval[J].JouralofSouthChinaUniversityTechnology,2000,28(2):60〜63.APiva,MBarni,FBartolini,VCappellini.DCT-basedwatermarkrecoveringwithoutresortingtotheuncorruptedoriginalimage[A].Proceedingsof4thIEEEInternationalConferenceonImageProcessingICIP'97[C].SantaBarbara,CA,U1S9A7CIPZhongwei.ImagewatermarkingusingLegendrearray[J].Journalofchinainstituteofcommunications,2001,22(1):1〜6.[4]王建卫.彩色图像的中值滤波算法的改进与应用 [J].哈尔滨商业大学大学报：自然科学版.2006,22(4):67〜69.[5]陈丹,张蜂,贺贵明.一种改进的文本图像二值化算法