CN109285159A - 一种基于新型纳什均衡的医学图像分割方法 - Google Patents
一种基于新型纳什均衡的医学图像分割方法 Download PDFInfo
- Publication number
- CN109285159A CN109285159A CN201810870190.0A CN201810870190A CN109285159A CN 109285159 A CN109285159 A CN 109285159A CN 201810870190 A CN201810870190 A CN 201810870190A CN 109285159 A CN109285159 A CN 109285159A
- Authority
- CN
- China
- Prior art keywords
- node
- entropy
- nash equilibrium
- nodes
- novel
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000003709 image segmentation Methods 0.000 claims abstract description 35
- 230000006870 function Effects 0.000 claims description 35
- 101150055297 SET1 gene Proteins 0.000 claims description 16
- 101150117538 Set2 gene Proteins 0.000 claims description 13
- 230000009977 dual effect Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000002776 aggregation Effects 0.000 claims 1
- 238000004220 aggregation Methods 0.000 claims 1
- 238000012706 support-vector machine Methods 0.000 description 35
- 230000011218 segmentation Effects 0.000 description 11
- 238000011160 research Methods 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 208000032612 Glial tumor Diseases 0.000 description 1
- 206010018338 Glioma Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 201000007983 brain glioma Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000012549 training 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/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
一种基于新型纳什均衡的医学图像分割方法,属于医学图像分割技术领域。本发明包括以下步骤:获取待分割图像的所有节点坐标和灰度值,初始化所有节点为两个节点集合,即对象集合和背景集合,按构建的近似熵和标准差双重约束下的纳什均衡确定每个节点实际所属集合,对于确定的两个节点集合给出每个节点的惩罚参数和核函数,使用新型纳什均衡修改后的SVM确定每个节点的决策函数,根据所有节点的决策函数对待分割图像进行分割。本发明在纳什均衡簇中引入双重分配约束,解决了决定节点簇的轮廓不明确的问题;基于所提出的纳什均衡来修改SVM,实现了更好的最大边缘,提高了医学图像分割的准确性。
Description
技术领域
本发明属于医学图像分割技术领域,具体涉及一种基于新型纳什均衡的医学图像分割方法。
背景技术
图像分割是将图像分成几个不同字符的子图像,以便提取图像中的一些感兴趣的对象。它是由图像处理到图像分析的关键步骤。现有的图像分割方法主要分以下几类:基于阈值的分割方法、基于区域的分割方法、基于边缘的分割方法以及基于特定理论的分割方法等。
随着影像医学的发展,图像分割在医学应用中具有了越来越特殊重要的意义。医学图像,如磁共振图像(MRI)通常是复杂的和不清楚的,需要准确无误地分割图像中的对象。随着新的理论和方法的频繁出现,图像分割领域已经提出了许多种方法,特别是一些基于支持向量机(SVM)的改进方法。
准确的图像分割是医学图像分析的一项重要任务,因为它可以帮助我们更好地区分肿瘤和正常组织。脑胶质瘤(一种脑肿瘤)的MRI图像的重要特征之一是肿瘤形状最常出现不规则,轮廓不明显。因此,轮廓上的节点集不易建立并聚集在一起。为了对节点集进行聚类分割,对神经胶质瘤图像进行分割,为了达到更高的分割精度,许多研究者一方面试图通过图像中的特征对节点集进行聚类,另一方面改进了SVM中的最大边缘参数。在节点集的聚类中,当前较为活跃的研究方向是使用马尔可夫随机场(Markov Random Field,MRF)通过先验特征对节点集进行聚类,通过先验特征支持向量机—马尔可夫随机场(pSVMRF)进行图像分割;也存在这样一种思路:通过对t-统计量和特征强度进行排序来训练每个集合,以找到聚类大脑图像的阈值。这些方法的使用都在一定程度上帮助了聚类集合,其目的旨在找到可决定最终聚类或图像分割准确性的阈值或平衡。
事实上,SVM中的最大值或阈值是实现高精度分割的核心和重要依据,因此寻找合理的SVM阈值或平衡成为图像分割当前研究重点。在这一研究过程中,一个又一个的难题相继被攻克。最初,研究者使用Arimoto熵作为阈值选择训练样本,提出了一种改进的最小二乘支持向量机(LS-SVM),并考虑了最大局部能量,最大梯度,为训练后的LS-SVM生成最大二次矩阵进行图像分割。为了解决传统的支持向量机不适合对不平衡任务进行聚类这一问题,加权余量准则应运而生,这一准则解决了不平衡学习任务,使得少数类比不重要类更加聚类。同时,通过改变拉普拉斯公式中的损失函数,避免了支持向量机中最大边缘集群的早熟收敛。最终我们可以得出结论:SVM方法为聚类或图像分割提供了重要的机会。图像分割的基本和重要部分是节点集聚类和提高SVM的最大容限。
旨在为所有元素找到最佳解决方案的纳什均衡实现了一个元素的最大值与所有元素的最大值之间的平衡。因此,研究图像分割的研究者们开始尝试将纳什均衡这一理论引入图像分割中。这一理论的引进使得研究者可以通过节点位置计算主要集合之间的最大对应关系,以找到重叠对象区域。纳什均衡在研究节点灰度分布的随机广义聚类上也起到了重要的作用。
然而,迄今为止的这些基于纳什均衡的图像分割研究都只是致力于具有节点灰度分布和位置的纳什均衡方法,其在不确定和不明确的聚类上比现有方法更好。但是,节点簇不仅由其位置决定,还由集合中的灰度分布确定。因此,本发明提出一种基于纳什均衡的新方法,同时考虑灰度分布和节点位置,并且构建了对应的新模型。
发明内容
本发明的目的在于提供一种基于新型纳什均衡的医学图像分割方法。本发明提出两大改进点:一个是使用具有熵和标准偏差的双重分配约束的改进的纳什均衡来分析聚类,另一个是基于所提出的纳什均衡来修改SVM,以便实现更好的最大边缘,从而提高医学图像分割的准确性。
本发明的目的是这样实现的,具体包括以下步骤:
(1)输入需要分割的目标图像,使用待分割图像中每个节点的坐标和灰度值两个特征值构建近似熵和标准差双重约束下的新型纳什均衡;
(2)根据每个节点构建出的新型纳什均衡确定该节点所属的节点集;
(3)使用提出的新型纳什均衡修改传统SVM;
(4)对于已确定的节点集合中的每个节点,使用修改后的SVM确定所有节点的决策函数,根据决策函数对待分割图像进行分割。
步骤(1)中构建近似熵和标准差双重约束下的新型纳什均衡具体方法如下:
(1.1)纳什均衡中的一个约束是熵,指示图像中灰度分布的聚集特征,使用简单的近似冯诺依曼熵来近似双重约束中的熵;熵的概念与香农信息理论的形成密切相关,此理论中负熵为:
其中Ep代表熵,k表示图像内的不同集合,Pi表示节点i的灰度发生概率;将归一化拉普拉斯熵谱相关的图的冯诺依曼熵定义为:
其中Ep代表熵,k表示图像内的不同集合,Pi表示节点i的灰度发生概率;使用简单的近似冯诺依曼熵来近似于熵,熵被二次熵替代为:则
熵Ep越大,区域内节点灰度值越均匀,反之亦然;k是集合P中的节点号;在医学图像中熵代表节点的灰度密度熵,因此熵Ep中的Pi为:
其中WS(i)表示集合中节点i的灰度值,WS表示集合中节点的灰度值的总值,因此图像中的近似熵为:
(1.2)约束条件标准差使用双重分配约束ap和bp表示节点位置密度值,则
其中N是集合P中的节点数量,M是集合Q中的节点数量;图像中的ap和bp定义为节点位置的特征,从图像中提取的特征是标准差;标准差是图像中节点与其邻域之间强度的位置测量,对于二维平面中的节点,标准差ap和bp为:
其中为计算的节点对应的节点集里所有节点的x,y坐标的平均值;
(1.3)近似熵和标准差双重约束下的新型纳什均衡模型为:
根据对熵和标准差的研究,对上式进行扩展为:
步骤(2)中确定该节点所属的节点集具体方法如下:
(2.1)输入需要分割的目标图像,获取并输入待分割目标图像所有节点位置(xi,yi)i=1,2,...,I;
(2.2)将输入的所有节点初始化为两个集合,即将所有节点随机分在不同的两个集合,一个命名为该对象的集合,记为set1,另一个命名为背景集合,记为set2;
(2.3)从对象节点集set1中选择一个节点node1(xi,yi);
(2.4)计算节点node1(xi,yi)属于set1时,近似熵和标准差双重约束下两个集合set1和set2的纳什均衡为
(2.5)将node1(xi,yi)从set1中剔除,放入set2中,计算此时近似熵和标准差双重约束下两个集合的纳什均衡
(2.6)比较步骤(2.4)和步骤(2.5)中所计算的纳什均衡的值,则node1(xi,yi)属于纳什均衡的值最大的集合;
(2.7)返回到步骤(2.3),直到计算出set1中初始化时的所有节点所属集合;
(2.8)返回步骤(2.3),对集合set2中初始化时的所有节点进行相似的处理,确定set2中初始化时的所有节点所属集合;
(2.9)输出计算完的节点集。
步骤(3)具体实现方法如下:
(3.1)通过新型纳什均衡中节点灰度和位置分布双重约束得到SVM中的修改惩罚参数ξi=apEp或ξi=bqEq;
(3.2)使用ci代表纳什均衡中ap和bp的最大值,将ξi定义为:
ξi=ciEi=maximum(apEp,bqEq)
(3.3)将修改后的SVM余量定义为:
限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I。
步骤(4)具体实现方法如下:
(4.1)从样本节点集中选取一点(xi,yi),对该点进行如下处理:
给定核函数k(xi,x)和惩罚参数ciEi,其中ci是该点的纳什均衡ap和bp的最大值,ciEi=maximum(apEp,bqEq),且样本函数的区间为限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I
(4.2)计算所选节点的步骤(4.1)中的参数w和参数b:
w=a1y1x1+a2y2x2+…+anynxn
b≥1-ciEi-yi(wxi)(i=1,2,…,I)
(4.3)计算所选节点的决策函数g(x):
(4.4)返回到步骤(4.1)并再次运行,直到计算出所有节点的决策函数;根据决策函数分割待分割图像。
本发明的有益效果在于:
(1)本发明在纳什均衡簇中引入双重分配约束,解决了决定节点簇的轮廓不明确的问题;
(2)本发明基于所提出的纳什均衡来修改SVM,实现了更好的最大边缘,从而提高医学图像分割的准确性。
说明书附图
图1为基于新型纳什均衡的图像分割流程图;
图2为新型纳什均衡流程图;
图3为修改的SVM分割目标图像流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
本发明提出了熵和标准偏差双重约束下的新型纳什均衡,并根据提出的纳什均衡来修改SVM,以便实现更好的最大边缘,从而提高医学图像分割的准确性。基于新型纳什均衡的图像分割流程如图1所示。本发明对目标图像进行分割主要经过两步,第一步:新型纳什均衡聚类节点;第二点:使用基于新型纳什均衡的SVM分割目标图像。新型纳什均衡聚类节点流程如图2所示;修改的SVM分割目标图像流程图如图3所示。
本发明公开了一种基于新型纳什均衡的医学图像分割模型与方法,该方法包括以下步骤:获取待分割图像的所有节点坐标和灰度值,初始化所有节点为两个节点集合:对象集合和背景集合,按构建的近似熵和标准差双重约束下的纳什均衡确定每个节点实际所属集合,对于确定的两个节点集合给出每个节点的惩罚参数和核函数,使用新型纳什均衡修改后的SVM确定每个节点的决策函数,根据所有节点的决策函数对待分割图像进行分割。本发明提出两大改进点:一个是使用具有熵和标准偏差的双重分配约束的改进的纳什均衡来分析聚类,另一个是基于所提出的纳什均衡来修改SVM,以便实现更好的最大边缘,从而提高医学图像分割的准确性。
第一个改进点:对于纳什均衡的改进。在Nash均衡簇中引入双重分配约束来解决决定节点簇的轮廓不明确的问题。纳什均衡中的一个约束是熵,指示图像中灰度分布的聚集特征,本文用Ep或Eq来表示;另一个是标准偏差,表示节点在聚类中的位置,用ap或bq表示。用w代表节点灰度值。考虑双重分配约束apEp或bqEq在纳什均衡聚类中的作用,提出双重约束下的新型纳什均衡模型:
当把Ep和ap放在一起,并且Eq和bq放在一起时,Ep和Eq有相同的意义,因此Eq可以被类似的定义为Ep。从该模型可以得出:对于想要获得的最佳聚类,应该寻找不同集合中节点的最大值。对于约束条件之一的熵进行了研究(以EP代表熵,Eq与EP相同,只是表示不同集合)。为追求更高效率,使用简单的近似冯诺依曼熵来近似于熵,熵被二次熵替代
熵EP越大,区域内节点灰度值越均匀,反之亦然。所以,根据熵值i,可以推断出哪个区域或某个节点可能属于哪个区域。k是集合P中的节点号。对于另一约束条件标准差,使用双重分配约束ap和bq表示节点位置密度值。
上面的模型包括双重分配,一个是考虑节点灰度值的密度和度量的近似熵,另一个是标准偏差,它约束了节点的位置密度值的限制。也就是说,本发明中纳什均衡的公式不仅包括节点的灰度值,还包括节点的位置密度值。
本发明中的纳什均衡寻求位置密度值和灰度值的最大值。因此,本发明提出了一种基于熵和标准差的双分配约束纳什均衡的聚类模型和方法,旨在寻求最大的平衡。这意味着它不仅在一个集合中最多获得一个节点,而且还获得所有集合中所有节点的所有最大值。这也意味着最大值不是绝对最大值,而是相对最大值或平衡最大值。如果没有新的最大值出现,所获得的余额就是纳什均衡。
第二个改进点:改进SVM用于图像分割。本发明通过新型纳什均衡中节点灰度和位置分布双重约束得到SVM中的修改惩罚参数。
用于图像分割的SVM方法旨在恢复由超平面分离的最大余量。余量越大,SVM中修改的惩罚参数区间就越小。
2.1传统SVM公式
超平面函数g(x)可以表述为:
g(x)的限制条件为:yi[(wxi+b)]≥1,(i=1,2,...I)
其中,i是样本数量,1表示样本函数的区间为1。样本函数的区间值越小,分类器的边际越大,准确度越高。有一个参数ξi,是一个惩罚参数,使间隔按下式变小:
yi=[(wxi+b)]≥1-ξi,i=1,2,...I
样本函数的间隔可以表述为:
样本函数限制条件为:yi(wxi)+b≥1-ξi,(i=1,2,...I)
2.2纳什均衡改进的SVM
每个样本都有相应的ξi值,用apEp或者bqEq代替ξi,即:
ξi=apEp orξi=bqEq
在SVM中,使用ci代表纳什均衡中ap和bq的最大值,故ξi可定义为:
ξi=ciEi=maximum(apEp,bqEq)
因此,修改后的SVM余量可以定义为:
限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
本发明采用如下步骤:
步骤001.输入需要分割的目标图像,获取并输入待分割目标图像所有节点位置(xi,yi),i=1,2,...,I。
步骤002.将输入的所有节点初始化为两个集合,即将所有节点随机分在不同的两个集合,一个命名为该对象的集合,记为set1,另一个命名为背景集合,记为set2。
步骤003.从对象节点集set1中选择一个节点node1(xi,yi)。
步骤004.计算节点node1(xi,yi)属于set1时,近似熵和标准差双重约束下两个集合set1和set2的纳什均衡
步骤005.将node1(xi,yi)从set1中剔除,放入set2中,计算此时近似熵和标准差双重约束下两个集合的纳什均衡
步骤006.比较步骤004和步骤005中所计算出的纳什均衡的值,node1(xi,yi)属于纳什均衡的值最大的集合。
步骤007.返回到步骤003,直到计算出set1中初始化时的所有节点所属集合。
步骤008.返回步骤003,对集合set2中初始化时的所有节点进行相似的处理,确定set2中初始化时的所有节点所属集合。
步骤009.输出计算完的节点集,作为使用修改后的SVM分割目标图像时的样本节点集。
步骤010.步骤009中得到的样本节点集记为{(xi,xj),i,j=1,2,...,I},其中xi,xj∈Rn。
步骤011.从样本节点集中选取一点(xi,yi),将该点进行如下处理:给定核函数k(xi,x)和惩罚参数ciEi,其中ci是该点的纳什均衡ap和bp的最大值,ciEi=maximum(apEp,bqEq),且样本函数的区间为:限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I
步骤012.计算所选节点在步骤011计算时所用等式描述的w的参数:
w=a1y1x1+a2y2x2+…+anynxn
步骤013.计算所选节点在步骤011计算时所用等式描述的b参数:
b≥1-ciEi-yi(wxi) (i=1,2,…,I)
步骤014.计算所选节点的决策函数g(x)(即超平面函数):
步骤015.返回到步骤011并再次运行,直到计算出所有节点的决策函数,然后退出。
步骤016.根据上述所有节点的决策函数,对应步骤009中根据新型纳什均衡所分得的两个节点集(步骤002中目标节点集和背景节点集)各个节点,对目标图像进行分割,对比传统方式,发现分割更为精确,证明本发明达到了最终目标。
本发明中用到的纳什均衡寻求的是节点灰度值和位置密度值的最大值,是在近似熵和标准差双重约束下的新型纳什均衡,旨在寻求最大的平衡。对于步骤004和步骤005中计算双重约束下的纳什均衡,具体的步骤如下(以步骤004为例,步骤005类似):
步骤00401.计算两个节点集的新型纳什均衡中的标准差ap(属于对象节点集)和bi(属于背景节点集):
其中,为计算的节点对应的节点集里所有节点的x,y坐标的平均值
步骤00402.计算两个节点集的新型纳什均衡中的近似熵ES1,ES2,已知式(6)为求节点集P近似熵公式,其中Pi可由式(5)得到,式(5)中WS(i)表示集合中节点i的灰度值,WS表示集合中节点的灰度值的总值。结合式(5)、式(6)可计算得到节点集set1、set2的近似熵。
步骤00403.将计算所得近似熵和标准差代入纳什均衡模型进行比较即可。
对一种基于新型纳什均衡的医学图像分割模型与方法,也可进行如下步骤描述:
(1)使用待分割图像中每个节点的两个特征值:坐标和灰度值构建近似熵和标准差双重约束下的新型纳什均衡,根据每个节点构建出的新型纳什均衡确定该节点所属的节点集。
使用简单的近似冯诺依曼熵来近似与双重约束中的熵,将归一化拉普拉斯熵谱相关的图的冯诺依曼熵定义为:其中,k表示图像内的不同集合,并且Pi是节点i的灰度发生概率。即此时的熵熵EP越大,区域内节点灰度值越均匀,反之亦然。所以,根据熵值i,可以推断出哪个区域或某个节点可能属于哪个区域。k是集合P中的节点号。
对于另一约束条件标准差,使用双重分配约束ap和bq表示节点位置密度值,即图像中的ap和bq可以定义为节点位置的特征,即从图像中提取的特征是标准差。标准偏差是图像中节点与其邻域之间强度的位置测量。对于二维平面中的节点,将其标准偏差定义为:
根据对熵和标准差的研究,提出近似熵和标准偏差双重约束下的新型纳什均衡模型:
(2)使用提出的新型纳什均衡修改传统SVM。
步骤(2)中所述的传统SVM公式为:超平面函数g(x)的限制条件为yi[(wxi+b)]≥1(i=1,2,...I),其中i是样本数量,1表示样本函数的区间为1,样本函数的区间值越小,分类器的边际越大,准确度越高;参数ξi是惩罚参数,使间隔按yi=[(wxi+b)]≥1-ξi i=1,2,...I变小,则样本函数的间隔为样本函数限制条件为yi(wxi)+b≥1-ξi i=1,2,...I;
(3)对于已确定的节点集合中每个节点,使用修改后的SVM确定所有节点的决策函数,根据决策函数对待分割图像进行分割。
Claims (5)
1.一种基于新型纳什均衡的医学图像分割方法,其特征在于,包括:
(1)输入需要分割的目标图像,使用待分割图像中每个节点的坐标和灰度值两个特征值构建近似熵和标准差双重约束下的新型纳什均衡;
(2)根据每个节点构建出的新型纳什均衡确定该节点所属的节点集;
(3)使用提出的新型纳什均衡修改传统SVM;
(4)对于已确定的节点集合中的每个节点,使用修改后的SVM确定所有节点的决策函数,根据决策函数对待分割图像进行分割。
2.根据权利要求1所述的一种基于新型纳什均衡的医学图像分割方法,其特征在于,所述的步骤(1)中构建近似熵和标准差双重约束下的新型纳什均衡具体方法如下:
(1.1)纳什均衡中的一个约束是熵,指示图像中灰度分布的聚集特征,熵的概念与香农信息理论的形成密切相关,此理论中负熵为:
其中Ep代表熵,k表示图像内的不同集合,Pi表示节点i的灰度发生概率;将归一化拉普拉斯熵谱相关的图的冯诺依曼熵定义为:
其中Ep代表熵,k表示图像内的不同集合,Pi表示节点i的灰度发生概率;使用简单的近似冯诺依曼熵来近似于熵,熵被二次熵替代为:则
熵Ep越大,区域内节点灰度值越均匀,反之亦然;k是集合P中的节点号;在医学图像中熵代表节点的灰度密度熵,因此熵Ep中的Pi为:
其中WS(i)表示集合中节点i的灰度值,WS表示集合中节点的灰度值的总值,因此图像中的近似熵为:
(1.2)约束条件标准差使用双重分配约束ap和bp表示节点位置密度值,则
其中N是集合P中的节点数量,M是集合Q中的节点数量;图像中的ap和bp定义为节点位置的特征,从图像中提取的特征是标准差;标准差是图像中节点与其邻域之间强度的位置测量,对于二维平面中的节点,标准差ap和bp为:
其中为计算的节点对应的节点集里所有节点的x,y坐标的平均值;
(1.3)近似熵和标准差双重约束下的新型纳什均衡模型为:
根据对熵和标准差的研究,对上式进行扩展为:
3.根据权利要求1所述的一种基于新型纳什均衡的医学图像分割方法,其特征在于:所述的步骤(2)中确定该节点所属的节点集具体方法如下:
(2.1)输入需要分割的目标图像,获取并输入待分割目标图像所有节点位置(xi,yi)i=1,2,...,I;
(2.2)将输入的所有节点初始化为两个集合,将所有节点随机分在不同的两个集合,一个命名为该对象的集合,记为set1,另一个命名为背景集合,记为set2;
(2.3)从对象节点集set1中选择一个节点node1(xi,yi);
(2.4)计算节点node1(xi,yi)属于set1时,近似熵和标准差双重约束下两个集合set1和set2的纳什均衡为
(2.5)将node1(xi,yi)从set1中剔除,放入set2中,计算此时近似熵和标准差双重约束下两个集合的纳什均衡
(2.6)比较步骤(2.4)和步骤(2.5)中所计算的纳什均衡的值,则node1(xi,yi)属于纳什均衡的值最大的集合;
(2.7)返回到步骤(2.3),直到计算出set1中初始化时的所有节点所属集合;
(2.8)返回步骤(2.3),对集合set2中初始化时的所有节点进行相似的处理,确定set2中初始化时的所有节点所属集合;
(2.9)输出计算完的节点集。
4.根据权利要求1所述的一种基于新型纳什均衡的医学图像分割方法,其特征在于:所述的步骤(3)的具体实现方法如下:
(3.1)通过新型纳什均衡中节点灰度和位置分布双重约束得到SVM中的修改惩罚参数ξi=apEp或ξi=bqEq;
(3.2)使用ci代表纳什均衡中ap和bp的最大值,将ξi定义为:
ξi=ciEi=maximum(apEp,bqEq);
(3.3)将修改后的SVM余量定义为:
限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I。
5.根据权利要求1所述的一种基于新型纳什均衡的医学图像分割方法,其特征在于:所述的步骤(4)的具体实现方法如下:
(4.1)从样本节点集中选取一点(xi,yi),对该点进行如下处理:
给定核函数k(xi,x)和惩罚参数ciEi,其中ci是该点的纳什均衡ap和bp的最大值,ciEi=maximum(apEp,bqEq),且样本函数的区间为限制条件:yi(wxi)+b≥1-ciEi,i=1,2,...I;
(4.2)计算所选节点的步骤(4.1)中的参数w和参数b:
(4.3)计算所选节点的决策函数g(x):
(4.4)返回到步骤(4.1)并再次运行,直到计算出所有节点的决策函数;根据决策函数分割待分割图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810870190.0A CN109285159A (zh) | 2018-08-02 | 2018-08-02 | 一种基于新型纳什均衡的医学图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810870190.0A CN109285159A (zh) | 2018-08-02 | 2018-08-02 | 一种基于新型纳什均衡的医学图像分割方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109285159A true CN109285159A (zh) | 2019-01-29 |
Family
ID=65182937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810870190.0A Pending CN109285159A (zh) | 2018-08-02 | 2018-08-02 | 一种基于新型纳什均衡的医学图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109285159A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766624A (zh) * | 2019-10-14 | 2020-02-07 | 中国科学院光电技术研究所 | 一种基于迭代修复的点目标和暗斑图像背景均衡方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105261017A (zh) * | 2015-10-14 | 2016-01-20 | 长春工业大学 | 基于路面约束的图像分割法提取行人感兴趣区域的方法 |
-
2018
- 2018-08-02 CN CN201810870190.0A patent/CN109285159A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105261017A (zh) * | 2015-10-14 | 2016-01-20 | 长春工业大学 | 基于路面约束的图像分割法提取行人感兴趣区域的方法 |
Non-Patent Citations (1)
Title |
---|
TIAN CHI ZHANG等: ""A Novel Model and Method Based on Nash Equilibrium for Medical Image Segmentation"", 《JOURNAL OF MEDICAL IMAGING AND HEALTH INFORMATICS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766624A (zh) * | 2019-10-14 | 2020-02-07 | 中国科学院光电技术研究所 | 一种基于迭代修复的点目标和暗斑图像背景均衡方法 |
CN110766624B (zh) * | 2019-10-14 | 2022-08-23 | 中国科学院光电技术研究所 | 一种基于迭代修复的点目标和暗斑图像背景均衡方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Machairas et al. | Waterpixels | |
Zhang et al. | Abnormal breast identification by nine-layer convolutional neural network with parametric rectified linear unit and rank-based stochastic pooling | |
Zhou et al. | Fuzziness parameter selection in fuzzy c-means: the perspective of cluster validation | |
CN108629783B (zh) | 基于图像特征密度峰值搜索的图像分割方法、系统及介质 | |
Hooda et al. | Brain tumor segmentation: A performance analysis using K-Means, Fuzzy C-Means and Region growing algorithm | |
Ashraf et al. | A multichannel Markov random field framework for tumor segmentation with an application to classification of gene expression-based breast cancer recurrence risk | |
Veenman et al. | A cellular coevolutionary algorithm for image segmentation | |
CN110188763B (zh) | 一种基于改进图模型的图像显著性检测方法 | |
CN110276764A (zh) | 基于K值估计的K-Means水下图像背景分割改进算法 | |
Kapse et al. | Literature survey on detection of brain tumor from MRI images | |
Yang et al. | Color texture segmentation based on image pixel classification | |
CN110176021A (zh) | 结合亮度校正的显著性信息的水平集图像分割方法及系统 | |
Devi et al. | Performance comparison of various clustering techniques for diagnosis of breast cancer | |
CN110473224B (zh) | 一种基于kl熵的rsf水平集图像自动分割方法 | |
CN106022359A (zh) | 基于有序信息熵的模糊熵空间聚类分析方法 | |
Song et al. | Automated segmentation of overlapping cytoplasm in cervical smear images via contour fragments | |
Xiang et al. | Interactive natural image segmentation via spline regression | |
CN109285159A (zh) | 一种基于新型纳什均衡的医学图像分割方法 | |
Afshar et al. | Lung tumor area recognition in CT images based on Gustafson-Kessel clustering | |
CN113052242B (zh) | 图像处理网络的训练方法及装置、图像处理方法及装置 | |
Gu et al. | A novel method for breast mass segmentation: from superpixel to subpixel segmentation | |
Li et al. | An Improved Image Segmentation Algorithm Based on GPU Parallel Computing. | |
CN111161282B (zh) | 基于深度种子的图像多层次分割的目标尺度选择方法 | |
CN110969639B (zh) | 一种基于lfmvo优化算法的图像分割方法 | |
CN107492101B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190129 |