CN115359056A - 分裂细胞检测方法、装置和计算机设备 - Google Patents

分裂细胞检测方法、装置和计算机设备 Download PDF

Info

Publication number
CN115359056A
CN115359056A CN202211279428.5A CN202211279428A CN115359056A CN 115359056 A CN115359056 A CN 115359056A CN 202211279428 A CN202211279428 A CN 202211279428A CN 115359056 A CN115359056 A CN 115359056A
Authority
CN
China
Prior art keywords
image
dividing cells
cell
candidate regions
candidate
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
CN202211279428.5A
Other languages
English (en)
Other versions
CN115359056B (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.)
Zhejiang Huanuokang Technology Co ltd
Original Assignee
Zhejiang Huanuokang Technology Co ltd
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 Zhejiang Huanuokang Technology Co ltd filed Critical Zhejiang Huanuokang Technology Co ltd
Priority to CN202211279428.5A priority Critical patent/CN115359056B/zh
Publication of CN115359056A publication Critical patent/CN115359056A/zh
Application granted granted Critical
Publication of CN115359056B publication Critical patent/CN115359056B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/56Extraction of image or video features relating to colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/764Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
    • G06V10/765Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects using rules for classification or partitioning the feature space
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/82Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/695Preprocessing, e.g. image segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698Matching; Classification
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Abstract

本申请涉及一种分裂细胞检测方法、装置和计算机设备,通过获取细胞染色图像,在细胞染色图像中确定分裂细胞的多个候选区域;提取各候选区域的图像特征、轮廓特征和结构特征;根据多个候选区域的图像特征、轮廓特征和结构特征,对多个候选区域进行分类以得到分裂细胞所在区域,解决了相关技术中存在分裂细胞检出率低的问题,提高了分裂细胞的检出率。

Description

