一种细胞形态扫描区域选择方法与流程

本发明涉及一种图像处理方法。特别是涉及一种细胞形态扫描区域选择方法。
背景技术：
：细胞形态学作为各种血细胞检验的基础，对大多数血液病的诊断与治疗起着至关重要的作用。很多疾病通过在显微镜下观察各种血细胞的形态与数量可以得到初步印象或诊断及分类，并为进一步检查提供依据或思路。传统的检验方法需要人工选择涂片中具有代表性的细胞镜检区域，然后从低倍镜转换到高倍镜进一步观察细胞形态、比例等。随着全自动显微镜扫描平台的出现，通过图像处理及人工智能技术可实现对涂片的全自动扫描及分析，彻底改变了传统的人工镜下对血细胞的识别及手工计数的工作方式，极大地提高了工作效率，并且改善了检验结果的精密度和准确度。而在全自动扫描流程中涉及到低倍物镜(通常是10x)和高倍物镜(通常是100x)的转换，此过程需要在低倍镜下寻找合适的扫描区域，然后再切换到高倍镜下进行扫描。高分辨率细胞形态学自动扫描流程中涉及到低倍物镜和高倍物镜的转换过程，首先需要从大量的低倍物镜扫描结果图中寻找合适的扫描区域(一个低倍视野＝100个高倍视野)，然后切换到高倍镜下进行详细扫描。因此，扫描区域的选择是直接影响扫描效果的关键因素。技术实现要素：本发明所要解决的技术问题是，提供一种能够从大量的低倍物镜扫描图中选择出合适的区域以进行高倍镜扫描的细胞形态扫描区域选择方法。本发明所采用的技术方案是：一种细胞形态扫描区域选择方法，在高分辨率细胞形态学自动扫描流程中，通过10倍物镜预扫得到若干张十倍图之后，首先估算每张十倍图中的有核细胞数和成熟红细胞数，计算区域面积大于设定阈值的空白区域轮廓，计算染色深度大于设定阈值且连通区域面积大于设定阈值的轮廓，然后按照成熟红细胞数从大到小的顺序选择区域，同时要避开所述的空白区域和染色深度大于设定阈值的区域，直到所选区域中的有核细胞数达到设定值，完成最优扫描区域的选择。本发明的一种细胞形态扫描区域选择方法，根据低倍图(清晰度不如高倍图)选择合适有效的高倍镜扫描区域的算法，其扫描区域选择的结果符合检验的要求，与人工选择基本一致。本发明的方法能够模拟人工在低倍镜下选择代表性区域的过程，能够避免高倍镜下对大量无效区域的扫描，直接对选择的有效区域进行高倍镜扫描，可以减少扫描时间，提高扫描效率，从而确保自动扫描流程高效顺利地进行。附图说明图1是本发明中有核细胞数和成熟红细胞数计算流程图；图2是本发明中较大空白区域轮廓计算流程图；图3是本发明中染色过深的大连通区域轮廓计算流程图。具体实施方式下面结合实施例和附图对本发明的基于丝杠和直线导轨的竖直方向波浪浮标检测装置及方法做出详细说明。本发明的一种细胞形态扫描区域选择方法，是在高分辨率细胞形态学自动扫描流程中，通过10倍物镜预扫得到若干张十倍图之后，包括：根据图1所示的流程估算每张十倍图中的有核细胞数、成熟红细胞数；根据图2所示的流程计算区域面积大于设定阈值的空白区域轮廓；根据图3所示的流程计算染色深度大于设定阈值且连通区域面积大于设定阈值的轮廓，然后按照成熟红细胞数从大到小的顺序选择区域，同时要避开所述的空白区域和染色深度大于设定阈值的区域，直到所选区域中的有核细胞数达到设定值，完成最优扫描区域的选择。如图1、图2、图3所示，具体包括如下步骤：1)对预扫得到的每一张十倍图的各参数进行计算，包括：1.1)计算有核细胞数有核细胞是需在高倍镜下进一步观察的目标，其总数需达到检验的要求。有核细胞数的计算方法具体如下：(1.1.1)对十倍图进行通道分离，得到r、g、b三个通道的灰度图；(1.1.2)根据有核细胞的特点，将b通道与g通道相减得到突出有核细胞的灰度图；(1.1.3)通过三类别的最大类间方差法及一维最大熵方法对得到的灰度图计算阈值，取其中最大阈值对灰度图进行二值化；(1.1.4)对二值化结果图进行设定尺寸的腐蚀操作，并去掉小于设定面积的连通区域；(1.1.5)统计最终图中的连通区域个数即为估算的有核细胞数；1.2)计算成熟红细胞数在成熟红细胞基本为单层、分布均匀的视野中有核细胞结构较清晰，适于细胞形态的观察，成熟红细胞数在一定程度上能够反应其分布情况，单个十倍图视野中识别出的红细胞数越多说明其分布越相对均匀。成熟红细胞数的计算方法具体如下：(1.2.1)将十倍图的r通道反色；(1.2.2)对第(1.1)步得到的二值化结果图进行设定尺寸的开运算，结果作为掩膜(mask)；(1.2.3)将r通道反色图中对应掩膜的非零位置赋值为零，以排除有核细胞对成熟红细胞检测的干扰；(1.2.4)进行霍夫(hough)圆检测，检测出圆的个数即识别出的成熟红细胞个数；1.3)计算区域面积大于设定阈值的空白区域轮廓(油滴或划痕)较大的空白区域在高倍镜扫描时是应尽量避开的，其计算方法如下：(1.3.1)将十倍图转为灰度图；(1.3.2)将灰度图的宽高都缩小为原图的1/10；(1.3.3)利用大津阈值方法对缩小后的灰度图做二值化；(1.3.4)对二值化结果图进行设定尺寸的腐蚀，并去掉小于设定面积的连通区域(1.3.5)对去掉小于设定面积的连通区域的结果图进行相同尺寸的膨胀，并恢复原图尺寸；(1.3.6)寻找恢复原图尺寸的结果图中的连通区域的外轮廓，即为大于设定阈值的空白区域轮廓；1.4)计算染色深度大于设定阈值且连通区域面积大于设定阈值的轮廓染色过深的大连通区域(或较大杂质)在高倍镜扫描时也是应尽量避开的，其计算方法如下：(1.4.1)将十倍图转为灰度图并反色；(1.4.2)将灰度图的宽高都缩小为原图的1/10；(1.4.3)通过三类别的最大类间方差方法计算阈值，取其中最大值对缩小后的灰度图做二值化；(1.4.4)对二值化结果图进行设定尺寸的腐蚀，并去掉小于设定面积的连通区域；(1.4.5)对上去掉小于设定面积的连通区域的结果图进行相同尺寸的膨胀，并恢复原图尺寸；(1.4.6)寻找恢复原图尺寸结果图中的连通区域的外轮廓，即为染色深度大于设定阈值且连通区域面积大于设定阈值的轮廓；2)寻找最优扫描区域(1)将所有十倍图按坐标位置组成m列、n行的矩阵，如表1所示；表10，01，02，03，00，11，12，13，10，21，22，23，2(2)将计算得到的成熟红细胞数、区域面积大于设定阈值的空白区域轮廓的面积和染色深度大于设定阈值且连通区域面积大于设定阈值的轮廓面积对应取值为b、c、d，计算每张十倍图评估值q的函数为：q＝k*b/(m*c+n*d)其中，k、m、n为大于零小于1的常量；q值取值越高，图像质量越好；(3)选出矩阵中，q值最大的图作为当前选择区域，如表1中的[2,1]；(4)判断当前选择区域有核细胞总数t是否满足扫描需求的有核细胞个数c，其中t>＝c；如果满足需求或者已经覆盖矩阵中的所有十倍图区域，则当前选择区域就是最优扫描区域，结束；否则进入下一步；(5)选出与当前选择区域相邻的上下左右区域中q值最大的区域，将所述的当前选择区域合并到q值最大的区域中，形成新的当前选择区域(如表1中的[2,1]合并到[3,1])，返回第(4)步。当前第1页12