[发明专利]基于形状先验和相似性约束的颈动脉内外壁自动分割方法

摘要

一种基于形状先验和相似性约束的颈动脉内外壁自动分割方法，所述分割方法包括如下步骤1)通过SVM方法检测颈动脉内腔，得到内腔所在位置；2)使用各向异性的分割方法分割内壁；3)在分割内壁的基础上，结合外壁的特征学习回归模型，估算外壁大致位置，并构造概率模型；4)使用相似性约束算法分割外壁。本发明提供了一种自动方式、有效实现内外壁分割、精确度较高的基于形状先验和相似性约束的颈动脉内外壁自动分割方法。

主权项

一种基于形状先验和相似性约束的颈动脉内外壁自动分割方法，其特征在于：所述分割方法包括如下步骤：1)通过SVM方法检测颈动脉内腔，内腔所在区域，记为S0区域；2)使用各向异性的分割方法分割内壁，过程如下：把步骤1)检测得到的S0区域作为分割内壁的先验；首先，在S0区域内提取前景种子，在S0区域外ROI区域内提取背景种子；计算S0的最小直径dmin，然后用直径为1/2dmin的菱形形状腐蚀S0直到收缩至区域里面的像素点被用作前景种子；使用同样的形状因子把S0膨胀至把在区域外但在ROI内的像素点作为背景种子；接下来，把前景和背景种子的灰度值作为特征来构建高斯混合模型，分别对应和通过下面的后验概率公式来计算：PLIB(yp=f|C)=P(fp|GfgLIB)P(fp|GfgLIB)+P(fp|GbgLIB)---(2)]]>PLIB(yp＝b|C)＝1‑PLIB(yp＝f|C)   (3)其中，fp表示像素p的特征向量；PLIB(yp＝f|C)和PLIB(yp＝b|C)分别表示像素p属于前景和背景的后验概率；接着，使用各向异性的分割算法来分割内壁，能量函数如下：ELIB(y|C)＝∑p∈yDLIB(yp)+∑(p，q)∈NVLIB(yp，yq)+∑(p，q)∈NSLIB(yp，yq)   (4)其中，DLIB(yp)是数据项被定义为DLIB(yp)＝‑log PLIB(yp|C)，VLIB(yp，yq)是平滑项被定义为VLIB(yp，yq)＝exp(‑(Ip‑Iq)2/2σ2)；(p，q)∈N表示像素p和q是在一个邻域内，最后SLIB(yp，yp)是一个约束项确保中心点C到像素p的路径上的点q和p拥有同样的标签，表示为：SLIB(yp,yq)=0if(p,q)∈S*({c∈C})∞otherwise---(5)]]>其中，像素p和q属于不同的邻域M，q是p和中心点集的最近测地距离的路径上的相邻像素；最后，最小化能量函数(4)，得到精确的内壁分割结果；3)在分割内壁的基础上，结合外壁的特征学习回归模型，估算外壁大致位置，并构造概率模型；首先，学习一个回归模型，将内壁的最大直径，平均直径、最小直径和内壁轮廓的力矩作为回归模型的特征，第(u+v)个力矩的计算公式为其中，i和j分别表示内壁轮廓中的某个像素的行和列；通过这个回归模型估计颈动脉管壁的最大最小厚度，回归模型用如下公式表示：F是特征向量t＝[tmax，tmin]和W是从训练样本中得到的参数集，tmax是测量得到的最大管壁厚度，同理，tmin是测量得到的最小管壁厚度，是S曲线函数；通过回归模型估算出最大和最小外壁轮廓后，分别用表示，构造前景和背景两个高斯混合模型；介于内壁和之间的像素被看作是前景，在之外和内的像素被看作是背景；此外，还构造了对应训练样本的高斯混合模型，结合两部分的高斯混合模型，通过以下公式计算一个像素点p属于前景还是背景的概率：PAB(yp=f|C)=αPAB1(yp=f|C)+(1-α)PAB2(yp=f|C)=αP(fp|GfgAB1)P(fp|GfgAB1)+P(fp|GbgAB1)+(1-α)P(fp|GfgAB2)P(fp|GfgAB2)+P(fp|GbgAB2)---(7)]]>PAB(yp＝b|C)＝α(1‑PAB1(yp＝f|C|))+(1‑α)(1‑PAB2(yp＝f||C))   (8)其中，和i∈{1，2}，分别表示第i部分的前景和背景高斯混合模型，和i∈{1，2}，分别表示第i部分高斯混合模型对应的前景或背景概率，α表示第一部分所占的权重，通过公式(7)(8)计算得到的概率的负的对数即数据项DAB；4)使用相似性约束算法分割外壁，过程如下；以分割得到的内壁Cl为基础，以每次一个像素的频率逐渐往外扩张；当第k次扩张时，如果Ck曲线上的像素点都被标记为前景，那么之前扩展的曲线C0到Ck‑1的轮廓都看作是前景；对于第k次曲线上像素p，计算p的8邻域内的点与Cl的距离，并分为两种情况，第一种情况：邻域点与Cl的距离不大于p到Cl的距离，第二种情况：邻域点与Cl的距离大于p到Cl的距离；把符合第一种情况的邻域点标为前景区域，所以，如果当前曲线上某一点属于前景，那么整个曲线也属于前景；相似性约束算法采用以下公式定义：SAB(yp,yq)=0ifdq≤dp∞otherwise]]>其中，dp和dq表示像素p和q到Cl的距离，dq是指像素q和Cl所有点的最小欧氏距离；得到的SAB和常规平滑项VAB，还有步骤3)中的得到的数据项DAB共同构成能量函数，能量函数最小时就得到了外壁的分割结果。