一种细胞显微图像感兴趣区域自动查找方法技术

本发明专利技术涉及一种细胞显微图像感兴趣区域自动查找方法。现有的方法处理速度慢、细胞边缘提取不准确、重叠细胞分割不精确。本发明专利技术方法首先利用兰索斯算法把原始图像Ａ压缩成图像Ｂ；其次对图像Ｂ进行灰度化处理，然后再采用多阈值融合算法提取各个细胞的边缘，生成二值图Ｃ；然后对二值图Ｃ进行形态学膨胀操作，计算图像边缘形状改变的速度，直到离散程度小于设定值；最后提取细胞图像的边缘并还原到原始图像Ａ中，即得到感兴趣区域。本发明专利技术方法处理速度快，提取细胞边缘的准确度更高。

【技术实现步骤摘要】

本专利技术属于显微图像处理
，具体涉及。
技术介绍
细胞显微图像计算机自动分析技术是一种用计算机代替人眼进行的显微图像快 速全自动处理技术，在加快图像的处理速度和节省人力物力资源等方面有着重要的意义。 细胞显微图像计算机自动分析技术虽然发展迅速，但是目前国内外的细胞显微图像计算机 自动分析技术仍然存在着处理速度慢、细胞边缘提取不准确、重叠细胞分割不精确及阴性 细胞提取不准确等诸多问题。一个微小的细胞切片在显微镜下被放大几百倍之后，经过图像采集和拼接得到的 图片是相当大的(甚至可以达到GByte)。如果不压缩直接处理，将耗费大量的时间和空间。研究发现组织的异常细胞显微图像的表面积、周长、等效圆直径等形态参数与正 常细胞相比有显著差异或高度显著差异。细胞显微图像中异常细胞往往呈多个聚集区分 布，聚集区内细胞相互距离近，而且重叠细胞比较多，这些异常细胞聚集区就是要提取的感 兴趣区域。如能找出这些感兴趣关键区域将会大大提高图像的处理和分析效率。假如对整 张图做全局的精细处理将耗费大量的时间，而且由于背景明暗的差异将会对分析的结果造 成一定的不良影响，而只对局部感兴趣区域进行处理会提高细胞边界提取的精度。因此对 异常细胞聚集的局部区域进行提取显得相当关键也尤为重要。
技术实现思路
本专利技术针对现有技术不足，提出。本专利技术方法的具体步骤为步骤(1)利用兰索斯(Lanczos)算法把原始图像A压缩成图像B，压缩倍数为采集 原始图像时放大倍数的1/10 ；步骤(2)对图像B进行灰度化处理，然后再采用多阈值融合算法提取各个细胞的 边缘，生成二值图C ；其中灰度化处理和多阈值融合算法均为现有技术。步骤(3)确定细胞图像的边缘形状；用细胞边缘点分布的离散程度S来表示细胞的边缘形状<formula>formula see original document page 3</formula>,否则 其中<formula>formula see original document page 3</formula>是边缘点横坐标的标准方差，<formula>formula see original document page 3</formula>是边缘点纵坐标的标准方差，Xk是细胞边缘点的横坐标，yk是细胞边缘点的纵坐标，η为细胞边缘点数。步骤(4)对二值图C进行形态学膨胀操作，计算图像边缘形状改变的速度<formula>formula see original document page 4</formula>其中Si为第i次膨胀后的S值。步骤(5)重复步骤⑷，一直到Δ Si小于设定值；步骤(6)提取细胞图像的边缘，并且滤去细胞边缘长度小于指定值的细胞；步骤(7)在原图像A中确定在步骤(6)中提取的细胞图像边缘，即得到感兴趣区 域。本专利技术相对于现有技术具有以下有益效果本方法采用形态学膨胀的方法，显然比传统方法的处理速度要快的多，而且本方 法更加适合于细胞显微图像中异常细胞往往呈多个聚集区分布，聚集区内细胞相互距离 近，而且重叠细胞比较多的特点。用本专利技术方法提取感兴趣区域后，再用阈值分割等算法提取细胞的边缘再提取细 胞的边缘，可以省去一些没有用的区域的处理，比直接提取细胞的边缘速度要快。找出感兴 趣区域后，再用阈值分割等算法提取细胞的边缘，这样可以提高阈值的准确度，所以比整张 显微图像直接提取边缘，提取细胞边缘的准确度更高。附图说明图1为膨胀操作后边缘扩展与内缩示意；图2为经过Lanczos算法压缩并灰度化后生成的图B ；图3为用多阈值融合算法提取各个细胞边缘后生成的二值图C ；图4为细胞边缘提取结果图；图5为自动膨胀结果图；图6为提取膨胀后的图像边缘图；图7为去除内部边缘点后的最终区域轮廓线图；图8为提取的细胞图像边缘图；图9为滤去细胞边缘长度小于指定值的细胞后的结果图。 具体实施例方式以下结合附图对本专利技术作进一步说明。步骤(1)利用Lanczos算法把原始图像A压缩成图像B，如图2所示；Lanczos算 法使用基于sine函数的滤镜。这是理论上最准确的滤镜，对于不具有锐度转换的照片图 像能产生最佳输出。压缩后图像质量会有所损失，只要缩小后的异常细胞不消失就不会影 响后面的工作，所以可以根据异常细胞的平均面积及采集图象时的放大倍数来确定压缩倍 数，通常压缩倍数为采集原始图像时放大倍数的1/10。步骤(2)对图像B进行灰度化处理，然后再采用多阈值融合算法提取各个细胞的 边界，最后形成二值图，如图3所示；步骤(3)确定细胞图像的边缘形状；提取二值图中的细胞边缘，如图4所示。然后用细胞边缘点分布的离散程度S来 表示细胞的边缘形状。步骤(4)对图像C进行膨胀，边缘的形状特征S就会发生改变，把这种改变定义 为Δ Si = I Si-Si^1(1)把它称作图像形状特征改变的速度。当图像膨胀完了之后，图像的边缘有的会向外扩展，有的则向内缩小，也行的会相 交重叠起来，当然重叠之后边缘也就不成其为边缘了，称这些边缘消亡了。随着图像的膨 胀，边缘也在逐步消亡，这种消亡的速度时快时慢，刚开始是比较慢的，随后逐渐加快，到了 一定程度之后出现一个拐点，速度会突然降下来。这是由于每张细胞显微图像都会有许 多细胞相对聚集的区域，而且同一区域内，细胞之间的间隔总是相差不大的；随着图像的膨 胀，细胞之间的间隔也在缩小，缩小到一定程度之后细胞开始重叠，边缘消亡也就发生了， 随着图像的继续膨胀，近距离内的细胞逐渐粘连在一起，当形成一个大的连通区域后，必然 会出现一个短暂的稳定期，细胞的重叠一下子不再发生或很少发生，表现在边缘的消亡上 就是消亡的速度突然变慢了。把边缘的扩展对应为Ei,内缩定义为Ci,消亡对应为\。则有Si = Ej+Ci+Vi(2)而边缘的扩展、内缩和消亡所引起的形状特征的变化依次设为ΔΕ” ACi和 Δ\。对它们有如下定义AEi = Ei-Ep1(3)ACi = Ci-C^1(4)AVi = Vi-Vp1(5)ASi = I Δ Ei+ Δ Ci+ Δ Vi |(6)图1示意了任意形状的图像(阴影部分)，在发生膨胀时边缘的走向，其中阴影外 围的边缘向外扩展，而阴影内部的边缘则向内收缩了。举扩展为例，根据膨胀的定义可知， 无论对任何形状的图像而言，扩展变化都是对称的，一个方向上的边缘点发生了改变意味 着在与它相反的方向上必有某些边缘点发生了与它几乎相反的改变，因此边缘的重心位置_ η— η并没有发生重大的改变，即I=Σ &/ w和f=Σα “几乎不变，从而边缘点横纵坐标的标准方差也是基本保持不变的。那么就可以做出如下推导Ei ^ Eh(7)AEi = Ei-Ep1 ^ 0(8)内缩变化与扩展变化类似，同理可以得到如下式子Ci ^ Ci^1(9)Δ Ci = Ci-Ci^1 ^ 0(10)另外，扩展变化和内缩变化总是朝着两个相反的方向进行的，假如定扩展变化的方向为正，那么内缩变化的方向则必然为负，即假如Δ Ei > 0则必有Δ Ci < 0，再结合式(8) 和式(10)可得如下式子Δ Ei+Δ Ci I ^ 0(11)再由式(本文档来自技高网...

【技术保护点】
一种细胞显微图像感兴趣区域自动查找方法，其特征在于该方法包括如下步骤：步骤（１）利用兰索斯算法把原始图像Ａ压缩成图像Ｂ，压缩倍数为采集原始图像时放大倍数的１／１０；步骤（２）对图像Ｂ进行灰度化处理，然后采用多阈值融合算法提取各个细胞的边缘，生成二值图Ｃ；步骤（３）确定细胞图像的边缘形状；用细胞边缘点分布的离散程度Ｓ来表示细胞的边缘形状：＊＊＊其中：＊＊＊是边缘点横坐标的标准方差，＊＊＊是边缘点纵坐标的标准方差，ｘ↓［ｋ］是细胞边缘点的横坐标，ｙ↓［ｋ］是细胞边缘点的纵坐标，ｎ为细胞边缘点数；步骤（４）对二值图Ｃ进行形态学膨胀操作，计算图像边缘形状改变的速度ΔＳ↓［ｉ］＝｜Ｓ↓［ｉ］－Ｓ↓［ｉ－１］｜其中Ｓ↓［ｉ］为第ｉ次形态学膨胀后的Ｓ值；步骤（５）重复步骤（４），一直到ΔＳ↓［ｉ］小于设定值；步骤（６）提取细胞图像的边缘，并且滤去细胞边缘长度小于指定值的细胞；步骤（７）在原图像Ａ中确定在步骤（６）中提取的细胞图像边缘，即得到感兴趣区域。

【技术特征摘要】

【专利技术属性】
技术研发人员：李训根，孙玲玲，周磊，叶向春，黄良孟，孔霆，林陆君，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：86[中国|杭州]