CN108172278A - 一种he染色病理图像颜色归一化方法 - Google Patents
一种he染色病理图像颜色归一化方法 Download PDFInfo
- Publication number
- CN108172278A CN108172278A CN201810067454.9A CN201810067454A CN108172278A CN 108172278 A CN108172278 A CN 108172278A CN 201810067454 A CN201810067454 A CN 201810067454A CN 108172278 A CN108172278 A CN 108172278A
- Authority
- CN
- China
- Prior art keywords
- matrix
- image
- pathological
- pathological image
- staining
- 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
- 230000001575 pathological effect Effects 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000000975 dye Substances 0.000 title description 2
- 239000011159 matrix material Substances 0.000 claims abstract description 59
- 238000010606 normalization Methods 0.000 claims abstract description 41
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000002835 absorbance Methods 0.000 claims description 25
- 238000010186 staining Methods 0.000 claims description 23
- 238000005457 optimization Methods 0.000 claims description 22
- 230000007170 pathology Effects 0.000 claims description 14
- 238000007490 hematoxylin and eosin (H&E) staining Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 6
- 239000003795 chemical substances by application Substances 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 4
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 3
- 238000004043 dyeing Methods 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 23
- 230000000694 effects Effects 0.000 description 9
- 210000003855 cell nucleus Anatomy 0.000 description 8
- 239000003086 colorant Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 210000004940 nucleus Anatomy 0.000 description 3
- BVRIUXYMUSKBHG-UHFFFAOYSA-N 4-[[4-(dimethylamino)phenyl]diazenyl]aniline Chemical compound C1=CC(N(C)C)=CC=C1N=NC1=CC=C(N)C=C1 BVRIUXYMUSKBHG-UHFFFAOYSA-N 0.000 description 2
- 238000004195 computer-aided diagnosis Methods 0.000 description 2
- 210000000805 cytoplasm Anatomy 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 208000030776 invasive breast carcinoma Diseases 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- 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/10024—Color 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/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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30024—Cell structures in vitro; Tissue sections in vitro
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种HE染色病理图像颜色归一化方法，先对病理切片图像进行通道分解，得到病理图像，通过建立优化模型，计算最优反卷积矩阵，并利用最优反卷积矩阵对病理图像进行颜色归一化，归一化后不同染色条件下的病理图像的颜色趋于一致，且计算速度快。
Description
技术领域
本发明涉及图像处理技术领域，更具体的说是涉及一种HE染色病理图像颜色归一化方法。
背景技术
数字病理全切片是通过专用的扫描成像系统将病理切片转换为可供计算机显示、传输及处理的高倍率大规模数字图像。随着数字病理图像成像技术不断成熟，基于数字病理图像的计算机辅助诊断方法发展迅速。随着研究的深入，辅助诊断方法已向适应多染色平台、扫描平台的病理全切片(Whole Slide Image，WSI)方向发展。然而，受染色剂配比、染色平台、成像平台的影响，来自不同平台的WSI在颜色上的差异很大。这一差异极大的影响了计算机辅助诊断方法的准确度。为了减少这一差异带来的影响，国内外学者在近年提出了多种病理图像颜色归一化方法。其中，苏木精-伊红(HE)染色剂作为病理图像最常用的颜色方法。
目前，有学者直接采用自然场景图像颜色增强的方法来对数字病理图像进行归一化，如直方图均衡，直方图标准化，Retinex增强等。这些方法没有考虑病理图像颜色的本质特点，固不能取得满意的病理图像归一化效果。大多数有效的病理图像归一化方法以基于颜色反卷积的染色分离方法作为基础，颜色反卷积算法利用“病理图像的颜色是由有限种染色剂混合而成的”这一先验知识，在光密度空间对病理图像进行线性变换，以分离出独立的颜色成分。然而，颜色反卷积算法需要通过实验的手段测得反卷积参数，且同一组参数无法适应来自不同平台的病理图像。近年，有学者利用细胞核检测算法，在病理图像中检测出细胞核，从而将病理图像中细胞核与细胞质区域进行分类，并以此作为先验，计算得到自适应的颜色反卷积参数。这种方法使颜色反卷积方算法不再依赖实验获得反卷积参数，极大的提高该算法的适应性，从而得到了更好的病理图像归一化效果。但是，该类方法的效果依赖于细胞核检测与分割算法，一旦细胞核检测算法失效，则无法得到合理的归一化效果，且细胞核检测与分割算法的计算量较大，拖慢了整个归一化过程。
因此，如何提供一种快速HE染色病理图像颜色归一化方法是本领域技术人员亟需解决的问题。
发明内容
有鉴于此，本发明提供了一种HE染色病理图像颜色归一化方法，通过优化算法得到最优反卷积矩阵，并利用最优反卷积矩阵对图像进行归一化，归一化后，不同染色条件下的病理图像的颜色趋于一致，且计算速度快。
为了实现上述目的，本发明采用如下技术方案：
一种HE染色病理图像颜色归一化方法，包括：
S1：对预先存储好的数字病理切片图像进行通道分解，得到病理图像I(x,y)；
其中，IR(x,y)、IG(x,y)、IB(x,y)分别表示病理图像在RGB三个颜色通道内的像素值；
S2：建立优化模型，计算最优反卷积矩阵，其具体步骤如下：
S21：定义优化变量
吸光率矩阵中的每个吸光率向量ms均为一个单位向量，用两个独立的角度变量(αc,βc)表示吸光率向量：
则反卷积矩阵由6个独立的角度变量决定，由此定义待优化变量：
S22：定义目标函数
将待优化变量决定的吸光率矩阵表示为通过计算得到反卷积矩阵
利用反卷积矩阵对所述病理图像I(x,y)进行染色分离，公式如下：
其中
其中，Imax为单通道最大值；表示病理图像I(x,y)的染色成分，分别代表病理图像I(x,y)的H、E、DAB染色剂对应的染色强度；
定义损失函数：
优化目标函数为：
其中，为基本损失函数：
为平衡损失函数：
为能量损失函数：
其中，P(I)表述病理图像I(x,y)中参与上述优化的像素点，在病理图像的非背景区域随机抽取，P(I)表示P(I)中包含像素点的数量，λ为权重系数，γ为平衡系数，ε为能量系数；
S23：利用梯度下降算法对目标函数进行求解，得到求解后的最优参数
其中，求解时的初始条件为
λ,γ,ε为超参数，在优化前预先设定；
S24：基于所述最优参数，计算最优的反卷积矩阵，具体包括：
根据得到的计算得到优化后的吸光率矩阵：
计算最优反卷积矩阵：
S3:根据所述最优反卷积矩阵对所述病理图像进行颜色归一化，具体步骤如下：
S31：计算归一化矩阵：
其中，为理想成像条件下的吸光率矩阵
S32：利用归一化矩阵T对所述病理图像进行归一化：
即为归一化后的病理图像，其中Oc(x,y)表示病理图像在RGB通道的光密度：Oc(x,y)＝-log10(Ic(x,y)/Imax),c＝R,G,B，Imax为单通道最大值。
优选的，所述数字病理切片图像通过切片扫描仪将病理切片扫描到计算机中而得到的。
优选的，对目标函数进行求解的方法包括但不限于SGD，SGBM，AdaGrad，AdaDelta或Adam。
优选的，对目标函数进行求解的方法推荐使用AdaGrad。
优选的，λ＝0.02,γ＝0.82,ε＝0.40
经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种HE染色病理图像颜色归一化方法，通过定义优化变量以及损失函数并设定优化目标计算最优参数基于最优参数得到优化后的吸光率矩阵从而得到优化后的反卷积矩阵基于优化后的反卷积矩阵计算归一化矩阵T，并基于归一化矩阵T对图像进行归一化，得到归一化的病理图像。利用本发明提供的HE染色病理图像颜色归一化方法对病理图像进行颜色归一化，得到的病理图像的颜色趋于一致，且计算速度快。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的一种HE染色病理图像颜色归一化方法的流程图；
图2附图为本发明提供的建立优化模型计算最优反卷积矩阵的流程图；
图3附图为本发明提供的根据最优反卷积矩阵对病理图像进行颜色归一化的流程图；
图4附图为本发明提供的HE染色病理图像颜色归一化方法的效果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见附图1，本发明实施例公开了一种HE染色病理图像颜色归一化方法，包括：
S1：对预先存储好的数字病理切片图像进行通道分解，得到病理图像I(x,y)；
其中，IR(x,y)、IG(x,y)、IB(x,y)分别表示病理图像在RGB三个颜色通道内的像素值；
S2：建立优化模型，计算最优反卷积矩阵，其具体步骤如下：
S21：定义优化变量
吸光率矩阵中的每个吸光率向量ms均为一个单位向量，用两个独立的角度变量(αc,βc)表示吸光率向量：
则反卷积矩阵由6个独立的角度变量决定，由此定义待优化变量：
S22：定义目标函数
为了既便于表达，将待优化变量决定的吸光率矩阵表示为通过计算得到反卷积矩阵
利用反卷积矩阵对所述病理图像I(x,y)进行染色分离，公式如下：
其中，
其中，Imax为单通道最大值；表示病理图像I(x,y)的染色成分，分别代表病理图像I(x,y)的H、E、DAB染色剂对应的染色强度；
定义损失函数：
优化目标函数为：
其中，为基本损失函数：
HE染色的病理图像没有DAB染色成分，因此，在对HE染色的病理图像进行染色分离时，DAB分量的理想值为0；另一方面，希望H和E两种染色成分的差异尽可能大，由此定义了基本损失函数；
为平衡损失函数：
引入平衡损失函数的目的是为了对分解后的H和E染色成分比率进行控制；
为能量损失函数：
引入能量损失函数的目的是为了对分解后的染色成分总能量进行控制；
其中，P(I)表述病理图像I(x,y)中参与上述优化的像素点，在病理图像的非背景区域随机抽取，|P(I)|表示P(I)中包含像素点的数量，λ为权重系数，γ为平衡系数，ε为能量系数；
为了达到稳定的优化效果，像素点的个数不少于10000。
S23：利用梯度下降算法对目标函数进行求解，得到求解后的最优参数
由于损失函数是关于的连续函数，则可以用基于梯度下降的优化算法进行求解。为了得到合理的优化结果，须使用理想条件下测得的参数进行初始化，本发明使用进行初始化
在求解时的初始条件为
λ＝0.02,γ＝0.82,ε＝0.40；
上述λ＝0.02,γ＝0.82,ε＝0.40为超参数，需要在优化前人为设定，不同的参数组合会取得不同的归一化效果，此处给出一组较优的参数组合。
S24：计算最优的反卷积矩阵；
根据得到的计算得到优化后的吸光率矩阵：
计算最优反卷积矩阵：
S3:根据所述最优反卷积矩阵对所述病理图像进行颜色归一化，具体步骤如下：
S31：计算归一化矩阵：
其中，为理想成像条件下的吸光率矩阵
S32：利用归一化矩阵T对所述病理图像进行归一化：
即为归一化后的病理图像，其中Oc(x,y)表示病理图像在RGB通道的光密度：Oc(x,y)＝-log10(Ic(x,y)/Imax),c＝R,G,B，Imax为单通道最大值。
为了进一步优化上述技术方案，所述病理图像通过切片扫描仪将病理切片扫描到计算机中而得到的。
为了进一步优化上述技术方案，对目标函数进行求解的方法包括但不限于SGD，SGBM，AdaGrad，AdaDelta或Adam。
为了进一步优化上述技术方案，对目标函数进行求解的方法推荐使用AdaGrad，该方法明显优于其他几个对目标函数的求解方法。
下面结合现有的算法对本发明提供的HE染色病理图像颜色归一化方法作进一步说明。
现有技术中HE染色的病理图像一般使用基于H-E-DAB染色测得的卷积参数进行染色分离，在此给出使用传统颜色反卷积算法对HE染色的病理图像进行染色成分分离的步骤，以进一步帮助理解本发明的发明思路以及算法。
在得到R、G、B三通道的像素值后，计算RGB通道对应的光密度：
Oc(x,y)＝-log10(Ic(x,y)/Imax),c＝R,G,B
Imax为单通道最大值，一般情况下，数字图像取值范围是Ic(x,y)∈[0,255]，即Imax＝255。
计算单独染色剂的着色强度As(x,y),s＝H,E,DAB
其中，矩阵称为反卷积矩阵；
令O＝[OR(x,y),OG(x,y),OB(x,y)]T,A＝[AH(x,y),AE(x,y),ADAB(x,y)]T，则上述公式可以简写为:A＝D·O，A＝[AH(x,y),AE(x,y),ADAB(x,y)]T即为分解后的染色强度。
在上述算法中，反卷积矩阵D是由实验测量并进行计算得到的，具体方法如下：
a)计算吸光率矩阵：通过使用单独染色剂对切片染色，可以测得染色剂s(s＝H,E,DAB)对通道R,G,B的吸光率，用向量表示为：
其中c＝R,G,B，表示染色剂s对通道c的吸光率；ms是一个单位向量，即||ms||＝1，对染色剂H,E,DAB进行单独测量后得到吸光率矩阵；
M＝(mH,mE,mDAB)T，对于H-E-DAB染色的数字病理图像，利用上述方法计算得到的理想情况下的吸光率矩阵为：
b)计算反卷矩阵：反卷积矩阵即为吸光率矩阵的逆，则计算得到的理想情况下的反卷积矩阵(即反卷积参数)为：
然而，上述反卷积参数只适用于理想的染色情况，不能适应不同染色条件和成像条件，本发明则针对这一问题，提出了自适应染色分离方法,通过对吸光率矩阵M进行优化，以达到改善反卷积矩阵D的目的。
下面结合具体实验结果对本发明提出的颜色归一化方法做进一步说明：
图4中示出在使用超参数λ＝0.02,γ＝0.82,ε＝0.40进行归一化时，不同染色条件下一组病理图像的归一化效果。从图4(a)中为来自不同染色平台和成像平台的病理图像，图4(b)为经过本发明提供的颜色归一化方法处理得到的病理图像，从中可以看出，病理图像的颜色趋于一致。
为了进一步说明本发明的优势，将本发明中提供的HE染色病理图像颜色归一化方法和文献[1](Bejnordi B E,Litjens G,Timofeeva N,et al.Stain SpecificStandardization of Whole-Slide Histopathological Images[J].IEEE Transactionson Medical Imaging,2016,35(2):404.)和文献[2](Khan A M,Rajpoot N,Treanor D,etal.A nonlinear mapping approach to stain normalization in digitalhistopathology images using image-specific color deconvolution[J].IEEE TransBiomedical Engineering,2014,61(6):1729-38.)中的两种方法进行比较；以上两种方法均需要病理图像细胞核分割算法的支持，归一化效果受到细胞核分割算法性能的影响。相比之下，本发明算法不需要对细胞核进行分割，只依赖病理图像中的像素信息就可以完成归一化，算法的稳定性更强。另一方面，由于本发明提供的方法不需要提取病理图像中的细胞核，计算速度有很大提升，对比计算归一化参数所用时间，文献[1]需要123.73秒，文献[2]需要10.1分钟，本发明方法只需9.23秒。
该方法能够有效提高基于深度学习的病理图像分类方法的性能，对浸润性乳腺癌图像分类实验中，选定深度学习模型为DenseNet-121的情况下，归一化前分类精度为78.7％，归一化后分类精度提高到了84.7％，ROC曲线下面积(AUC)从0.868提高到了0.914。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (5)
1.一种HE染色病理图像颜色归一化方法，其特征在于，包括：
S1：对预先存储好的数字病理切片图像进行通道分解，得到病理图像I(x,y)；
其中，IR(x,y)、IG(x,y)、IB(x,y)分别表示病理图像在RGB三个颜色通道内的像素值；
S2：建立优化模型，计算最优反卷积矩阵，其具体步骤如下：
S21：定义优化变量
吸光率矩阵中的每个吸光率向量ms均为一个单位向量，用两个独立的角度变量(αc,βc)表示吸光率向量：
则反卷积矩阵由6个独立的角度变量决定，由此定义待优化变量：
S22：定义目标函数
将待优化变量决定的吸光率矩阵表示为通过计算得到反卷积矩阵
利用反卷积矩阵对所述病理图像I(x,y)进行染色分离，公式如下：
其中
其中，Imax为单通道最大值；表示病理图像I(x,y)的染色成分，分别代表病理图像I(x,y)的H、E、DAB染色剂对应的染色强度；
定义损失函数：
定义优化目标函数为：
其中，为基本损失函数：
为平衡损失函数：
为能量损失函数：
其中，P(I)表述病理图像I(x,y)中参与上述优化的像素点，在病理图像的非背景区域随机抽取，|P(I)|表示P(I)中包含像素点的数量，λ为权重系数，γ为平衡系数，ε为能量系数；
S23：利用梯度下降算法对目标函数进行求解，得到求解后的最优参数
其中，求解时的初始条件为
λ,γ,ε为超参数，在优化前预先设定；
S24：基于所述最优参数，计算最优的反卷积矩阵，具体包括：
根据得到的计算得到优化后的吸光率矩阵：
计算最优反卷积矩阵：
S3:根据所述最优反卷积矩阵对所述病理图像进行颜色归一化，具体步骤如下：
S31：计算归一化矩阵：
其中，为理想成像条件下的吸光率矩阵
S32：利用归一化矩阵T对所述病理图像进行归一化：
即为归一化后的病理图像，其中Oc(x,y)表示病理图像在RGB通道的光密度：Oc(x,y)＝-log10(Ic(x,y)/Imax),c＝R,G,B，Imax为单通道最大值。
2.根据权利要求1所述的HE染色病理图像颜色归一化方法，其特征在于，所述数字病理切片图像通过切片扫描仪将病理切片扫描到计算机中进行存储。
3.根据权利要求1所述的HE染色病理图像颜色归一化方法，其特征在于，对目标函数进行求解的方法包括但不限于SGD，SGBM，AdaGrad，AdaDelta或Adam。
4.根据权利要求3所述的HE染色病理图像颜色归一化方法，其特征在于，对目标函数进行求解的方法推荐使用AdaGrad。
5.根据权利要求1所述的HE染色病理图像颜色归一化方法，其特征在于，λ＝0.02,γ＝0.82,ε＝0.40。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810067454.9A CN108172278B (zh) | 2018-01-24 | 2018-01-24 | 一种he染色病理图像颜色归一化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810067454.9A CN108172278B (zh) | 2018-01-24 | 2018-01-24 | 一种he染色病理图像颜色归一化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108172278A true CN108172278A (zh) | 2018-06-15 |
| CN108172278B CN108172278B (zh) | 2021-09-10 |
Family
ID=62515307
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810067454.9A Active CN108172278B (zh) | 2018-01-24 | 2018-01-24 | 一种he染色病理图像颜色归一化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108172278B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109345601A (zh) * | 2018-09-21 | 2019-02-15 | 合肥工业大学 | 基于低秩嵌入非负矩阵分解的病理图像颜色标准化的方法 |
| CN111028923A (zh) * | 2019-10-18 | 2020-04-17 | 平安科技（深圳）有限公司 | 数字病理图像染色归一化方法、电子装置及存储介质 |
| CN111353987A (zh) * | 2020-03-02 | 2020-06-30 | 中国科学技术大学 | 一种细胞核分割方法及装置 |
| CN112435243A (zh) * | 2020-11-26 | 2021-03-02 | 山东第一医科大学附属省立医院(山东省立医院) | 一种全切片数字病理图像的自动分析系统及方法 |
| CN115994874A (zh) * | 2023-03-22 | 2023-04-21 | 赛维森(广州)医疗科技服务有限公司 | 玻片图像处理方法、装置、玻片、计算机设备和存储介质 |
| WO2023228229A1 (ja) * | 2022-05-23 | 2023-11-30 | 日本電気株式会社 | 画像処理装置、画像処理方法、およびプログラム |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090324126A1 (en) * | 2008-06-26 | 2009-12-31 | Microsoft Corporation | Image deconvolution using color priors |
| US20140016853A1 (en) * | 2012-07-13 | 2014-01-16 | Sony Corporation | Method and apparatus for stain separation using vector analysis |
| CN104408717A (zh) * | 2014-11-24 | 2015-03-11 | 北京航空航天大学 | 一种基于着色分离的病理图像颜色质量综合评价方法 |
| CN105354581A (zh) * | 2015-11-10 | 2016-02-24 | 西安电子科技大学 | 融合颜色特征与卷积神经网络的彩色图像特征提取方法 |
| CN105741266A (zh) * | 2016-01-22 | 2016-07-06 | 北京航空航天大学 | 一种病理图像细胞核快速定位方法 |
-
2018
- 2018-01-24 CN CN201810067454.9A patent/CN108172278B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090324126A1 (en) * | 2008-06-26 | 2009-12-31 | Microsoft Corporation | Image deconvolution using color priors |
| US20140016853A1 (en) * | 2012-07-13 | 2014-01-16 | Sony Corporation | Method and apparatus for stain separation using vector analysis |
| CN104408717A (zh) * | 2014-11-24 | 2015-03-11 | 北京航空航天大学 | 一种基于着色分离的病理图像颜色质量综合评价方法 |
| CN105354581A (zh) * | 2015-11-10 | 2016-02-24 | 西安电子科技大学 | 融合颜色特征与卷积神经网络的彩色图像特征提取方法 |
| CN105741266A (zh) * | 2016-01-22 | 2016-07-06 | 北京航空航天大学 | 一种病理图像细胞核快速定位方法 |
Non-Patent Citations (1)
| Title |
|---|
| 孟军 等: "一种新的双重标记免疫组化图像自动分割方法", 《电子显微学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109345601A (zh) * | 2018-09-21 | 2019-02-15 | 合肥工业大学 | 基于低秩嵌入非负矩阵分解的病理图像颜色标准化的方法 |
| CN109345601B (zh) * | 2018-09-21 | 2021-09-17 | 合肥工业大学 | 基于低秩嵌入非负矩阵分解的病理图像颜色标准化的方法 |
| CN111028923A (zh) * | 2019-10-18 | 2020-04-17 | 平安科技（深圳）有限公司 | 数字病理图像染色归一化方法、电子装置及存储介质 |
| CN111028923B (zh) * | 2019-10-18 | 2024-01-30 | 平安科技（深圳）有限公司 | 数字病理图像染色归一化方法、电子装置及存储介质 |
| CN111353987A (zh) * | 2020-03-02 | 2020-06-30 | 中国科学技术大学 | 一种细胞核分割方法及装置 |
| CN112435243A (zh) * | 2020-11-26 | 2021-03-02 | 山东第一医科大学附属省立医院(山东省立医院) | 一种全切片数字病理图像的自动分析系统及方法 |
| WO2023228229A1 (ja) * | 2022-05-23 | 2023-11-30 | 日本電気株式会社 | 画像処理装置、画像処理方法、およびプログラム |
| CN115994874A (zh) * | 2023-03-22 | 2023-04-21 | 赛维森(广州)医疗科技服务有限公司 | 玻片图像处理方法、装置、玻片、计算机设备和存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108172278B (zh) | 2021-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108172278B (zh) | 一种he染色病理图像颜色归一化方法 | |
| Zanjani et al. | Stain normalization of histopathology images using generative adversarial networks | |
| Salehi et al. | Pix2pix-based stain-to-stain translation: A solution for robust stain normalization in histopathology images analysis | |
| CN109872285A (zh) | 一种基于变分约束的Retinex低照度彩色图像增强方法 | |
| CN109325938B (zh) | 一种基于光密度空间的中心化图像通道差的色素分离方法 | |
| CN107705336A (zh) | 一种病理图像染色成分调节方法 | |
| CN111223110B (zh) | 一种显微图像增强方法、装置及计算机设备 | |
| CN110298792A (zh) | 低光照图像增强与去噪方法、系统及计算机设备 | |
| CN116664462B (zh) | 一种基于ms-dsc和i_cbam的红外和可见光图像融合方法 | |
| CN107203980B (zh) | 自适应多尺度暗通道先验的水下目标探测图像增强方法 | |
| Peng et al. | Multi-algorithm fusion of RGB and HSV color spaces for image enhancement | |
| CN105761292A (zh) | 一种基于颜色转移和修正的图像着色方法 | |
| CN111882555B (zh) | 基于深度学习的网衣检测方法、装置、设备及存储介质 | |
| CN115731146B (zh) | 基于色彩梯度直方图特征光流估计多曝光图像融合方法 | |
| CN112102928B (zh) | 病理图像染色风格归一化方法及装置 | |
| Parihar et al. | A comprehensive analysis of fusion-based image enhancement techniques | |
| CN105118076A (zh) | 基于过度分割和局部与全局一致性的图像彩色化方法 | |
| Chen et al. | Semantic attention and relative scene depth-guided network for underwater image enhancement | |
| CN112330790A (zh) | 基于对抗学习与自适应色度修正的ct图像自动着色方法 | |
| Yan et al. | A novel biologically-inspired method for underwater image enhancement | |
| Xiao et al. | USLN: A statistically guided lightweight network for underwater image enhancement via dual-statistic white balance and multi-color space stretch | |
| CN117593235A (zh) | 基于深度CNN去噪先验的Retinex变分水下图像增强方法和装置 | |
| Pierre et al. | Exemplar-based colorization in RGB color space | |
| Song et al. | Dual-model: Revised imaging network and visual perception correction for underwater image enhancement | |
| CN113810683A (zh) | 一种客观评估水下视频质量的无参考评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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 |