分裂细胞检测方法、装置和计算机设备
技术领域
本申请涉及图像处理技术领域,特别是涉及一种分裂细胞检测方法、装置和计算机设备。
背景技术
传统的分裂细胞检出方法,通常由病理学家利用显微镜在高放大倍率下对HE(Hematoxylin and eosin,苏木精伊红)染色的组织学载玻片制剂中的分裂细胞进行检查和计数。然而人工检查和计数存在主观性,非常耗时且容易出错,为此,相关技术提出了一种细胞有丝分裂检测方法,其设计要点是利用图像矩阵的列向量构造成新矩阵,再对新矩阵低秩分解,分解为大小相同的低秩矩阵和稀疏矩阵。分解的稀疏矩阵的列向量被用来构造新矩阵,迭代新矩阵的地址分解。最终稀疏部分认为是分裂细胞部分,低秩部分认为是非分裂细胞部分。但是,这种方法容易受组织染色差异影响,不适用于检测所有时期的分裂细胞,导致分裂细胞检出率较低。
针对相关技术中存在分裂细胞检出率低的问题,目前还没有提出有效的解决方案。
发明内容
基于此,有必要针对上述技术问题,提供一种能够提高分裂细胞检出率的分裂细胞检测方法、装置和计算机设备。
第一方面,本申请提供了一种分裂细胞检测方法,所述方法包括:
获取细胞染色图像,在所述细胞染色图像中确定分裂细胞的多个候选区域;
提取各所述候选区域的图像特征、轮廓特征和结构特征;
根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域。
在其中一些实施例中,在所述细胞染色图像中确定分裂细胞的多个候选区域包括:
提取所述细胞染色图像中的最佳染色向量,将所述最佳染色向量输入至YOLOv5模型进行检测,得到所述分裂细胞的多个候选区域。
在其中一些实施例中,提取所述细胞染色图像中的最佳染色向量包括:
将所述细胞染色图像的多种颜色分量转换成对应的光密度值;
根据奇异值对所述多种颜色分量进行分解;
根据各所述光密度值与奇异值主分解方向之间的角度确定所述最佳染色向量。
在其中一些实施例中,提取各所述候选区域的图像特征包括:
将所述多个候选区域输入至Resnet-18网络,基于所述Resnet-18网络,对各所述候选区域的细胞核的颜色进行特征提取,得到所述图像特征。
在其中一些实施例中,提取各所述候选区域的轮廓特征包括:
对各所述候选区域进行边缘检测,得到边缘检测结果;
将所述多个候选区域输入至Resnet-18网络,基于所述Resnet-18网络,对所述边缘检测结果的形态进行特征提取,得到所述轮廓特征。
在其中一些实施例中,提取各所述候选区域的结构特征包括:
根据Hessian矩阵获取各所述候选区域的特征值;
根据所述特征值计算与细胞结构相关的形状指数和归一化曲率指数;
根据所述形状指数和归一化曲率指数,得到所述结构特征。
在其中一些实施例中,根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域包括:
将各所述候选区域的图像特征、轮廓特征和结构特征,输入至随机森林模型进行决策,根据决策结果确定所述多个候选区域中分裂细胞所在区域。
在其中一些实施例中,获取细胞染色图像包括:
获取拍摄染色的组织学载玻片得到的第一图像;
将所述第一图像裁剪为1024*1024像素的所述细胞染色图像。
第二方面,本申请提供了一种分裂细胞检测装置,所述装置包括:
获取模块,用于获取细胞染色图像,在所述细胞染色图像中确定分裂细胞的多个候选区域;
提取模块,用于提取各所述候选区域的图像特征、轮廓特征和结构特征;
决策模块,用于根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域。
第三方面,本申请提供了一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述第一方面所述分裂细胞检测方法的步骤。
第四方面,本申请提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述第一方面所述分裂细胞检测方法的步骤。
上述分裂细胞检测方法、装置和计算机设备,通过获取细胞染色图像,在细胞染色图像中确定分裂细胞的多个候选区域;提取各候选区域的图像特征、轮廓特征和结构特征;根据多个候选区域的图像特征、轮廓特征和结构特征,对多个候选区域进行分类以得到分裂细胞所在区域,解决了相关技术中存在分裂细胞检出率低的问题,提高了分裂细胞的检出率。
附图说明
图1是本申请一实施例中分裂细胞检测方法的终端的硬件结构框图;
图2是本申请一实施例中分裂细胞检测方法的流程图;
图3是本申请一实施例中宫颈鳞状上皮组织非分裂细胞检测方法的流程图;
图4是本申请一实施例中分裂细胞检测装置的结构框图;
图5是本申请一实施例中计算机设备的内部结构图。
具体实施方式
为更清楚地理解本申请的目的、技术方案和优点,下面结合附图和实施例,对本申请进行了描述和说明。
除另作定义外,本申请所涉及的技术术语或者科学术语应具有本申请所属技术领域具备一般技能的人所理解的一般含义。在本申请中的“一”、“一个”、“一种”、“该”、“这些”等类似的词并不表示数量上的限制,它们可以是单数或者复数。在本申请中所涉及的术语“包括”、“包含”、“具有”及其任何变体,其目的是涵盖不排他的包含;例如,包含一系列步骤或模块(单元)的过程、方法和系统、产品或设备并未限定于列出的步骤或模块(单元),而可包括未列出的步骤或模块(单元),或者可包括这些过程、方法、产品或设备固有的其他步骤或模块(单元)。在本申请中所涉及的“连接”、“相连”、“耦接”等类似的词语并不限定于物理的或机械连接,而可以包括电气连接,无论是直接连接还是间接连接。在本申请中所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系,表示可以存在三种关系,例如,“A和/或B”可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。通常情况下,字符“/”表示前后关联的对象是一种“或”的关系。在本申请中所涉及的术语“第一”、“第二”、“第三”等,只是对相似对象进行区分,并不代表针对对象的特定排序。
在本实施例中提供的方法实施例可以在终端、计算机或者类似的运算装置中执行。比如在终端上运行,图1是本申请一实施例的分裂细胞检测方法的终端的硬件结构框图。如图1所示,终端可以包括一个或多个(图1中仅示出一个)处理器102和用于存储数据的存储器104,其中,处理器102可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置。上述终端还可以包括用于通信功能的传输设备106以及输入输出设备108。本领域普通技术人员可以理解,图1所示的结构仅为示意,其并不对上述终端的结构造成限制。例如,终端还可包括比图1中所示更多或者更少的组件,或者具有与图1所示出的不同配置。
存储器104可用于存储计算机程序,例如,应用软件的软件程序以及模块,如在本实施例中的分裂细胞检测方法对应的计算机程序,处理器102通过运行存储在存储器104内的计算机程序,从而执行各种功能应用以及数据处理,即实现上述的方法。存储器104可包括高速随机存储器,还可包括非易失性存储器,如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中,存储器104可进一步包括相对于处理器102远程设置的存储器,这些远程存储器可以通过网络连接至终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
传输设备106用于经由一个网络接收或者发送数据。上述的网络包括终端的通信供应商提供的无线网络。在一个实例中,传输设备106包括一个网络适配器(NetworkInterface Controller,简称为NIC),其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中,传输设备106可以为射频(Radio Frequency,简称为RF)模块,其用于通过无线方式与互联网进行通讯。
在一个实施例中,如图2所示,提供了一种分裂细胞检测方法,以该方法应用于图1中的终端为例进行说明,包括以下步骤:
步骤S201,获取细胞染色图像,在细胞染色图像中确定分裂细胞的多个候选区域。
细胞染色图像,指的是对染色的组织学载玻片进行拍摄后得到的图像。例如,获取经HE(Hematoxylin and eosin,苏木精伊红)染色的组织学载玻片,采用高放大倍率摄像机对组织学载玻片进行拍摄,得到细胞染色图像。候选区域,指的是在细胞染色图像中对分裂细胞进行初步识别所得到的区域。
步骤S202,提取各候选区域的图像特征、轮廓特征和结构特征。
图像特征,指的是细胞核颜色,具体可以是基于像素级的RGB值。轮廓特征,指的是细胞核是否分裂成两部分的分裂信息、细胞核所处位置的排列信息。结构特征,指的是细胞核的形状和不规则度。
步骤S203,根据多个候选区域的图像特征、轮廓特征和结构特征,对多个候选区域进行分类以得到分裂细胞所在区域。
将各候选区域的图像特征、轮廓特征和结构特征,输入至随机森林模型进行决策,根据决策结果确定多个候选区域中分裂细胞所在区域,即将多个候选区域分为细胞分裂区域和非细胞分裂区域。一种可选的决策过程如下:利用训练好的21棵决策树对输入的图像特征、轮廓特征和结构特征进行决策,每棵树都将得到是否为细胞分裂区域的决策结果;采取投票的方式,获得最多决策树投票的类别为最终类别。
在上述步骤S201至S203中,在细胞染色图像中初步识别分裂细胞的候选区域,利用细胞染色图像中多个维度的特征,对初步识别的候选区域进行精细分类,进而得到分裂细胞所在区域,如此设置,能够充分利用像素信息和形状信息,有利于检测处于不同时期的分裂细胞,提高分裂细胞检出率和检出准确度。
相关技术利用细胞分裂区域显著特征和时空关联特性提取细胞分裂候选区域,再用方向梯度直方图对候选区域进行分类,其提取和确定候选区域的特征较为简单,容易造成漏检,以及不能很好地区分较为相似的分裂细胞和非分裂细胞。针对该问题,在一个实施例中,在细胞染色图像中确定分裂细胞的多个候选区域包括:提取细胞染色图像中的最佳染色向量,将最佳染色向量输入至YOLOv5模型进行检测,得到分裂细胞的多个候选区域。具体地,YOLOv5模型的相关参数可以作如下设置:conf-thres(置信度)设置为0.05,iou-thres(交并比)设置为0.3,模型大小设置为YOLOv5m。如此设置,提升了分裂细胞的检出率。
可选的,YOLOv5模架构由BackBone、Neck和YOLO输出端构成;其中,BackBone用于特征提取,将最佳染色向量图像中的物体信息通过卷积网络进行提取,以用于后续进行目标检测;Neck用于组合特征,增强网络的鲁棒性,加强物体检测能力,并且将这些特征传递给Head层进行预测;YOLO输出端用于进行候选区域的预测输出。具体地,BackBone包括Focus模块、SPP模块和CSP_X模块;其中,Focus模块用于采用切片操作把高分辨率的最佳染色向量图像拆分成多个低分辨率的最佳染色向量,即隔列采样和拼接;SPP模块用于进行空间金字塔池化,将任意大小的特征图转换成固定大小的特征向量;CSP_X模块用于增加残差结构,例如增加层与层之间反向传播的梯度值,避免因为加深而带来的梯度消失,从而可以提取到更细粒度的特征并且不用担心网络退化。Neck包括FPN和PAN;其中,FPN用于以自顶向下的方式,将高层特征通过上采样和低层特征做融合得到可预测的特征图;PAN用于FPN的基础上增加从下至上的特征融合,把浅层的定位信息传导到深层,增强多个尺度上的定位能力。YOLO输出端设置有Bounding Box损失函数和NMS非极大值抑制;其中,BoundingBox损失函数,根据真实检测框和模型预测输出框的吻合程度,以反向传播优化模型;NMS非极大值抑制,用于判断相邻网格识别的是否是同一物体,以消除掉多余检测框。
相关技术利用对图像矩阵的低秩分解检测分裂细胞,其容易受组织染色差异影响且不适用于检测所有时期的分裂细胞。针对该问题,在一个实施例中,提取细胞染色图像中的最佳染色向量包括:将细胞染色图像的多种颜色分量转换成对应的光密度值;根据奇异值对多种颜色分量进行分解;根据各光密度值与奇异值主分解方向之间的角度确定最佳染色向量。
如此设置,将RGB通道转换为光密度值,减少了组织染色差异对检测的影响。具体地,将细胞染色图像的RGB颜色转换为其对应的光密度值,并用奇异值对RGB颜色进行分解,利用每个光密度值与奇异值主分解方向之间的角度来确定最佳染色向量。RGB转换为光密度值的公式如下:
Figure DEST_PATH_IMAGE001
其中,OD代表光密度值,I代表经过归一化后的RGB颜色向量。
在一个实施例中,提取各候选区域的图像特征包括:将多个候选区域输入至Resnet-18网络,基于Resnet-18网络,对各候选区域的细胞核的颜色进行特征提取,得到图像特征。
Resnet-18网络,用于提取候选区域的空间特征。可选的,ResNet-18网络结构由17个卷积层和1个全联接层组成,由堆叠的残差块表示,每个不同卷积结构的残差块之间通过stride为2的下采样联接;ResNet-18网络采用统一的卷积内核大小和池化层内核大小,具体为由尺寸大小为3×3的滤波器构成卷积层,尺寸为2×2的滤波器构成池化层,卷积之后的平均池化操作之后是一个全联接层,ResNet-18网络最后一层的FC层包含1000个通道,模型的最后一部分由SoftMax分类器执行分类。
在一个实施例中,提取各候选区域的轮廓特征包括:对各候选区域进行边缘检测,得到边缘检测结果;将多个候选区域输入至Resnet-18网络,基于Resnet-18网络,对边缘检测结果的形态进行特征提取,得到轮廓特征。
具体地,将候选区域进行灰度化和二值化,利用Canny边缘检测算法提取细胞边缘特征,将细胞边缘特征作为Resnet-18网络的输入,利用Resnet-18网络提取分裂细胞的形态相关的特征,得到轮廓特征。
在一个实施例中,提取各候选区域的结构特征包括:根据Hessian矩阵获取各候选区域的特征值;根据特征值计算与细胞结构相关的形状指数和归一化曲率指数;根据形状指数和归一化曲率指数,得到结构特征。
Hessian矩阵,是一个多元函数的二阶偏导数构成的方阵,描述了函数的局部曲率。
形态指数θ计算公式如下:
Figure DEST_PATH_IMAGE002
归一化曲率指数Φ计算公式如下:
Figure DEST_PATH_IMAGE003
其中,
Figure DEST_PATH_IMAGE005
Figure DEST_PATH_IMAGE007
是Hessian矩阵的两个特征值,I是经过归一化后的RGB颜色向量。
相关技术的分裂细胞检测方法,通常是将图像裁剪成非重叠的图像块,再对其进行核分裂区域和非核分裂区域的分类。这种方法不能确保分裂细胞位于图像块的中心。因此,一个分裂细胞容易被分割到不同的图像块,从而影响分裂细胞的检出率。针对该问题,在一个实施例中,获取细胞染色图像包括:获取拍摄染色的组织学载玻片得到的第一图像;将第一图像裁剪为1024*1024像素的细胞染色图像。如此设置,能够确保分裂细胞位于图像块的中心,避免一个分裂细胞被分割到不同的图像块,从而提升分裂细胞的检出准确度。
在一个实施例中,请参阅图3,给出了一种宫颈鳞状上皮组织非分裂细胞检测方法的流程图,该流程包括如下步骤:
步骤S301,裁剪图像;将宫颈鳞状上皮组织的HE图像裁剪为1024*1024像素的图像;
步骤S302,转换通道;将RGB颜色转换为其对应的光密度值,并用奇异值对其进行分解,利用每个光密度值与奇异值主分解方向之间的角度来确定最佳染色向量;
步骤S303,基于YOLOv5模型检测候选区域;将步骤S302提取的最佳染色向量作为输入,使用YOLOv5模型检测图像中分裂细胞的候选区域
步骤S304,基于Resnet-18网络提取候选区域的图像特征;
步骤S305,基于Resnet-18网络提取候选区域的轮廓特征;利用Canny边缘检测算法提取细胞边缘特征,将细胞边缘特征作为Resnet-18网络的输入,利用Resnet-18网络提取分裂细胞的形态相关的特征,得到轮廓特征;
步骤S306,提取候选区域的结构特征;使用Hessian矩阵获取候选区域的特征值,利用特征值计算与细胞结构相关的形状指数和归一化曲率指数;
步骤S307,将多维特征输入至随机森林进行决策,根据决策结果对候选区域进行分类。
通过上述步骤S301至S307,将图像裁剪为1024*1024像素的细胞染色图像,能够确保分裂细胞位于图像块的中心,避免一个分裂细胞被分割到不同的图像块,使用YOLOv5模型提取分裂细胞的候选区域,以及充分利用了多维度的特征,可以较好地区分相似地分裂细胞和非分裂细胞,检出率高且更为准确。
需要说明的是,在上述流程中或者附图的流程图中示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。例如,步骤S303至步骤S305之间的顺序可以任意互换。
结合上述实施例的分裂细胞检测方法,在本实施例中还提供了一种分裂细胞检测装置,用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。以下所使用的术语“模块”、“单元”、“子单元”等可以实现预定功能的软件和/或硬件的组合。尽管在以下实施例中所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。如图4所示,该装置包括:获取模块,用于获取细胞染色图像,在细胞染色图像中确定分裂细胞的多个候选区域;提取模块,用于提取各候选区域的图像特征、轮廓特征和结构特征;决策模块,用于根据多个候选区域的图像特征、轮廓特征和结构特征,对多个候选区域进行分类以得到分裂细胞所在区域。
上述分裂细胞检测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端,其内部结构图可以如图5所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信,无线方式可通过WIFI、移动蜂窝网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种分裂细胞检测方法的步骤。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图5中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
此外,结合上述实施例中提供的分裂细胞检测方法,在本实施例中还可以提供一种存储介质来实现。该存储介质上存储有计算机程序;该计算机程序被处理器执行时实现上述实施例中的任意一种分裂细胞检测方法的步骤。
需要说明的是,本申请所涉及的用户信息(包括但不限于用户设备信息、用户个人信息等)和数据(包括但不限于用于分析的数据、存储的数据、展示的数据等),均为经用户授权或者经过各方充分授权的信息和数据。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-OnlyMemory,ROM)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、磁变存储器(Magnetoresistive Random Access Memory,MRAM)、铁电存储器(Ferroelectric Random Access Memory,FRAM)、相变存储器(Phase Change Memory,PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器等。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic RandomAccess Memory,DRAM)等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等,不限于此。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等,不限于此。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利要求为准。

