基于改进的距离调整水平集算法的医学图像分割方法技术

本发明专利技术公开了一种基于改进的距离调整水平集算法的医学图像分割方法，包括以下步骤：输入二维医学图像；利用Matlab软件采用高斯核函数对输入的二维医学图像进行二维空间内的高斯滤波平滑，平滑滤波后的图像的梯度信息为从而得到该二维医学图像的轮廓曲线，对DRLSE模型的隐函数进行改进得到ADRLSE模型，利用ADRLSE模型设置诺埃曼条件的边界停止函数；轮廓曲线的演化使轮廓曲线自适应调整演化方向，直至满足停止条件为止。本发明专利技术方法不仅保持了符号距离函数的属性，而且避免了重新初始化，减少了计算量，解决了符号函数的周期性初始化问题，提高了算法的效率，避免了图像分割中边缘失真。

【技术实现步骤摘要】

本专利技术涉及一种医学图像分割方法。特别是涉及一种改进的距离调整水平集方法。
技术介绍
图像分割是图像处理、图像分析和计算机视觉等领域经典的研究课题之一，也是难点之一，其目的是将图像分成各具特性的区域并将目标从背景和噪声中分离出来，从而为后续的定量、定性分析提供基础支持。由于各种图像，尤其是医学图像的方式各异，图像处理过程中又不可避免地引入一定的随机噪声，再加上图像中感兴趣目标的形态和强度差别很大，至今还没有一种通用的方法适合于所有图像的分割。近年来发展起来的水平集方法综合利用了区域与边界信息，以其精确性、自动性和最终分割结果的连续性被广泛应用于图像分割和计算机视觉等领域。该类方法易于结合分割物体的形状等先验知识约束，逐渐成为医学图像的分割领域的研究热点。水平集方法是欧拉法求解隐式偏微分方程的一种方法，水平集方法是最先由Osher and Sethian在上世纪80年代提出的，用于流体运动和火焰传播界面的追踪，随后又被引入到计算机视觉和图像处理领域。其主要思想是将被动曲线C(p,t)嵌入到高一维空间的水平集函数中，通过计算水平集函数零水平的位置来追踪曲线的演化。其中，p是任意参数化变量，t表示时间。在传统的水平集方法中，初始水平集函数通常取为符号距离函数。水平集函数在迭代的过程中可能发生退化，使它不再保持符号函数，因此要进行重新初始化操作，以保证水平集函数接近一个符号函数。在演化过程中周期性地对水平集函数进行校正，即重新初始化为符号函数，这一操作计算量非常大。Li等引入了符号距离函数来惩罚能量泛函，即无需初始化的水平集方法，解决了符号函数的周期性初始化问题，提高了算法的效率，但分割边界模糊或不连续的图像时往往出现边界泄露现象。2005年Li等引入距离调整函数的水平集方法(Advanced Distance Regularized Level Set Evolution)，在符号距离函数的基础上对该函数进行了调整。
技术实现思路
针对上述现有技术，本专利技术提供一种基于改进的距离调整水平集算法的医学图像分割方法，主要是在距离调整水平集算法(DRLSE，Distance Regularized Level Set Evolution Mode,)基础上提出了改进的距离调整水平集算法(ADRLSE，Advanced Distance Regularized Level Set Evolution)，在曲线演化即水平集函数演化过程中，不仅保持符号距离函数的属性，避免重新初始化，提供了更有效的计算，减少了经过DRLSE处理后的二维医学图像边缘锐化，利用本专利技术中建立的ADRLSE模型处理后的图像边缘更光滑。为了解决上述技术问题，本专利技术一种基于改进的距离调整水平集算法的医学图像分割方法予以实现的技术方案是，该方法包括以下步骤：步骤一、输入二维医学图像；步骤二、利用Matlab软件采用高斯核函数对输入的二维医学图像进行二维空间内的高斯滤波平滑，二维医学图像函数为I(x,y)；高斯核函数为其中σ为标准方差，平滑滤波后的图像为I′，I′＝Gσ*I(x,y)，平滑滤波后的图像的梯度信息为从而得到该二维医学图像的轮廓曲线；步骤三、对DRLSE模型的隐函数进行改进，使水平集函数接近于符号距离函数，改进后的隐函数为：p4(s)=0.001s,s<112(s-1)2,s≥1]]>式中，s为平面区域。基于改进后的隐函数得到ADRLSE模型，利用ADRLSE模型设置诺埃曼条件的边界停止函数；步骤四、开始轮廓曲线的演化，即使轮廓曲线自适应调整演化方向；步骤五、在轮廓曲线的演化过程中，当轮廓曲线满足诺埃曼条件的边界停止函数时，轮廓曲线收敛到目标边界，则停止演化，最终得到目标图像轮廓；否则，返回步骤四。进一步讲，步骤三中，利用ADRLSE模型设置诺埃曼条件的边界停止函数为：其中，Gσ*I表示用高斯核函数Gσ对图像I进行平滑滤波，表示对平滑后的图像求梯度。与现有技术相比，本专利技术的有益效果是：本专利技术中提出了一种改进的距离调整水平集算法，从而建立用于曲线演化的ADRLSE模型；利用ADRLSE模型进行医学图像分割可以避免图像分割中边缘失真，是目标图像轮廓更为平滑；利用ADRLSE模型相对于DRLSE模型可以大幅度提高计算效率。本专利技术方法主要应用临床影像工程、数字医学技术、放射治疗、手术规划、手术导航等领域中的图像处理。附图说明图1是本专利技术基于改进的距离调整水平集算法的医学图像分割方法流程图；图2是一合成葫芦原始图像；图3是本专利技术ADRLSE模型和现有技术DRLSE模型应用于合成图像分割实验效果图对比结果，其中,(a)是本专利技术ADRLSE模型23次迭代计算处理结果，(b)是DRLSE模型23次迭代计算处理结果；图4是图3中方框对应区域的局部放大图；图5是细胞显微原始图像；图6是细胞图像处理效果对比，其中，(a)是ADRLSE模型，(b)是DRLSE模型80次迭代计算处理结果，(c)是ADRLSE模型85次迭代计算处理结果，(d)是DRLSE模型85次迭代计算处理结果。具体实施方式下面结合附图和具体实施例对本专利技术技术方案作进一步详细描述。本专利技术基于改进的距离调整水平集算法的医学图像分割方法，在图像处理中，减少了图像边缘锐化，使其更光滑。在距离调整水平集方法的基础上，提出了改进的方法，根据引入改进的隐函数，不仅保持了符号距离函数的属性，而且避免了重新初始化，减少了计算量，解决了符号函数的周期性初始化问题，提高了算法的效率，避免了图像分割中边缘失真。如图1所示，本专利技术一种基于改进的距离调整水平集算法的医学图像分割方法，包括以下步骤：输入二维医学图像，高斯滤波平滑：利用Matlab软件采用高斯核函数对输入的二维医学图像进行二维空间内的高斯滤波平滑，以去掉图像背景中的高频噪声，二维医学图像函数为I(x,y)；高斯核函数为其中σ为标准方差，平滑滤波后的图像为I′，I′＝Gσ*I(x,y)，平滑滤波后的图像的梯度信息为从而得到该二维医学图像的轮廓曲线；建立ADRLSE模型：本专利技术中通过设置改进的距离调整水平集函数，使演化曲线能够自适应调整演化方向，对DRLSE模型的隐函数进行改进，改进的距离调整水平集函数使水平集函数接近于符号距离函数，不仅减少重新初始化过程，而且相比现有技术DRLSE模型和传统水平集方法实现了更有效的计算，具体说明如下：传统水平集方法中变化的本文档来自技高网...

【技术保护点】
一种基于改进的距离调整水平集算法的医学图像分割方法，其特征在于：包括以下步骤：步骤一、输入二维医学图像；步骤二、利用Matlab软件采用高斯核函数对输入的二维医学图像进行二维空间内的高斯滤波平滑，二维医学图像函数为I(x,y)；高斯核函数为其中σ为标准方差，平滑滤波后的图像为I′，I′＝Gσ*I(x,y)，平滑滤波后的图像的梯度信息为从而得到该二维医学图像的轮廓曲线；步骤三、对DRLSE模型的隐函数进行改进，使水平集函数接近于符号距离函数，改进后的隐函数为：p4(s)=0.001s,s<112(s-1)2,s≥1]]>式中，p为隐函数，s为平面区域；基于改进后的隐函数得到ADRLSE模型，利用ADRLSE模型设置诺埃曼条件的边界停止函数；步骤四、开始轮廓曲线的演化，即使轮廓曲线自适应调整演化方向；步骤五、在轮廓曲线的演化过程中，当轮廓曲线满足诺埃曼条件的边界停止函数时，轮廓曲线收敛到目标边界，则停止演化，最终得到目标图像轮廓；否则，返回步骤四。

【技术特征摘要】
1.一种基于改进的距离调整水平集算法的医学图像分割方法，其特征在于：包括以下
步骤：
步骤一、输入二维医学图像；
步骤二、利用Matlab软件采用高斯核函数对输入的二维医学图像进行二维空间内的高
斯滤波平滑，二维医学图像函数为I(x,y)；高斯核函数为其中
σ为标准方差，平滑滤波后的图像为I′，I′＝Gσ*I(x,y)，平滑滤波后的图像的梯度信息为
从而得到该二维医学图像的轮廓曲线；
步骤三、对DRLSE模型的隐函数进行改进，使水平集函数接近于符号距离函数，改进
后的隐函数为：
p4(s)=0.001s,s<112(s-1)2,s≥...

【专利技术属性】
技术研发人员：倪爱娟，田丰，谢新武，孙秋明，杨健，刘长军，杜振杰，赵杰，
申请(专利权)人：中国人民解放军军事医学科学院卫生装备研究所，中国人民解放军总后勤部东北军用物资采购局，
类型：发明
国别省市：天津;12