一种基于先验局部统计与动态轮廓模型的人体器官重建方法与流程

技术特征：

1.一种基于先验局部统计与动态轮廓模型的人体器官重建方法，其特征在于，包括以下步骤：

(1)按下式初始化水平集函数φ：

$\mathrm{\φ}(x,t=0)=\left\{\begin{array}{cc}1,& x\∈{\mathrm{\Ω}}_1\\ -1,& x\∈{\mathrm{\Ω}}_2\end{array}\right.$

其中Ω₁是手工标注的轮廓内部区域,Ω₂是余下的轮廓外部区域；

(2)根据下式使用高斯滤波约束水平集函数：

φ＝φ*G_σ，k

其中G_σ，k是标准差为σ的高斯滤波函数，并且只考虑高斯中心处的部分以减少计算量；

(3)根据下式分别计算c₁和c₂：

$c_1=\frac{{\∫}_{\mathrm{\Ω}}H\left(\mathrm{\φ}\right)I\left(x\right)dx}{{\∫}_{\mathrm{\Ω}}H\left(\mathrm{\φ}\right)dx},$

$c_2=\frac{{\∫}_{\mathrm{\Ω}}(1-H(\mathrm{\φ}\left)\right)I\left(x\right)dx}{{\∫}_{\mathrm{\Ω}}(1-H(\mathrm{\φ}\left)\right)dx}$

其中c₁和c₂都是常数，它们分别是轮廓内部区域灰度的平均值和轮廓外部区域灰度的平均值；

(4)根据下式分别计算和：

$C_1^{Loc}\left(x\right)=\frac{K*H\left(\mathrm{\φ}\right)I\left(x\right)}{K*H\left(\mathrm{\φ}\right)}$

$C_2^{Loc}\left(x\right)=\frac{K*(1-H(\mathrm{\φ}\left)\right)I\left(x\right)}{K*(1-H(\mathrm{\φ}\left)\right)}$

其中，和是两个常量，是轮廓内部区域的局部权重灰度平局值，是轮廓外部区域的局部灰度加权平局值；

(5)根据下式进化水平集函数：

${\mathrm{\φ}}_t=\mathrm{\λ}-SPF\left(x\right)\·LIDF\left(x\right)|\▿\mathrm{\φ}|,x\∈\mathrm{\Ω}.;$

(6)如果φ(x)＞0，则设置φ(x)＝1，其余设置φ(x)＝-1；

(7)检查水平集函数的进化过程是否收敛，如果收敛则结束；否则返回步骤(2)。

2.根据权利要求1所述的基于先验局部统计与动态轮廓模型的人体器官重建方法，其特征在于，所述动态轮廓模型由λ-SPF函数和局部灰度差值力(LIDF)构成。