基于自适应形态学的细胞边缘分割方法和装置与流程

1.本发明涉及医学图像处理技术领域，尤其涉及一种基于自适应形态学的细胞边缘分割方法和装置。


背景技术：

2.近年来，循环遗传异常细胞(cacs, circulating genetically abnormal cells) 被报道在非小细胞肺癌(nsclc ,non small cell lung cancer)病人的外周血中被发现。针对于外周血中cacs的检测可以更早的预测到肿瘤存在，具有广阔的应用前景。目前，对cacs进行检测的方法，是通过采集外周血液的暗视野(df，dark field)显微镜图像，利用cacs的形态学信息，例如，形状、大小，并基于人工干预，标记df显微镜图像中的cacs，对标记的cacs进行计数，依据计数结果确定cacs在血液中的含量。但依据cacs的形态学信息并结合人工干预进行cacs的检测、确认和计数的方式，检测的主观性较强、可靠性不高，且由于需要人工参与，检测的效率较低、检测成本较高。如果采用计算机辅助分析，首先需要对df图像中的细胞进行精准分割，这是进一步对各个细胞进行自动化分析的基础。
3.对于df图像中细胞的分割，传统的方法是基于形态学或者分水岭等方法来实现，这些方法对于拥挤的细胞效果较差，而df图像中细胞拥挤的情况非常普遍。近年来，人们越来越多地采用基于深度卷积神经网络的实例分割算法来实现对于df图像中细胞的分割。基于深度学习的方法，例如实例分割网络如mask-rcnn，虽然目前能够获取对于拥挤情况下的细胞的检测与分割，但是往往细胞间边界容易产生大量重叠，这种情况对于下游任务，如循环细胞的各路通道进行分析和检测时，需要依据各细胞的边缘对于其中的内容物进行归属的判定，叠加的情形不利于判别内容物的归属。此外，还有部分工作将基于深度学习的分割方式和分水岭分割方式进行了简单结合，以解决细胞重叠的问题。然而，此种方式仍然存在细胞边缘分割不准确的问题，包括多个细胞边缘之间仍存在重叠、细胞边缘不平滑，以及部分细胞信息点被排除在细胞边缘之外从而产生信息点“孤岛”等问题，难以适用于边缘敏感的cacs细胞分割任务中。可见，现有的细胞分割方法难以保证细胞边缘的精确分割，难以适用于cacs细胞的分割任务。


技术实现要素：