Claims (11)

1.一种分裂细胞检测方法,其特征在于,所述方法包括:
获取细胞染色图像,在所述细胞染色图像中确定分裂细胞的多个候选区域;
提取各所述候选区域的图像特征、轮廓特征和结构特征;
根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域。
2.根据权利要求1所述的分裂细胞检测方法,其特征在于,在所述细胞染色图像中确定分裂细胞的多个候选区域包括:
提取所述细胞染色图像中的最佳染色向量,将所述最佳染色向量输入至YOLOv5模型进行检测,得到所述分裂细胞的多个候选区域。
3.根据权利要求2所述的分裂细胞检测方法,其特征在于,提取所述细胞染色图像中的最佳染色向量包括:
将所述细胞染色图像的多种颜色分量转换成对应的光密度值;
根据奇异值对所述多种颜色分量进行分解;
根据各所述光密度值与奇异值主分解方向之间的角度确定所述最佳染色向量。
4.根据权利要求1所述的分裂细胞检测方法,其特征在于,提取各所述候选区域的图像特征包括:
将所述多个候选区域输入至Resnet-18网络,基于所述Resnet-18网络,对各所述候选区域的细胞核的颜色进行特征提取,得到所述图像特征。
5.根据权利要求1所述的分裂细胞检测方法,其特征在于,提取各所述候选区域的轮廓特征包括:
对各所述候选区域进行边缘检测,得到边缘检测结果;
将所述多个候选区域输入至Resnet-18网络,基于所述Resnet-18网络,对所述边缘检测结果的形态进行特征提取,得到所述轮廓特征。
6.根据权利要求1所述的分裂细胞检测方法,其特征在于,提取各所述候选区域的结构特征包括:
根据Hessian矩阵获取各所述候选区域的特征值;
根据所述特征值计算与细胞结构相关的形状指数和归一化曲率指数;
根据所述形状指数和归一化曲率指数,得到所述结构特征。
7.根据权利要求1所述的分裂细胞检测方法,其特征在于,根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域包括:
将各所述候选区域的图像特征、轮廓特征和结构特征,输入至随机森林模型进行决策,根据决策结果确定所述多个候选区域中分裂细胞所在区域。
8.根据权利要求1至权利要求7中任一项所述的分裂细胞检测方法,其特征在于,获取细胞染色图像包括:
获取拍摄染色的组织学载玻片得到的第一图像;
将所述第一图像裁剪为1024*1024像素的所述细胞染色图像。
9.一种分裂细胞检测装置,其特征在于,所述装置包括:
获取模块,用于获取细胞染色图像,在所述细胞染色图像中确定分裂细胞的多个候选区域;
提取模块,用于提取各所述候选区域的图像特征、轮廓特征和结构特征;
决策模块,用于根据所述多个候选区域的图像特征、轮廓特征和结构特征,对所述多个候选区域进行分类以得到分裂细胞所在区域。
10.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至权利要求7中任一项所述分裂细胞检测方法的步骤。
11.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至权利要求7中任一项所述分裂细胞检测方法的步骤。
CN202211279428.5A 2022-10-19 2022-10-19 分裂细胞检测方法、装置和计算机设备 Active CN115359056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211279428.5A CN115359056B (zh) 2022-10-19 2022-10-19 分裂细胞检测方法、装置和计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211279428.5A CN115359056B (zh) 2022-10-19 2022-10-19 分裂细胞检测方法、装置和计算机设备

Publications (2)

Publication Number Publication Date
CN115359056A true CN115359056A (zh) 2022-11-18
CN115359056B CN115359056B (zh) 2023-03-24

Family

ID=84008803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211279428.5A Active CN115359056B (zh) 2022-10-19 2022-10-19 分裂细胞检测方法、装置和计算机设备

Country Status (1)

Country Link
CN (1) CN115359056B (zh)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005076197A2 (en) * 2004-01-31 2005-08-18 Bioimagene, Inc. Method and system for morphology based mitosis identification and classification of digital images
WO2010003041A2 (en) * 2008-07-03 2010-01-07 Nec Laboratories America, Inc. Mitotic figure detector and counter system and method for detecting and counting mitotic figures
GB201212090D0 (en) * 2012-07-06 2012-08-22 Univ Warwick Method and apparatus
CN102737232A (zh) * 2012-06-01 2012-10-17 天津大学 一种分裂细胞识别方法
CN108601509A (zh) * 2015-11-26 2018-09-28 奥林巴斯株式会社 图像处理装置、图像处理方法以及程序
CN110111354A (zh) * 2019-03-28 2019-08-09 青岛大学附属医院 肝脏医学图像动脉区域分割方法及系统
CN110580699A (zh) * 2019-05-15 2019-12-17 徐州医科大学 基于改进Faster RCNN算法的病理图像细胞核检测方法
CN111886630A (zh) * 2018-02-20 2020-11-03 密歇根大学董事会 用于细胞和亚细胞形态学建模和分类的三维细胞和组织图像分析
CN112365440A (zh) * 2020-10-10 2021-02-12 天津大学 基于yolov3的乳腺病理学图像有丝分裂检测装置和方法
CN112884737A (zh) * 2021-02-08 2021-06-01 武汉大学 基于多级迭代的乳腺癌病理图像中有丝分裂自动检测方法
CN113066061A (zh) * 2021-03-24 2021-07-02 同心医联科技(北京)有限公司 一种基于mra的动脉瘤检测方法、系统、终端及介质
CN114937032A (zh) * 2022-06-28 2022-08-23 山东省计算中心(国家超级计算济南中心) 基于染色分离的病理切片颜色相关质量评估方法
CN114973244A (zh) * 2022-06-12 2022-08-30 桂林电子科技大学 一种乳腺癌h&e染色病理图像有丝分裂自动识别系统和方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005076197A2 (en) * 2004-01-31 2005-08-18 Bioimagene, Inc. Method and system for morphology based mitosis identification and classification of digital images
WO2010003041A2 (en) * 2008-07-03 2010-01-07 Nec Laboratories America, Inc. Mitotic figure detector and counter system and method for detecting and counting mitotic figures
CN102737232A (zh) * 2012-06-01 2012-10-17 天津大学 一种分裂细胞识别方法
GB201212090D0 (en) * 2012-07-06 2012-08-22 Univ Warwick Method and apparatus
CN108601509A (zh) * 2015-11-26 2018-09-28 奥林巴斯株式会社 图像处理装置、图像处理方法以及程序
CN111886630A (zh) * 2018-02-20 2020-11-03 密歇根大学董事会 用于细胞和亚细胞形态学建模和分类的三维细胞和组织图像分析
CN110111354A (zh) * 2019-03-28 2019-08-09 青岛大学附属医院 肝脏医学图像动脉区域分割方法及系统
CN110580699A (zh) * 2019-05-15 2019-12-17 徐州医科大学 基于改进Faster RCNN算法的病理图像细胞核检测方法
CN112365440A (zh) * 2020-10-10 2021-02-12 天津大学 基于yolov3的乳腺病理学图像有丝分裂检测装置和方法
CN112884737A (zh) * 2021-02-08 2021-06-01 武汉大学 基于多级迭代的乳腺癌病理图像中有丝分裂自动检测方法
CN113066061A (zh) * 2021-03-24 2021-07-02 同心医联科技(北京)有限公司 一种基于mra的动脉瘤检测方法、系统、终端及介质
CN114973244A (zh) * 2022-06-12 2022-08-30 桂林电子科技大学 一种乳腺癌h&e染色病理图像有丝分裂自动识别系统和方法
CN114937032A (zh) * 2022-06-28 2022-08-23 山东省计算中心(国家超级计算济南中心) 基于染色分离的病理切片颜色相关质量评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
IRSHAD H. ET AL: "Automated mitosis detection in histopathology using morphological and multi-channel statistics features", 《JOURNAL OF PATHOLOGY INFORMATICS》 *
SOMMER C. ET AL: "Learning-based mitotic cell detection in histopathological images", 《ICPR2012》 *
单晓雯: "乳腺病理图像中有丝分裂检测算法研究", 《中国优秀硕士学位论文全文数据库》 *
马明林: "基于深度学习的乳腺病理图像分析技术研究", 《中国优秀硕士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN115359056B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
CN107358260B (zh) 一种基于表面波cnn的多光谱图像分类方法
US6404920B1 (en) System for generalizing objects and features in an image
Aytekın et al. Unsupervised building detection in complex urban environments from multispectral satellite imagery
WO2022000862A1 (zh) 鱼眼图像中的对象检测方法、装置及存储介质
CN110781756A (zh) 基于遥感图像的城市道路提取方法及装置
WO2020062360A1 (zh) 一种图像融合分类的方法及装置
CN106295613A (zh) 一种无人机目标定位方法及系统
CN102737243A (zh) 获取多幅图像的描述信息的方法及装置与图像匹配方法
GB2534903A (en) Method and apparatus for processing signal data
CN114241478A (zh) 宫颈细胞图像中异常细胞图像识别方法、装置
CN113159300A (zh) 图像检测神经网络模型及其训练方法、图像检测方法
CN114155365A (zh) 模型训练方法、图像处理方法及相关装置
CN113901900A (zh) 一种同源或异源遥感图像用无监督变化检测方法及系统
CN111192239A (zh) 遥感影像变化区域检测方法、装置、存储介质及电子设备
Fan et al. Region-merging method with texture pattern attention for SAR image segmentation
CN113870196A (zh) 一种基于锚点切图的图像处理方法、装置、设备和介质
CN111382638A (zh) 一种图像检测方法、装置、设备和存储介质
CN115359056B (zh) 分裂细胞检测方法、装置和计算机设备
CN113284153A (zh) 卫星云层图像处理方法、装置、计算机设备和存储介质
CN115880362B (zh) 码区定位方法、装置、计算机设备及计算机可读存储介质
CN116310899A (zh) 基于YOLOv5改进的目标检测方法及装置、训练方法
CN115035313B (zh) 黑颈鹤识别方法、装置、设备及存储介质
CN115719416A (zh) 刀闸状态识别方法、装置、计算机设备和存储介质
CN113284066B (zh) 遥感影像自动云检测方法和装置
CN110633627A (zh) 在视频中定位对象的方法、装置、计算机设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant