一种图像中的粘连细胞分割方法技术

本发明专利技术提出一种图像中的粘连细胞分割方法，首先对图像进行二值化，接着对二值图做距离变换，得到各点的距离值。再在距离图中检测局部极值点，筛选过后得到代表细胞的中心点。然后对距离图做水域分割，根据各子块中心与细胞中心的距离确定子块与细胞的隶属关系。最后对毗邻多个子块的中间块进行再分割，使分割结果更为精细。本方法有效地解决了传统分水岭分割中易出现的过分割及误分割现象，对多细胞粘连有较好的适用性。

【技术实现步骤摘要】

本专利技术涉及图像中的粘连细胞的分割方法，应用于分离图像中粘连的多个细 胞。
技术介绍
在血液的染色图片中，常常会遇到细胞粘连的情况，为准确获取每个细胞的形 态信息，我们必须对这些粘连的细胞进行分离。国内外在粘连细胞的分离方面已有一些 研究，目前采用的方法主要有以下几种1.基于形态学的分水岭分割方法分水岭分割方法，是一种基于拓扑理论的数学形态学的分割方法，其基本思想 是把图像看作是测地学上的拓扑地貌，图像中每一点像素的灰度值表示该点的海拔高 度，每一个局部极小值及其影响区域称为集水盆，而集水盆的边界则形成分水岭。分水 岭的概念和形成可以通过模拟浸入过程来说明。在每一个局部极小值表面，刺穿一个小 孔，然后把整个模型慢慢浸入水中，随着浸入的加深，每一个局部极小值的影响域慢慢 向外扩展，在两个集水盆汇合处构筑大坝，即形成分水岭。形态学分水岭算法的优点是边界连续，在种子点选取较好的情况下，分割结果 稳定。由于种子点的选取是通过连续腐蚀或者通过先做距离变换再通过某个阈值来选取 的，该算法对微弱边缘响应良好，但易受噪声影响，存在过分割和误分割问题，迭代膨 胀计算量较大，适用于形态无明显凹陷，噪声较小，细胞外形规则的情况。2.基于边界角点配对的分割方法边界角点配对法，是将图像中的物体边界用封闭的连续像点表示，通过分析边 界上各点与相邻点的关系得到该点所在边界区段的走向，即凹入，凸起或趋于平滑，以 一阈值限定得到所有凹点。将所有凹点按照一定规则进行匹配后得到分割结果。该方法直观，简单，且实现快速。但凹点的选择基于轮廓，对轮廓上毛刺和噪 声敏感，且存在凹点不能精确求取，凹点配对的困难。3.基于类圆的分割方法基于类圆的分割方法以血细胞形状呈类圆为前提，首先寻找物体边界上距离最 远的两点a和b，以a和b距离的一半作为外接圆的半径，a与b连线的中点作为外接圆的 圆心画圆，检测物体边界上距离圆心最近的两点，如果该两点的连线满足一定的条件， 则连接两点分割细胞，并继续递归分割，否则分割结束。该方法的优点是充分利用了血细胞的先验信息，对于圆形细胞分割效果好。但 在实际操作中，由于细胞形状呈多样性，且切割条件判定复杂，常使分割出现较大偏差。
技术实现思路
本专利技术的目的在于，针对现有方法的不足提出一种粘连细胞分割方法，通过对过分割区域进行重分割，避免传统分水岭分割带来的过分割与误分割。本专利技术提出的方案如下步骤一，对图像进行二值化，将图像中的粘连细胞作为前景色进行标记，其余 的为背景色，得到二值图。对上述二值图做距离变换，在距离变换的过程中记录下各个 点到最近背景点的距离值，保存为距离图。步骤二，确定真实细胞的中心点首先，在距离图中搜索局部极大值点，即距离值大于其八邻域点的距离值的 点，以这些局部极大值点为起始点对二值图作水域分割，得到分割后的区域图。其次，对所有的局部极大值点进行筛选，去除符合筛选条件的局部极大值点， 剩余的极大值点每个均代表单个血细胞的中心点。其中，筛选的条件为(1)该局部极 大值点的距离值低于门限阈值D ； (2)该局部极大值点的距离值小于任一阈值点的距离 值，其中所述阈值点指与当前局部极大值点的距离低于门限阈值T的局部极大值点。步骤三，计算区域图中各区域中心与各中心点的距离，获得各区域与各单个血 细胞的距离值，使每个区域归属于与之距离值最小的血细胞，完成血细胞的基本分割。 其中，所述的区域中心指区域中距离值最大的点。步骤四，对分割边界进一步修正。检测与两个或两个以上区域相邻的区域，在 该区域的边界上搜索代表点，连接代表点对当前区域重新分割，分割得到的区域分属于 与各自相邻的区域。本专利技术能够达到的有益效果如下由于血细胞形态各异，其距离图中存在较多伪极值点，这些点易造成传统分水 岭分割中的过分割现象，解决过分割的基本方法是简单地进行区域合并，但原本由伪极 值点分割得到的区域不论归属于哪个区域都将造成误分割。而本专利技术提出的方法通过筛 选极值点可以有效地避免过分割。消除过分割后，对过分割区域进行重分割，分割结果 更符合人眼观测的结果。附图说明图1 3x3倒角距离变换模板示意2 八连通链码模板示意图具体实施例方式下面结合附图和具体实施例对本专利技术作进一步说明。步骤一，将含有粘连细胞的图像进行阈值分割，将图像中的粘连细胞作为前景 色进行标记，其余的为背景色，得到二值图f(x，y)，对该二值图f(x，y)做倒角距离变 换，得到图中各点的距离值，保存为距离图d(x，y)。二值图数据为宽w，高h的0-1型矩阵，其中0表示背景点，1表示前景点。步骤二，确定真实细胞的中心点(1)在距离图d(x，y)中搜索 局部极大值点(局部极大值点指该点的距离值大于 其八邻域各点的距离值)，将其坐标值保存于数组p，其中N为检测到的局部极大值点 的个数。以这些局部极大值点为起始点对二值图f(x，y)进行水域分割，得到区域块图area, M为水域分割得到的总区域个数。 (2)根据以下两个原则筛选局部极大值点依次检测数组p中的各点，当满足式(1)时，去除当前点。d(p) <D, 0<t1 t2 < N (1)其中，D是阈值，为常数，本实施例中取为常见细胞半径值的一半。然后对数组p中剩余的点进行两两配对，如果两点距离小于门限阈值T则继 续进行判断，如果满足式(2)则去除距离值较小的点，保留距离值较大的点。其中，T 为常数，本实施例中取为常见细胞直径的三分之一。UpttJ-pMII <min ((KpCt1]), d(p)), 0<t1 t2 < N，I1 ^ t2 (2)p中最终保留的为真实细胞的中心点。步骤三，在区域块图area中寻找每个区域块中距离值最大的点作为该区域的 中心点。对每个区域中心点依次进行处理在p中寻找与当前区域中心点距离最近的 细胞中心点，并将当前区域标记为该细胞的一部分。所有区域计算完毕后，得到初步的 粘连细胞分割图。步骤四，修正分割边界。首先，寻找毗邻两个或多个细胞的区域。根据区域块图area统计各区域的边 界点所毗邻的细胞编号，当细胞编号多于一个且区域内不包含细胞中心点时，即为毗邻 两个或多个细胞的区域，标记该区域。其次，对于已标记的区域，计算其边界链码差，计算步骤如下(1)用链码表示区域边界按逆时针方向搜索区域边界，起始点链码值标记为0，其余各点i的链码值由当 前点C与其前一点c的相对关系决定。当C相对于c作逆时针旋转时，按 偏角的大小取值0 7，分别对应逆时针旋转0°、45°、90°、135°、180°、225°、 270°、315°。如图2所示，箭头起点为前一点C，箭头终点为当前点C。将各 点链码值保存于数组C中，Q为边界像素点数。(2)计算相对链码对于边界上的任意一点i都有前一点指向此点的链码C和此点指向后一点的链 码C，相对链码是指C和C的相互关系。C相对于C作逆时针旋转时， 按偏角的大小取值1 3，分别对应逆时针旋转45°、90°、135°作顺时针旋转时，相 对链码取值-1 -4，分别对应顺时针旋转45°、90°、135°、180°。式(3)为相对 链码R的计算公式。R = (C-C+8)MOD8(3)R = R_8，如果 R > 3(3)计算绝对链码绝对链码则是指从本文档来自技高网...