4.本发明提供一种基于自适应形态学的细胞边缘分割方法和装置，用以解决现有技术中难以保证细胞边缘的精确分割，难以适用于cacs细胞的分割任务。
5.本发明提供一种基于自适应形态学的细胞边缘分割方法，包括：确定待分割细胞图像；利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理
过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
6.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述卷积参数包括卷积核尺寸；所述基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果，具体包括：确定以初始卷积核尺寸为首的多个尺寸递减的卷积核尺寸；依次基于各个卷积核尺寸，对上一腐蚀处理结果进行迭代腐蚀，得到当前腐蚀处理结果；其中，每次迭代腐蚀会在达到对应卷积核尺寸的迭代停止条件时终止；将当前腐蚀处理结果作为所述最优腐蚀结果；其中，首次腐蚀时，基于所述初始卷积核尺寸对所述初始分割结果进行腐蚀；任一卷积核尺寸的迭代停止条件为当前腐蚀处理结果中的细胞数量与所述初始分割结果中的细胞数量之间的差值小于所述任一卷积核尺寸对应的预设差值阈值，且所述任一卷积核尺寸越小，所述任一卷积核尺寸对应的预设差异阈值越小。
7.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述初始卷积核尺寸是基于如下步骤确定的：基于所述初始分割结果中各像素所属细胞的标签信息，统计所述初始分割结果中各细胞的面积；基于所述初始分割结果中各细胞的面积，计算各细胞的平均面积；基于所述各细胞的平均面积，确定所述初始卷积核尺寸。
8.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果，具体包括：精细分割步骤：利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到当前的精细分割结果；反馈调节步骤：若所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征不匹配，则调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，得到当前的精细分割结果；循环步骤：重复执行所述反馈调节步骤，直至所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征相匹配；结果确定步骤：确定当前的精细分割结果为所述最终分割结果。
9.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，具体包括：对所述待分割细胞图像进行子图划分，得到所述待分割细胞图像的多个子图区域；
获取所述初始分割结果和所述精细分割结果中对应各个子图区域的细胞数量；若所述初始分割结果和所述精细分割结果中对应任一子图区域的细胞数量不匹配，则调整所述初始卷积核尺寸，并对所述初始分割结果中对应所述任一子图区域的部分重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果；基于所述新的最优腐蚀结果，对所述初始分割结果中对应所述任一子图区域的部分重新进行精细分割。
10.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述调整所述初始卷积核尺寸，具体包括：基于多个预设卷积核尺寸，缩小所述初始卷积核尺寸。
11.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，具体包括：对所述最优腐蚀结果进行连通域分析，得到所述最优腐蚀结果中的多个连通域；基于所述最优腐蚀结果中的多个连通域，利用分水岭算法对所述初始分割结果进行精细分割。
12.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述细胞分割模型是基于如下步骤确定的：获取训练样本，并将所述训练样本划分为训练集和测试集；所述训练样本包括样本细胞图像，或，包括样本细胞图像和所述样本细胞图像中各细胞的掩膜标注信息；基于所述训练集对分割模型进行训练，并在训练过程中利用所述测试集对所述分割模型进行阶段性的性能测试，保存性能测试结果最优时的分割模型作为所述细胞分割模型。
13.根据本发明提供的一种基于自适应形态学的细胞边缘分割方法，所述获取训练样本，具体包括：获取细胞显微图像；若所述细胞显微图像的通道数大于所述细胞分割模型的输入条件对应的预设通道数，则将所述细胞显微图像分解为多个低通道图像；对所述多个低通道图像进行图像质量评估，得到所述多个低通道图像的质量评估得分；基于所述多个低通道图像的质量评估得分，获取质量评估得分最高的多个低通道图像；将所述质量评估得分最高的多个低通道图像组合为所述样本细胞图像；其中，所述质量评估得分最高的多个低通道图像组合后的通道数等于所述预设通道数。
14.本发明还提供一种基于自适应形态学的细胞边缘分割装置，包括：图像确定单元，用于确定待分割细胞图像；实例分割单元，用于利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；自适应腐蚀单元，用于基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与
当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；精细分割单元，用于利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
15.本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述基于自适应形态学的细胞边缘分割方法的步骤。
16.本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述基于自适应形态学的细胞边缘分割方法的步骤。
17.本发明还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述基于自适应形态学的细胞边缘分割方法的步骤。
18.本发明提供的基于自适应形态学的细胞边缘分割方法和装置，在利用细胞分割模型对待分割细胞图像进行实例分割，得到初始分割结果后，对初始分割结果进行自适应迭代腐蚀处理，得到最优腐蚀结果，其中，自适应迭代腐蚀处理过程中的卷积参数会随着初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，以较大程度地将粘连重叠的细胞分离开，同时尽可能地保留原有细胞的形态，再在最优腐蚀结果的指引下利用分水岭算法进行精细分割，可以分割出更精确且更符合细胞原有形态的平滑边缘，提高了细胞边缘分割的准确性。
附图说明
19.为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
20.图1是本发明提供的细胞边缘分割方法的流程示意图之一；图2是本发明提供的细胞边缘分割方法与现有方法的效果对比图；图3是本发明提供的自适应迭代腐蚀方法的流程示意图；图4是本发明提供的初始卷积核确定方法的流程示意图；图5是本发明提供的精细分割方法的流程示意图之一；图6是本发明提供的反馈调节步骤的流程示意图；图7是本发明提供的精细分割方法的流程示意图之二；图8是本发明提供的细胞边缘分割方法的详细示意图；图9是本发明提供的细胞边缘分割装置的结构示意图；图10是本发明提供的电子设备的结构示意图。
具体实施方式
21.为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，
而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
22.图1为本发明实施例提供的细胞分割方法的流程示意图，如图1所示，该方法包括：步骤110，确定待分割细胞图像。
23.具体地，可以获取包含细胞的显微镜图像，作为待分割细胞图像。其中，可以选取dapi视野图像作为上述待分割细胞图像。此处，待分割图像包含的细胞可以是循环肿瘤细胞中的cacs细胞，也可以是其他需要进行精确边缘分割的细胞，本发明实施例对此不作具体限定。
24.步骤120，利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果。
25.具体地，预先训练细胞分割模型，该细胞分割模型为实例分割模型，例如mask-rcnn网络等。其中，细胞分割模型可以采用resnet 50 特征金字塔fpn作为特征提取网络，进一步采用roi align模块来实现小目标的感兴趣区域提取，最终通过多分支网络，同步实现对于感兴趣区域内的前景与背景的分割，以及前景区域的分类。由于前景中仅有细胞，因而分类任务是一个二分类任务。在细胞分割模型的训练过程中，将所述样本图像和样本图像中每个细胞的掩膜二值图像等作为元数据输入到细胞分割模型中进行多次epoches训练，从而得到训练好的细胞分割模型。将待分割细胞图像输入到细胞分割模型中进行实例分割，可以得到待分割图像的初始分割结果。其中，初始分割结果中包含待分割图像中每个细胞像素点以及各个细胞像素点所属的细胞标签。
26.步骤130，基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配。
27.具体地，如图2左侧所示，由于待分割细胞图像中细胞之间可能存在重叠，因此细胞分割模型输出的初始分割结果中各个细胞边界之间容易产生大量重叠，细胞边缘分割不准确。因此，需要对重叠的细胞边缘进行精细分割，从而得到精准的细胞边缘。为了解决重叠细胞的问题，部分工作会利用分水岭算法等对实例分割得到的初始分割结果进行进一步的分割操作。然而，针对细胞拥挤、细胞重叠情况较为严重的待分割图像，上述方式的边缘分割效果欠佳。如图2中间所示，利用分水岭算法对初始分割结果进行精细分割后得到的分割结果中，细胞边缘之间仍然存在大量重叠，且细胞边缘不平滑，另外，部分细胞像素点被划分到所有细胞边缘之外，产生了大量信息孤岛，在进行下游任务时，这些处于信息孤岛中的细胞像素点将没有办法判别其归属。
28.对此，本发明实施例在对初始分割结果进行精细分割之前，基于初始分割结果中细胞的统计特征，对初始分割结果对应的二值图像进行自适应迭代腐蚀处理，实现重叠细胞的分离，从而得到待分割细胞图像对应的最优腐蚀结果。此处，自适应迭代腐蚀处理过程中的卷积参数会随着初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小。其中，细胞的统计特征可以包括初始分割结果或当前腐蚀结果中细胞的总数，还可以包括初始分割结果或当前腐蚀结果中细胞的分布情况，例如各个子图
区域中细胞的数量。
29.虽然实例分割过程中对于重叠细胞的细胞边缘分割效果欠佳，但是实例分割所识别出的细胞数量或细胞分布（即统计特征）是可靠的，因此当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征之间的差异可以表征当前腐蚀结果中重叠细胞的分割情况与实际情况之间的差距。此处，当前腐蚀结果中重叠细胞的分割情况与实际情况之间的差距越大，越需要进行效果更明显的腐蚀操作，以更大程度地分离重叠细胞，相应地，当前腐蚀结果中重叠细胞的分割情况与实际情况之间的差距越小，表明大部分的重叠细胞已被分离，仅有少部分重叠部分较不明显的重叠细胞，故需要进行更精细的细胞分离操作，同时尽可能保留细胞的原有形态，以保证细胞边缘的平滑性。
30.具体而言，在腐蚀过程中，根据当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征之间的差异调整腐蚀过程中使用的卷积参数，例如卷积核尺寸和卷积步长。其中，当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征之间的差异越大，可以使用更大的卷积参数以进行更快速、效果更明显的迭代腐蚀操作，以将大部分重叠细胞迅速分离。随着当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征之间的差异变小，可以缩小卷积参数以进行更精细、幅度更小的迭代腐蚀操作，在分离重叠细胞的基础上尽量保留细胞的原有形态从而保证细胞边缘的平滑。在当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征相匹配，即，当前腐蚀结果中的细胞总数与初始分割结果中的细胞总数相同，或当前腐蚀结果中的细胞分布与初始分割结果中的细胞分布相同时，即可停止迭代腐蚀，得到最优腐蚀结果。其中，最优腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征相匹配。
31.可见，通过对初始分割结果对应的二值图像进行自适应迭代腐蚀处理，得到待分割细胞图像对应的最优腐蚀结果后，再进行后续的精细分割操作，一方面可以对初始分割结果中重叠的细胞进行分割，较大程度地将粘连重叠的细胞分离开，从而优化后续分水岭算法的效果，避免由于细胞重叠导致分水岭算法失效；另一方面，由于自适应迭代腐蚀过程中使用的卷积参数能够随着当前腐蚀结果中细胞的统计特征与初始分割结果中细胞的统计特征之间差异的减少而减小，因此可以避免腐蚀过度，从而尽可能保留原有细胞的形态，后续分水岭算法在最优腐蚀结果的指引下进行精细分割时，可以分割出更符合细胞原有形态的平滑边缘。
32.步骤140，利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
33.具体地，基于上述最优腐蚀结果，利用分水岭算法，对初始分割结果进行精细分割。其中，可以将最优腐蚀结果中得到的细胞区域的标识作为函数参数传递给分水岭算法，分水岭算法在进行精细分割时，可以以最优腐蚀结果中得到的细胞区域作为初始值进行精细分割，得到待分割细胞图像的最终分割结果。由于自适应迭代腐蚀操作较大程度地将粘连重叠的细胞分离开，且尽可能地保留了原有细胞的形态，因此，在自适应迭代腐蚀后得到的最优腐蚀结果的指引下，分水岭算法可以分割出更精确且更符合细胞原有形态的平滑边缘，如图2右侧所示。
34.本发明实施例提供的方法，在利用细胞分割模型对待分割细胞图像进行实例分割，得到初始分割结果后，对初始分割结果进行自适应迭代腐蚀处理，得到最优腐蚀结果，
其中，自适应迭代腐蚀处理过程中的卷积参数会随着初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，以较大程度地将粘连重叠的细胞分离开，同时尽可能地保留原有细胞的形态，再在最优腐蚀结果的指引下利用分水岭算法进行精细分割，可以分割出更精确且更符合细胞原有形态的平滑边缘，提高了细胞边缘分割的准确性。
35.基于上述实施例，所述卷积参数包括卷积核尺寸；图3为本发明实施例提供的自适应迭代腐蚀方法的流程示意图，如图3所示，步骤130具体包括：步骤131，确定以初始卷积核尺寸为首的多个尺寸递减的卷积核尺寸；步骤132，依次基于各个卷积核尺寸，对上一腐蚀处理结果进行迭代腐蚀，得到当前腐蚀处理结果；其中，每次迭代腐蚀会在达到对应卷积核尺寸的迭代停止条件时终止；步骤133，将当前腐蚀处理结果作为所述最优腐蚀结果；其中，首次腐蚀时，基于所述初始卷积核尺寸对所述初始分割结果进行腐蚀；任一卷积核尺寸的迭代停止条件为当前腐蚀处理结果中的细胞数量与所述初始分割结果中的细胞数量之间的差值小于所述任一卷积核尺寸对应的预设差值阈值，且所述任一卷积核尺寸越小，所述任一卷积核尺寸对应的预设差异阈值越小。
36.具体地，可以设置多个以初始卷积核尺寸为首的尺寸递减的卷积核尺寸，例如7*7、5*5和3*3，其中7*7为初始卷积核尺寸。其中，初始卷积核尺寸可以基于待分割细胞图像中细胞的形态确定得到，细胞越大，初始卷积核尺寸可以设置得越大。
37.依次基于各个卷积核尺寸，对上一腐蚀处理结果进行迭代腐蚀，得到当前腐蚀处理结果，其中，每次迭代腐蚀会在达到对应卷积核尺寸的迭代停止条件时终止。即，按照尺寸递减的顺序，基于当前卷积核尺寸对上一腐蚀处理结果进行一次迭代腐蚀，此次迭代腐蚀会在达到当前卷积核尺寸的迭代停止条件后终止，并将下一个卷积核尺寸作为当前卷积核尺寸进行又一次迭代腐蚀。以7*7、5*5和3*3三个卷积核尺寸为例，首先基于7*7卷积核，迭代对上一腐蚀处理结果进行腐蚀，直至达到7*7卷积核的迭代停止条件，其中，首次腐蚀时的腐蚀对象为实例分割得到的初始分割结果；随后，基于5*5卷积核，接着迭代对上一腐蚀处理结果进行腐蚀，直至达到5*5卷积核的迭代停止条件；最后，基于3*3卷积核，继续迭代对上一腐蚀处理结果进行腐蚀，直至达到3*3卷积核的迭代停止条件。当基于最后一个卷积核尺寸进行的迭代腐蚀终止后，将得到的当前腐蚀处理结果作为最优腐蚀结果。
38.此处，任一卷积核尺寸的迭代停止条件为当前腐蚀处理结果中的细胞数量与初始分割结果中的细胞数量之间的差值小于该卷积核尺寸对应的预设差值阈值，且该卷积核尺寸越小，该卷积核尺寸对应的预设差异阈值越小。另外，最后一个卷积核尺寸对应的预设差异阈值可以设置为0，表明需要在当前腐蚀处理结果中细胞的统计特征与初始分割结果中细胞的统计特征相匹配时停止此次迭代。例如，基于7*7卷积核进行迭代腐蚀时，迭代停止条件可以为当前腐蚀处理结果中的细胞数量与初始分割结果中的细胞数量之间的差值小于初始分割结果中细胞数量的10%；基于5*5卷积核进行迭代腐蚀时，迭代停止条件可以为当前腐蚀处理结果中的细胞数量与初始分割结果中的细胞数量之间的差值小于初始分割结果中细胞数量的5%；基于3*3卷积核进行迭代腐蚀时，迭代停止条件可以为当前腐蚀处理结果中的细胞数量等于初始分割结果中的细胞数量。
39.基于上述任一实施例，图4为本发明实施例提供的初始卷积核确定方法的流程示意图，如图4所示，所述初始卷积核尺寸是基于如下步骤确定的：步骤410，基于所述初始分割结果中各像素所属细胞的标签信息，统计所述初始分割结果中各细胞的面积；步骤420，基于所述初始分割结果中各细胞的面积，计算各细胞的平均面积；步骤430，基于所述各细胞的平均面积，确定所述初始卷积核尺寸。
40.具体地，基于初始分割结果中各像素所属细胞的标签信息，可以确定初始分割结果中的各个细胞包含的像素数量，从而统计得到各个细胞的面积。根据各个细胞的面积，计算各个细胞的平均面积。基于各个细胞的平均面积，可以确定初始卷积核尺寸。其中，各个细胞的平均面积越大，确定的初始卷积核尺寸越大。由于初始卷积核尺寸与各个细胞的平均面积正相关，因此可以选择与各个细胞大小相适应的初始卷积核尺寸，避免选取过大的初始卷积核尺寸导致过度腐蚀，同时还可以在保证腐蚀效果的同时选择尽量大的初始卷积核尺寸，以提高自适应迭代腐蚀的效率。
41.基于上述任一实施例，步骤140具体包括：精细分割步骤：利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到当前的精细分割结果；反馈调节步骤：若所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征不匹配，则调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，得到当前的精细分割结果；循环步骤：重复执行所述反馈调节步骤，直至所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征相匹配；结果确定步骤：确定当前的精细分割结果为所述最终分割结果。
42.具体地，如图5所示，在进行自适应迭代腐蚀操作得到最优腐蚀结果后，执行精细分割步骤，利用分水岭算法，基于最优腐蚀结果对初始分割结果进行精细分割，得到当前的精细分割结果。随后，获取当前的精细分割结果中细胞的统计特征。此处，可以利用初始分割结果中细胞的统计特征作为调整标准，根据初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征之间的差异，判断是否需要重新进行自适应迭代腐蚀处理和精细分割操作，以进一步提高细胞边缘分割的准确性。
43.具体而言，若初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征不匹配，例如初始分割结果中的细胞总数和当前的精细分割结果中的细胞总数不同，或初始分割结果中的细胞分布和当前的精细分割结果中的细胞分布不同，则表明此次精细分割结果不准确，需要对最优腐蚀结果进行形态学调整，执行反馈调节步骤。其中，可以调整自适应迭代腐蚀处理过程中的初始卷积核尺寸，然后基于调整后的初始卷积核尺寸，对初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于新的最优腐蚀结果，再执行上述精细分割步骤，得到当前的精细分割结果。
44.若当前的精细分割结果中细胞的统计特征与初始分割结果中细胞的统计特征仍不匹配，则执行循环步骤，重复执行上述反馈调节步骤，继续调整初始卷积核尺寸，重新执行自适应迭代腐蚀处理和精细分割步骤，直至初始分割结果中细胞的统计特征和当前的精
细分割结果中细胞的统计特征相匹配。随后，将当前的精细分割结果确定为最终分割结果。
45.基于上述任一实施例，图6为本发明实施例提供的反馈调节步骤的流程示意图，如图6所示，所述调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，具体包括：步骤610，对所述待分割细胞图像进行子图划分，得到所述待分割细胞图像的多个子图区域；步骤620，获取所述初始分割结果和所述精细分割结果中对应各个子图区域的细胞数量；步骤630，若所述初始分割结果和所述精细分割结果中对应任一子图区域的细胞数量不匹配，则调整所述初始卷积核尺寸，并对所述初始分割结果中对应所述任一子图区域的部分重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果；步骤640，基于所述新的最优腐蚀结果，对所述初始分割结果中对应所述任一子图区域的部分重新进行精细分割。
46.具体地，可以对待分割细胞图像进行子图划分，得到待分割细胞图像的多个子图区域。随后，获取初始分割结果和精细分割结果中对应各个子图区域的细胞数量。此处，初始分割结果和精细分割结果均与待分割细胞图像相对应，因此，将待分割细胞图像划分为多个子图区域后，也可以获取初始分割结果和精细分割结果中对应各个子图区域的部分，并获取初始分割结果和精细分割结果中对应各个子图区域的部分包含的细胞数量。
47.若初始分割结果和精细分割结果中对应任一子图区域i的细胞数量不匹配，表明针对该子图区域i的精细分割不准确，需要进行调整。因此，调整自适应迭代腐蚀过程中的初始卷积核尺寸，并对初始分割结果中对应该子图区域i的部分重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，再基于上述新的最优腐蚀结果，对初始分割结果中对应该子图区域i的部分重新进行精细分割。此处，将待分割细胞图像划分为多个更小的子图区域后，判断各个子图区域的精细分割是否准确，从而进队精细分割不准确的部分子图区域重新进行自适应迭代腐蚀和精细分割，减少了自适应迭代腐蚀和精细分割的计算量，提高了细胞边缘分割的效率。
48.基于上述任一实施例，所述调整所述初始卷积核尺寸，具体包括：基于多个预设卷积核尺寸，缩小所述初始卷积核尺寸。
49.具体地，当初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征不匹配，或，初始分割结果和精细分割结果中对应任一子图区域的细胞数量不匹配时，表明精细分割不准确。导致精细分割不准确的原因可能在于自适应迭代腐蚀处理时腐蚀过度，使得最优腐蚀结果中的细胞状态和待分割细胞图像中细胞的实际状态有所差别，故而使得精细分割结果不准确。因此，为了提高细胞边缘分割的准确性，可以缩小自适应迭代腐蚀处理过程中的初始卷积核尺寸，从较小的卷积核尺寸开始迭代腐蚀，以避免腐蚀过度。此处，可以预先设置多个预设卷积核尺寸，然后从中选取比当前的初始卷积核尺寸更小的卷积核尺寸，作为调整后的初始卷积核尺寸。
50.基于上述任一实施例，图7为本发明实施例提供的精细分割方法的流程示意图之二，如图7所示，所述利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精
细分割，具体包括：步骤710，对所述最优腐蚀结果进行连通域分析，得到所述最优腐蚀结果中的多个连通域；步骤720，基于所述最优腐蚀结果中的多个连通域，利用分水岭算法对所述初始分割结果进行精细分割。
51.具体地，为了提高细胞边缘分割的精准性，分水岭算法可以根据细胞的形态特征考虑分割形状的规则性，从而提高分割的最终性能。此处，可以通过连通域分割方法，对最优腐蚀结果进行连通域分析，得到多个连通域，并给予上述最优腐蚀结果每个独立连通域一个标识。其中，单个独立连通域对应自适应腐蚀处理得到的一个细胞。
52.随后，基于上述最优腐蚀结果中的多个连通域所指示的细胞形态特征，利用分水岭算法对初始分割结果进行精细分割。其中，可以将最优腐蚀结果中每个独立连通域的标识将作为函数参数传递给分水岭分割算法进行精细分割，得到相应的细胞分割结果。
53.基于上述任一实施例，所述细胞分割模型是基于如下步骤确定的：获取训练样本，并将所述训练样本划分为训练集和测试集；所述训练样本包括样本细胞图像，或，包括样本细胞图像和所述样本细胞图像中各细胞的掩膜标注信息；基于所述训练集对分割模型进行训练，并在训练过程中利用所述测试集对所述分割模型进行阶段性的性能测试，保存性能测试结果最优时的分割模型作为所述细胞分割模型。
54.具体地，首先需要获取大量的训练样本。训练样本可以包括样本细胞图像，也可以既包括样本细胞图像又包括样本细胞图像中各细胞的掩膜标注信息。其中，训练样本中是否需要样本细胞图像中各细胞的掩膜标注信息，可以根据细胞分割模型是采用全监督方式训练还是采用半监督方式训练而定。此处，可以获取大量细胞显微图像作为样本细胞图像以供训练细胞分割模型。掩膜标注信息可以由具有专业资质的标注人员通过标注软件来形成。例如，可以采用标注软件如labelme等开源软件，标注出样本细胞图像中每个细胞的掩膜图，掩膜图对于细胞核边缘的描绘应尽可能的准确。
55.在获取训练样本之后，还可以对训练样本进行清理与预处理，去除一些成像不合格或者标注不完整的数据。随后，根据上述掩膜标注信息，确定样本细胞图像中各个细胞的位置及边缘以生成细胞核区域的掩膜二值图，并合理划分数据集，将训练样本划分为训练集和测试集。
56.然后，利用上述训练集和测试集对选取的分割模型进行训练及调优。其中，分割模型可以采用mask-rcnn网络来对细胞进行实例分割，具体可以采用resnet 50 特征金字塔fpn作为特征提取网络，采用了roi align模块来实现小目标的感兴趣区域提取，最终通过多分支网络，同步实现对于感兴趣区域内的前景与背景的分割，以及前景区域的分类。在分割模型的训练过程中，可以将训练集中的样本细胞图像、上述样本细胞图像中各细胞的掩膜二值图像等作为元数据输入到分割模型中，进行多次（例如100次）epoches训练。
57.在训练过程中，可以按照一定的程序在测试集上进行阶段性的性能测试（这里可以采用map），最终根据性能表现结果合理选择和存储最优模型，将性能最优的分割模型作为细胞分割模型，以提高细胞分割模型的分割性能。然后，将细胞分割模型封装并部署到应用服务器或者工作站中，并建立数据推断流程；数据推断流程一般应包括上述的数据预处
理及组织过程，以及模型的调用和推断，以及模型推断结果的输出和存储等。其中性能最优的分割模型应该在测试集上的表现大于某些阈值，例如在map上表现大于0.9，在dice上大于0.95等。同时，性能最优的分割模型还应具有较好的泛化能力，例如，在测试集和训练集的性能表现上相差不大于10%。
58.基于上述任一实施例，所述获取训练样本，具体包括：获取细胞显微图像；若所述细胞显微图像的通道数大于所述细胞分割模型的输入条件对应的预设通道数，则将所述细胞显微图像分解为多个低通道图像；对所述多个低通道图像进行图像质量评估，得到所述多个低通道图像的质量评估得分；基于所述多个低通道图像的质量评估得分，获取质量评估得分最高的多个低通道图像；将所述质量评估得分最高的多个低通道图像组合为所述样本细胞图像；其中，所述质量评估得分最高的多个低通道图像组合后的通道数等于所述预设通道数。
59.具体地，根据所选择的细胞分割模型的网络结构，其输入条件中要求的输入图像的通道数有所不同，例如目前大部分的实例分割网络通常会要求三通道图像输入。然而，细胞图像一般会多于三通道，因而需要将多通道的细胞图像分解为满足细胞分割模型输入条件要求的较低通道的图像，得到细胞分割模型的输入数据。
60.因此，可以获取细胞显微图像，然后判断上述细胞显微图像的通道数是否大于细胞分割模型的输入条件对应的预设通道数（例如三通道）。若上述细胞显微图像的通道数大于上述预设通道数，则将细胞显微图像逐层分解为多个低通道图像。其中，低通道图像可以为单通道图像。随后，对上述多个低通道图像进行图像质量评估，得到上述多个低通道图像的质量评估得分。此处，可以采用自动聚焦评判算法对上述多个低通道图像进行图像质量评估，例如采用无参图像清晰度评价函数，如能量梯度函数、tengengrad梯度函数、lapacian梯度函数以及一些统计函数，如信息熵及range函数等对多个低通道图像进行图像质量评估。
61.根据上述多个低通道图像的质量评估得分，获取质量评估得分最高的多个低通道图像，再将上述质量评估得分最高的多个低通道图像进行组合，得到通道数满足细胞分割模型输入条件的样本细胞图像。其中，上述质量评估得分最高的多个低通道图像组合后的通道数应当等于预设通道数。
62.基于上述任一实施例，图8为本发明实施例提供的细胞边缘分割方法的详细示意图，该方法可以用于循环肿瘤细胞（ctc）的精准分割，可应用于ctc细胞中cacs细胞的检出，但具体的检测目标可以根据实际应用场景确定，本发明实施例对此不作具体限定。如图8所示，该方法包括：步骤810，获取原始训练样本，上述原始训练样本包括细胞显微图像和细胞显微图像中各细胞的掩膜标注信息。
63.其中，可以选取dapi视野图像作为上述细胞显微图像，基于此实现对细胞的精准检测与边缘分割。掩膜标注信息由具有专业资质的标注人员通过标注软件来形成，例如可以采用标注软件如labelme等开源软件，掩膜标注信息中应至少包括细胞显微图像中每个
细胞的掩膜图，掩膜图对于细胞核边缘的描绘应尽可能的准确。根据上述掩膜标注信息，确定细胞显微图像中细胞核的区域，从图像中提取细胞核区域可以生成细胞核的掩膜二值图。获取原始训练样本之后，可以对原始训练样本进行清理与预处理，去除一些成像不合格或者标注不完整的数据。
64.步骤820，根据掩膜标注信息，确定细胞显微图像中各个细胞的位置及边缘以生成细胞核区域的掩膜二值图，并合理划分数据集。
65.此步骤的目的在于组织细胞分割模型的输入数据。为了提高细胞分割模型的性能，细胞分割模型的图像输入为经过质量控制的图像，送入模型进行训练以及进行实际分割的图像皆为一次扫描中质量最优的图像数据。此处，细胞分割模型一般会要求三通道图像输入，而细胞显微图像一般会多于三通道，因而需要采用自动图像质量评估的方式来选取质量最佳的多通道图像来组成模型的数据输入。其中，自动图像评估方式可以包括无参图像清晰度评价函数，如能量梯度函数，tengengrad梯度函数，lapacian梯度函数 以及一些统计函数，如信息熵及range函数等。另一方面，需要把上述掩膜标注信息转化为细胞坐标或者掩膜二值图，并将数据整理成适合于进一步网络训练时需要的元数据格式。另外，在样本收集的过程中，需要注意训练集和测试集的分离，尽量做到测试集和实际样品分布与形态相一致。
66.步骤830，利用样本细胞图像及其掩膜标注信息对细胞分割模型进行训练和调优。
67.其中，利用训练集中的样本细胞图像及其掩膜标注信息对细胞分割模型进行训练，直至达到预设终止条件，以完成细胞分割模型的训练。上述细胞分割网络是一种端到端的实例分割网络，可以为单阶段实例分割或者多阶段实例分割网络，其要点在于能够有效分割出各个独立细胞核，并给予分别的标识。
68.此外，在模型训练过程中，可以建立合理的评估指标，在测试集上对模型进行阶段性性能评估，选择出性能最优的模型作为投入使用的细胞分割模型。其中，性能最优的模型指的是在测试集上的表现大于某些阈值，如在map上表现大于0.9，在dice上大于0.95等的模型。同时，该模型也应具有较好的泛化能力，例如在测试集和训练集的性能表现相差应不大于10%。
69.步骤840，对最优模型进行封装和部署，并对新的待分割细胞图像进行实例分割。
70.具体的，将上述最优模型进行封装并部署到应用服务器或者工作站中，并建立数据推断流程；数据推断流程一般应包括上述的数据预处理及组织过程，以及模型的调用和推断，以及模型推断结果的输出和存储等。通过训练完成的细胞分割模型对待分割细胞图像进行实例分割，可以获得待分割细胞图像视野中每个细胞的初始分割结果。
71.步骤850，对细胞分割模型输出的初始分割结果进行形态学后处理，获取待分割细胞图像中每个细胞核的腐蚀结果与标识。
72.此处，首先根据初始分割结果统计出待分割细胞图像中细胞的平均面积，根据细胞的平均面积，设定自适应迭代腐蚀的初始卷积核尺寸。设定以初始卷积核尺寸为首的多个递减的卷积核尺寸，并以此基于上述多个卷积核尺寸进行迭代腐蚀。
73.以5*5和3*3两个递减的卷积核尺寸为例，首先设定一个while循环，使用第一个卷积核尺寸（即初始卷积核尺寸5*5）对初始分割结果的二值图或上一腐蚀结果（首轮腐蚀时的腐蚀对象为初始分割结果的二值图，第二轮及以后腐蚀时的腐蚀对象为上一轮腐蚀得到
的腐蚀结果，即上一腐蚀结果）进行迭代腐蚀处理。当腐蚀后统计出的当前腐蚀结果中连通域的数量与初始分割结果中细胞的数量间的差值小于5*5卷积核尺寸对应的阈值（例如初始分割结果中细胞数量的5%）时，停止迭代，终止该while循环。
74.随后，设定另一个while循环，选择下一卷积核尺寸3*3，利用该卷积核尺寸对上一腐蚀结果进行迭代腐蚀处理，直至腐蚀后统计出的当前腐蚀结果中连通域的数量等于初始分割结果中细胞的数量，停止迭代，终止该while循环，得到最优腐蚀结果。
75.完成自适应迭代腐蚀处理后，对最优腐蚀结果进行连通域分析，得到多个连通域，并给每一个独立连通域赋予一个唯一标识，从而完成对于最优腐蚀结果中各细胞的标识。
76.步骤860，基于最优腐蚀结果中各细胞的标识，通过分水岭算法对每个细胞进行进一步精细分割，获取待分割细胞图像的最终分割结果。
77.此处，分水岭算法可以根据细胞的形态特征考虑分割形状的规则性以提高分割的最终性能，因此将最优腐蚀结果中各细胞的标识作为函数参数传递给分水岭分割算法以获取精细分割结果。精细分割后，还可以对精细分割结果进行细胞数量的统计，如果统计得到的细胞数量与初始分割结果中的细胞数量存在差异，则将对最优腐蚀结果进行形态学调整，并重新进行上述分水岭算法，直到精细分割结果中的细胞数量与初始分割结果中的细胞数量一致。其中，对最优腐蚀结果进行形态学调整时，可以缩小自适应迭代腐蚀过程中的初始卷积核尺寸，随后再重复进行自适应迭代腐蚀，得到新的最优腐蚀结果。
78.通过上述方式，可以分割出更精确且更符合细胞原有形态的平滑边缘，提高了细胞边缘分割的准确性，因此更适用于cacs细胞分割等边缘敏感的分割任务。
79.下面对本发明提供的细胞边缘分割装置进行描述，下文描述的细胞边缘分割装置与上文描述的细胞边缘分割方法可相互对应参照。
80.基于上述任一实施例，图9为本发明实施例提供的细胞边缘分割装置的结构示意图，如图9所示，该装置包括：图像确定单元910、实例分割单元920、自适应腐蚀单元930和精细分割单元940。
81.其中，图像确定单元910用于确定待分割细胞图像；实例分割单元920用于利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；自适应腐蚀单元930用于基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；精细分割单元940用于利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
82.本发明实施例提供的装置，在利用细胞分割模型对待分割细胞图像进行实例分割，得到初始分割结果后，对初始分割结果进行自适应迭代腐蚀处理，得到最优腐蚀结果，其中，自适应迭代腐蚀处理过程中的卷积参数会随着初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，以较大程度地将粘连重叠的细胞分离开，同时尽可能地保留原有细胞的形态，再在最优腐蚀结果的指引下利用分水岭算法进
行精细分割，可以分割出更精确且更符合细胞原有形态的平滑边缘，提高了细胞边缘分割的准确性。
83.基于上述实施例，所述卷积参数包括卷积核尺寸；自适应腐蚀单元930具体用于：确定以初始卷积核尺寸为首的多个尺寸递减的卷积核尺寸；依次基于各个卷积核尺寸，对上一腐蚀处理结果进行迭代腐蚀，得到当前腐蚀处理结果；其中，每次迭代腐蚀会在达到对应卷积核尺寸的迭代停止条件时终止；将当前腐蚀处理结果作为所述最优腐蚀结果；其中，首次腐蚀时，基于所述初始卷积核尺寸对所述初始分割结果进行腐蚀；任一卷积核尺寸的迭代停止条件为当前腐蚀处理结果中的细胞数量与所述初始分割结果中的细胞数量之间的差值小于所述任一卷积核尺寸对应的预设差值阈值，且所述任一卷积核尺寸越小，所述任一卷积核尺寸对应的预设差异阈值越小。
84.基于上述任一实施例，所述初始卷积核尺寸是基于如下步骤确定的：基于所述初始分割结果中各像素所属细胞的标签信息，统计所述初始分割结果中各细胞的面积；基于所述初始分割结果中各细胞的面积，计算各细胞的平均面积；基于所述各细胞的平均面积，确定所述初始卷积核尺寸。
85.基于上述任一实施例，精细分割单元940具体用于执行：精细分割步骤：利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到当前的精细分割结果；反馈调节步骤：若所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征不匹配，则调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，得到当前的精细分割结果；循环步骤：重复执行所述反馈调节步骤，直至所述初始分割结果中细胞的统计特征和当前的精细分割结果中细胞的统计特征相匹配；结果确定步骤：确定当前的精细分割结果为所述最终分割结果。
86.基于上述任一实施例，所述调整所述初始卷积核尺寸，对所述初始分割结果重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果，并基于所述新的最优腐蚀结果，执行所述精细分割步骤，具体包括：对所述待分割细胞图像进行子图划分，得到所述待分割细胞图像的多个子图区域；获取所述初始分割结果和所述精细分割结果中对应各个子图区域的细胞数量；若所述初始分割结果和所述精细分割结果中对应任一子图区域的细胞数量不匹配，则调整所述初始卷积核尺寸，并对所述初始分割结果中对应所述任一子图区域的部分重新进行自适应迭代腐蚀处理，得到新的最优腐蚀结果；基于所述新的最优腐蚀结果，对所述初始分割结果中对应所述任一子图区域的部分重新进行精细分割。
87.基于上述任一实施例，所述调整所述初始卷积核尺寸，具体包括：
基于多个预设卷积核尺寸，缩小所述初始卷积核尺寸。
88.基于上述任一实施例，所述利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，具体包括：对所述最优腐蚀结果进行连通域分析，得到所述最优腐蚀结果中的多个连通域；基于所述最优腐蚀结果中的多个连通域，利用分水岭算法对所述初始分割结果进行精细分割。
89.基于上述任一实施例，所述细胞分割模型是基于如下步骤确定的：获取训练样本，并将所述训练样本划分为训练集和测试集；所述训练样本包括样本细胞图像，或，包括样本细胞图像和所述样本细胞图像中各细胞的掩膜标注信息；基于所述训练集对分割模型进行训练，并在训练过程中利用所述测试集对所述分割模型进行阶段性的性能测试，保存性能测试结果最优时的分割模型作为所述细胞分割模型。
90.基于上述任一实施例，所述获取训练样本，具体包括：获取细胞显微图像；若所述细胞显微图像的通道数大于所述细胞分割模型的输入条件对应的预设通道数，则将所述细胞显微图像分解为多个低通道图像；对所述多个低通道图像进行图像质量评估，得到所述多个低通道图像的质量评估得分；基于所述多个低通道图像的质量评估得分，获取质量评估得分最高的多个低通道图像；将所述质量评估得分最高的多个低通道图像组合为所述样本细胞图像；其中，所述质量评估得分最高的多个低通道图像组合后的通道数等于所述预设通道数。
91.图10示例了一种电子设备的实体结构示意图，如图10所示，该电子设备可以包括：处理器(processor)1010、通信接口(communications interface)1020、存储器(memory)1030和通信总线1040，其中，处理器1010，通信接口1020，存储器1030通过通信总线1040完成相互间的通信。处理器1010可以调用存储器1030中的逻辑指令，以执行细胞边缘分割方法，该方法包括：确定待分割细胞图像；利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
92.此外，上述的存储器1030中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施
例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器（rom，read-only memory）、随机存取存储器（ram，random access memory）、磁碟或者光盘等各种可以存储程序代码的介质。
93.另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的细胞边缘分割方法，该方法包括：确定待分割细胞图像；利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
94.又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的细胞边缘分割方法，该方法包括：确定待分割细胞图像；利用细胞分割模型，对所述待分割细胞图像进行实例分割，得到所述待分割细胞图像的初始分割结果；基于所述初始分割结果中细胞的统计特征，对所述初始分割结果进行自适应迭代腐蚀处理，得到所述待分割细胞图像对应的最优腐蚀结果；其中，所述自适应迭代腐蚀处理过程中的卷积参数会随着所述初始分割结果中细胞的统计特征与当前腐蚀结果中细胞的统计特征之间差异的减少而减小，所述最优腐蚀结果中细胞的统计特征与所述初始分割结果中细胞的统计特征相匹配；利用分水岭算法，基于所述最优腐蚀结果，对所述初始分割结果进行精细分割，得到所述待分割细胞图像的最终分割结果。
95.以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
96.通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施例或者实施例的某些部分所述的方法。
97.最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。