一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法与流程

本发明涉及一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法，具体涉及医学超声图像分割，属于图像处理
技术领域：
。
背景技术：
：超声图像作为医学影像中的一种，在医学临床诊断中有着重要作用。随着影像医学的发展，超声图像分割在医学应用中越来越具有重要意义。超声图像分割是计算机辅助诊断(computeraideddiagnosis,cad)系统中的一个重要组成部分，其作用是对超声图像进行分割，标定图像中的病灶部位，为对病灶部位的特征提取、诊断和分类提供必要条件。目前的研究中，所采用的方法大致可以分为：阈值法、聚类法还有活动轮廓。其中，活动轮廓方法是一种常见的图像分割方法，经常用于医学图像的分割，包括ct，mri和超声图像。活动轮廓(activecontourmodel，acm)是基于形变模型的分割技术。2001年chan和vese提出了chan-vese模型，简称cv模型。cv是一种几何形变活动轮廓模型，可以较好地分割具有较大噪声图像，经常用于超声图像分割。然而cv模型在处理同质的图像可以得到较好的分割结果，但在处理非均匀，异质的图像得到的边缘不太准确。为此，许多学者在这方面提出许多改进，有提出局域二元拟合能量驱动的隐式活动轮廓，lbf(localbinaryfitting)模型。也有提出局部的活动轮廓模型，如lcv模型。这些在处理非均匀的图像能取得一定的效果，但lbf和lcv模型对初始轮廓敏感，特别对于噪声大，同时又有其他伪影的超声图像，容易因为初始轮廓选取不恰当而导致分割效果不理想。技术实现要素：为克服上述现有技术的不足，本发明提供了一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法，本发明能自动获取病灶区域，解决靠人为选择初始轮廓的问题以及选择不恰当的问题。进行边界暗区域去除，克服了超声图像不均匀的特点，以及结合全局和局部的cv模型，达到更准确分割的效果。本发明的技术方案是：一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法，对图像进行预处理，首先利用双边滤波器进行去噪，用于解决图像噪声的存在，影响分割准确性的问题；接着对图像进行模糊聚类，图像二值化处理和形态学处理后，进行病灶区域选择，同时获得初始轮廓；再对边界暗区域去除，用于解决超声图像病灶区域外不均匀的特点，来提高分割的准确性；最后利用全局和局部cv活动轮廓模型的结合分割，得到分割结果。进一步地，所述方法的具体步骤如下：step1、图像滤波；采用双边滤波器对超声图像进行滤波，对超声图像存在的斑点噪声，进行滤波；step2、模糊聚类分割；对滤波后的图像进行模糊c均值聚类分割，聚类数目设定为3，得到将灰度值信息分成三类的图像；step3、图像二值化；选用一个阈值，这个阈值为根据病灶区域的灰度均值设定的，将其中包含病灶区域的一类与其他两类再分类，得到二值图像，以便后续形态处理；step4、形态学处理；得到的二值图像，先对二值图像进行膨胀处理，接着对图像取反，然后对边界的亮区域抑制，再对图像取反，再进行对图像进行腐蚀处理，得到多个连通区域的二值图像；step5、病灶区域选择；首先，统计连通区域的像素个数，选择最多像素的连通区域，接着比较其他的连通区域，如果其他区域的像素点大于最大连通区域的25％，则将视其为病灶区域roi，否则将其视为非病灶区域；step6、初始轮廓输入；将选择的病灶区域作为cv模型分割的初始轮廓；step7、边界暗区域去除；在进行分割时，轮廓将病灶区域作为内部，除此以外的区域作为外部区域，而外部边界其他组织区域的存在会使外部区域更加不均匀，分割时将图像边缘的暗区域去除，那么得到较均匀的外部区域，使cv活动轮廓分割得到的轮廓更加准确；step8、cv模型分割；cv模型将会根据内外能量的变化，不断将轮廓趋近病灶区域的边缘，得到病灶区域的大致轮廓；step9、局部cv模型lcv分割；对cv模型分割得到的大致轮廓再作为lcv模型的初始轮廓，利用lcv模型，对病灶区域的大致轮廓在边缘局部上进一步微调，从而得到更准确的超声病灶区域的边缘轮廓，并获得病灶区域；step10、获得边缘；利用边缘算子canny获得病灶区域的边缘，并标到原图像上。所述步骤step1中，采用的滤波是双边滤波；双边滤波是一种非线性滤波器，它可以达到保持边缘、降噪平滑的效果。双边滤波是采用加权平均的方法，用周边像素亮度值的加权平均代表某个像素的强度，所用的加权平均基于高斯分布。同时，双边滤波的权重不仅考虑了像素的欧氏距离，还考虑了像素范围域中的灰度值差异，在计算中心像素的时候同时考虑这两个权重。cv模型是chan和vese提出的一种基于区域的水平集活动轮廓模型。所述步骤step7中，边界暗区域一般是其他非病灶组织，但往往跟病灶区域有相同的灰度。边界暗区域的去除可以较好解决超声图像外部灰度不均匀而导致cv模型分割不准确的现象。所述步骤step8、9中，cv模型分割可以得到病灶区域的大致轮廓，lcv模型能对将cv模型得到的大致轮廓加以微调，得到更准确的边缘。所述步骤step11中，若得到的tp的值高，fp，fn低，说明分割结果的准确性高，效果好。本发明的有益效果是：本发明采用双边滤波，模糊聚类分割，形态学处理，边界暗区域去除，全局和局部cv活动轮廓模型结合的方法，有效实现病灶区域的选择，较为准确实现提取出病灶区域的轮廓，解决超声图像病灶区域自动选择，分割不精确的问题。另外，本发明在实验得到数据tp在实例中均达到最高，fp，fn也是最低。这在定量分析上也说明提取的病灶区域轮廓也较为准确。附图说明图1为本发明的总体流程图；图2为本发明实施例1超声病灶区域分割的实验图；图2中：(a)带病灶区域超声原图；(b)模糊聚类分割及二值化后的图；(c)形态学处理后的图；(d)病灶区域选择后的初始轮廓；(e)边界暗区域去除及cv分割后的图；(f)局部lcv模型分割后的图；(g)分割结果及获得边缘后的图；(h)人工分割结果图；图3为本发明实施例2超声病灶区域分割的实验图。图3中：(a)带病灶区域超声原图；(b)模糊聚类分割及二值化后的图；(c)形态学处理后的图；(d)病灶区域选择后的初始轮廓；(e)边界暗区域去除及cv分割后的图；(f)局部lcv模型分割后的图；(g)分割结果及获得边缘后的图；(h)人工分割结果图。具体实施方式实施例1：如图1-2所示，一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法，对图像如图2(a)进行预处理，首先利用双边滤波器进行去噪，用于解决图像噪声的存在，影响分割准确性的问题；接着对图像进行模糊聚类，图像二值化处理和形态学处理后，进行病灶区域选择，同时获得初始轮廓；再对边界暗区域去除，用于解决超声图像病灶区域外不均匀的特点，来提高分割的准确性；最后利用全局和局部cv活动轮廓模型的结合分割，得到分割结果。进一步地，所述方法的具体步骤如下：step1、图像滤波；采用双边滤波器对超声图像进行滤波，对超声图像存在的斑点噪声，进行滤波；step2、模糊聚类分割；对滤波后的图像进行模糊c均值聚类分割，聚类数目设定为3，得到将灰度值信息分成三类的图像；step3、图像二值化；选用一个阈值，这个阈值为根据病灶区域的灰度均值设定的，将其中包含病灶区域的一类与其他两类再分类，得到二值图像，以便后续形态处理；如图2(b)所示；step4、形态学处理；得到的二值图像，先对二值图像进行膨胀处理，接着对图像取反，然后对边界的亮区域抑制，再对图像取反，再进行对图像进行腐蚀处理，得到多个连通区域的二值图像；如图2(c)所示；step5、病灶区域选择；首先，统计连通区域的像素个数，选择最多像素的连通区域，接着比较其他的连通区域，如果其他区域的像素点大于最大连通区域的25％，则将视其为病灶区域roi，否则将其视为非病灶区域；step6、初始轮廓输入；初始轮廓如图2(d)所示，将选择的病灶区域作为cv模型分割的初始轮廓；step7、边界暗区域去除；在进行分割时，轮廓将病灶区域作为内部，除此以外的区域作为外部区域，而外部边界其他组织区域的存在会使外部区域更加不均匀，分割时将图像边缘的暗区域去除，那么得到较均匀的外部区域，使cv活动轮廓分割得到的轮廓更加准确；step8、cv模型分割；cv模型将会根据内外能量的变化，不断将轮廓趋近病灶区域的边缘，得到病灶区域的大致轮廓；如图2(e)所示；step9、局部cv模型lcv分割；对cv模型分割得到的大致轮廓再作为lcv模型的初始轮廓，利用lcv模型，对病灶区域的大致轮廓在边缘局部上进一步微调，从而得到更准确的超声病灶区域的边缘轮廓，并获得病灶区域；如图2(f)所示；step10、获得边缘；利用边缘算子canny获得病灶区域的边缘，并标到原图像上。如图2(g)所示。效果评价。将三种不同方法分割后的结果与人工分割的结果进行比较，人工分割结果如图2(h)，采取基于分类评判的三种指标，包括tp(truepositive)，fp(falsepositive),fn(falsenegative)。分割方法比较：比较方法包括仅用的模糊c均值聚类和cv模型方法。先对每个方法分割效果进行分析对比，再通过定量分析的对比，采用tp,fp,fn三种指标。表1为本实施例分割结果的定量分析及方法比较。表1实施1分割结果的定量分析及方法比较methodtpfpfnfmc0.86783.98170.1322cv0.87710.00650.1229本方法0.92510.00520.0749实施例2：如图1、图3所示，一种基于全局和局部活动轮廓模型的超声图像病灶自动分割的方法，实例采取的方案是先对带有病灶的超声图像，图3(a)进行预处理，首先利用双边滤波器进行去噪，用于解决图像噪声的存在，影响分割准确性的问题；接着对图像进行模糊聚类，图像二值化处理和形态学处理后，进行病灶区域选择，同时获得初始轮廓；再对边界暗区域去除，用于解决超声图像病灶区域外不均匀的特点，来提高分割的准确性；最后利用全局和局部cv活动轮廓模型的结合分割，得到较准确的分割结果。进一步地，所述方法的具体步骤如下：step1、图像滤波；采用双边滤波器对超声图像进行滤波，对超声图像存在的斑点噪声，进行滤波；step2、模糊聚类分割；对滤波后的图像进行模糊c均值聚类分割，聚类数目设定为3，得到将灰度值信息分成三类的图像；step3、图像二值化；选用一个阈值，这个阈值为根据病灶区域的灰度均值设定的，将其中包含病灶区域的一类与其他两类再分类，得到二值图像，以便后续形态处理；如图3(b)所示；step4、形态学处理；得到的二值图像，先对二值图像进行膨胀处理，接着对图像取反，然后对边界的亮区域抑制，再对图像取反，再进行对图像进行腐蚀处理，得到多个连通区域的二值图像；如图3(c)所示；step5、病灶区域选择；首先，统计连通区域的像素个数，选择最多像素的连通区域，接着比较其他的连通区域，如果其他区域的像素点大于最大连通区域的25％，则将视其为病灶区域roi，否则将其视为非病灶区域；step6、初始轮廓输入；初始轮廓如图3(d)所示，将选择的病灶区域作为cv模型分割的初始轮廓；step7、边界暗区域去除；在进行分割时，轮廓将病灶区域作为内部，除此以外的区域作为外部区域，而外部边界其他组织区域的存在会使外部区域更加不均匀，分割时将图像边缘的暗区域去除，那么得到较均匀的外部区域，使cv活动轮廓分割得到的轮廓更加准确；step8、cv模型分割；cv模型将会根据内外能量的变化，不断将轮廓趋近病灶区域的边缘，得到病灶区域的大致轮廓；如图3(e)所示；step9、局部cv模型lcv分割；对cv模型分割得到的大致轮廓再作为lcv模型的初始轮廓，利用lcv模型，对病灶区域的大致轮廓在边缘局部上进一步微调，从而得到更准确的超声病灶区域的边缘轮廓，并获得病灶区域；如图3(f)所示；step10、获得边缘；利用边缘算子canny获得病灶区域的边缘，并标到原图像上。如图3(g)所示。效果评价。将三种不同方法分割后的结果与人工分割的结果进行比较，人工分割结果如图3(h)，采取基于分类评判的三种指标，包括tp(truepositive)，fp(falsepositive),fn(falsenegative)。分割方法比较：比较方法包括仅用的模糊c均值聚类和cv模型方法。先对每个方法分割效果进行分析对比，再通过定量分析的对比，采用tp,fp,fn三种指标。表2为本实施例分割结果的定量分析及方法比较。表2实施例2分割结果的定量分析及方法比较methodtpfpfnfmc0.819614.96320.1804cv0.85030.05520.1497所用方法0.93100.01720.0690上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。当前第1页12