CN104463892A - 基于水平集及GVF Snake精确定位的菌落图像分割方法 - Google Patents

基于水平集及GVF Snake精确定位的菌落图像分割方法 Download PDF

Info

Publication number
CN104463892A
CN104463892A CN201410813281.2A CN201410813281A CN104463892A CN 104463892 A CN104463892 A CN 104463892A CN 201410813281 A CN201410813281 A CN 201410813281A CN 104463892 A CN104463892 A CN 104463892A
Authority
CN
China
Prior art keywords
bacterium colony
image
area
bacterial colony
threshold
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
CN201410813281.2A
Other languages
English (en)
Other versions
CN104463892B (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.)
Fuzhou University
Original Assignee
Fuzhou 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 Fuzhou University filed Critical Fuzhou University
Priority to CN201410813281.2A priority Critical patent/CN104463892B/zh
Publication of CN104463892A publication Critical patent/CN104463892A/zh
Application granted granted Critical
Publication of CN104463892B publication Critical patent/CN104463892B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30242—Counting objects in image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种基于水平集及GVF Snake精确定位的菌落图像分割方法，根据菌落图像灰度变化快慢，区域灰度变化相似特性，采用基于图像灰度等高线及灰度变化加速度进行图像分割，即采用基于图像灰度等高线及按区域顺序合并来进行图像初分割，进而基于菌落的边界扫描和边界线连接确定菌落初步轮廓，再用GVF Snake算法对菌落轮廓进行精确地定位，然后对菌落中的孔洞进行充填和对菌落边界进行拟合，得到最终图像分割结果。本发明所提出的一种基于水平集及GVF Snake精确定位的菌落图像分割方法，不同于常规的直接利用相邻像素间灰度值的相似性进行图像分割的方法，有效的解决了菌落图像分割困难的问题。

Description

基于水平集及GVF Snake精确定位的菌落图像分割方法
技术领域
本发明涉及图像处理，特别是一种基于水平集及GVF Snake精确定位的菌落图像分割方法。
背景技术
传统评判水和食物质量的传统和标准的方法是将试样均匀分布地散放在一标准的玻璃盘内，放入一定温度的冷藏室内，放置一定时间后，将玻璃盘进行人工计数细菌群的个数及分析其它有关参数。而人工计数及分析方法存在以下几方面的不利因数：(1)计数的局限性：如在一直径为九厘米的玻璃盘内，细菌群的数目超过一百时，人工计数的准确性将大大降低，一般都是取二分之一到三十二分之一盘来计数，然后进行整盘的细菌群数目估算；(2)无法对细菌群参数进行定量地分析，如形状尺才表面颜色等；(3)由于细菌群是随时间而变化的，无法隔天甚至隔时复查；(4) 沉长和重复性的计数影响计数员的视力；(5)刺痹的气味影响操作员的身体健康；(6)计数速度也有待提高；(7)计数的不稳定性：例如对同盘细菌群的计数，不同的计数员往往给出不同的计数结果。
为了克服以上不足，提高计数的精度，增强定量分析的能力和提高自动化程度，近年来出现的最佳方法之一是基于图像分析的方法。如北美和欧洲一些国家相继地基于计算机图像的技术，研究和发展了计算机图像计数分析系统，但该类系统仍未被广泛应用。当然有其各种各样的原因，但其主要原因还是系统本身的稳定性和精度问题。
发明内容
本发明的目的在于提供一种基于图像灰度等高线、区域膨胀、一阶微分边界扫描及GVF Snake算法对菌落轮廓进行精确定位的图像分割方法，以解决菌落图像分割困难的问题。
为实现上述目的，本发明的技术方案是：一种基于水平集及GVF Snake精确定位的菌落图像分割方法，其特征在于，按照如下步骤实现：
S1：采集并输入菌落图像，通过高斯算子对该菌落图像进行平滑滤波操作，获取平滑图像；
采用大津法对该平滑图像进行阈值分割，获得第一阈值T1，并获取由该第一阈值T1确定的第一等高线和第一区域；
S2：统计所述平滑图像的灰度直方图，检测出该灰度直方图种主波峰的左谷点灰度值，并在该灰度值与所述第一阈值T1的差值区域中，采用大津法进行阈值分割，获得第二阈值T2，并获取由该第二阈值T2确定的第二等高线和第二区域；
S3：用一阶微分边界扫描或山脊边界扫描算法获取菌落边界；
S4：对所述第一等高线向所述第二等高线方向进行加速度判定操作，对所述第一区域向所述第二区域进行扩展操作；将扩展操作后处于所述第一区域与所述第二区域交集区域中，且与所述第一区域对应的菌落目标种子点进行合并，扩大种子区域；继续对所述第一区域向所述第二区域进行扩展操作，进一步扩大所述种子区域；
S5：对所述菌落边界进行细线化操作和端点检测操作，并根据该菌落边界中端点的方向及端点之间的距离对应连接端点，以使菌落边界闭合，消除该菌落边界中存在的断点和空隙；
S6：以所述步骤S5所得的菌落边界作为约束条件，对所述种子区域进行膨胀操作，再通过GVF Snake算法对该菌落边界进行精确地定位；
S7：对所述步骤S6中完成菌落边界精确定位后的菌落图像中的孔洞进行充填，并用二次曲线对菌落边界进行拟合操作，得到最终图像分割结果。
进一步的，所述步骤S3还包括：对所述平滑图像进行一阶微分获得梯度图像，并将梯度图像转换成二值图像。
相较于现有技术，本发明具有以下有益效果：本发明所提出的一种基于水平集及GVF Snake精确定位的菌落图像分割方法，对于从中央到边缘亮度逐渐减弱的菌落目标以及对应作图像反值运算后的亮度逐渐增强的菌落目标，该方法能够快速准确地提取菌落目标，从而提高菌落分析的精度。
附图说明
图1为本发明中基于水平集及GVF Snake精确定位的菌落图像分割方法的流程图。
图2（a）为本发明一实施例中理想的椭圆形状的菌落示意图。
图2（b）为本发明一实施例中将理想的椭圆形状的菌落经处理后得到的灰度值等高线示意图。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
本发明提供一种基于水平集及GVF Snake精确定位的菌落图像分割方法，如图1所示，其特征在于，按照如下步骤实现：
S1：采集并输入菌落图像，通过高斯算子对该菌落图像进行平滑滤波操作，去除所述图像中的噪声，获取平滑图像；在本实施例中，将该过程表示为：                                                ，其中h表示平滑滤波器，为原始图像，平滑滤波器采用高斯滤波器：；采用大津法或迭代法对该平滑图像进行阈值分割，获得第一阈值T1，并获取由该第一阈值T1确定的第一等高线和第一区域；且在本实施例中，将大津法作如下定义：
二维直方图中存在两类和，分别代表物体和背景，具有两个不同的概率密度函数分布设阈值为，那么两类出现的概率分别为：
；
且该两类对应的均值矢量为：
其中：  。
二维直方图上总的均值矢量：
，
 而，，使用其迹作为类间的离散度测度，则有：，最佳阈值，满足：；
S2：统计图像灰度直方图，该灰度直方图一般符合正态分布，检测出该灰度直方图种主波峰的左谷点灰度值a，并在该灰度值与所述第一阈值T1的差值区域（a-T1）中，并以该差值区域作为新的背景图像，采用大津法或迭代法进行阈值分割，获得第二阈值T2，并获取由该第一阈值T2确定的第二等高线和第二区域；
S3：用一阶微分边界扫描或山脊边界扫描算法获取菌落边界或菌落部分边界；对所述平滑图像进行一阶微分获得梯度图像，并将梯度图像转换成二值图像。，并将梯度图像转换成二值图像；在本实施例中，用表示g的梯度向量，，是平滑图像g的梯度图像；
S4：对所述第一等高线向所述第二等高线方向进行加速度判定操作，对所述第一区域向所述第二区域进行扩展操作，以去除掉菌落图像背景产生的误识别曲线及去掉菌落图像中的黑洞和其它噪声的影响，将扩展操作后处于所述第一区域与所述第二区域交集区域中，且与所述第一区域对应的菌落目标种子点进行合并，扩大种子区域；继续对所述第一区域向所述第二区域进行扩展操作，进一步去除掉菌落图像背景产生的误识别第一等高线曲线及去掉黑洞及其它噪声的影响，进一步扩大所述种子区域；
S5：对所述菌落边界或所述菌落的部分边界进行细线化操作和端点检测操作，并根据该菌落边界中端点的方向及端点之间的距离对应连接端点，尽可能使菌落边界闭合，消除该菌落边界中存在的断点和空隙；
S6：以所述步骤S5所得的菌落边界作为约束条件，对所述种子区域进行膨胀操作，再通过GVF Snake算法对菌落轮廓进行精确地定位；
在本实施例中，所采用的Snake方法较为稳定，但也有定位误差等问题，为了克服这些问题，采用梯度矢量流(Gradient Vector Flow, GVF) 代替传统外力场，它是通过扩散边缘图的负梯度矢量得到，在扩大边缘势能作用范围的同时，保持了边界区域梯度矢量流的性质。令为图像的边缘图，则有，设的梯度场为，把梯度场向图像的边缘扩散，则形成了扩散的梯度矢量流场，其对能量函数进行极小化，即，而求解能量最小的矢量场可以通过欧拉方程来求解：，其中为拉普拉斯算子；、为边缘图对的偏导；
S7：由于受到噪声的影响及其黑洞，菌落中会有一些孔洞，为此，在上述图像分割的基础上，再对步骤S6中完成菌落边界精确定位后的菌落图像中的孔洞进行充填，如果需要，还可以用二次曲线对菌落边界进行拟合或平滑操作，得到最终图像分割结果。
为了让本领域的技术人员进一步了解本发明中所提出的一种基于水平集及GVF Snake轮廓精确地定位的的菌落图像分割方法，下面结合具体的实施例进行说明。
如图2所示，根据该菌落中的白色光点和黑色光点及菌落边界弱的特点，通过以如图1所示的方法获得所有菌落的闭合轮廓。本发明所提出的菌落图像分割方法是基于基于图像灰度等高线、区域膨胀、一阶微分边界扫描及GVF Snake算法对菌落轮廓进行精确定位的图像分割算法，是建立在大津阈值和三条等高线及一阶微分边界扫描及GVF Snake算法的基础之上，是一种新的图像分割方法，该方法首先检测出白色光区域的边界，再按一定的阈值增长检测出第二阈值，按等高线的顺序逐步膨胀区域面积到菌落边界，使等高线内部区域不断向外扩展，解决区域由多余黑、白点及其它噪声引起的过分分割问题，然后用一阶微分边界算子或山脊扫描算子检测出图像中菌落的边界和部分边界，再进行菌落边界闭合的后处理。再以该菌落边界作为约束条件，膨胀上述的种子区域，再用GVF Snake算法对菌落轮廓进行精确地定位，这里主要采用了梯度矢量流(Gradient Vector Flow, GVF) 代替传统外力场，它是通过扩散边缘图的负梯度矢量得到，在扩大边缘势能作用范围的同时，保持了边界区域梯度矢量流的性质。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (2)

1.一种基于水平集及GVF Snake精确定位的菌落图像分割方法，其特征在于，按照如下步骤实现：
S1：采集并输入菌落图像，通过高斯算子对该菌落图像进行平滑滤波操作，获取平滑图像；
采用大津法对该平滑图像进行阈值分割，获得第一阈值T1，并获取由该第一阈值T1确定的第一等高线和第一区域；
S2：统计所述平滑图像的灰度直方图，检测出该灰度直方图种主波峰的左谷点灰度值，并在该灰度值与所述第一阈值T1的差值区域中，采用大津法进行阈值分割，获得第二阈值T2，并获取由该第二阈值T2确定的第二等高线和第二区域；
S3：用一阶微分边界扫描或山脊边界扫描算法获取菌落边界；
S4：对所述第一等高线向所述第二等高线方向进行加速度判定操作，对所述第一区域向所述第二区域进行扩展操作；将扩展操作后处于所述第一区域与所述第二区域交集区域中，且与所述第一区域对应的菌落目标种子点进行合并，扩大种子区域；继续对所述第一区域向所述第二区域进行扩展操作，进一步扩大所述种子区域；
S5：对所述菌落边界进行细线化操作和端点检测操作，并根据该菌落边界中端点的方向及端点之间的距离对应连接端点，以使菌落边界闭合，消除该菌落边界中存在的断点和空隙；
S6：以所述步骤S5所得的菌落边界作为约束条件，对所述种子区域进行膨胀操作，再通过GVF Snake算法对该菌落边界进行精确地定位；
S7：对所述步骤S6中完成菌落边界精确定位后的菌落图像中的孔洞进行充填，并用二次曲线对菌落边界进行拟合操作，得到最终图像分割结果。
2.根据权利要求1所述的基于水平集及GVF Snake精确定位的菌落图像分割方法，其特征在于，所述步骤S3还包括：对所述平滑图像进行一阶微分获得梯度图像，并将梯度图像转换成二值图像。
CN201410813281.2A 2014-12-24 2014-12-24 基于水平集及GVF Snake精确定位的菌落图像分割方法 Expired - Fee Related CN104463892B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410813281.2A CN104463892B (zh) 2014-12-24 2014-12-24 基于水平集及GVF Snake精确定位的菌落图像分割方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410813281.2A CN104463892B (zh) 2014-12-24 2014-12-24 基于水平集及GVF Snake精确定位的菌落图像分割方法

