CN110702615B - 一种彩色数码透射显微镜颜色校正方法 - Google Patents
一种彩色数码透射显微镜颜色校正方法 Download PDFInfo
- Publication number
- CN110702615B CN110702615B CN201910907878.6A CN201910907878A CN110702615B CN 110702615 B CN110702615 B CN 110702615B CN 201910907878 A CN201910907878 A CN 201910907878A CN 110702615 B CN110702615 B CN 110702615B
- Authority
- CN
- China
- Prior art keywords
- color
- rgb
- values
- matrix
- value
- 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.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005286 illumination Methods 0.000 claims abstract description 22
- 238000012512 characterization method Methods 0.000 claims abstract description 9
- 230000004044 response Effects 0.000 claims description 87
- 239000011159 matrix material Substances 0.000 claims description 49
- 239000011521 glass Substances 0.000 claims description 19
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 230000007935 neutral effect Effects 0.000 claims description 11
- 238000002834 transmittance Methods 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000004456 color vision Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004438 eyesight Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种彩色数码透射显微镜颜色校正方法,利用微型透射色卡校正彩色数码透射显微镜所采集的显微图像的颜色,由于数码显微镜的照明光源通常是固定的,固定条件下的数码相机色度特征化特别适用于数码显微镜的颜色校正,因此本发明利用微型透射色卡中的色块作为训练样本建立彩色数码透射显微镜的正向色度特征化模型,利用该模型可以实现拍摄目标的颜色精确校正,提升彩色数码透射显微镜的颜色精度,符合人眼的颜色感知,实现所见即所得的效果,显著提高显微检测中目标对象的辨识度,实现目标对象的准确识别与定位。
Description
技术领域
本发明涉及一种显微镜颜色校正方法,尤其是涉及一种彩色数码透射显微镜颜色校正方法。
背景技术
在彩色数字成像***中,人们对图像的颜色质量提出了越来越高的要求。数字显微成像***的颜色形成有其特殊性,一方面要符合人眼的感观特性,另一方面还要尽可能真实地反映物体的原色。现有的彩色数码透射显微镜很大一部分应用于医学检测,许多病理特征依靠细胞的颜色变化来区分,对色彩真实性要求比较高。显微图像处理软件是数码显微成像***不可分割的一部分,它将图像传感器采集到的数据转换成和人眼视觉一致的彩色图像,并且提供强大的图像后期处理与分析功能。但处理软件对数字显微成像***的图像颜色质量要求也非常高,如果图像颜色质量不好,不仅影响视觉效果,还会对图像后期的处理、分析造成很大的负面影响,因此需要对彩色数码透射显微镜的图像进行颜色校正。目前颜色校正的解决方法大多采用白平衡算法,然而白平衡算法只是一种粗略的颜色校正方法,无法实现数字显微成像***颜色的精确校正。
发明内容
本发明所要解决的技术问题是提供一种能够提高图像颜色质量的彩色数码透射显微镜颜色校正方法。
本发明解决上述技术问题所采用的技术方案为:一种彩色数码透射显微镜颜色校正方法,包括以下步骤:
(1)在彩色数码透射显微镜使用彩色相机获得微型透射色卡中所有色块的RGB值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的RGB矩阵,其中,每一行表示一个色块对应的RGB 3通道响应值,每一列表示一个通道输出的所有色块的响应值;
(2)关闭显微镜照明光源后继续拍摄,获得彩色相机的暗电流响应值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的暗电流响应矩阵,其中,每一行表示一个色块对应的RGB的3通道的暗电流响应值,然后利用分光光度计测量微型透射色卡中所有色块的XYZ三刺激值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的三刺激矩阵,其中,每一行表示一个色块对应的XYZ三刺激值;
(3)将透过率为70%的玻璃片置于载物台上,使用彩色相机拍摄玻璃片的RGB图像,对RGB响应值进行光照空间不均匀性校正得,得到光照空间不均匀性校正后微型透射色卡的色块响应矩阵;
(4)利用微型透射色卡中的中性色块的三刺激值中的Y值及其对应的RGB值对RGB响应值进行非线性校正;
(5)建立彩色数码透射显微镜的正向色度特征化模型,即RGB值到XYZ三刺激值之间的转换关系;
(6)拍摄并采集目标对象的RGB图像,逐像素扣除暗电流并校正RGB响应值的光照空间不均匀性和非线性,得到目标对象校正后的响应值矩阵;
(7)利用目标对象的三刺激值计算其sRGB响应值,得到sRGB图像即为颜色校正后的显微图像。利用目标对象的三刺激值St计算其sRGB响应值,sRGB图像即为颜色校正后的显微图像。
具体的颜色校正方法为:
(1)打开彩色数码透射显微镜照明光源,预热一段时间,将微型透射色卡置于载物台目镜下方,拍摄微型透射色卡,获得该色卡图像中所有色块的RGB值,记为P为n行3列的矩阵,其中n为微型透射色卡上的色块数,3为彩色相机通道的数量,每一行表示一个色块对应的RGB 3通道响应值,每一列表示一个通道输出的n个色块的响应值。
(2)关闭显微镜照明光源并使用彩色相机拍摄,获得彩色相机的暗电流响应值,并记为表示D为n行3列的矩阵,其中n为微型透射色卡上的色块数,每一行表示一个色块对应的RGB 3通道的暗电流响应值,扣除暗电流后的色块响应值为PD=P-D,然后利用分光光度计测量微型透射色卡中所有色块的XYZ三刺激值,记为表示S为n行3列的矩阵,n为微型透射色卡上的色块数,每一行表示一个色块对应的XYZ三刺激值;
(3)将透过率为70%的玻璃片置于载物台上,使用彩色相机拍摄玻璃片的RGB图像,利用玻璃片的RGB值逐通道校正微型透射色卡中色块响应值的光照不均匀性,
其中,为色卡中第i个色块第k个通道校正后的响应值,k∈{R,G,B},为玻璃片图像第k个通道的平均响应值,为玻璃片图像第k个通道对应微型透射色卡图像中第i个色块处的响应值,为色卡中第i个色块第k个通道的响应值,得到光照空间不均匀性校正后色卡的色块响应矩阵,记为表示PL为n行3列的矩阵,n为微型透射色卡上的色块数,每一行表示一个色块校正后的RGB值,每一列表示所有色块某一颜色通道下的值;
(4)利用微型透射色卡中的中性色块的三刺激值中的Y值及其对应的RGB值,分别建立RGB三通道的非线性校正关系,
yi=c1pik 3+c2pik 2+c3pik+c4
式中,yi表示微型透射色卡上第i个中性色块的归一化Y刺激值,pik表示第i个中性色块在第k个通道下的归一化响应值,k∈{R,G,B},c1,c2,c3,c4分别为三次多项式的拟合系数,四个拟合系数的取值分别在0-3之间,并且不能同时为0,利用上式校正微型透射色卡中所有色块的RGB响应值PL,并将校正后的色块响应值记为表示Pc为n行3列的矩阵;
(5)建立彩色数码显微镜的正向色度特征化模型,即建立RGB值到XYZ三刺激值之间的转换关系,利用校正后的色块响应值Pc及其三刺激值S建立两者之间的转换数学模型,由最小二乘法计算两者之间的转换矩阵Z,
Z=STPc(Pc TPc)-1
其中,上标T表示矩阵的转置,上标-1表示矩阵的逆矩阵;
(6)使用彩色相机拍摄并采集目标对象的RGB图像,逐像素扣除暗电流并按公式和公式yi=c1pik 3+c2pik 2+c3pik+c4分别校正目标对象RGB响应值的光照空间不均匀性和非线性,得到目标对象校正后的响应值矩阵,记为表示Pt为m行3列的矩阵,m为图像中的像素数,利用转换矩阵Z计算目标对象的XYZ三刺激值St,
St=ZPt;
与现有技术相比,本发明的优点在于利用微型透射色卡校正彩色数码透射显微镜所采集的显微图像的颜色,由于数码显微镜的照明光源通常是固定的,固定条件下的数码相机色度特征化特别适用于数码显微镜的颜色校正,因此本发明利用微型透射色卡中的色块作为训练样本建立彩色数码透射显微镜的正向色度特征化模型,利用该模型可以实现拍摄目标的颜色精确校正,提升彩色数码透射显微镜的颜色精度,符合人眼的颜色感知,实现所见即所得的效果,显著提高显微检测中目标对象的辨识度,实现目标对象的准确识别与定位。
附图说明
图1是利用本发明的方法校正彩色数码透射显微镜颜色的流程图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
以一台生物显微镜为例,阐述基于微型透射色卡校正彩色数码透射显微镜颜色的方法,微型透射色卡采用包含45个色块的透射色卡,色卡尺寸为3mm×4mm。
本实施例基于微型透射色卡校正彩色数码透射显微镜颜色的方法,具体包括以下步骤:
(1)打开彩色数码透射显微镜照明光源,预热一段时间,将微型透射色卡置于载物台目镜下方,拍摄微型透射色卡,获得该色卡图像中所有色块的RGB值,记为P为45行3列的矩阵,其中45为微型透射色卡上的色块数,3为彩色相机通道的数量,每一行表示一个色块对应的RGB三通道响应值,每一列表示一个通道输出的45个色块的响应值。
(2)关闭显微镜照明光源并拍摄,获得色卡的暗电流响应值45为微型透射色卡上的色块数,每一行表示一个色块对应的RGB 3通道的暗电流响应值,扣除暗电流后的色块响应值为PD=P-D。利用爱色丽公司生产的分光光度计CE7000A测量微型透射色卡中所有色块的XYZ三刺激值45为微型透射色卡上的色块数,每一行表示一个色块对应的XYZ三刺激值。
(3)对RGB响应值进行光照空间不均匀性校正。将透过率为70%的玻璃片置于载物台上,拍摄玻璃片的RGB图像,利用玻璃片的RGB值逐通道校正色卡中色块响应值的光照不均匀性,校正公式如下,
其中,为玻璃片图像第k(k∈{R,G,B})个通道的平均响应值,为玻璃片图像第k个通道对应色卡图像中第i个色块处的响应值,为色卡中第i个色块第k个通道的响应值,为色卡中第i个色块第k个通道校正后的响应值。光照空间不均匀性校正后色卡的色块响应值矩阵为PL。
(4)对RGB响应值进行非线性校正。利用微型透射色卡中的中性色块的三刺激值中的Y值及其对应的RGB值,分别建立RGB三通道的非线性校正关系,如下式所示,
yi=c1pik 3+c2pik 2+c3pik+c4 (2)
上式中,yi表示第i个中性色块的归一化Y刺激值,pik表示第i个中性色块在第k(k∈{R,G,B})个通道下的归一化响应值,c1,c2,c3,c4分别为三次多项式的拟合系数,四个拟合系数的取值分别在0-3之间,并且不能同时为0,本实施例中四个拟合系数分别取值c1=1.637,c2=-1.997,c3=1.502,c4=-0.136。利用上式校正微型透射色卡中所有色块的RGB响应值PL,校正后的色块响应值记为
(5)建立彩色数码显微镜的正向色度特征化模型,即建立RGB值到XYZ三刺激值之间的转换关系。利用校正后的色块响应值Pc及其三刺激值S建立两者之间的转换数学模型,由最小二乘法计算两者之间的转换矩阵Z,即
Z=STPc(Pc TPc)-1 (3)
其中,上标T表示矩阵的转置,上标-1表示矩阵的逆矩阵。
(6)拍摄并采集生物组织切片的RGB图像,逐像素扣除暗电流并按式(1)和式(2)分别校正生物组织切片RGB响应值的光照空间不均匀性和非线性,得到生物组织切片校正后的响应值矩阵m为图像中的像素数。利用式(3)中的转换矩阵Z计算生物组织切片的XYZ三刺激值St,即
St=ZPt (4)
(7)利用生物组织切片的三刺激值St计算其sRGB响应值,如式(5)所示。sRGB图像即为颜色校正后的显微图像。
sRGB=MSt (5)
Claims (2)
1.一种彩色数码透射显微镜颜色校正方法,其特征在于,包括以下步骤:
(1)在彩色数码透射显微镜使用彩色相机获得微型透射色卡中所有色块的RGB值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的RGB矩阵,其中,每一行表示一个色块对应的RGB 3通道响应值,每一列表示一个通道输出的所有色块的响应值;
(2)关闭显微镜照明光源后继续拍摄,获得彩色相机的暗电流响应值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的暗电流响应矩阵,其中,每一行表示一个色块对应的RGB的3通道的暗电流响应值,然后利用分光光度计测量微型透射色卡中所有色块的XYZ三刺激值,并定义为一个以微型透射色卡上的色块数为行数,彩色相机通道数量为列数的三刺激矩阵,其中,每一行表示一个色块对应的XYZ三刺激值;
(3)将透过率为70%的玻璃片置于载物台上,使用彩色相机拍摄玻璃片的RGB图像,对RGB响应值进行光照空间不均匀性校正得,得到光照空间不均匀性校正后微型透射色卡的色块响应矩阵;
(4)利用微型透射色卡中的中性色块的三刺激值中的Y值及其对应的RGB值对RGB响应值进行非线性校正;
(5)建立彩色数码透射显微镜的正向色度特征化模型,即RGB值到XYZ三刺激值之间的转换关系;
(6)拍摄并采集目标对象的RGB图像,逐像素扣除暗电流并校正RGB响应值的光照空间不均匀性和非线性,得到目标对象校正后的响应值矩阵;
(7)利用目标对象的三刺激值计算其sRGB响应值,得到sRGB图像即为颜色校正后的显微图像。
2.如权利要求1所述的一种彩色数码透射显微镜颜色校正方法,其特征在于,具体的校正方法为:
(1)打开彩色数码透射显微镜照明光源,预热一段时间,将微型透射色卡置于载物台目镜下方,拍摄微型透射色卡,获得该色卡图像中所有色块的RGB值,记为 P为n行3列的矩阵,其中n为微型透射色卡上的色块数,3为彩色相机通道的数量,每一行表示一个色块对应的RGB 3通道响应值,每一列表示一个通道输出的n个色块的响应值;
(2)关闭显微镜照明光源并使用彩色相机拍摄,获得彩色相机的暗电流响应值,并记为表示D为n行3列的矩阵,其中n为微型透射色卡上的色块数,每一行表示一个色块对应的RGB 3通道的暗电流响应值,扣除暗电流后的色块响应值为PD=P-D,然后利用分光光度计测量微型透射色卡中所有色块的XYZ三刺激值,记为表示S为n行3列的矩阵,n为微型透射色卡上的色块数,每一行表示一个色块对应的XYZ三刺激值;
(3)将透过率为70%的玻璃片置于载物台上,使用彩色相机拍摄玻璃片的RGB图像,利用玻璃片的RGB值逐通道校正微型透射色卡中色块响应值的光照不均匀性,
其中,为色卡中第i个色块第k个通道校正后的响应值,k∈{R,G,B},为玻璃片图像第k个通道的平均响应值,为玻璃片图像第k个通道对应微型透射色卡图像中第i个色块处的响应值,为色卡中第i个色块第k个通道的响应值,得到光照空间不均匀性校正后色卡的色块响应矩阵,记为表示PL为n行3列的矩阵,n为微型透射色卡上的色块数,每一行表示一个色块校正后的RGB值,每一列表示所有色块某一颜色通道下的值;
(4)利用微型透射色卡中的中性色块的三刺激值中的Y值及其对应的RGB值,分别建立RGB三通道的非线性校正关系,
yi=c1pik 3+c2pik 2+c3pik+c4
式中,yi表示微型透射色卡上第i个中性色块的归一化Y刺激值,pik表示第i个中性色块在第k个通道下的归一化响应值,k∈{R,G,B},c1,c2,c3,c4分别为三次多项式的拟合系数,四个拟合系数的取值分别在0-3之间,并且不能同时为0,利用上式校正微型透射色卡中所有色块的RGB响应值PL,并将校正后的色块响应值记为表示Pc为n行3列的矩阵;
(5)建立彩色数码显微镜的正向色度特征化模型,即建立RGB值到XYZ三刺激值之间的转换关系,利用校正后的色块响应值Pc及其三刺激值S建立两者之间的转换数学模型,由最小二乘法计算两者之间的转换矩阵Z,
Z=STPc(Pc TPc)-1
其中,上标T表示矩阵的转置,上标-1表示矩阵的逆矩阵;
(6)使用彩色相机拍摄并采集目标对象的RGB图像,逐像素扣除暗电流并按公式和公式yi=c1pik 3+c2pik 2+c3pik+c4分别校正目标对象RGB响应值的光照空间不均匀性和非线性,得到目标对象校正后的响应值矩阵,记为表示Pt为m行3列的矩阵,m为图像中的像素数,利用转换矩阵Z计算目标对象的XYZ三刺激值St,
St=ZPt;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910907878.6A CN110702615B (zh) | 2019-09-25 | 2019-09-25 | 一种彩色数码透射显微镜颜色校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910907878.6A CN110702615B (zh) | 2019-09-25 | 2019-09-25 | 一种彩色数码透射显微镜颜色校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110702615A CN110702615A (zh) | 2020-01-17 |
CN110702615B true CN110702615B (zh) | 2022-06-17 |
Family
ID=69196049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910907878.6A Active CN110702615B (zh) | 2019-09-25 | 2019-09-25 | 一种彩色数码透射显微镜颜色校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110702615B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113033550B (zh) * | 2021-03-15 | 2022-04-15 | 合肥联宝信息技术有限公司 | 一种图像的检测方法、装置及计算机可读介质 |
CN113628135A (zh) * | 2021-07-28 | 2021-11-09 | Oppo广东移动通信有限公司 | 图像颜色校正方法、装置、计算机装置及存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101959072B (zh) * | 2010-07-30 | 2012-06-27 | 浙江大学 | 一种彩色表面投影显示颜色再现方法 |
CN103905803B (zh) * | 2014-03-18 | 2016-05-04 | 中国科学院国家天文台 | 一种图像的颜色校正方法及装置 |
CN105651386B (zh) * | 2016-03-04 | 2018-12-11 | 温州佳易仪器有限公司 | 一种利用标准光源对色观察箱进行颜色测试的方法 |
CN107389194B (zh) * | 2017-07-12 | 2018-10-12 | 蔚彩信息科技(上海)有限公司 | 基于色卡的开放光环境下颜色高保真数字成像*** |
CN109525749B (zh) * | 2018-11-07 | 2020-02-21 | 厦门亿联网络技术股份有限公司 | 一种校正图像色彩偏差的方法及装置 |
-
2019
- 2019-09-25 CN CN201910907878.6A patent/CN110702615B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110702615A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105960586B (zh) | 用于使用监督式高品质成像的图像采集的***和方法 | |
JP2005516188A5 (zh) | ||
CN103402117A (zh) | 基于Lab色度空间的视频图像偏色检测方法 | |
CN110702615B (zh) | 一种彩色数码透射显微镜颜色校正方法 | |
CN101933321A (zh) | 用于估测场景光源的图像传感器装置及方法 | |
CN107534727A (zh) | 图像处理装置、摄像装置、图像处理方法及图像处理程序 | |
CN110175969B (zh) | 图像处理方法及图像处理装置 | |
CN104113743A (zh) | 低照度下彩色摄像机自动白平衡处理方法及装置 | |
CN109242792B (zh) | 一种基于白色物体的白平衡校对方法 | |
CN115665565A (zh) | 在线烟叶图像颜色校正方法、***及装置 | |
CN101947101B (zh) | 一种舌色还原色卡的制作方法 | |
CN106296658B (zh) | 一种基于相机响应函数的场景光源估计准确性提升方法 | |
CN110726536B (zh) | 一种彩色数码反射显微镜颜色校正方法 | |
JP2015194567A (ja) | 表示装置 | |
CN107635124B (zh) | 人脸拍摄的白平衡处理方法、装置及设备 | |
EP2205002A1 (en) | Spectral characteristic correction device and spectral characteristic correction method | |
US8169475B2 (en) | Image processing system, imaging system, and microscope imaging system | |
CN106023238A (zh) | 一种相机模组颜色数据标定方法 | |
US7777786B2 (en) | Method for calculating color correction | |
CN111010557B (zh) | 白点范围确定方法、白平衡校正方法及装置、存储介质 | |
CN108735010A (zh) | 一种用于英语教学的智能型英语教学*** | |
CN109561292A (zh) | 一种用于相机***的白平衡处理方法及装置 | |
EP2204980B1 (en) | Image processing system, imaging system, and microscope imaging system | |
CN111896109B (zh) | 一种基于数码相机原始响应值的光谱重建方法 | |
CN101500074A (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 |