CN107977682A - 基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 - Google Patents
基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 Download PDFInfo
- Publication number
- CN107977682A CN107977682A CN201711376577.2A CN201711376577A CN107977682A CN 107977682 A CN107977682 A CN 107977682A CN 201711376577 A CN201711376577 A CN 201711376577A CN 107977682 A CN107977682 A CN 107977682A
- Authority
- CN
- China
- Prior art keywords
- image
- cell
- polar coordinate
- data
- lymphocyte
- 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
Links
- 210000004027 cell Anatomy 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000002708 enhancing effect Effects 0.000 title claims abstract 3
- 210000002751 lymph Anatomy 0.000 title abstract description 5
- 210000004698 lymphocyte Anatomy 0.000 claims abstract description 71
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 230000000877 morphologic effect Effects 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims description 22
- 238000012549 training Methods 0.000 claims description 21
- 210000003743 erythrocyte Anatomy 0.000 claims description 12
- 238000013528 artificial neural network Methods 0.000 claims description 8
- 210000000265 leukocyte Anatomy 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 210000000601 blood cell Anatomy 0.000 claims description 5
- 210000003855 cell nucleus Anatomy 0.000 claims description 5
- 210000000805 cytoplasm Anatomy 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims 1
- 238000004891 communication Methods 0.000 claims 1
- 238000003672 processing method Methods 0.000 claims 1
- 210000004369 blood Anatomy 0.000 description 7
- 239000008280 blood Substances 0.000 description 7
- 238000010186 staining Methods 0.000 description 5
- 210000001772 blood platelet Anatomy 0.000 description 4
- 238000013135 deep learning Methods 0.000 description 3
- 238000000513 principal component analysis Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 210000004940 nucleus Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于极坐标变换数据增强的淋巴类细胞分类方法及其装置,属于计算机视觉领域。分类方法包括:S1获取细胞图像;S2计算图像的R通道与G通道的比值矩阵;S3对比值矩阵进行直方图统计并二值化;S4对二值图像做形态学处理;S5提取细胞的边缘像素点并分割细胞;S6选取细胞边缘像素点为极点,建立极坐标系,把图片用极坐标变换映射到直角坐标系中;S7遍历边缘像素点,得到变换后的细胞图像;S8用ResNet训练细胞图像,并测试分类结果。本发明提出了一种新的数据增强方法,可以用极少的标记数据进行分类,而且有效地分开了难以识别的大颗粒淋巴细胞与异型淋巴细胞。
Description
技术领域
本发明涉及计算机视觉领域,特别涉及一种基于极坐标变换数据增强的淋巴类细胞分类方法及其装置。
背景技术
近年来,随着深度学习在各个领域的迅猛发展,在医学病理检测中,如何用深度学习的方法检测并分类各种细胞成为一个广泛关注的热点问题。基于显微图像的细胞检测和分类过程如下:首先获取细胞的显微图像,再运用一系列算法,将细胞与背景分割开,放入网络中训练。在训练过程中,为了出现过拟合的现象,需要输入充足的数据量,即要做数据增强。
目前,深度学习中常用的数据增强方法大致有四种:
第一种方法是对颜色的数据增强,包括调整色彩的饱和度、亮度和对比度等方面。具体操作方法为,在图像的HSV颜色空间,改变饱和度S和亮度V分量,保持色调H不变。对每个像素的S和V分量进行指数运算(指数因子在0.25到4之间),增加光照变化。
第二类是PCA Jittering,最早是由Alex在他2012年赢得ImageNet竞赛的那篇NIPS中提出来的。首先按照RGB三个颜色通道计算均值和标准差,对网络的输入数据做normalization,随后在整个训练集上计算了协方差矩阵,进行特征分解,得到特征向量和特征值,用来做PCA Jittering。
第三类是Random Crop,该方法采用随机图像差值方式,对图像进行裁剪、缩放。包括Scale Jittering方法(VGG及ResNet模型使用)或者尺度和长宽比增强变换。缩放(zoom)就是按照一定的比例放大或缩小图像。翻转(flip)是沿着水平或者垂直方向翻转图像。旋转(rotation)是随机旋转图像一定角度。
第四类方法是噪声扰动。对图像的每个像素点RGB进行随机扰动,常用的噪声模式是椒盐噪声和高斯噪声。
在以上这些算法中,最常用的数据增强算法是对图像进行缩放、翻转和旋转。但是该方法至少具有以下缺点:在处理医学细胞图像时,普通的旋转、翻转和缩放并不能提高网络的泛化能力,而裁剪和特征分解等方法,在没有人工干预的条件下,很难准确找到细胞之间的差异部分,从而有效判断细胞分类。另外,由于人体血液中各类细胞的含量不同,已有的大颗粒淋巴细胞数据量不到异型淋巴数据量的十分之一,在极大的数据量差别下,上述的几种方法都不能准确分类。
发明内容
针对上述现有技术中存在的缺陷,为了增加淋巴类细胞分类的准确度,本发明的目的是提供一种基于极坐标变换的数据增强方法,能准确高效地针对细胞间的差异做数据增强,并且能够处理数据量相差非常大的情况。本发明的另外一个目的是提供实现该方法的装置。
为了实现上述发明目的,本发明方法采用的技术方案如下:
基于极坐标变换数据增强的淋巴类细胞分类方法,包括如下步骤:
S1,获取染色血细胞的原始图像;
S2,分离原始图像的RGB通道,计算R通道与G通道的比值矩阵;
S3,根据R通道与G通道的比值绘制比值矩阵的直方图,根据阈值范围将原始图像二值化,得到淋巴细胞和异型淋巴细胞的二值图像;
S4,将二值图像进行一系列的形态学处理;
S5,用连通域的方法提取经形态学处理后的细胞图像的边缘像素点,找到细胞上下左右的边缘像素点,然后分割细胞;
S6,取细胞边缘的任一像素点为极点,建立极坐标系,将所有像素点用极坐标变换一一映射到直角坐标系中;
S7,遍历细胞的边缘像素点,每个像素点产生一张变换后的图像;
S8,将步骤S7得到的变换后的图像作为神经网络的输入,训练网络,测试分类结果,分出淋巴细胞和异型淋巴细胞。
本发明装置采用的技术方案如下:
基于极坐标变换数据增强的淋巴类细胞分类装置,包括:显微图像获取模块,用于获取淋巴类细胞的显微图像;图像预处理模块,用于选取血细胞显微图像中的淋巴类细胞;细胞分割模块,用于把淋巴类细胞从背景中分割出来,去除周围的红细胞,保留细胞浆和细胞核;细胞边缘提取模块,用于根据淋巴细胞和异型淋巴细胞的特征,提取淋巴类细胞的边缘像素点信息;图像数据增强模块,用于增大训练样本数,放大细胞的边缘特征,平衡已标记细胞的数据量;神经网络训练模块,用于学习已标记样本的特征;淋巴类细胞分类模块,用于对新的淋巴类细胞数据做分类。
本发明提出了一种新的数据增强方法,通过极坐标变换做数据增强,将淋巴类细胞的边缘像素点投射到极坐标系上,再用神经网络训练后分类。该方法不仅可以用极少的标记数据进行分类,而且有效地分开了难以识别的大颗粒淋巴细胞与异型淋巴细胞,把原本需要用非线性分类器分类的细胞,简单地用线性分类器就能分出良好的效果,减少了过拟合现象,在各类数据量相差极大的情况下,大大增加了细胞分类的准确率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性的劳动前提下,还可以根据这些附图获得其他的附图。
图1是根据极坐标变换推导出对应空间点坐标的示意图。
图2是极坐标变换实现数据增强的具体算法流程表。
图3是进行基于极坐标数据增强前后细胞图片对比图,左图为数据增强前,右图为数据增强后。
图4是本发明极坐标变换数据增强方法的流程图。
图5是本发明淋巴细胞分类装置的结构示意图。
具体实施方式
为使本发明的目的,技术方案和优点更加清楚,下面将结合附图对本发明实施方法作进一步地详细描述。
实施例1
本实施例提供了一种基于极坐标变换数据增强的淋巴细胞检测与分类方法,参见图4,具体包括:
S1:对原始图像分离通道,计算R通道与G通道数值比矩阵。
将样本血液染色后制作成血涂片。在显微镜下观看,血涂片中包含白细胞、红细胞和血小板等。经过染色,白细胞核与血小板呈紫色,白细胞浆呈淡紫色,红细胞呈淡粉色。用DP27相机拍摄血细胞图片传到计算机中,获取原始图像。
每张原始图像中有若干个白细胞和无数红细胞。为了提取白细胞,分离原始图像的RGB三通道,并绘制三个通道的直方图。计算R通道与G通道的比值。
S2:绘制矩阵的直方图,根据直方图信息确定阈值范围,得到二值化的图像。
其中,淋巴类细胞属于白细胞,为前景。其余细胞为背景。根据直方图,找到分离背景与前景的临界点为1.19。用for循环遍历矩阵中的每一个像素,访问像素值,大于1.19的区域为细胞浆和细胞核,将像素值设为255;小于1.19的为背景区域,即红细胞和其他背景,将像素值设为0。由此得到淋巴类细胞的二值图像。
S3:对上述图像进行一系列形态学运算,提取细胞边缘像素坐标信息。
具体地,先对淋巴类二值图像进行用imfill函数进行填充运算,填补漏洞。然后用10×10的模板对图像进行闭运算处理。在血涂片的制作和染色过程中,由于人工操作的失误,会将一些细胞弄破,或是染色不均衡,这对上述步骤S2中的二值化结果产生影响。用闭运算处理可以将一些染色不当的背景区域滤除。最后,因为血小板和细胞核的染色结果一样,去除连通区域像素点小于1000的区域就可以去除血小板,得到淋巴类细胞。
S4:对上述图像进行分割,得到带有边缘红细胞的细胞图像,图像大小为n*n。
具体地,淋巴类细胞分为大颗粒淋巴细胞、淋巴细胞、异型淋巴细胞。在这三种细胞中,大颗粒淋巴细胞和异型淋巴细胞最难区分,目前的血液分析仪等检测器不能有效分类这两种淋巴类细胞,造成很大误诊率,需要医生的人工判断。而上述两种细胞的差别在于细胞浆与红细胞的接触边缘有无变形和颜色加深。因此,单纯分割出淋巴类细胞并不能进行有效分类。本方法保留淋巴类细胞和红细胞接触的边缘区域。
用连通域的方法提取处理后的细胞图像的边缘像素点,找到细胞上下左右的边缘像素点,分割细胞,这里采用8邻域连通。
首先从上往下、从左往右对图像进行扫描,找到连通区域的第一个目标段,标记该段并且压入堆栈,作为“区域增长”的种子段。检查当前段的上下两行是否有重叠且未标记的目标段。如果不存在重叠的目标段就把当前段弹出堆栈;如果存在重叠且未标记的目标段,则标记该段并且压入堆栈,作为新的“种子段”。后续操作不断从堆栈中取出种子段,重复上述操作直到堆栈为空(标记完一个连通区域)。接着搜索图像中下一个未标记的连通区域,重复上述操作直到图像中所有的连通区域标记完毕。
本行目标段(Xs1,Xe1)和上下两行目标段(Xs2,Xe2)重叠的准则是:
Xs1-1≤Xe2并且Xe1+1≥Xs2
其中,Xs1表示本行目标段的起始坐标,Xe1表示本行目标段的结束坐标。Xs2表示上下两行目标段的起始坐标,Xe2表示上下两行目标段的结束坐标。
设连通区域总数为num,则对每一个区域inum,找出行的最小值与最大值:x1,x2,找出列的最小值与最大值:y1,y2。考虑到神经网络的输入,将图片大小统一设置成224×224。这样既能保留边缘的红细胞信息,又能利于网络的计算。
S5:取细胞边缘任一坐标点为极点,建立极坐标系,将其余像素点通过极坐标变换得到直角坐标,得到新图像。
在一张M×N大小的图像I中,取直角坐标系中细胞边缘的任一像素点(m,n)为极点,建立极坐标系。用M条射线将2π划分成M个角度,单位角度为2π/M。比如,在一张M×N=224×224的图像I中,选取(90,100)这个点为极点,则图像I以(90,100)为中心,被224条射线划分成224个单位角度。那么在极坐标射线上的每个像素都可以用直角坐标表示为(90+x,100+y),也可以用极坐标表示为(θu,v)。其中,
x=v·cos(θu)
y=v·sin(θu)
设变换后的图像为如果满足条件0≤m+x<M且0≤n+y<V,即确保点(m+x,n+y)落在图像I内,则
即可得到变换后的新图像作为数据增强的结果。
S6:顺着边缘移动像素点,以新像素点为极点再次将图像映射到直角坐标系中。循环此步骤,直到遍历整个细胞边缘。
S7:将所有变换后的图片作为ResNet的输入,经过网络的训练及测试,得到分类的淋巴细胞与异型淋巴细胞。
具体地,本方法采用的神经网络为深度残差网络,即ResNet。这个网络的优点是,随着层数的增加,网络的精度不断上升,不会出现梯度爆炸和梯度弥散的现象。采用残差作为网络的输出,加快了收敛速度,对细胞间微小的差异更加敏感。
将所有原始细胞数据(大颗粒淋巴细胞、淋巴细胞和异型淋巴细胞)随机分为两类,一类为训练数据,占总数据的2/3,一类为测试数据,占总数据的1/3。其中大颗粒淋巴数据量为异型淋巴数据量的十分之一。将训练数据做极坐标数据增强,得到所有的训练数据。把所有训练数据输入ResNet中,经过网络的学习,得到checkpoint文件,该文件存储了网络模型的各类参数。用测试图片进行测试,得到分类准确率。
实施例2
参见图5,本发明实施例提供了一种基于极坐标变换数据增强的淋巴类细胞检测与分类的装置,包括:
显微图像获取模块201,用于获取淋巴类细胞的原始图像;
显微图像采集工具可以包括普通相机、扫片机等可以拍摄图像的工具。在医学检验科,最好的细胞拍摄仪器为扫片机,该仪器带有电动载物台,可以自动锁定白细胞并聚焦拍摄。拍摄效果清晰。其次可用显微镜和相机组合拍摄,该方法需要人工对焦和移动载玻片。为节省成本,本装置采用显微镜和相机的组合获取原始图像。
图像预处理模块202,对图像做一系列形态学运算,分离背景和前景;
具体地,以400倍显微镜拍摄的血涂片为例,一张图片有数个白细胞和无数红细胞。对图像做二值化、填充、开运算、去除小块等运算,定位淋巴类细胞。
细胞分割模块203,用于从原图中分割出淋巴类细胞,将无用的背景和其他细胞剔除,并保留对分类结果有影响的边缘红细胞;
细胞边缘提取模块204,用于保留淋巴类细胞的细胞浆与红细胞接触的边缘像素点,这些像素点的值对大颗粒淋巴细胞和异型淋巴细胞的分类至关重要;
图像数据增强模块205,将所有图像对应的深度图转化为数据矩阵,具体方法为:
遍历保留的每一个边缘像素点,以该像素点为极点,做极坐标变换,将图像投射到极坐标系中,由于极坐标变换的特点,靠近极点的区域过采样,远离极点的区域欠采样。这样放大了细胞边缘的特征。
神经网络训练模块206,用于学习细胞分类的各种参数。
具体地,本实例采用深度残差网络ResNet进行学习。将输入图片向量化,每张图片默认大小为224×224,输入50层的ResNet中,迭代次数为1000次,初始学习率为0.001。网络末端用softmax分类器连接。训练参数保存在checkpoint文件中。
淋巴类细胞分类模块207,将未标记的淋巴类细胞输入网络中,网络用test对淋巴类细胞分类。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改,等同替换,改进等,均应包含在本发明的保护范围之中。
Claims (8)
1.基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,包括如下步骤:
S1,获取染色血细胞的原始图像;
S2,分离原始图像的RGB通道,计算R通道与G通道的比值矩阵;
S3,根据R通道与G通道的比值绘制比值矩阵的直方图,根据阈值范围将原始图像二值化,得到淋巴细胞和异型淋巴细胞的二值图像;
S4,将二值图像进行一系列的形态学处理;
S5,用连通域的方法提取经形态学处理后的细胞图像的边缘像素点,找到细胞上下左右的边缘像素点,然后分割细胞;
S6,取细胞边缘的任一像素点为极点,建立极坐标系,将所有像素点用极坐标变换一一映射到直角坐标系中;
S7,遍历细胞的边缘像素点,每个像素点产生一张变换后的图像;
S8,将步骤S7得到的变换后的图像作为神经网络的输入,训练网络,测试分类结果,分出淋巴细胞和异型淋巴细胞。
2.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S3中,阈值范围的确定方法具体为:淋巴类细胞属于白细胞为前景,其余细胞为背景,找到分离背景与前景的临界点为1.19,大于1.19的区域为细胞浆和细胞核,将像素值设为255;小于1.19的区域为背景区域,将像素值设为0。
3.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S4中,形态学处理的方法为:先对二值图像进行填充运算,填补漏洞;然后用10×10的模板进行闭运算处理,将一些颜色较深的背景细胞滤除;最后去除连通区域像素点小于1000的区域。
4.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S5中,连通域的方法采用8邻域连通,具体为:
(1)首先从上往下、从左往右对图像进行扫描,找到连通区域的第一个目标段,标记该段并且压入堆栈,作为“区域增长”的种子段;
(2)检查当前段的上下两行是否有重叠且未标记的目标段;如果不存在重叠的目标段,则把当前段弹出堆栈;如果存在重叠且未标记的目标段,则标记该段并且压入堆栈,作为新的“种子段”;本行目标段(Xs1,Xe1)和上下两行目标段(Xs2,Xe2)重叠的准则是:
Xs1-1≤Xe2并且Xe1+1≥Xs2;
其中,Xs1表示本行目标段的起始坐标,Xe1表示本行目标段的结束坐标。Xs2表示上下两行目标段的起始坐标,Xe2表示上下两行目标段的结束坐标;
(3)后续操作不断从堆栈中取出种子段,重复步骤(2)直到堆栈为空,即标记完一个连通区域;
(4)接着搜索图像中下一个未标记的连通区域,重复上述步骤(1)-(3)直到图像中所有的连通区域标记完毕。
5.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S6中,具体实现方法为:在一张M×N大小的图像I中,取直角坐标系中细胞边缘的任一像素点(m,n)为极点,建立极坐标系;用M条射线将2π划分成M个角度,单位角度为2π/M;那么在极坐标射线上的每个像素都可以用直角坐标表示为(m+x,n+y),也可以用极坐标表示为(θu,v),其中,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>&CenterDot;</mo>
<mi>u</mi>
<mo>/</mo>
<mi>M</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>u</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>0</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
x=v·cos(θu)
y=v·sin(θu)
设变换后的图像为如果满足条件0≤m+x<M且0≤n+y<V,即确保点(m+x,n+y)落在图像I内,则
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>x</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
即可得到变换后的图像作为数据增强的结果。
6.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S7中,具体实现方法为:存储每个细胞边缘像素点的坐标,以这些坐标为极点,建立极坐标系,重复步骤S6,直到所有边缘像素都被变换;如果有n个像素点,就会产生n张变换的图像。
7.根据权利要求1所述的基于极坐标变换数据增强的淋巴类细胞分类方法,其特征在于,所述步骤S8中,具体实现方法为:将所有原始细胞数据,包括淋巴细胞和异型淋巴细胞,随机分为两类:一类为训练数据,占总数据的2/3;一类为测试数据,占总数据的1/3;将训练数据做极坐标数据增强,得到所有的训练数据;然后把所有训练数据输入深度残差网络中,经过网络的学习,得到checkpoint文件,该文件存储了网络模型的各类参数;用测试图片进行测试,得到分类准确率。
8.基于极坐标变换数据增强的淋巴类细胞分类装置,其特征在于,所述装置包括:
显微图像获取模块,用于获取淋巴类细胞的显微图像;
图像预处理模块,用于选取血细胞显微图像中的淋巴类细胞;
细胞分割模块,用于把淋巴类细胞从背景中分割出来,去除周围的红细胞,保留细胞浆和细胞核;
细胞边缘提取模块,用于根据淋巴细胞和异型淋巴细胞的特征,提取淋巴类细胞的边缘像素点信息;
图像数据增强模块,用于增大训练样本数,放大细胞的边缘特征,平衡已标记细胞的数据量;
神经网络训练模块,用于学习已标记样本的特征;
淋巴类细胞分类模块,用于对新的淋巴类细胞数据做分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711376577.2A CN107977682B (zh) | 2017-12-19 | 2017-12-19 | 基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711376577.2A CN107977682B (zh) | 2017-12-19 | 2017-12-19 | 基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107977682A true CN107977682A (zh) | 2018-05-01 |
CN107977682B CN107977682B (zh) | 2021-09-28 |
Family
ID=62006760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711376577.2A Active CN107977682B (zh) | 2017-12-19 | 2017-12-19 | 基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107977682B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109034208A (zh) * | 2018-07-03 | 2018-12-18 | 怀光智能科技(武汉)有限公司 | 一种高低分辨率组合的宫颈细胞病理切片分类方法 |
CN109087283A (zh) * | 2018-07-03 | 2018-12-25 | 怀光智能科技(武汉)有限公司 | 基于细胞团的宫颈细胞病理切片病变细胞识别方法及系统 |
CN109145941A (zh) * | 2018-07-03 | 2019-01-04 | 怀光智能科技(武汉)有限公司 | 一种非规则宫颈细胞团图像分类方法及系统 |
CN110705403A (zh) * | 2019-09-19 | 2020-01-17 | 平安科技(深圳)有限公司 | 细胞分类方法、装置、介质及电子设备 |
CN110827259A (zh) * | 2019-11-04 | 2020-02-21 | 上海芯翌智能科技有限公司 | 一种用于目标关键点检测的样本增广的方法与设备 |
CN111127434A (zh) * | 2019-12-25 | 2020-05-08 | 志诺维思(北京)基因科技有限公司 | 染色分离方法及装置 |
CN111458269A (zh) * | 2020-05-07 | 2020-07-28 | 厦门汉舒捷医疗科技有限公司 | 一种外周血淋巴微核细胞图像人工智能识别方法 |
CN112132253A (zh) * | 2019-06-24 | 2020-12-25 | 北京眼神智能科技有限公司 | 3d动作识别方法、装置、计算机可读存储介质及设备 |
CN112184696A (zh) * | 2020-10-14 | 2021-01-05 | 中国科学院近代物理研究所 | 一种细胞核和细胞器计数及其面积计算方法与系统 |
CN112330701A (zh) * | 2020-11-26 | 2021-02-05 | 山东师范大学 | 基于极坐标表示的组织病理图像细胞核分割方法及系统 |
CN112446427A (zh) * | 2020-11-27 | 2021-03-05 | 王伟佳 | 一种髓系血细胞识别方法、装置、存储介质和电子设备 |
CN114300099A (zh) * | 2021-11-24 | 2022-04-08 | 大连工业大学 | 一种基于YOLOv5和显微高光谱图像的异型淋巴细胞分型方法 |
CN114387196A (zh) * | 2021-12-24 | 2022-04-22 | 汉姆德(宁波)智能医疗科技有限公司 | 超分辨显微镜欠采样图像的生成方法、装置 |
CN115050024A (zh) * | 2022-06-21 | 2022-09-13 | 广州医科大学 | 一种可解释性的粒细胞智能实时识别方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101923648A (zh) * | 2009-06-15 | 2010-12-22 | 深圳迈瑞生物医疗电子股份有限公司 | 支持向量机的聚类方法与装置 |
US20110243417A1 (en) * | 2008-09-03 | 2011-10-06 | Rutgers, The State University Of New Jersey | System and method for accurate and rapid identification of diseased regions on biological images with applications to disease diagnosis and prognosis |
CN103473739A (zh) * | 2013-08-15 | 2013-12-25 | 华中科技大学 | 一种基于支持向量机的白细胞图像精确分割方法与系统 |
CN103745210A (zh) * | 2014-01-28 | 2014-04-23 | 爱威科技股份有限公司 | 一种白细胞分类方法及装置 |
CN104732524A (zh) * | 2015-02-07 | 2015-06-24 | 嘉善加斯戴克医疗器械有限公司 | 血液白细胞显微图像的随机权网络分割方法 |
CN104834914A (zh) * | 2015-05-15 | 2015-08-12 | 广西师范大学 | 宫颈细胞图像特征识别方法及宫颈细胞特征识别装置 |
CN106248559A (zh) * | 2016-07-14 | 2016-12-21 | 中国计量大学 | 一种基于深度学习的白细胞五分类方法 |
CN106940889A (zh) * | 2017-03-30 | 2017-07-11 | 福建师范大学 | 基于像素邻域特征聚类的淋巴结he染色病理图像分割方法 |
-
2017
- 2017-12-19 CN CN201711376577.2A patent/CN107977682B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110243417A1 (en) * | 2008-09-03 | 2011-10-06 | Rutgers, The State University Of New Jersey | System and method for accurate and rapid identification of diseased regions on biological images with applications to disease diagnosis and prognosis |
CN101923648A (zh) * | 2009-06-15 | 2010-12-22 | 深圳迈瑞生物医疗电子股份有限公司 | 支持向量机的聚类方法与装置 |
CN103473739A (zh) * | 2013-08-15 | 2013-12-25 | 华中科技大学 | 一种基于支持向量机的白细胞图像精确分割方法与系统 |
CN103745210A (zh) * | 2014-01-28 | 2014-04-23 | 爱威科技股份有限公司 | 一种白细胞分类方法及装置 |
CN104732524A (zh) * | 2015-02-07 | 2015-06-24 | 嘉善加斯戴克医疗器械有限公司 | 血液白细胞显微图像的随机权网络分割方法 |
CN104834914A (zh) * | 2015-05-15 | 2015-08-12 | 广西师范大学 | 宫颈细胞图像特征识别方法及宫颈细胞特征识别装置 |
CN106248559A (zh) * | 2016-07-14 | 2016-12-21 | 中国计量大学 | 一种基于深度学习的白细胞五分类方法 |
CN106940889A (zh) * | 2017-03-30 | 2017-07-11 | 福建师范大学 | 基于像素邻域特征聚类的淋巴结he染色病理图像分割方法 |
Non-Patent Citations (7)
Title |
---|
KUSE, M.,ET.AL: "A classification scheme for lymphocyte segmentation in H&E stained histology images", 《IN INTERNATIONAL CONFERENCE ON PATTERN RECOGNITION》 * |
SAJI, H.,ET.AL: "A proposal for combination of total number and anatomical location of involved lymph nodes for nodal classification in non-small cell lung cancer", 《CHEST》 * |
YU, W.,ET.AL: "Automatic classification of leukocytes using deep neural network", 《 IN 2017 IEEE 12TH INTERNATIONAL CONFERENCE ON ASIC (ASICON)》 * |
刘晶: "自动化血细胞形态学分析及分类关键技术研究", 《中国优秀硕士学位论文全文数据库》 * |
周颖颖: "彩色白细胞图像的特征选择与分类识别", 《中国优秀硕士学位论文全文数据库》 * |
尹晓红: "骨髓淋巴细胞自动识别技术的研究", 《中国优秀硕士学位论文全文数据库》 * |
李一安: "全血细胞成像五分类算法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109087283B (zh) * | 2018-07-03 | 2021-03-09 | 怀光智能科技(武汉)有限公司 | 基于细胞团的宫颈细胞病理切片病变细胞识别方法及系统 |
CN109087283A (zh) * | 2018-07-03 | 2018-12-25 | 怀光智能科技(武汉)有限公司 | 基于细胞团的宫颈细胞病理切片病变细胞识别方法及系统 |
CN109145941A (zh) * | 2018-07-03 | 2019-01-04 | 怀光智能科技(武汉)有限公司 | 一种非规则宫颈细胞团图像分类方法及系统 |
CN109034208B (zh) * | 2018-07-03 | 2020-10-23 | 怀光智能科技(武汉)有限公司 | 一种高低分辨率组合的宫颈细胞切片图像分类系统 |
CN109145941B (zh) * | 2018-07-03 | 2021-03-09 | 怀光智能科技(武汉)有限公司 | 一种非规则宫颈细胞团图像分类方法及系统 |
CN109034208A (zh) * | 2018-07-03 | 2018-12-18 | 怀光智能科技(武汉)有限公司 | 一种高低分辨率组合的宫颈细胞病理切片分类方法 |
CN112132253A (zh) * | 2019-06-24 | 2020-12-25 | 北京眼神智能科技有限公司 | 3d动作识别方法、装置、计算机可读存储介质及设备 |
CN112132253B (zh) * | 2019-06-24 | 2024-06-11 | 北京眼神智能科技有限公司 | 3d动作识别方法、装置、计算机可读存储介质及设备 |
CN110705403A (zh) * | 2019-09-19 | 2020-01-17 | 平安科技(深圳)有限公司 | 细胞分类方法、装置、介质及电子设备 |
CN110827259A (zh) * | 2019-11-04 | 2020-02-21 | 上海芯翌智能科技有限公司 | 一种用于目标关键点检测的样本增广的方法与设备 |
CN111127434A (zh) * | 2019-12-25 | 2020-05-08 | 志诺维思(北京)基因科技有限公司 | 染色分离方法及装置 |
CN111127434B (zh) * | 2019-12-25 | 2024-03-08 | 志诺维思(北京)基因科技有限公司 | 染色分离方法及装置 |
CN111458269A (zh) * | 2020-05-07 | 2020-07-28 | 厦门汉舒捷医疗科技有限公司 | 一种外周血淋巴微核细胞图像人工智能识别方法 |
CN112184696A (zh) * | 2020-10-14 | 2021-01-05 | 中国科学院近代物理研究所 | 一种细胞核和细胞器计数及其面积计算方法与系统 |
CN112184696B (zh) * | 2020-10-14 | 2023-12-29 | 中国科学院近代物理研究所 | 一种细胞核和细胞器计数及其面积计算方法与系统 |
CN112330701A (zh) * | 2020-11-26 | 2021-02-05 | 山东师范大学 | 基于极坐标表示的组织病理图像细胞核分割方法及系统 |
CN112446427A (zh) * | 2020-11-27 | 2021-03-05 | 王伟佳 | 一种髓系血细胞识别方法、装置、存储介质和电子设备 |
CN114300099A (zh) * | 2021-11-24 | 2022-04-08 | 大连工业大学 | 一种基于YOLOv5和显微高光谱图像的异型淋巴细胞分型方法 |
CN114387196A (zh) * | 2021-12-24 | 2022-04-22 | 汉姆德(宁波)智能医疗科技有限公司 | 超分辨显微镜欠采样图像的生成方法、装置 |
CN114387196B (zh) * | 2021-12-24 | 2024-08-27 | 汉姆德(宁波)智能医疗科技有限公司 | 超分辨显微镜欠采样图像的生成方法、装置 |
CN115050024A (zh) * | 2022-06-21 | 2022-09-13 | 广州医科大学 | 一种可解释性的粒细胞智能实时识别方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107977682B (zh) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107977682B (zh) | 基于极坐标变换数据增强的淋巴类细胞分类方法及其装置 | |
CN108021903B (zh) | 基于神经网络的人工标注白细胞的误差校准方法及装置 | |
US10565479B1 (en) | Identifying and excluding blurred areas of images of stained tissue to improve cancer scoring | |
Li et al. | Segmentation of white blood cell from acute lymphoblastic leukemia images using dual‐threshold method | |
CN106815604B (zh) | 基于多层信息融合的注视点检测方法 | |
CN109978822B (zh) | 一种基于机器视觉的香蕉成熟度评判建模方法和评判方法 | |
Agaian et al. | Automated screening system for acute myelogenous leukemia detection in blood microscopic images | |
Ramesh et al. | Isolation and two-step classification of normal white blood cells in peripheral blood smears | |
JP5394485B2 (ja) | 印環細胞検出器及び関連する方法 | |
US9239281B2 (en) | Method and device for dividing area of image of particle in urine | |
Salem et al. | A comparative study of white blood cells segmentation using otsu threshold and watershed transformation | |
CN106610969A (zh) | 基于多模态信息的视频内容审查系统及方法 | |
CN104392240A (zh) | 一种基于多特征融合的寄生虫虫卵识别方法 | |
Poomcokrak et al. | Red blood cells extraction and counting | |
CN104636754B (zh) | 基于舌体分区颜色特征的智能图像分类方法 | |
Soltanzadeh et al. | Extraction of nucleolus candidate zone in white blood cells of peripheral blood smear images using curvelet transform | |
Galsgaard et al. | Circular hough transform and local circularity measure for weight estimation of a graph-cut based wood stack measurement | |
US11790673B2 (en) | Method for detection of cells in a cytological sample having at least one anomaly | |
Yingthawornsuk et al. | Automatic Thai Coin Calculation System by Using SIFT | |
CN108073940A (zh) | 一种非结构化环境中的3d目标实例物体检测的方法 | |
Nguyen et al. | Using contextual information to classify nuclei in histology images | |
CN109741351A (zh) | 一种基于深度学习的类别敏感型边缘检测方法 | |
JP3653804B2 (ja) | 粒子画像の領域分割方法及び装置 | |
CN109615630A (zh) | 基于图像处理技术的半连续铸造铝硅合金微观组织分析方法 | |
WO2014181024A1 (es) | Método implementado por ordenador para reconocimiento y clasificación de células sanguíneas anormales y programas informáticos para llevar a cabo el método |
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 |