一种显微图像的复杂粘连细胞个数统计方法技术

本发明专利技术涉及一种显微图像的复杂粘连细胞个数统计方法，该方法首先将原始彩色显微图像进行预处理，然后再进行二值化；再用链码差方法找出所有细胞轮廓的凹点，连接相互配对的凹点后，将略微粘连的细胞进行分割；最后用采用区域面积统计法计算细胞个数。与现有技术相比，本发明专利技术的优点在于：由于本方法首先利用彩色显微图像的颜色特点对原始彩色显微图像进行预处理，然后再二值化，在利用略微粘连细胞存在明显凹点的特点用链码差检测出细胞轮廓并进行分割，最后用面积统计法计算粘连严重的细胞和单细胞的个数，该方法误差小，精确度高。

【技术实现步骤摘要】

本专利技术涉及显微镜细胞统计领域，特别是涉及一种显微图像的复杂粘连细胞个数统计方法。
技术介绍
显微细胞图像分析包括对显微细胞图像的识别测定、个数统计以及面积、周长、形 状因子等各种参数的测量。在实际的显微细胞图像中，由于切片和细胞本身的原因，经常会 出现多个细胞重叠在一起形成一个较大的区域，如血细胞、淋巴细胞、癌细胞等。对团聚在 一起的重叠粘连细胞计数是显微细胞图像分析的难点所在。目前，已提出许多分割算法处 理细胞重叠粘连，如基于数学形态学的腐蚀膨胀法、二值图像多边形近似找凹点切割、测地 重建法和流域(亦称分水岭或水坝)变换等。由于简单且可以并行计算，目前基于数学形 态学的方法是主流，但是它们在应用中表现出局限性。如腐蚀膨胀法对重叠区域采用先腐 蚀分离再膨胀恢复的策略，获得具有几何意义的分界线。但是该法比较简单，腐蚀和膨胀的 次数须人工经验性确定，且腐蚀和膨胀的不可逆性会造成一定误差。凹点分割对于少量细 胞且粘连细胞的面积少于三分之一效果比较好，但对于粘连一半或者过半的细胞，效果不 是很明显；流域分割种子点比较难找而且对于粘连厉害的种子点错误的概率很大。
技术实现思路
本专利技术所要解决的技术问题是针对上述现有技术提供一种误差小、精确度较高的 显微图像的复杂粘连细胞个数统计方法。 本专利技术解决上述技术问题所采用的技术方案为该显微图像的复杂粘连细胞个数 统计方法，其特征在于包括以下步骤 步骤一、利用色度分量直方图分别获得原始彩色显微图像中细胞、背景、和除了细 胞、背景之外的杂质噪声的颜色取值区间； 步骤二、将除了细胞、背景之外的杂质噪声的颜色取值区间内的像素点的颜色值 赋值为背景的颜色取值区间内的一个值； 步骤三、将经过步骤二处理后的色度分量直方图进行canny边缘检测，将两个细 胞的边缘线进行增粗，并将边缘点的颜色值赋值为背景的颜色取值区间内的一个值； 步骤四、将经过步骤三处理后的色度分量直方图进行二值化，得到二值化显微图 像； 步骤五、在二值化后的显微图像上，用链码差方法找出所有细胞轮廓的凹点，连接 相互配对的凹点后，将略微粘连的细胞进行分割； 步骤六、采用区域面积统计法统计经过步骤五处理后的二值化显微图像中所有细 胞轮廓的面积，将所有细胞轮廓中面积最小值设定为单位细胞面积a，设某一细胞轮廓的面积为13，该细胞轮廓中细胞个数为"=^"，其中k为粘连系数，k的取值范围为0. 9 < k < 1 ;将所有细胞轮廓中细胞个数相加后的和即为该显微图像中细胞总个数。 所述步骤五中，首先采用链码差方法寻找二值化显微图像上所有细胞轮廓的凹 点，然后对所有凹点的凹陷角的平分线进行比较，若两个凹点的凹陷角的平分线在一条直 线上，则这两个凹点为相互配对的凹点，将这两个凹点直接相连，就可以将略微粘连的细胞 进行分割。 所述步骤六中，采用基于Freeman 8_链码的区域面积统计法统计经过步骤五处 理后的二值化显微图像中所有细胞轮廓的面积。 当显微图像显示的是血细胞时，所述步骤六中k取值较好为0. 94。 与现有技术相比，本专利技术的优点在于由于本方法首先利用彩色显微图像的颜色特点对原始彩色显微图像进行预处理，然后再二值化，在利用略微粘连细胞存在明显凹点的特点用链码差检测出细胞轮廓并进行分割，最后用面积统计法计算粘连严重的细胞和单细胞的个数，该方法误差小，精确度高。附图说明 图1为本专利技术实施例方法的流程图。 具体实施例方式以下结合附图实施例对本专利技术作进一步详细描述。 本专利技术提供了，其包括以下步骤如 图1所示， 步骤一、利用色度分量直方图分别获得原始彩色显微图像中细胞、背景、和除了细 胞、背景之外的杂质噪声的颜色取值区间； 通过对彩色显微图像特别是血细胞图像的分析可以看出，色度信息对于细胞的判 别分析有着很重要的作用，如果仅仅将彩色图像转化为灰度图，无疑将会损失一些对分类 判决有效的信息，所以本步骤直接针对彩色图像进行处理；色度特征是在图像识别中应用 最为广泛的视觉特征，颜色特征对图像本身的尺寸、方向、视角的依赖性较小；一幅彩色图 像中色度分量直方图最能体现显微图像中细胞与背景之间的颜色区别，因此本步骤采用色 度分量图进行相关处理； 步骤二、将除了细胞、背景之外的杂质噪声的颜色取值区间内的像素点的颜色值 赋值为背景的颜色取值区间内的一个值； 步骤三、将经过步骤二处理后的色度分量直方图进行canny边缘检测，将两个细 胞的边缘线进行增粗，并将边缘点的颜色值赋值为背景的颜色取值区间内的一个值； 这是因为显微图像中大部分细胞的粘连处都会有一些轻微的颜色上的突变，而直 接二值化对这些细微变化不敏感，因此，我们将这些颜色突变强化，使其在二值化过程中能 被区分； 步骤四、将经过步骤三处理后的色度分量直方图进行二值化，得到二值化显微图 像； 二值化是一个常规处理方法；例如同样以血细胞为例，将细胞的颜色取值区间内的像素点的颜色值赋值为o，将背景和除了细胞、背景之外的杂质噪声的颜色取值区间的像素点的颜色值赋值为l，就实现了二值化； 步骤五、在二值化后的显微图像上，用链码差方法找出所有细胞轮廓的凹点，连接 相互配对的凹点后，将略微粘连的细胞进行分割； 本步骤中，首先采用链码差方法寻找二值化显微图像上所有细胞轮廓的凹点，然 后对所有凹点的凹陷角的平分线进行比较，若两个凹点的凹陷角的平分线在一条直线上， 则这两个凹点为相互配对的凹点，将这两个凹点直接相连，就可以将略微粘连的细胞进行 分割。链码差方法可以采用《科学技术与工程》杂志(2007. 7印刷的第五期中69-72页)中 公布的《基于链码差的边界凸凹性判别》中采用的方法。 步骤六、采用基于Freeman 8_链码的区域面积统计法统计经过步骤五处理后的 二值化显微图像中所有细胞轮廓的面积，将所有细胞轮廓中面积最小值设定为单位细胞面积a，设某一细胞轮廓的面积为b，该细胞轮廓中细胞个数为n = ^ ，其中k为粘连系数，k的取值范围为0. 9 < k < 1 ;将所有细胞轮廓中细胞个数相加后的和即为该显微图像中细胞总个数。同样以血细胞为例，k取值为0.94，计算的误差小于3X。权利要求，其特征在于包括以下步骤步骤一、利用色度分量直方图分别获得原始彩色显微图像中细胞、背景和除了细胞、背景之外的杂质噪声的颜色取值区间；步骤二、将除了细胞、背景之外的杂质噪声的颜色取值区间内的像素点的颜色值赋值为背景的颜色取值区间内的一个值；步骤三、将经过步骤二处理后的色度分量直方图进行canny边缘检测，将两个细胞的边缘线进行增粗，并将边缘点的颜色值赋值为背景的颜色取值区间内的一个值；步骤四、将经过步骤三处理后的色度分量直方图进行二值化，得到二值化显微图像；步骤五、在二值化后的显微图像上，用链码差方法找出所有细胞轮廓的凹点，连接相互配对的凹点后，将略微粘连的细胞进行分割；步骤六、采用区域面积统计法统计经过步骤五处理后的二值化显微图像中所有细胞轮廓的面积，将所有细胞轮廓中面积最小值设定为单位细胞面积a，设某一细胞轮廓的面积为b，该细胞轮廓中细胞个数为 <mrow><mi>n</mi><mo>=</mo><mfrac> &本文档来自技高网...

【技术保护点】
一种显微图像的复杂粘连细胞个数统计方法，其特征在于：包括以下步骤步骤一、利用色度分量直方图分别获得原始彩色显微图像中细胞、背景和除了细胞、背景之外的杂质噪声的颜色取值区间；步骤二、将除了细胞、背景之外的杂质噪声的颜色取值区间内的像素点的颜色值赋值为背景的颜色取值区间内的一个值；步骤三、将经过步骤二处理后的色度分量直方图进行ｃａｎｎｙ边缘检测，将两个细胞的边缘线进行增粗，并将边缘点的颜色值赋值为背景的颜色取值区间内的一个值；步骤四、将经过步骤三处理后的色度分量直方图进行二值化，得到二值化显微图像；步骤五、在二值化后的显微图像上，用链码差方法找出所有细胞轮廓的凹点，连接相互配对的凹点后，将略微粘连的细胞进行分割；步骤六、采用区域面积统计法统计经过步骤五处理后的二值化显微图像中所有细胞轮廓的面积，将所有细胞轮廓中面积最小值设定为单位细胞面积ａ，设某一细胞轮廓的面积为ｂ，该细胞轮廓中细胞个数为ｎ＝ｂ／ｋａ，其中ｋ为粘连系数，ｋ的取值范围为０．９＜ｋ＜１；将所有细胞轮廓中细胞个数相加后的和即为该显微图像中细胞总个数。

【技术特征摘要】

【专利技术属性】
技术研发人员：杨鸣，杨翠翠，蒋金涛，邵华，章湖，王建宏，
申请(专利权)人：宁波大学，
类型：发明
国别省市：97[中国|宁波]