CN107833235A - 一种构建心肌细胞图像灰度层级结构信息的方法及系统 - Google Patents
一种构建心肌细胞图像灰度层级结构信息的方法及系统 Download PDFInfo
- Publication number
- CN107833235A CN107833235A CN201710843667.1A CN201710843667A CN107833235A CN 107833235 A CN107833235 A CN 107833235A CN 201710843667 A CN201710843667 A CN 201710843667A CN 107833235 A CN107833235 A CN 107833235A
- Authority
- CN
- China
- Prior art keywords
- image
- cardiac muscle
- muscle cell
- layers
- information
- 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
Classifications
-
- 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
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种构建心肌细胞图像灰度层级结构信息的方法及系统,用以解决现有技术不能很好地检测与分析心肌细胞染色图像中细胞核独立情况与粘连情况的问题。该方法包括:S1、将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;S2、采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;S3、采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。本发明通过构建染色图像灰度层级结构信息,并基于层级结构信息关系检测和分析心肌细胞核粘连情况与独立情况,从而理解心肌细胞免疫组化图像内部的场景信息,有利于心肌细胞免疫组化图像的后处理。
Description
技术领域
本发明涉及医学图像处理技术领域,尤其涉及一种构建心肌细胞图像灰度层级结构信息的方法及系统。
背景技术
心肌细胞又称心肌纤维,有横纹,受植物性神经支配,属于有横纹的不随意肌,具有兴奋收缩的能力。呈短圆柱形,有分支,其细胞核位于细胞中央,一般只有一个。各心肌纤维分支的末端可相互连接构成肌纤维网。广义的心肌细胞包括组成窦房结、房内束、房室交界部、房室束(即希斯束)和浦肯野纤维等的特殊分化了的心肌细胞,以及一般的心房肌和心室肌工作细胞。
细胞运动的研究一直是细胞学和生物学研究的重要组成部分,但是传统的研究方法在技术日益革新的现在已经渐渐变得不这么适用了,传统的在显微镜下利用细胞计数板用人眼进行染色、分类、计数、跟踪等这类不但需要大量繁琐的人为操作而且容易使得操作者变得疲劳从而影响结果的正确性,且其可重用性比较低。国内外的医学专家经过长期地实践与研究取得一致共识认为应该在细胞运动研究中引入数字视频技术和数字图像处理技术,从而极大地提高研究效率,减轻研究人员负担。用计算机来跟踪细胞运动,部分代替人类始终盯着显微镜来观察细胞的眼睛,尽量地把人从繁重的重复劳动中解脱出来进行更有创造性的工作。
因此,如何利用计算机图像处理、视频分析等相关技术手段来提高生物研究过程中的自动化程度已成为当前急需解决的非常有意义的难题。
公开号为CN103559724A的专利提供了一种高粘连度细胞环境下的多细胞同步跟踪方法。细胞序列图像中,多细胞的分割和同步跟踪是一个尚未解决的难题,尤其在高粘连度情况下多细胞检测与分割,更加迫切需要解决。该发明首先提出了一种改进的基于分水岭和多特征匹配的分割算法实现细胞分割,然后,建立适用于卡尔曼滤波的运动模型并加入多特征匹配实现细胞的预测和跟踪。该发明不能很好地处理心肌细胞核粘连情况与独立情况的分析。
发明内容
本发明要解决的技术问题目的在于提供一种构建心肌细胞图像灰度层级结构信息的方法及系统,用以检测心肌细胞核染色响应位点与响应强度,从而用以解决现有技术不能很好地检测与分析心肌细胞免疫组化染色图像中细胞核独立情况与粘连情况等问题。
为了实现上述目的,本发明采用的技术方案为:
一种构建心肌细胞图像灰度层级结构信息的方法,包括步骤:
S1、将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
S2、采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
S3、采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
进一步地,其特征在于,所述步骤S1具体包括:
基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
进一步地,其特征在于,所述步骤S2具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤K≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
那么,调整参数C的值便能获得不同的图像灰度层级结构信息结果。
进一步地,所述最佳参数C的确定方法为:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合的信息都称为p 层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
…
…
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
…
…
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3) 的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
一种构建心肌细胞图像灰度层级结构信息的系统,具体包括:
染色模块,用于将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
算法模块,用于采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
评估模块,用于采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
进一步地,其特征在于,所述染色模块具体用于基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
进一步地,其特征在于,所述算法模块具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤K≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
那么,调整参数C的值便能获得不同的图像灰度层级结构信息结果。但是,并不是任意参数都能得到较好的图像划分结果,需要确定最佳算法参数进一步地,所述算法模块包括确定最佳参数C单元,具体用于:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合的信息都称为p 层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与 p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
…
…
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
…
…
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3) 的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
本发明与传统的技术相比,有如下优点:
本发明通过构建染色图像灰度层级结构信息,并基于层级结构信息关系检测和分析心肌细胞核粘连情况与独立情况,从而理解心肌细胞免疫组化图像内部的场景信息,有利于心肌细胞免疫组化图像的后处理。
附图说明
图1是本发明实施例提供的一种构建心肌细胞图像灰度层级结构信息的方法流程图;
图2是本发明实施例提供的一种构建心肌细胞图像灰度层级结构信息的系统结构图;
图3是本发明实施例提供的三幅心肌细胞染色图像;
图4是图3中每幅心肌细胞染色图像对应的红色、绿色以及蓝色染色信息;
图5是本发明实施例提供的数值整体关系;
图6是本发明实施例提供的数值局部关系;
图7是图3中每幅心肌细胞染色图像中蓝色染色信息的最佳灰度层级结构信息。
图8是图3中Pos_0001和Pos_0002的红色染色信息的最佳灰度层级结构信息。
具体实施方式
以下是本发明的具体实施例,并结合附图对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。
本发明实施例提供了一种构建心肌细胞图像灰度层级结构信息的方法,如图1所示,包括步骤:
S11:将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
S12:采用多阈值最大类间方差算法依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
S13:确定最佳参数以及最佳灰度层级结构信息。
现有技术不能很好地检测和分析心肌细胞核粘连情况与独立情况,为了更好地基于层级结构信息关系检测心肌细胞核粘连情况与独立情况,需要构建心肌细胞染色图像的灰度层级结构信息。
本实施例中,步骤S11具体包括:
基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
具体的,基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,用于定位不同蛋白的表达位点和表达强度。
如图3所示,步骤S11具体为大多数心肌细胞核内部被染成蓝色,但由于某些细胞核内部同时存在不同特异性的蛋白,这些心肌细胞核区域内包含了两种颜色,即红色与蓝色。通过区分红色、绿色以及蓝色的染色信息,能将三幅原图分别分为对应的三幅具有不同颜色的染色信息图像,如图4所示。
此时假设x={1,2,3}。其中,图4(ax)为原图的红色染色信息,图4(bx)为原图的绿色染色信息,图4(cx)为原图的蓝色染色信息。针对每一幅图像,因其仅反映某一种特异性蛋白在心肌细胞组织上的分布情况,其也能单独地被看作是一幅灰度图像。由于绿色染色信息响应区域较大且不是需要被重点关注研究的细胞核区域,所以在之后的处理过程中不考虑。那么仅需要针对原图中的红色与蓝色的染色响应信息进行处理。由于Neg_0001中不存在红色染色信息,所以之后处理中仅考虑Neg_0001_B的蓝色染色信息。
首先,需要构建针对待处理染色信息的灰度层级结构信息,对不同的响应位点与响应强度进行一定地评估。此处,步骤S12 采用多阈值最大类间方差算法构建图像灰度层级结构信息。
本实施例中,步骤S12具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤K≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
具体的,假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤K≤tC-1<L;
那么,此时便能得到C个灰度集合,且这些集合之间存在包含与被包含的关系。假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
通常,多阈值最大类间方差算法的参数C都是针对特定的图像通过人为设置的。为了能够更为自动化地设定最佳的参数C,针对蓝色染色信息和红色染色信息评价了在[1,r)区间内的最佳数值,通常r=10。当C数值过大的时候,相当于对灰度图像细节进行了过度地划分,而当C数值过小的时候,图片中的结构信息不是非常突出,不利于之后的处理。以下阐述了确定多阈值最大类间方差算法最佳参数C的评价方法。
本实施例中,步骤S13中所述最佳参数C的确定方法为:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合的信息都称为p 层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与 p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
…
…
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
…
…
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3) 的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
具体的,假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合的信息都称为p层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与 p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
…
…
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
…
…
为了检验条件是否满足,需要求解此2p-1个等式,从而求得这2p-2个未知数。但是求解等式的方式确定2p-2个未知数的数值比较繁琐,且求得的2p-2个数值仅用于判定是否为正整数。
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(1)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
那么,从第1层一直不断地拓展,直至当p层存在不满足(2) 或者(3)的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
本发明实施例还提供了一种构建心肌细胞图像灰度层级结构信息的系统,如图2所示,包括:
染色模块21,用于将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
算法模块22,用于采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
评估模块23,用于采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
具体的,基于不同特异性蛋白的结合特性,将具有不同蛋白表达的细胞核区域染色为具有不同颜色,用于定位不同蛋白的表达位点和表达强度。染色模块21中大多数细胞核内部被染成蓝色,但由于某些细胞核内部同时存在不同特异性的蛋白,这些细胞核区域内包含了两种颜色,即红色与蓝色。通过区分红色、绿色以及蓝色的染色信息,能将三幅原图分别分为对应的三幅不同染色信息图像。算法模块22采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息。评估模块23采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
本实施例中,染色模块21具体用于基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
本实施例中,算法模块22具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤...≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
那么,调整参数C的值能获得不同的灰度图像层级结构信息结果。但是,并不是任意参数都能得到较好的图像划分结果,需要确定最佳算法参数。
本实施例中,算法模块22包括确定最佳参数C单元,具体用于:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合相关的数值信息都称为p层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与 p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
…
…
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
…
…
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3) 的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.一种构建心肌细胞图像灰度层级结构信息的方法,其特征在于,包括步骤:
S1、将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
S2、采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
S3、采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
2.根据权利要求1所述的一种构建心肌细胞图像灰度层级结构信息的方法,其特征在于,所述步骤S1具体包括:
基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
3.根据权利要求1所述的一种构建心肌细胞图像灰度层级结构信息的方法,其特征在于,所述步骤S2具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤...≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
4.根据权利要求1所述的一种构建心肌细胞图像灰度层级结构信息的方法,其特征在于,步骤S3所述多阈值最大类间方差算法的最佳参数C的确定方法具体包括:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合相关的数值信息都称为p层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3)的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
5.一种构建心肌细胞图像灰度层级结构信息的系统,其特征在于,包括:
染色模块,用于将心肌细胞区域染色为具有不同颜色的响应位点和响应强度;
算法模块,用于采用多阈值最大类间方差算法,依据不同参数构建针对不同染色信息的不同灰度层级结构信息;
评估模块,用于采用最佳参数选取方法确定最佳参数,并将其所对应的响应位点和响应强度的评估结果作为最佳灰度层级结构信息。
6.根据权利要求5所述的一种构建心肌细胞图像灰度层级结构信息的系统,其特征在于,所述染色模块具体用于基于不同特异性蛋白的结合特性,将具有不同蛋白表达的所述心肌细胞区域染色为红色、绿色以及蓝色,且心肌细胞核区域主要呈现为红色与蓝色。
7.根据权利要求5所述的一种构建心肌细胞图像灰度层级结构信息的系统,其特征在于,所述算法模块具体包括:
针对任意一种颜色的染色信息图像,视为是灰度信息图像。此时假设一幅灰度信息图像I中存在灰度级别在0到L之间的N个像素,且处于灰度级别i的像素个数标记为ni,那么灰度级别i的像素在此图片中的概率为:
pi=ni/N;
假设图像中存在C-1个阈值,记为{t1,t2,...,tC-1};第一个集合G1由灰度范围为[0,...,t1]的像素确定;第i个集合由灰度范围为(ti-1,...,ti]的像素确定,标记为Gi;最后一个集合由灰度范围为(tC-1,...,L]的像素确定,标记为GC;
假设μT是此幅图像的灰度均值,则多阈值最大类间方差算法的目标函数定义为:
其中,
0≤t1≤K≤tC-1<L;
假设此时存在两个相邻集合,记为Gi与Gj,且|i-j|=1;那么此两集合的某些局部区域能直接融合为一个局部连通区域,且融合边界无间隙。
那么,调整参数C的值便能获得不同的图像灰度层级结构信息结果。
8.根据权利要求5所述的一种构建心肌细胞图像灰度层级结构信息的系统,其特征在于,所述评估模块具体包括:
假设图像I已被划分为m个集合,则第k个集合中的像素个数记为其中1≤k≤m;随着值m的变化,相应地也会得到关于像素个数的数值集合,那么将图像划分为p个集合的信息都称为p层信息;如此,便能形成层间数值关系,即处于p层的第q个数值记录仅与p+1层的第q个数值与第q+1个数值相互关联;
假设p层的第q个集合中的像素个数(记为)是由与中的一部分像素构成的,那么此时存在连续两层(即p层与p+1层)的数值集合信息,即与在p+1层内部,存在如下的p-1个数值等式关系:
其中,需要保证其中k={2,...,p};
如此,便能得到如下p个层间等式关系:
为了简化等式求解与判别,随着类别个数m的增加,为了满足条件,存在符合如下标准的数值关系:
判别时,除了需要满足如公式(2)所示的数值关系,也需要考虑集合灰度均值。此时,假设针对存在其所对应的灰度均值此时需要满足如下的数值不等式关系:
从第1层一直不断地拓展,直至当p层存在不满足(2)或者(3)的时候,将参数p作为标准参数数值,赋值给多阈值最大类间方差算法的参数C,并将其对应的结果作为最佳多层级结构信息结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710843667.1A CN107833235B (zh) | 2017-09-19 | 2017-09-19 | 一种构建心肌细胞图像灰度层级结构信息的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710843667.1A CN107833235B (zh) | 2017-09-19 | 2017-09-19 | 一种构建心肌细胞图像灰度层级结构信息的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107833235A true CN107833235A (zh) | 2018-03-23 |
CN107833235B CN107833235B (zh) | 2020-08-21 |
Family
ID=61643413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710843667.1A Active CN107833235B (zh) | 2017-09-19 | 2017-09-19 | 一种构建心肌细胞图像灰度层级结构信息的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107833235B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020585A (zh) * | 2012-11-06 | 2013-04-03 | 华南师范大学 | 一种免疫组织阳性细胞和阴性细胞识别方法 |
CN104331893A (zh) * | 2014-11-14 | 2015-02-04 | 东南大学 | 一种复杂图像多阈值分割方法 |
CN104835149A (zh) * | 2015-04-20 | 2015-08-12 | 浙江工业大学 | 一种牛奶体细胞图像分割方法 |
US20170060348A1 (en) * | 2015-08-26 | 2017-03-02 | Sap Se | Compact display of hierarchical structure on user interface |
CN106875394A (zh) * | 2016-12-22 | 2017-06-20 | 浙江工业大学 | 适用于无帧间参考信息的相差显微镜细胞主要信息粘连情况的分离与分组方法 |
-
2017
- 2017-09-19 CN CN201710843667.1A patent/CN107833235B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020585A (zh) * | 2012-11-06 | 2013-04-03 | 华南师范大学 | 一种免疫组织阳性细胞和阴性细胞识别方法 |
CN104331893A (zh) * | 2014-11-14 | 2015-02-04 | 东南大学 | 一种复杂图像多阈值分割方法 |
CN104835149A (zh) * | 2015-04-20 | 2015-08-12 | 浙江工业大学 | 一种牛奶体细胞图像分割方法 |
US20170060348A1 (en) * | 2015-08-26 | 2017-03-02 | Sap Se | Compact display of hierarchical structure on user interface |
CN106875394A (zh) * | 2016-12-22 | 2017-06-20 | 浙江工业大学 | 适用于无帧间参考信息的相差显微镜细胞主要信息粘连情况的分离与分组方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107833235B (zh) | 2020-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112070772B (zh) | 基于UNet++和ResNet的血液白细胞图像分割方法 | |
CN105205804B (zh) | 血细胞图像中白细胞的核浆分离方法、分类方法及装置 | |
CN110853009B (zh) | 基于机器学习的视网膜病理图像分析系统 | |
CN104680524B (zh) | 一种叶类蔬菜病害诊断方法 | |
CN109118485A (zh) | 基于多任务神经网络的消化道内镜图像分类及早癌检测系统 | |
CN108921825A (zh) | 基于深度学习的检测面部皮肤点状缺陷的方法及装置 | |
CN107330892A (zh) | 一种基于随机森林法的向日葵病害识别方法 | |
CN107845085A (zh) | 一种心肌细胞核粘连区域分离与分组的方法及系统 | |
CN105844621A (zh) | 一种印刷品质量检测方法 | |
CN102214306B (zh) | 叶片病斑识别方法及装置 | |
CN102930249A (zh) | 基于颜色和模型的农田害虫识别和计数方法 | |
CN109684967A (zh) | 一种基于ssd卷积网络的大豆植株茎荚识别方法 | |
CN105095870A (zh) | 基于迁移学习的行人重识别方法 | |
CN109117877A (zh) | 一种基于深度卷积生成对抗网络的黄颡鱼及其套养种识别方法 | |
CN109035254A (zh) | 基于改进K-means聚类的运动鱼体阴影去除及图像分割方法 | |
CN107480607A (zh) | 一种智能录播系统中站立人脸检测定位的方法 | |
CN103034872A (zh) | 基于颜色和模糊聚类算法的农田害虫识别方法 | |
CN109615616A (zh) | 一种基于abc-pcnn的裂缝识别方法及系统 | |
CN105095867A (zh) | 基于深度学习的快速动态人脸提取、识别方法 | |
Hasan et al. | Fish diseases detection using convolutional neural network (CNN) | |
CN109886346A (zh) | 一种心肌mri图像分类系统 | |
CN108288048A (zh) | 基于改进头脑风暴优化算法的面部情绪识别特征选择方法 | |
CN110288574A (zh) | 一种超声辅助诊断肝肿块系统及方法 | |
CN103020198B (zh) | 一种在视频中搜索目标的方法 | |
CN109087310A (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 |