CN105913434A - 一种白细胞定位和迭代分割方法 - Google Patents

一种白细胞定位和迭代分割方法 Download PDF

Info

Publication number
CN105913434A
CN105913434A CN201610227867.XA CN201610227867A CN105913434A CN 105913434 A CN105913434 A CN 105913434A CN 201610227867 A CN201610227867 A CN 201610227867A CN 105913434 A CN105913434 A CN 105913434A
Authority
CN
China
Prior art keywords
window
leukocyte
location
anchor
pixel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610227867.XA
Other languages
English (en)
Other versions
CN105913434B (zh
Inventor
曹飞龙
刘月华
黄震
楚建军
赵建伟
周正华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Jiliang University
Maccura Medical Electronics Co Ltd
Original Assignee
Jasdaq Medical Device Co ltd
China Jiliang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jasdaq Medical Device Co ltd, China Jiliang University filed Critical Jasdaq Medical Device Co ltd
Priority to CN201610227867.XA priority Critical patent/CN105913434B/zh
Publication of CN105913434A publication Critical patent/CN105913434A/zh
Application granted granted Critical
Publication of CN105913434B publication Critical patent/CN105913434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Landscapes

  • Image Analysis (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明属于医学图像处理领域,特别涉及一种白细胞定位和迭代分割方法。根据白细胞图片的边界密度信息和颜色信息,计算每一个滑动窗口的边界密度特征因子和颜色特征因子的得分,进而选取得分较高的一部分窗口,之后将这些定位窗口进行处理,整合为最终的白细胞定位窗口,截取白细胞图片中定位窗口的区域即为定位子图。而白细胞的分割是基于定位进行的,用定位子图初始化GrabCut算法,通过膨胀迭代的方法逐步实现白细胞的准确分割。同时,对GrabCut算法的初始图也做了替换像素的预处理,使得能够在一定程度上提高GrabCut的分割精确度。本发明能够有效地避免白细胞直接在原图上分割受到红细胞等其他因素的影响,而且能够产生准确且有效的分割结果。

Description

一种白细胞定位和迭代分割方法
技术领域
本发明属于医学图像处理领域,特别涉及一种白细胞定位和迭代分割方法。
背景技术
白细胞图像分割是白细胞自动识别的一个重要的步骤,是细胞分类计数的基础,甚至于会影响后期的医疗诊断。血细胞图像形状复杂,大小不一,其染色往往不均匀,甚至粘连重叠,加上显微系统的不稳定性和不同的光照条件等,这些因素使得细胞分割成为医学图像处理研究领域中的一个富有挑战性的工作。到目前为止,尽管已有许许多多的学者讨论并提出了多种分割方法,但是,还没有一个足够有效并且快速的算法能很好地适用于白细胞图像的分割。
目前,关于分割的方法主要有几个大类:基于阈值的分割、基于边缘检测的分割、基于区域的分割、基于分类思想的分割、基于活动轮廓模型的分割以及其他混合分割方法。早期的Otsu阈值方法已被广泛应用于白细胞的分割算法中,Cseke在Otsu方法的基础上提出了自动阈值方法,但此方法仅当目标与背景灰度对比度较为明显时有效。基于边缘检测的分割方法,例如Roberts算子,Sobel算子,Prewitt算子,Canny算子等,此类微分算子虽然具有边缘定位准确的优点,但却对噪声敏感,并且容易产生不闭合曲线。Kass等人在1987年提出Snake算法,1988年Osher和Sethian提出了水平集算法,这两种方法都是基 于活动轮廓模型的经典算法,但其迭代运算偏大,且计算结果极易受轮廓线初始位置的影响,对模糊边界和胞浆颗粒大的细胞分割效果不佳。另外,Ko等人在2011年提出了均值漂移和边缘去除逐步融合的算法,但是该算法对细胞质的分割还有待改进,同时该算法复杂,运行时间较长。Mohapatra等人在2012年提出用泛函链接型网络进行基于像素的分割,并利用神经网络训练使相近颜色划分为一类。但此方法需要学习,因此计算量大,同时白细胞染色不同也导致其检测效果泛化性不好等不良效果。
发明内容
本发明的目的是提供一种白细胞定位和迭代分割算法,先根据白细胞特征定位出白细胞子图,然后利用膨胀迭代的GrabCut算法分割出最终的白细胞。
为了实现上述目的,本发明提供了如下技术方案:
本发明提供一种白细胞定位和迭代分割方法,该方法包括如下步骤:
a、计算白细胞图片的边界密度图,设定边界密度滑动窗口w1,在整张白细胞边界密度图中从任意像素点开始,逐像素滑动,直至边界密度滑动窗口w1遍历整张白细胞边界密度图;根据图片的边界密度信息,计算所有边界密度滑动窗口w1的边界密度特征因子的得分值;
b、设定颜色对比度滑动窗口w3,在整张白细胞图中从任意像素点开始,逐像素滑动,直至颜色对比度滑动窗口w3遍历整张白细胞图;根据颜色对比度信息,计算所有颜色对比度滑动窗口w3的颜色对比度特征因子的得分值;
c、根据步骤a和步骤b中边界密度特征因子和颜色对比度特征因子的得分值,对边界密度特征因子和颜色对比度特征因子,分别选取分 值较高的一部分窗口作为定位窗口;
d、将步骤c中得到的所有定位窗口叠加在一起,整合成一个覆盖所有窗口的最小定位窗口,截取白细胞图片中定位窗口的区域作为最终定位子图;
e、对步骤d得到的最终定位子图,用膨胀作差的方法提取细胞核周围的像素值作为细胞质的代表,来替换定位子图中的细胞核像素,进行预处理;
f、选取预处理或非处理的定位子图,来初始化GrabCut算法,然后以膨胀迭代的方式得到白细胞的分割结果。
所述步骤a中,对于白细胞图片,通过canny边缘检测算子检测得到其边界密度;其中,边界密度滑动窗口w1的窗口尺寸为40×40~80×80像素;
边界密度特征因子SED(w1)的计算表示为:
其中,In(w1)是定位窗口w1的内围区域,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θED是缩放参数,L(In(w1))表示In(w1)的周长, 表示In(w1)窗口中所有像素B(p)之和,其中,B为通过canny边缘检测算子得到的边界图,其在每个像素p上的值为0或1。
根据边界密度信息,选取两个具有不同窗口尺寸的边界密度滑动窗口w1和w2,w1和w2的窗口尺寸为40×40~80×80像素,且w1≠w2
所述步骤b中,颜色对比度滑动窗口w3的窗口尺寸为60×60~90×90像素;
采用窗口w3和其外围区域Out(w3)在Lab空间上颜色信息直方图的卡方距离。颜色对比度特征因子SCC计算公式为:
其中,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θCC是缩放参数;
h和g分别是窗口w3和其外围区域Out(w3)的颜色直方图,h(i)是窗口w3的颜色直方图h在第i个组的频数,g(i)是外围区域Out(w3)的颜色直方图g在第i个组的频数,SCC就是两个区域颜色直方图的卡方距离X2(g,h)。
所述步骤c中,对边界密度特征因子和颜色对比度特征因子,均选取得分较高的前3%窗口作为定位窗口。
所述步骤d中,先将步骤c中得到的定位窗口进行去除离群点处理,然后再整合成一个覆盖所有窗口的最小定位窗口;
对分数较高的边界密度特征因子得到的定位窗口,根据窗口之间的欧氏距离的直方图,进行离群点去除处理:根据每个窗口之间的欧氏距离,去除距离较大的窗口;
对分数较高的颜色对比度特征因子得到的定位窗口,根据窗口之间的距离排序,进行离群点去除处理。
所述步骤d中,将所有定位窗口整合成最终的定位窗口的方法为:保留所有的颜色对比度特征因子定位窗口Ci,从所有的边界密度特征因子定位窗口Ei中移除一部分窗口Ek,则最终的定位窗口w*覆盖所有定位窗口w(j)。其过程可以表示为:
w(j)∈Ci∪Ei\Ek,j=1,2,3,...,;
其中,μ是阈值参数。
所述步骤e中,将细胞核像素替换为细胞核周围的细胞质代表像素,其中,细胞质是通过膨胀细胞核,并提取细胞核周围像素得到的;细胞质区域NS可以表示为:
NS=NUdil-NU;
其中NU、NUdil和NS分别代表细胞核、细胞核膨胀之后的区域和细胞质区域。
所述步骤f中,选取第一次GrabCut分割结果区域面积较大的预处理或非处理的定位子图:
其中Is代表用于初始化GrabCut算法的定位子图,Ip_s和Inp_s分别代表预处理和非处理的定位子图,代表迭代一次GrabCut的分割结果,分别表示的像素个数。
所述步骤f中,每迭代一次都对分割结果膨胀一次,再进行下一次迭代,直至迭代前后分割结果不变,则迭代停止。
与现有技术相比,本发明的有益效果在于:
由于白细胞自身复杂特性、红细胞和血小板的干扰以及染色不均与光照亮度等的影响,背景技术中提及的已有方法,总是不能完全考虑这些因素,也不能解决这些因素带来的不良影响。而本发明则考虑通过先定位的方法去除红细胞和血小板的干扰,同时使得染色不均和 光照亮度所带来的影响随之减弱,然后再嵌入GrabCut算法,成功地分割出白细胞。
综上所述,白细胞直接在原图上分割总是不可避免地受到红细胞等其他因素的影响,不能产生准确且有效的分割结果,而本发明的方法能够有效地避免了这些问题。
附图说明
图1a为白细胞图片;
图1b为白细胞边界密度图;
图1c为效果较好的边界密度特征因子分数较高定位窗口图;
图1d为效果较好的颜色对比度特征因子分数较高定位窗口图;
图2a为效果不好的边界密度特征因子分数较高定位窗口图;
图2b为窗口之间的欧氏距离的直方图;
图2c为边界密度特征因子定位窗口处理图;
图3a为效果不好的颜色对比度特征因子分数较高定位窗口图;
图3b为窗口之间的距离排序图;
图3c为颜色对比度特征因子定位窗口处理图;
图4a为所有定位窗口整合图;
图4b为最终定位子图;
图5为根据定位子图进行分割白细胞的流程图;
图6为本发明白细胞定位和迭代分割方法总体流程示意图。
其中的附图标记为:
(a)最终定位子图
(b)预处理定位图
(c)初始化GrabCut算法的初始图
(d)最终白细胞分割结果
具体实施方式
本发明基于白细胞先定位后分割,从而减少红细胞、血小板和其他背景对白细胞分割的影响。针对白细胞图片独有的特性,以密度特征和颜色特征来定位白细胞,又采用多尺度和多窗口的思想实现定位窗口的自适应调整。在分割过程中,对定位子图的预处理在一定程度上提高了GrabCut算法分割的精确度,同时以膨胀迭代的方式,分割出更加准确的结果。
本发明的基本思路为:
首先,根据白细胞图片的边界密度信息和颜色信息,计算边界密度特征因子和颜色对比度特征因子的得分,分别选取得分较高的一部分窗口;
其次,对由边界密度特征因子和颜色对比度特征因子得到的窗口进行去除离群点处理,并整合成覆盖所有窗口的最小定位窗口,窗口内部区域作为最终定位子图;
再次,用像素替换方法对最终定位子图进行预处理;
最后,引入GrabCut算法进行膨胀迭代,得到白细胞分割结果。
下面结合实施例对本发明进行进一步说明。
本发明提供一种白细胞定位和迭代分割方法,包括如下步骤:
a、计算白细胞图片的边界密度图,设定边界密度滑动窗口w1,在整张白细胞边界密度图中从任意像素点开始,逐像素滑动,直至边界密度滑动窗口w1遍历整张白细胞边界密度图;根据图片的边界密度信息,计算所有边界密度滑动窗口w1的边界密度特征因子的得分值。
优选地,在整张白细胞边界密度图中从左上角开始,逐像素滑动,直至边界密度滑动窗口w1遍历整张白细胞边界密度图。
如图1a所示,对于一张白细胞图片,通过canny边缘检测算子检测得到其边界密度图,如图1b所示;根据边界密度信息,计算所有边界密度滑动窗口w1的边界密度因子得分值,其中,边界密度滑动窗口w1的窗口尺寸为40×40~80×80像素。
边界密度特征因子SED(w1)的计算表示为:
其中,In(w1)是定位窗口w1的内围区域,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θED是缩放参数;
L(In(w1))表示In(w1)的周长,表示In(w1)窗口中所有像素B(p)之和,其中,B为通过canny边缘检测算子得到的边界图,其在每个像素p上的值为0或1。
图中的边界区域,即白色像素值为1;非边界区域,即黑色像素值为0。所以边界越密集的地方,也就是白细胞所在区域,SED(w1)得分值越大。因此得分大的区域也就是定位白细胞的窗口。
优选地,根据边界密度信息,选取两个具有不同窗口尺寸的边界密度滑动窗口w1和w2,w1和w2的窗口尺寸为40×40~80×80像素,且w1≠w2
b、设定颜色对比度滑动窗口w3,在整张白细胞图中从任意像素点开始,逐像素滑动,直至颜色对比度滑动窗口w3遍历整张白细胞图;根据颜色对比度信息,计算所有颜色对比度滑动窗口w3的颜色对比度特征因子得分值。
优选地,在整张白细胞图中从左上角开始,逐像素滑动,直至颜色对比度滑动窗口w3遍历整张白细胞图。
其中,颜色对比度滑动窗口w3的窗口尺寸为60×60~90×90像素。
采用窗口w3和其外围区域Out(w3)在Lab空间上颜色信息直方图的卡方距离。颜色对比度特征因子SCC计算公式为:
其中,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θCC是缩放参数;
h和g分别是窗口w3和其外围区域Out(w3)的颜色直方图,h(i)是窗口w3的颜色直方图h在第i个组的频数,g(i)是外围区域Out(w3)的颜色直方图g在第i个组的频数,SCC就是两个区域颜色直方图的卡方距离X2(g,h)。
c、根据步骤a和步骤b中边界密度特征因子和颜色对比度特征因子的得分值,对边界密度特征因子和颜色对比度特征因子,分别选取分值较高的一部分窗口作为定位窗口。
选取边界密度特征因子得分值较高的窗口作为定位窗口;如图1c和图2a所示。其中,图1c为效果较好的边界密度特征因子分数较高定位窗口图;图2a为效果不好的边界密度特征因子分数较高定位窗口图。
选取颜色对比度特征因子得分值较高的窗口作为定位窗口;如图1d和图3a所示。其中,图1d为效果较好的颜色对比度特征因子分数较高定位窗口图;图3a为效果不好的颜色对比度特征因子分数较高定位窗口图。
优选地,针对边界密度特征因子和颜色对比度特征因子,均选取得分较高的前3%窗口作为定位窗口。
d、将步骤c中得到的所有定位窗口叠加在一起,整合成一个覆盖 所有窗口的最小定位窗口,截取白细胞图片中定位窗口的区域作为最终定位子图。
优选地,考虑到步骤c中存在个别白细胞定位效果不好的情况,先将步骤c中得到的定位窗口进行去除离群点处理,然后再整合成一个覆盖所有窗口的最小定位窗口。
对效果不好的分数较高的边界密度特征因子得到的定位窗口,如图2a所示,根据窗口之间的欧氏距离的直方图,如图2b所示,进行离群点去除处理,其中图中的虚线即为阈值,从而得到处理之后的窗口,如图2c所示。
所述离群点去除处理为根据每个窗口之间的欧氏距离,去除距离较大的窗口。
类似地,对效果不好的分数较高的颜色对比度特征因子得到的定位窗口,如图3a所示,根据窗口之间的距离排序,如图3b所示,得到处理之后的窗口,如图3c所示。
整合如图2c和图3c所示处理之后的所有定位窗口,如图4a,得到一个覆盖所有窗口的最小的最终定位窗口,在白细胞图片中截取其区域即为如图4b所示的最终定位子图。整合过程为:
保留所有的颜色对比度特征因子定位窗口Ci,从所有的边界密度特征因子定位窗口Ei中移除一部分窗口Ek,则最终的定位窗口w*覆盖所有定位窗口w(j)。其过程可以表示为:
w(j)∈Ci∪Ei\Ek,j=1,2,3,...,
其中μ是阈值参数。
e、对步骤d得到的最终定位子图,用膨胀作差的方法提取细胞核周围的像素值作为细胞质的代表,来替换其中的细胞核像素,进行定位子图的预处理。
根据步骤d中得到的最终定位子图,进行替换像素的预处理,得到如图5中(b)所示的预处理定位子图。即将细胞核像素替换为细胞核周围的细胞质代表像素,其中,细胞质是通过膨胀细胞核,并提取细胞核周围像素得到的。这样能进一步提高GrabCut算法的准确性,其中细胞质代表区域NS可以表示为:
NS=NUdil-NU
其中NU,NUdil和NS分别代表细胞核、细胞核膨胀之后的区域和细胞质区域;
f、选取预处理或非处理的定位子图,来初始化GrabCut算法,然后以膨胀迭代的方式得到白细胞的分割结果。
选取首次GrabCut算法分割区域面积较大的所对应的预处理或非处理的定位子图,如图5中(c)所示:
其中Is代表用于初始化GrabCut算法的定位子图,Ip_s和Inp_s分别代表预处理和非处理定位子图,代表迭代一次GrabCut的分割结果,分别表示的像素个数。
将所选定位子图作为GrabCut算法的初始图进行分割。为了分割结果的精确性,每一次GrabCut运行之后,将分割结果膨胀一次,然后在输入GrabCut的下一次迭代分割中,直至迭代前后分割结果不变,则迭代停止,如图5中(d)所示。
如图6所示,为本发明的白细胞定位和迭代分割方法总体流程示 意图。作为一个实施例,对于一张白细胞图片,设定边界密度滑动窗口w1的窗口尺寸为50×50像素,边界密度滑动窗口w2的窗口尺寸为70×70像素,颜色对比度滑动窗口w3的窗口尺寸为90×90像素;然后计算所有边界密度特征因子和颜色对比度特征因子的得分值,均选取得分较高的前3%窗口作为定位窗口,即得到图6中的50×50大小的ED定位窗口ED、90×90大小的CC定位窗口CC和70×70大小的ED定位窗口ED;然后将得到的所有定位窗口叠加在一起,整合成一个覆盖所有窗口的最小定位窗口,如图6中④所示;所有定位窗口叠加在一起,如图6中⑤所示,多窗口融合成最终定位窗口。对最终定位子图,用膨胀作差的方法提取细胞核周围的像素值作为细胞质的代表,来替换其中的细胞核像素,进行定位子图的预处理;选取预处理或非处理的定位子图,来初始化GrabCut算法,然后以膨胀迭代的方式得到白细胞的分割结果,每一次GrabCut运行之后,将分割结果膨胀一次,然后在输入GrabCut的下一次迭代分割中,直至迭代前后分割结果不变,则迭代停止。

Claims (10)

1.一种白细胞定位和迭代分割方法,其特征在于:该方法包括如下步骤:
a、计算白细胞图片的边界密度图,设定边界密度滑动窗口w1,在整张白细胞边界密度图中从任意像素点开始,逐像素滑动,直至边界密度滑动窗口w1遍历整张白细胞边界密度图;根据图片的边界密度信息,计算所有边界密度滑动窗口w1的边界密度特征因子的得分值;
b、设定颜色对比度滑动窗口w3,在整张白细胞图中从任意像素点开始,逐像素滑动,直至颜色对比度滑动窗口w3遍历整张白细胞图;根据颜色对比度信息,计算所有颜色对比度滑动窗口w3的颜色对比度特征因子的得分值;
c、根据步骤a和步骤b中边界密度特征因子和颜色对比度特征因子的得分值,对边界密度特征因子和颜色对比度特征因子,分别选取分值较高的一部分窗口作为定位窗口;
d、将步骤c中得到的所有定位窗口叠加在一起,整合成一个覆盖所有窗口的最小定位窗口,截取白细胞图片中定位窗口的区域作为最终定位子图;
e、对步骤d得到的最终定位子图,用膨胀作差的方法提取细胞核周围的像素值作为细胞质的代表,来替换定位子图中的细胞核像素,进行预处理;
f、选取预处理或非处理的定位子图,来初始化GrabCut算法,然后以膨胀迭代的方式得到白细胞的分割结果。
2.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤a中,对于白细胞图片,通过canny边缘检测算子检测得到其边界密度;其中,边界密度滑动窗口w1的窗口尺寸为40×40~80×80像素;
边界密度特征因子SED(w1)的计算表示为:
S E D ( w 1 ) = Σ p ∈ I n ( w 1 ) B ( p ) L ( I n ( w 1 ) ) ;
| I n ( w 1 ) | | w 1 | = 1 θ E D 2 ;
其中,In(w1)是定位窗口w1的内围区域,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θED是缩放参数,L(In(w1))表示In(w1)的周长,表示In(w1)窗口中所有像素B(p)之和,其中,B为通过canny边缘检测算子得到的边界图,其在每个像素p上的值为0或1。
3.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:根据边界密度信息,选取两个具有不同窗口尺寸的边界密度滑动窗口w1和w2,w1和w2的窗口尺寸为40×40~80×80像素,且w1≠w2
4.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤b中,颜色对比度滑动窗口w3的窗口尺寸为60×60~90×90像素;
采用窗口w3和其外围区域Out(w3)在Lab空间上颜色信息直方图的卡方距离。颜色对比度特征因子SCC计算公式为:
S C C = χ 2 ( g , h ) = 1 2 Σ i ( g ( i ) - h ( i ) ) 2 g ( i ) + h ( i ) ;
| O u t ( w 3 ) | | w 3 | = θ C C 2 - 1 ;
其中,|In(w1)|和|w1|分别表示In(w1)和w1的像素个数,θCC是缩放参数;
h和g分别是窗口w3和其外围区域Out(w3)的颜色直方图,h(i)是窗口w3的颜色直方图h在第i个组的频数,g(i)是外围区域Out(w3)的颜色直方图g在第i个组的频数,SCC就是两个区域颜色直方图的卡方距离X2(g,h)。
5.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤c中,对边界密度特征因子和颜色对比度特征因子,均选取得分较高的前3%窗口作为定位窗口。
6.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤d中,先将步骤c中得到的定位窗口进行去除离群点处理,然后再整合成一个覆盖所有窗口的最小定位窗口;
对分数较高的边界密度特征因子得到的定位窗口,根据窗口之间的欧氏距离的直方图,进行离群点去除处理:根据每个窗口之间的欧氏距离,去除距离较大的窗口;
对分数较高的颜色对比度特征因子得到的定位窗口,根据窗口之间的距离排序,进行离群点去除处理。
7.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤d中,将所有定位窗口整合成最终的定位窗口的方法为:保留所有的颜色对比度特征因子定位窗口Ci,从所有的边界密度特征因子定位窗口Ei中移除一部分窗口Ek,则最终的定位窗口w*覆盖所有定位窗口w(j)。其过程可以表示为:
I = ∪ i C i ;
{ E k &Element; E i , ( i = 1 , 2 , 3 , ... , ) : | E k &cap; I | | I | < &mu; } ;
w(j)∈Ci∪Ei\Ek,j=1,2,3,...,;
其中,μ是阈值参数。
8.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤e中,将细胞核像素替换为细胞核周围的细胞质代表像素,其中,细胞质是通过膨胀细胞核,并提取细胞核周围像素得到的;细胞质区域NS可以表示为:
NS=NUdil-NU;
其中NU、NUdil和NS分别代表细胞核、细胞核膨胀之后的区域和细胞质区域。
9.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤f中,选取第一次GrabCut分割结果区域面积较大的预处理或非处理的定位子图:
I s = I n p _ s , i f | I n p _ s ( 1 ) | &GreaterEqual; | I p _ s ( 1 ) | I p _ s , i f | I n p _ s ( 1 ) | < | I p _ s ( 1 ) | ;
其中Is代表用于初始化GrabCut算法的定位子图,Ip_s和Inp_s分别代表预处理和非处理的定位子图,代表迭代一次GrabCut的分割结果,分别表示的像素个数。
10.根据权利要求1所述的白细胞定位和迭代分割方法,其特征在于:所述步骤f中,每迭代一次都对分割结果膨胀一次,再进行下一次迭代,直至迭代前后分割结果不变,则迭代停止。
CN201610227867.XA 2016-04-13 2016-04-13 一种白细胞定位和迭代分割方法 Active CN105913434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610227867.XA CN105913434B (zh) 2016-04-13 2016-04-13 一种白细胞定位和迭代分割方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610227867.XA CN105913434B (zh) 2016-04-13 2016-04-13 一种白细胞定位和迭代分割方法

Publications (2)

Publication Number Publication Date
CN105913434A true CN105913434A (zh) 2016-08-31
CN105913434B CN105913434B (zh) 2018-07-17

Family

ID=56745926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610227867.XA Active CN105913434B (zh) 2016-04-13 2016-04-13 一种白细胞定位和迭代分割方法

Country Status (1)

Country Link
CN (1) CN105913434B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106780522A (zh) * 2016-12-23 2017-05-31 杭州华卓信息科技有限公司 一种基于深度学习的骨髓液细胞分割方法
CN110059656A (zh) * 2019-04-25 2019-07-26 山东师范大学 基于卷积对抗生成神经网络的白细胞分类方法及系统
CN112330690A (zh) * 2020-11-02 2021-02-05 腾讯科技(深圳)有限公司 基于显微图像的细胞分割方法、装置、设备及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732524A (zh) * 2015-02-07 2015-06-24 嘉善加斯戴克医疗器械有限公司 血液白细胞显微图像的随机权网络分割方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732524A (zh) * 2015-02-07 2015-06-24 嘉善加斯戴克医疗器械有限公司 血液白细胞显微图像的随机权网络分割方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ALEXE, B 等: "Measuring the Objectness of Image Windows", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》 *
贺锦鹏: "基于图论的彩色纹理图像分割技术研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106780522A (zh) * 2016-12-23 2017-05-31 杭州华卓信息科技有限公司 一种基于深度学习的骨髓液细胞分割方法
CN106780522B (zh) * 2016-12-23 2019-08-30 杭州华卓信息科技有限公司 一种基于深度学习的骨髓液细胞分割方法
CN110059656A (zh) * 2019-04-25 2019-07-26 山东师范大学 基于卷积对抗生成神经网络的白细胞分类方法及系统
CN112330690A (zh) * 2020-11-02 2021-02-05 腾讯科技(深圳)有限公司 基于显微图像的细胞分割方法、装置、设备及存储介质
CN112330690B (zh) * 2020-11-02 2022-08-05 腾讯医疗健康(深圳)有限公司 基于显微图像的细胞分割方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN105913434B (zh) 2018-07-17

Similar Documents

Publication Publication Date Title
CN110533684B (zh) 一种染色体核型图像切割方法
US10438096B2 (en) Identifying and excluding blurred areas of images of stained tissue to improve cancer scoring
CN103942794B (zh) 一种基于置信度的图像协同抠图方法
CN103473739B (zh) 一种基于支持向量机的白细胞图像精确分割方法与系统
CN106651872B (zh) 基于Prewitt算子的路面裂缝识别方法及系统
CN106056118B (zh) 一种用于细胞的识别计数方法
Nosrati et al. Segmentation of overlapping cervical cells: A variational method with star-shape prior
CN103077521B (zh) 一种用于视频监控的感兴趣区域提取方法
CN103810491B (zh) 融合深度和灰度图像特征点的头姿估计兴趣点检测方法
CN103198315B (zh) 基于字符轮廓和模板匹配的车牌字符分割算法
CN104851085B (zh) 自动获取图像中检测区的方法及系统
CN104574418A (zh) 基于神经网络的压力容器焊缝缺陷识别方法及其装置
CN106650553A (zh) 车牌识别方法及系统
CN105844252A (zh) 一种面部关键部位的疲劳检测方法
CN104050471A (zh) 一种自然场景文字检测方法及系统
CN105761260B (zh) 一种皮肤图像患处分割方法
CN110400287B (zh) 结直肠癌ihc染色图像肿瘤侵袭边缘和中心的检测系统及方法
Chen et al. Facial expression recognition based on edge detection
CN104463138A (zh) 基于视觉结构属性的文本定位方法及系统
CN104658003A (zh) 舌体图像的分割方法和装置
CN104598907A (zh) 一种基于笔画宽度图的图像中文字数据提取方法
CN105913434A (zh) 一种白细胞定位和迭代分割方法
CN103049756A (zh) 基于cema算法和纹理块匹配修复技术的自动提取与去除彩色图像内文字的方法
CN109508720A (zh) 一种车辆颜色识别方法及装置
CN104834890A (zh) 一种对书法作品中文字神采信息的提取方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20171102

Address after: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park source Street No. 258

Applicant after: CHINA JILIANG UNIVERSITY

Applicant after: Mike Medical Electronics Ltd

Address before: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park source Street No. 258

Applicant before: China Jiliang University

Applicant before: Jiashan Just-Tech Medical Equipment Co., Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant