软组织核磁图像自适应外力水平集自动分割及实现方法

摘要

一种软组织核磁图像自适应外力水平集自动分割及实现方法：在交互式操作情况下选取符合灰度阈值的目标区域组织；计算平滑后图像的梯度信息；设置曲线演化引导函数；设置边界停止函数；如果演化曲线越过目标区域边界，则对边界停止函数参数进行调整，然后返回曲线演化，重新迭代，直到演化曲线收敛到目标区域边界为止。本发明具有实时性强、运算效率高、能对多个离散区域同时进行分割、能够准确识别人体软组织的模糊边界等特点，且分割精度高、图像细节特征清晰、智能化程度高、无需人工干预、运行稳定可靠。

主权项

一种软组织核磁图像的自适应外力水平集自动分割方法，其特征在于：包括如下步骤： 1)在交互式操作情况下，对人体软组织核磁共振序列二维医学图像按照图像序列和输入的灰度阈值进行初始轮廓框选，选取符合灰度阈值的目标区域组织； 2)对二维软组织核磁共振图像进行二维空间内的高斯滤波平滑，并计算平滑后图像的梯度信息，首先采用高斯核函数对原始图像进行平滑滤波，去掉图像背景中的高频噪声，然后计算平滑滤波后图像的梯度信息； 3)设置曲线演化引导函数，使演化曲线根据初始框选轮廓的形状自适应调整演化方向，所述的曲线演化引导函数是根据实际图像中物体边缘处两侧二阶微分符号相反，一阶微分的极值点对应于二阶微分中的零交叉点的特性而提出的，根据这一特性，定义水平集曲线演化引导函数： f_indicator(I)＝k·sign(I″)·(const+|I′|) 其中，sign(I″)表示原始图像I的二阶微分符号，(const+|I′|)用于加快曲线在目标边界附近的演化速度，常数项const≥1用于保证即使在绝对平滑的区域演化速度仍大于零，I′表示平滑滤波后的图像； 4)设置边界停止函数，开始曲线演化运算，如果演化曲线收敛到目标边界，则分割结束；否则返回步骤3)的曲线演化，继续进行迭代，如果演化曲线越过目标区域边界进入下一步骤； 5)如果演化曲线越过目标区域边界，则对边界停止函数参数进行调整，然后返回步骤3)曲线演化，重新迭代，直到演化曲线收敛到目标区域边界为止， 所述的对边界停止函数参数进行调整，是在快速变分水平集模型框架下，将曲线演化引导函数加权到演化外力项，并引入新的边界停止函数，得到一种新的自适应外力水平集自动分割模型，通过最小化下述能量泛函实现： ε(φ)＝α·g_ξ·Length(C)+β·f_indicator·g_ξ·S(C)+γ·P_SDF(φ) 其中φ为高一维的水平集函数，g_ξ表示边界停止函数； 利用变分法计算能量泛函ε相对于φ的一阶变分，则有： 若泛函ε取极值，则函数φ满足欧拉‑拉格朗日方程添加时间变量t，最小化能量泛函ε的梯度下降流为：其中，第一项为加权弧长项，保持迭代过程中曲线的平滑性；第二项为自适应外力项，根据演化曲线与目标边界的相互位置关系自动调整曲线的推进方向，保证演化曲线始终朝着 目标的边界逼近；第三项用于保持演化函数的符号距离函数特性；α,β,γ为各项的加权系数，δ(φ)为德拉克函数，用以检测水平集函数φ的零水平集。