【技术保护点】
一种图像中的粘连细胞分割方法，用于分离图像中粘连的多个细胞，该方法包括如下步骤：步骤一，对含有粘连细胞的图像进行二值化，将图像中的粘连细胞作为前景色进行标记，其余的为背景色，得到二值图，再对上述二值图做距离变换，在距离变换的过程中记录下各个点到最近背景点的距离值，保存为距离图；步骤二，确定真实细胞的中心点首先，在所述距离图中搜索局部极大值点，即距离值大于其八邻域点的距离值的点，以这些局部极大值点为起始点对所述二值图作水域分割，得到分割后的区域图；其次，对所有的局部极大值点进行筛选，去除符合筛选条件的局部极大值点，剩余的极大值点每个均代表单个血细胞的中心点，其中，所述筛选的条件为：（１）该局部极大值点的距离值低于预设的阈值Ｄ；（２）该局部极大值点的距离值小于任一阈值点的距离值，其中所述阈值点指与当前局部极大值点的距离低于门限阈值Ｔ的局部极大值点；步骤三，计算区域图中各区域中心与各上述中心点的距离，从而获得各区域与各单个血细胞的距离值，使每个区域归属于与之距离值最小的血细胞，完成血细胞的初步分割，其中，所述的区域中心指区域中距离值最大的点；步骤四，对初步分割的边界进一步修正检测获取与两个或两个以上区域毗邻的区域，在上述获取的毗邻区域的边界上搜索用于分割的连接角点，利用该连接角点对当前区域重新分割，分割得到的区域分属于与各自相邻的区域。通过上述步骤，即完成粘连细胞的分割。...

【技术特征摘要】

【专利技术属性】
技术研发人员：汪国有，林晨，郑馨，刘建国，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：83