一种基于分水岭变换的细胞图像分割方法
- 格式:pdf
- 大小:322.06 KB
- 文档页数:3
基于分水岭算法的遥感图像分割方法研究
遥感图像分割是将遥感图像分成不同区域的过程,用于提取地物信息和地图制图等应用。
基于分水岭算法的遥感图像分割方法是一种常用的分割方法。
分水岭算法基于图像中明显的梯度梯度变化,将图像分割成不同的区域。
本文将讨论基于分水岭算法的遥感图像分割方法。
分水岭算法的思想是将图像中的像素点看作梯度值,然后从梯度最小值的位置开始,让梯度值不断上升,直到与相邻像素梯度发生冲突。
这样,图像就可以被分割为相邻的区域,并且每个区域内的像素梯度值相同。
在分水岭算法中,一些重要的操作,如距离变换和标记等,都是用于解决边界问题的。
1. 对遥感图像进行预处理,将图像平滑化、去噪等操作。
2. 使用梯度检测算法对图像进行梯度计算,得到梯度图像。
3. 对梯度图像进行距离变换,得到距离图像。
5. 对分割图像进行形态学操作,如膨胀、腐蚀等操作,得到分水岭标记。
7. 对分水岭图像进行过滤和合并操作,得到最终分割结果。
基于分水岭算法的遥感图像分割方法可以有效地提取地物的空间信息,尤其在大尺度的遥感图像上,具有较好的效果。
但是,该方法容易出现过分割和欠分割的情况,这需要根据实际情况对算法进行优化和调整。
总之,基于分水岭算法的遥感图像分割方法是一种比较常用的算法,可以用于提取地物信息、地图制图等应用。
随着计算机算力的提高和算法的改进,分水岭算法的应用前景将会更加广泛。
基于改进分水岭方法的癌细胞图像分割技术(doc 31页)分类号:TP391.1 U D C:D10621-408-(2012)0670-0 密级:公开编号:2008023005成都信息工程学院学位论文基于改进分水岭方法的癌细胞图像分割技术论文作者姓名:张春香申请学位专业:生物医学工程申请学位类别:工学学士指导教师姓名(职称):赵旭论文提交日期:2012年06月01日独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。
据我所知,除了文中特别加以标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得成都信息工程学院或其它教育机构的学位或证书而使用过的材料。
与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。
签名:日期:2012年06月09日关于论文使用授权的说明本学位论文作者完全了解成都信息工程学院有关保留、使用学位论文的规定,有权保留并向国家有关部门或机构送交论文的复印件和磁盘,允许论文被查阅和借阅。
本人授权成都信息工程学院可以将学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。
(保密的学位论文在解密后应遵守此规定)签名:日期:2012年06月09日基于改进分水岭方法的癌细胞图像分割技术摘要细胞识别系统是针对癌变组织来分割粘连的细胞,用来进行细胞识别,获取它们的个数,面积,周长等特征。
但是在获取癌细胞图像的过程中,由于很多外在和人为的因素,会存在细胞粘连在一起、并且分布不均的情况。
那么就对后续的研究造成了很大的困难。
传统的分水岭算法在图像分割时,对噪声极为敏感和易于出现过分割。
量化误差、噪声及梯度纹理的局部不规则均会使在后面分割时许多的过小封闭区域, 容易丢失细胞的重要轮廓。
针对以上问题,本人在传统分水岭的基础上,实现了一种解决此问题的有效方法。
一、通过中值滤波器滤除部分噪声,再使用形态学滤波平滑图像。
基于分水岭算法的彩色细胞图像分割摘要随着影像医学的发展,通过对细胞涂片影像的分析,从而对细胞影像进行区分和识别成为重要的研究课题。
细胞图像分割是细胞图像分析和识别的重要步骤。
图像分割是将图像中具有特殊含义的不同区域区分开来,是图像处理的关键步骤。
分割后的子区域互不交叉,每一个区域满足特定性质的一致性。
人体细胞种类繁多、形态多样且图像质量也很不相同,而分析应用中对细胞图像分割的质量却要求较高,所以细胞图像的自动分割极为重要且困难很大。
彩色图像与灰度图像相比,信息量更为丰富,而且颜色的描述方法也,较多。
很多经典算法只能对二值图像或灰度图像进行运算。
为此,关于彩色细胞图像的分割研究成为一个非常活跃的研究领域。
本文针对彩色细胞图像经过染色处理的特点,提出了一种彩色细胞图像的分割方法。
以快速分水蛉算法为主要分割算法,为了较好地抑制彩色细胞图像背景噪声,选择更符合人类视觉感知的HSI颜色空间,结合自动阈值和色度提出去除图像背景的方法。
同时,使用中值滤波和均匀化处理,有效地克服了分水岭算法的过分割现象。
针对细胞图像特点改进了区域合并算法。
得到了较准确的分割结果。
本文首先概括介绍了图像分割的意义及发展现状,概述了当前主要的图像分割算法。
其次,介绍了彩色图像颜色空间和快速分水岭算法的基本思想及实现方法。
最后列出了实验流程和实验结果并进行了讨论。
关键词:图像分割,HSI颜色空间,分水岭二.分水岭算法本章从对分水岭算法的定义出发,对分水岭算法的发展过程中的不同实现方法进行比较,阐明快速分水岭算法的优越特性及实现方法。
(一)分水岭算法的定义分水岭分割的最初算法是针对地形数字高程模型提出的.目前分水岭算法在图像分割领域正得到广泛应用.分水岭算法的定义121J对一幅二维灰度图像,,Jr的定义域为Dr cZ2,,取离散灰度值【0,N】,将该值视为对应像素点的高度,Ⅳ为一正整数。
用G表示相应的数字格网(以四邻域为例)o图像I中点p和g之间一条长度为z的路程≯为由点,Pl,⋯Pt-l,P1)组成的(斗1)元组,有Po=P,Pl=q,且Vf∈【l,,】,(Pf-I,Pi)∈G (3-2)将路径P的长度标识为^纠,点p的邻域集标识为: (p):%Q)=p’∈Z2,(DP’)∈回(3‘3)图像f在高度矗的一个极小区膨(minimum)定义为由高度值为h的点组成的一连通区域,从该区域肘中的一点出发到达任一高度低于h的像素点。
毕业论文声明本人郑重声明:1.此毕业论文是本人在指导教师指导下独立进行研究取得的成果。
除了特别加以标注地方外,本文不包含他人或其它机构已经发表或撰写过的研究成果。
对本文研究做出重要贡献的个人与集体均已在文中作了明确标明。
本人完全意识到本声明的法律结果由本人承担。
2.本人完全了解学校、学院有关保留、使用学位论文的规定,同意学校与学院保留并向国家有关部门或机构送交此论文的复印件和电子版,允许此文被查阅和借阅。
本人授权大学学院可以将此文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本文。
3.若在大学学院毕业论文审查小组复审中,发现本文有抄袭,一切后果均由本人承担,与毕业论文指导老师无关。
4.本人所呈交的毕业论文,是在指导老师的指导下独立进行研究所取得的成果。
论文中凡引用他人已经发布或未发表的成果、数据、观点等,均已明确注明出处。
论文中已经注明引用的内容外,不包含任何其他个人或集体已经发表或撰写过的研究成果。
对本文的研究成果做出重要贡献的个人和集体,均已在论文中已明确的方式标明。
学位论文作者(签名):年月关于毕业论文使用授权的声明本人在指导老师的指导下所完成的论文及相关的资料(包括图纸、实验记录、原始数据、实物照片、图片、录音带、设计手稿等),知识产权归属华北电力大学。
本人完全了解大学有关保存,使用毕业论文的规定。
同意学校保存或向国家有关部门或机构送交论文的纸质版或电子版,允许论文被查阅或借阅。
本人授权大学可以将本毕业论文的全部或部分内容编入有关数据库进行检索,可以采用任何复制手段保存或编汇本毕业论文。
如果发表相关成果,一定征得指导教师同意,且第一署名单位为大学。
本人毕业后使用毕业论文或与该论文直接相关的学术论文或成果时,第一署名单位仍然为大学。
本人完全了解大学关于收集、保存、使用学位论文的规定,同意如下各项内容:按照学校要求提交学位论文的印刷本和电子版本;学校有权保存学位论文的印刷本和电子版,并采用影印、缩印、扫描、数字化或其它手段保存或汇编本学位论文;学校有权提供目录检索以及提供本学位论文全文或者部分的阅览服务;学校有权按有关规定向国家有关部门或者机构送交论文的复印件和电子版,允许论文被查阅和借阅。
基于分水岭算法的图像分割在图像分割领域中,分水岭算法是一种经典的方法。
分水岭算法最初是根据地理学领域的概念而来,类似于山脉与河流的分割。
在图像中,局部极小值和极大值处形成了类似于山脉和河流的结构,而分水岭算法就是通过泛洪法来模拟这个过程,将局部极小值连接起来,形成一个分水岭,从而完成图像的划分。
分水岭算法的主要思想是基于图像中的梯度信息进行分割。
具体而言,首先计算图像的梯度幅值,并根据阈值将其转化为二值图像。
然后使用距离变换的方法,计算图像中每个像素点到最近边界的距离,得到一副距离变换图像。
接下来,通过将距离变换图像中的局部极小值进行标记,并进行连通区域的合并,构建无向图。
最后,使用分水岭变换的方法将图像划分成不同的区域。
分水岭算法的优点是能够获得较为准确的图像分割结果,并且具有良好的鲁棒性。
然而,分水岭算法也存在一些问题。
首先,分水岭算法较为复杂,计算量较大,运行时间长。
其次,分水岭算法对图像中的噪声较为敏感,在噪声较多的图像上容易产生过度分割或欠分割的情况。
此外,分水岭算法对初始种子点的选择也较为敏感,不同的初始种子点可能得到不同的分割结果。
针对上述问题,研究者们提出了一些改进的方法。
例如,可以采用图割算法来替代分水岭算法的标记和分割阶段,以提高分割的效果。
同时,可以通过预处理步骤来降低图像中的噪声或平滑图像,以改善分割的质量。
此外,也可以采用自适应的分水岭算法来选择合适的初始种子点,以减少分割结果的不确定性。
总结来说,基于分水岭算法的图像分割是一种经典而有效的方法。
虽然该算法存在一些问题,但通过改进和优化,可以得到更好的分割结果。
随着计算机视觉和图像处理的不断发展,我们相信分水岭算法在未来将会得到更广泛的应用。
基于分水岭算法的遥感图像分割方法研究遥感图像分割是遥感图像处理中的一项重要任务,它的目的是将图像中的不同区域分割成具有相似特征的子区域。
这项工作在农业、环境监测、城市规划等领域有着重要的应用价值。
分水岭算法是一种常用的图像分割方法,它基于图像中的灰度梯度信息来实现分割。
本文将围绕基于分水岭算法的遥感图像分割方法展开研究,探讨其原理、实现过程以及应用效果。
一、分水岭算法原理及应用分水岭算法最初来源于地质学中的地下水分割理论,后来被引入到图像处理领域中。
它的原理是将图像看作地形地貌,图像中的亮度信息对应地形的高度,然后利用不同区域之间的梯度信息来确定分割线,实现图像的分割。
在遥感图像处理中,分水岭算法被广泛应用于不同类型的地物分割,包括植被、水域、建筑等。
分水岭算法的基本思路是从图像中的局部最小点(或者称为浸没点)出发,构建出一系列的水域,然后根据这些水域的相互关系来确定整个图像的分割线。
在遥感图像中,这些局部最小点往往对应着不同的地物或者地物边界,因此通过分水岭算法可以实现对图像中不同地物的精确分割。
分水岭算法还可以应用于图像的边缘检测、纹理分割等领域。
基于分水岭算法的遥感图像分割方法一般包括以下几个步骤:预处理、特征提取、分水岭算法实现和结果后处理。
在预处理阶段,需要对原始遥感图像进行几何校正、大气校正等操作,以保证图像的质量和准确性。
接下来进行特征提取,一般采用像元级的特征提取方法,包括灰度、颜色、纹理等特征。
然后利用这些特征信息构建图像的梯度信息,为后续的分水岭算法做准备。
分水岭算法的实现一般使用连通区域分割算法(Watershed segmentation algorithm),它是一种基于梯度信息的像素聚类算法,能够根据图像的梯度分布实现对图像的分割。
在算法实现过程中需要注意对梯度信息进行分析和处理,以保证分割结果的准确性和可靠性。
最后对得到的分割结果进行后处理,包括去除小面积的噪声点、填补分割线等操作,以得到最终的分割图像。
ComputerEngineeringandApplications计算机工程与应用2008,44(4)1引言基于浸水模型的分水岭(Watersheds)[1]变换将图像看成是一幅地形图。
在最低点处刺穿一个小孔,然后湖水开始由此小孔向对应的聚水盆地注水形成水库。
为了防止两个不同的局部最小点对应的水库汇聚到一起,在它们的相接处建立起一个水坝,当地形被完全浸没在湖面以下之后,所有的水坝就构成了分水岭。
分水岭变换具有边界定位准确,且分割结果连续的特点,被广泛的应用于医学图像分割中。
然而直接应用分水岭算法会产生过分割问题[2]:由于图像中存在很多极小值点,分割结果淹没在大量的不相关结果中。
目前主要有两种解决方式,一种方式是利用先验知识构造标记图像,然后用标记图像对原图进行重构,填充掉多余的盆地[3]。
基于此,HaiGao等[4]提出了一种依赖于不相交集的基于标记的快速分割方法;Roberto等[5]则针对血管分割问题提出了一种稳定的标记提取策略。
基于标记的方法分割结果准确,且速度快,但分割的结果的好坏直接依赖于标记,而标记的提取需要大量的领域知识。
另一种方式则是利用区域的相似性或差异性,采用分层多尺度的方法,对分割结果进行合并。
Jung[6]借助于小波变换,实现了多尺度的分水岭变换;赵晨光等[7]提出了基于分水岭差异性的多级图像分割方法。
分层多尺度方法相对稳定,且不依赖于先验知识,但该方法速度较慢,而且在有些情况下,由于在底层小尺度时,区域的特征容易受到噪声干扰而产生的错误,进而在后来高层大尺度合并过程中,错误被放大,导致最终分割结果的不准确。
本文综合运用上述两种方法:在取得了彩色细胞图像的饱和度分量的基础上,首先利用双阈值分割确定细胞的大体位置和范围,然后运用所提出的快速的h-dome提取算法和形态学运算得到细胞对象的粗糙标记,再以此为基础,对原图像的梯度图做形态学灰度重建,并借助于快速分水岭变换,得到粗糙的分割结果。
ComputerEngineeringandApplications计算机工程与应用2008,44(4)1引言基于浸水模型的分水岭(Watersheds)[1]变换将图像看成是一幅地形图。
在最低点处刺穿一个小孔,然后湖水开始由此小孔向对应的聚水盆地注水形成水库。
为了防止两个不同的局部最小点对应的水库汇聚到一起,在它们的相接处建立起一个水坝,当地形被完全浸没在湖面以下之后,所有的水坝就构成了分水岭。
分水岭变换具有边界定位准确,且分割结果连续的特点,被广泛的应用于医学图像分割中。
然而直接应用分水岭算法会产生过分割问题[2]:由于图像中存在很多极小值点,分割结果淹没在大量的不相关结果中。
目前主要有两种解决方式,一种方式是利用先验知识构造标记图像,然后用标记图像对原图进行重构,填充掉多余的盆地[3]。
基于此,HaiGao等[4]提出了一种依赖于不相交集的基于标记的快速分割方法;Roberto等[5]则针对血管分割问题提出了一种稳定的标记提取策略。
基于标记的方法分割结果准确,且速度快,但分割的结果的好坏直接依赖于标记,而标记的提取需要大量的领域知识。
另一种方式则是利用区域的相似性或差异性,采用分层多尺度的方法,对分割结果进行合并。
Jung[6]借助于小波变换,实现了多尺度的分水岭变换;赵晨光等[7]提出了基于分水岭差异性的多级图像分割方法。
分层多尺度方法相对稳定,且不依赖于先验知识,但该方法速度较慢,而且在有些情况下,由于在底层小尺度时,区域的特征容易受到噪声干扰而产生的错误,进而在后来高层大尺度合并过程中,错误被放大,导致最终分割结果的不准确。
本文综合运用上述两种方法:在取得了彩色细胞图像的饱和度分量的基础上,首先利用双阈值分割确定细胞的大体位置和范围,然后运用所提出的快速的h-dome提取算法和形态学运算得到细胞对象的粗糙标记,再以此为基础,对原图像的梯度图做形态学灰度重建,并借助于快速分水岭变换,得到粗糙的分割结果。
最后运用基于区域相似度的区域合并算法对分割结果进行合并,从而提取出细胞轮廓。
整个过程的流程图如图1。
2基于粗糙标记的分割2.1形态学灰度重建形态学灰度重建(GrayscaleMorphologicalReconstruction)[8]能够在保留外轮廓细节的同时将高于标记的山峰削平。
设I、J一种基于分水岭变换的细胞图像分割方法包振健,邸书灵BAOZhen-jian,DIShu-ling石家庄铁道学院计算机系,石家庄050043DepartmentofComputerScienceandTechnology,ShijiazhuangRailwayInstitute,Shijiazhuang050043,ChinaE-mail:heaven_bao@yahoo.com.cnBAOZhen-jian,DIShu-ling.Cellimagesegmentationmethodbasedonwatershedtransformation.ComputerEngineeringandApplications,2008,44(4):230-232.Abstract:Ahybridstrategybasedonwatershedtransformationisproposedandappliedtobonemarrowcellimagesegmentation.Throughtheintegrationofmarker-basedmorphologicalreconstructionandregional-similarity-basedregionmergingalgorithm,thedrawbackthattraditionalmarker-basedwatershedtransformisvulnerabletothemarkers,andtheproblemthatregionmergingalgorithmisslow,andinaccurateisovercome.Experimentalresultsshowthatthemethodisaccurateandefficientforbonemarrowcellimagesegmentation.Keywords:watersheds;marker;grayscalemorphologicalreconstruction;regionmerging摘要:提出了一种混合的分水岭变换策略并应用于骨髓细胞图像分割中。
该策略综合运用基于标记的形态学重建和基于区域相似度的区域合并算法,克服了传统的基于标记的分水岭变换对于标记提取敏感的缺点,及区域合并算法速度慢,结果不精确的问题。
实验结果表明,该方法能够准确高效地分割骨髓细胞图像。
关键词:分水岭;标记;形态学灰度重建;区域合并文章编号:1002-8331(2008)04-0230-03文献标识码:A中图分类号:TP317.4作者简介:包振健(1981-),男,硕士研究生,主要研究方向:图像处理、模式识别;邸书灵(1962-),女,教授,博士,主要研究方向:人工智能、数据挖掘。
收稿日期:2007-06-05修回日期:2007-08-122302008,44(4)是属于相同的图像空间的灰度图,且J≤I,则图像I经由标记图像J的形态学灰度重建ρI(J)定义为:ρI(J)=∧n≥1ε(n)I(J)(1)其中,ε(n)I(J)=(J⊙B)∨I,表示J被B腐蚀后再与I逐点取最大值。
图2演示了这一过程。
2.2粗糙标记提取彩色骨髓细胞图像在饱和度(Saturation)分量上具有较好的可分离性[9]。
根据这一特点,首先提取彩色骨髓细胞图像的S分量图S(I),然后利于双阈值分割S(I)得到二值化结果S1(I),用大圆盘结构体对S1(I)作膨胀运算再求补,得到背景标记B(I);对象标记O(I)按如下步骤进行提取:(1)对S1(I)做距离变换,得到距离图D(I);(2)取D1(I)=D(I)-h,其中h是常量;(3)以D1(I)做标记图像,对D(I)做形态学灰度重建,设得到的重建结果为R(I);(4)取R1(I)=D(I)-R(I);(5)以0为阈值,对R1(I)进行二值化,得到对象标记O(I)(见图4(b))。
上述步骤提取二值图像核心区域的算法又称为h-dome提取,当h=1时上述算法退化为条件膨胀。
设图像高为M,宽为N,借助于双向快速距离变换以及基于先进先出队列的快速形态学重建算法,上述步骤的时间复杂度为O(M*N),是一个高效的算法。
叠加背景标记和对象标记,令M(I)=B(I)+O(I),M(I)就是求得的标记图像(见图4(c))。
2.3基于粗糙标记的分水岭变换分水岭变换在图像的梯度图上进行,经过变换后,分水岭位于梯度最大的位置上。
利用公式(2)作高斯边缘映射,令σ=1.5,得到S(I)的梯度图G(I)(见图4(d))。
f(x,y)=|&[Gσ(x,y)*I(x,y)]|(2)本文采用Vincent[1]中提出的基于FIFO队列的快速分水岭算法,8-邻域搜索。
对G(I)直接运行分水岭变换,得到分割结果V(I),如图4(e)所示,过分割现象非常严重。
使用标记图像M(I)对梯度图G(I)进行反向形态学重建,得到重建结果F(I)。
对F(I)进行分水岭变换,得到初步的分割结果W(I)。
对比W(I)和V(I),不难发现,借助于粗糙标记的重建,过分割现象被显著的消除了,然而,由于某些标记提取的并不精确,过分割现象依然存在,还需要进一步对分水岭区域进行合并。
3基于相似度的区域合并3.1RAG区域合并区域合并算法将分水岭分割图像转换成区域邻接图(Re-gionAdjacencyGraphs,RAG)M,M=(V,E),V是顶点(Vertex)的集合,图的顶点对应被分水岭所划分出来的区域,E是边(Edge)的集合,图的边对应分水岭,它们之间的转换关系如图3所示。
传统的RAG区域合并过程是根据邻接区域的相异度,先合并相异度最小的一对邻接区域,用一个新的点取代,重新计算新生成点的特征,然后更新RAG图,再进行下一次合并,如此迭代,直到满足规定的终止条件。
RAG迭代是较费时的,而且迭代的终止条件很难确定。
3.2基于相似度的快速区域合并立足于借助粗糙标记所产生的分水岭区域的特点,提出了基于相似度的快速区域合并算法,与传统的RAG区域合并不同,该算法比较邻接区域之间的相似度,若相似则直接将小区域并入大区域,从而避免了频繁更新RAG图,提高了算法的执行效率。
设顶点v∈V,v=(labe1,size,gray),其中labe1是区域的编号,size表示区域的大小,gray表示区域的平均灰度值;边e∈E,e=(label1,label2,length,grad),其中label1、label2是与边e所邻接的区域的编号,且label1<label2,length是边界的长度,grad则记录了边e在梯度图G(I)上梯度平均值。
对于具有公共边界e两个邻接区域v1,v2,算法所采用的相似准则是:(1)v1.size<v2.size且v1.size<sizeThre(b)|v1.gray-v2.gray|<grayThre,且e.length<lenThre(c)|v1.gray-v2.gray|<grayThre,且e.grad<gradThre其中,sizeThre、lenThre、grayThre、gradThre是需实验确定的参数,分别表示区域面积最小值,边界长度最小值,灰度差值最小值和边界梯度最小值。
上述三个准则,满足任意一个就可以进行合并。
具体的合并过程如下:(1)初始化阶段:扫描W(I),初始化V和E;(2)排序阶段:按length从大到小,对E进行排序;按size从大到小,对V进行排序;(3)合并阶段:Foreache∈E(3.1)找到e.label1、e.label2所对应的区域v1、v2;(3.2)IFv1、v2相似:(3.2.1)所有与v2邻接的边都并入v1;(3.2.2)从集合V中删除v2;(3.3)从集合E中删除e。
Nexte设图像的宽为M,高为N,W(I)中区域数为R,边界数为Q,则初始化阶段所需的时间为M*N,排序阶段所需时间为R*R+Q*Q,合并阶段所需时间为R*(Q+R),整个合并过程需要的时间T为T=M*N+R*(2*Q+R)+Q*Q(3)由于采用了粗糙标记,R*R<<M*N,R*Q<<M*N,Q*Q<<M*N,故整个合并过程的时间复杂度为O(M*N)。
4实验过程及结果原彩色图像为I,如图4(a)。
令h=5,利用h-dome算法提取标记,提取出的对象标记如图4(b)所示。