Publications (2)

Publication Number Publication Date
CN104463892A true CN104463892A (zh) 2015-03-25
CN104463892B CN104463892B (zh) 2017-06-06

Family

ID=52909874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410813281.2A Expired - Fee Related CN104463892B (zh) 2014-12-24 2014-12-24 基于水平集及GVF Snake精确定位的菌落图像分割方法

Country Status (1)

Country Link
CN (1) CN104463892B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446917A (zh) * 2016-10-18 2017-02-22 电子科技大学中山学院 一种量体特征点提取方法及装置
CN106934785A (zh) * 2015-12-28 2017-07-07 哈尔滨工业大学 一种用于机器人虚拟训练系统中肝脏模型的医学图像分割方法
CN107203766A (zh) * 2017-04-19 2017-09-26 杭州泽火科技有限公司 精确定位图像中字符的方法、装置及系统
CN109064475A (zh) * 2018-09-11 2018-12-21 深圳辉煌耀强科技有限公司 针对宫颈脱落细胞图像的图像分割方法和装置
CN109214393A (zh) * 2018-07-25 2019-01-15 南华大学 一种多孔培养板板孔菌群生长情况图像检测方法及系统
CN110874733A (zh) * 2018-08-31 2020-03-10 北京意锐新创科技有限公司 支持外接设备的被动扫码支付方法和装置
CN110874732A (zh) * 2018-08-31 2020-03-10 北京意锐新创科技有限公司 基于移动支付设备的被扫支付方法和装置
CN111091668A (zh) * 2018-10-24 2020-05-01 北京意锐新创科技有限公司 基于二维码识读装置的零距离收银方法和装置
CN111091366A (zh) * 2018-10-24 2020-05-01 北京意锐新创科技有限公司 应用于远距离的扣扫式支付方法和装置
CN114317675A (zh) * 2022-01-06 2022-04-12 福州大学 基于机器学习定性和定量不同创面细菌的检测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101571419A (zh) * 2009-06-15 2009-11-04 浙江大学 采用图像分割的汽车仪表led指示灯自动检验方法
CN101901477A (zh) * 2010-07-27 2010-12-01 中国农业大学 植物叶片田间图像边缘提取方法及其系统
CN102068281A (zh) * 2011-01-20 2011-05-25 深圳大学 一种占位性病变超声图像的处理方法
CN102676633A (zh) * 2012-03-08 2012-09-19 天津大学 一种菌落自动计数方法
KR101212788B1 (ko) * 2012-07-12 2012-12-18 태양정보시스템(주) 통합 벡터 모델 기반 지도 제작을 위한 수치지도 시스템 및 이를 이용한 수치지도 작성 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101571419A (zh) * 2009-06-15 2009-11-04 浙江大学 采用图像分割的汽车仪表led指示灯自动检验方法
CN101901477A (zh) * 2010-07-27 2010-12-01 中国农业大学 植物叶片田间图像边缘提取方法及其系统
CN102068281A (zh) * 2011-01-20 2011-05-25 深圳大学 一种占位性病变超声图像的处理方法
CN102676633A (zh) * 2012-03-08 2012-09-19 天津大学 一种菌落自动计数方法
KR101212788B1 (ko) * 2012-07-12 2012-12-18 태양정보시스템(주) 통합 벡터 모델 기반 지도 제작을 위한 수치지도 시스템 및 이를 이용한 수치지도 작성 방법

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
EMRE CELEBI M ET AL.: "Unsupervised border detection in dermoscopy images", 《SKIN RESEARCH AND TECHNOLOGY》 *
丁茂实: "基于运动背景的运动目标检测技术研究", 《中国优秀硕士学位论文全文数据库_信息科技辑》 *
朱国普: "基于活动轮廓模型的图像分割", 《中国博士学位论文全文数据库_信息科技辑》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934785B (zh) * 2015-12-28 2020-06-09 哈尔滨工业大学 一种用于机器人虚拟训练系统中肝脏模型的医学图像分割方法
CN106934785A (zh) * 2015-12-28 2017-07-07 哈尔滨工业大学 一种用于机器人虚拟训练系统中肝脏模型的医学图像分割方法
CN106446917A (zh) * 2016-10-18 2017-02-22 电子科技大学中山学院 一种量体特征点提取方法及装置
CN106446917B (zh) * 2016-10-18 2019-08-16 电子科技大学中山学院 一种量体特征点提取方法及装置
CN107203766A (zh) * 2017-04-19 2017-09-26 杭州泽火科技有限公司 精确定位图像中字符的方法、装置及系统
CN107203766B (zh) * 2017-04-19 2019-08-20 杭州泽火科技有限公司 精确定位图像中字符的方法、装置及系统
CN109214393A (zh) * 2018-07-25 2019-01-15 南华大学 一种多孔培养板板孔菌群生长情况图像检测方法及系统
CN109214393B (zh) * 2018-07-25 2020-12-22 南华大学 一种多孔培养板板孔菌群生长情况图像检测方法及系统
CN110874733A (zh) * 2018-08-31 2020-03-10 北京意锐新创科技有限公司 支持外接设备的被动扫码支付方法和装置
CN110874732A (zh) * 2018-08-31 2020-03-10 北京意锐新创科技有限公司 基于移动支付设备的被扫支付方法和装置
CN109064475A (zh) * 2018-09-11 2018-12-21 深圳辉煌耀强科技有限公司 针对宫颈脱落细胞图像的图像分割方法和装置
CN111091366A (zh) * 2018-10-24 2020-05-01 北京意锐新创科技有限公司 应用于远距离的扣扫式支付方法和装置
CN111091668A (zh) * 2018-10-24 2020-05-01 北京意锐新创科技有限公司 基于二维码识读装置的零距离收银方法和装置
CN114317675A (zh) * 2022-01-06 2022-04-12 福州大学 基于机器学习定性和定量不同创面细菌的检测方法及系统

