[发明专利]血管超声图像中颈动脉血管内中膜的计算机自动分割和厚度均匀度分析方法

权利要求书

1.一种血管超声图像中颈动脉血管内中膜的计算机自动分割和厚度均匀度分析方法，其特征在于，包括以下步骤：

(1)在计算机中人工载入一幅采集的原始血管超声图像；

(2)人工选择包含血管远侧壁颈动脉内中膜的一个矩形感兴趣区域；

(3)在矩形感兴趣区域中检测演化得到颈动脉血管内中膜的LII和MAI两条边界的精确轮廓线：

(3.1)在矩形感兴趣区域中检测得到颈动脉血管内中膜的LII和MAI两条边界的初始轮廓线；

(3.2)将步骤(3.1)检测得到的LII和MAI两条边界的初始轮廓线代入下述双Snake模型推导得到的差分方程进行反复迭代直至稳定，得到LII的精确轮廓线Γ₁(q)＝(x₁(q)，y₁(q))和MAI的精确轮廓线Γ₂(q)＝(x₂(q)，y₂(q))，所述的差分方程的推导过程如下：

定义双Snake模型的能量泛函为：

E_d-snake＝E_int+E_ext                      (1)

Eint=12Σi=12∫01wi1|∂Γi∂q|2+wi2|∂2Γi∂q2|2dq---(2)]]>

Eext=12α∫0t|∂(Γ1-Γ2)∂q|2dq+(1-α)Σi=12∫01-f(Γi)dq---(3)]]>

其中，Γ₁(q)＝(x₁(q)，y₁(q))表示LII的轮廓线，Γ₂(q)＝(x₂(q)，y₂(q))表示MAI的轮廓线，q表示曲线的参数，q∈[0，1]，q为0表示轮廓线的左端点，q为1表示轮廓线的右端点，q为中间值则表示轮廓线中间的一点，w_i1，i＝1，2分别是LII轮廓线的张力系数和MAI轮廓线的张力系数，w_i1，i＝1，2∈[0，1]，w_i2，i＝1，2分别是LII轮廓线的硬度系数和MAI轮廓线的硬度系数，w_i2，i＝1，2∈[0，1]，α∈[0，1]，使E_d-snake达到极小，根据梯度下降流可以得到：

∂xi∂t=wi1∂2xi∂q2-wi2∂4xi∂q4+α(∂2xi∂q2-∂2x3-i∂q2)+(1-α)∂f∂x(xi,yi)---(4)]]>

∂yi∂t=wi1∂2yi∂q2-wi2∂4yi∂q4+α(∂2yi∂q2-∂2y3-i∂q2)+(1-α)∂f∂y(xi,yi)---(5)]]>

其中，i＝1，2，对(4)和(5)中的4个偏微分方程在空间上使用中心差分，在时间上使用向前差分，得到差分方程；

(4)厚度均匀度分析计算：根据步骤(3)分割得到的颈动脉血管内中膜LII和MAI两条边界的精确轮廓线，计算厚度均匀度。