Also Published As

Publication number Publication date
CN104463892B (zh) 2017-06-06

Similar Documents

Publication Publication Date Title
CN104463892A (zh) 基于水平集及GVF Snake精确定位的菌落图像分割方法
Liao et al. Automatic segmentation for cell images based on bottleneck detection and ellipse fitting
CN115249246B (zh) 一种光学玻璃表面缺陷检测方法
CN110838126B (zh) 细胞图像分割方法、装置、计算机设备和存储介质
Wang et al. Modified watershed technique and post-processing for segmentation of skin lesions in dermoscopy images
CN107256558A (zh) 一种无监督式的宫颈细胞图像自动分割方法及系统
CN110930427B (zh) 一种基于语义轮廓信息的图像分割方法、设备和存储介质
CN102982545B (zh) 一种图像深度估计方法
Liu et al. Active contour model driven by local histogram fitting energy
CN103035013A (zh) 一种基于多特征融合的精确运动阴影检测方法
CN105205832B (zh) 一种运动目标检测的方法
CN107220596A (zh) 基于级联纠错机制的人体姿态估计方法
CN106127205A (zh) 一种适用于室内轨道机器人的数字仪表图像的识别方法
CN106991686A (zh) 一种基于超像素光流场的水平集轮廓跟踪方法
CN105513053A (zh) 一种用于视频分析中背景建模方法
CN103679677A (zh) 一种基于模型互更新的双模图像决策级融合跟踪方法
Zhao et al. Separate degree based Otsu and signed similarity driven level set for segmenting and counting anthrax spores
CN102663723A (zh) 一种基于颜色样本与电场模型的图像分割方法
CN110910367A (zh) 生物反应器细胞培养质量评价方法
CN106447673A (zh) 一种非均匀光照条件下的芯片引脚提取方法
CN106327490A (zh) 一种基于白细胞检测的细胞核分割方法
CN107123102A (zh) 一种贴壁细胞生长融合度自动分析方法
CN105488475A (zh) 手机中人脸检测方法
CN107610148B (zh) 一种基于双目立体视觉系统的前景分割方法
CN106651860A (zh) 一种图像分割方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170606

Termination date: 20191224