专利名称：：一种增量式变分水平集快速医学图像分割方法
技术领域：
：本发明涉及医学图像分割领域，尤其涉及一种快速的医学图像分割方法。(二)
背景技术：
：在医学图像处理与分析应用中，图像分割技术起着关键的作用。医学图像分割的任务是从医学图像中提取包含重要诊断信息的感兴趣区域，为临床珍疗和病理学研究提供可靠的依据。由于医学图像的成像原理和人体组织结构本身的复杂性和差异性，医学图像与普通图像比较，不可避免地具有模糊、不均匀性等特点；同时，医学成像技术的快速发展使得获取各种复杂的海量医学图像数据成为可能，这些都对图像分割技术提出了更高的要求。近年來，基于曲线演化几何流的变形模型成为医学图像分割领域的研究热点之一。基于几何流的变形模型又分为两类一类是直接得到问题的几何流描述，以基于曲线演化的方法为代表，曲线演化方程即是一个几何流；另一类是来源T能量最小化的几何流，其思想是将所研究问题归结为一个泛函能量函数的极小问题，然后应用变分方法导出相应的几何流。近些年提出的水平集方法对基于几何流的图像分割方法的研究和应用具有深远的影响。水平集的基本思想是将运动的轮廓曲线(曲面)隐含地表达为高一维的曲面的零水平集，并通过曲面的运动来隐含地求解轮廓曲线(曲面)的演化。该方法是一种能自适应边界拓扑变化的稳定数值方法，不但提供了曲线演化的精确数值解法，而目.很好地解决了非常棘手的拓扑变化问题，因此已发展为当今医学图像分割中最为活跃和成功的研究领域之一。医学图像分割方法中，Chan-Vese模型是Chan和Vese在2001年提出的一种不依赖于梯度的分片常值简化Mumford-Shah模型，该模型根据区域的均值将图像仅分为两个区域(目标和背景)，即两相分割，并采用水平集方法进行数值求解。基于变分水平集方法的Chan-Vese模型综合利用了整个图像区域的全部信息，因而其显著特点是全局优化性。初始轮廓曲线无需完全定位于目标的内部或者外部，而且仅需要一条初始轮廓曲线，就可以分割即使内部含有空洞的目标。另外，该模型不依赖于边缘梯度，即使图像中的边缘比较模糊，仍然能够得到比较理想的结果。但该模型存在一个重要缺陷，即Chan-Vese模型的演化过程需要利用整个图像的全部信息，并且每次迭代都需要(且必须)在整个图像域h计算，因而无法采用窄带法或快速步进法等快速算法来提高分割效率。(三)
发明内容本发明的目的在于提供一种根据区域和区域中像素的动态变化，采用增量方式迭代求取区域的平均灰度，将其解析公式改进为递进迭代公式，从而可以采用窄带法等快速算法，能够较大地提高分割效率的增量式变分水平集快速医学图像分割方法。本发明的目的是这样实现的首先选取初始边界，根据初始边界计算的子区域及其平均灰度，采用窄带法等快速算法求解水平集曲线演化过程，提取零水平集即新的边界，判断是否满足停止条件，若是则得到分割结果，若否则利用边界的运动导致区域的改变，根据增量在窄带范围内计算新区域的平均灰度，再进行窄带法等快速算法求解水平集曲线演化过程，最终的零水平集即为分割结果。本发明的详细过程为首先根据实际分割问题选择合适的模型参数//,^^1,^2，进行如下分割过程设AQ,，AQ2分别表示区域QpQ2增加或减少的部分，符号lQl表示区域Q的面积，对于离散问题来说即区域上的点的个数；设动态区域QcW"为有界连通丌集，AQ为其改变的区域，xe9T，(1)给定初始边界Co'以此边界可得子区域Q，Q〖，分别计算其平均灰度^,W，并将《，《分别作为q，2;根据公式(2)、《(3)(4)(5)鸣=对轮廓曲线Co进行单歩演化，实现上可以采用窄带法或快速步进法来快速求解；提取上一歩演化结果的零水平集，并以此为新的边界，得到两个更新后的区域Qpf^,检查演化停止条件，满足则转到歩骤(7);分别计算区域Q,，A的变化区域AQPAQ2;采用递进方式，分别按照公式《l淵+.4|淵+JJ线/.#(一2)Cte勿W2计算区域Qp^的平均灰度^，《2;(6)转到歩骤(2);(7)区域QpQ2即为分割结果。本发明根据区域和区域中像素的动态变化，采用增量方式迭代求取区域的平均灰度，将其解析公式改进为递进迭代公式，从而可以采用窄带法等快速算法，较大地提高了分割效率。一种快速的医学图像分割方法，属于医学图像分割领域。本发明针对Chan-Vese模型在演化过程中每次迭代都必须在整个图像区域上计算而导致无法釆用快速算法的缺陷，从基本理论公式的推导出发，提出一种增量式变分水平集快速医学图像分割方法。该方法根据区域和区域中像素的动态变化，采用增量方式迭代求取区域的平均灰度，将其解析公式改进为递进迭4代公式，从而可以采用窄带法等快速算法，较大地提高了分割效率，使得该模型更具实际意义。本发明在医学图像分割领域有着广阔的应用前景。本发明针对Chan-Vese模型无法采用快速算法的缺陷进行改进，提出一种增量式变分水平集快速医学图像分割方法。该方法给出了区域的平均灰度的增量式迭代求取公式，不再同时需要全部图像区域上的信息，因而能够利用窄带法或快速歩进法來加快求解过程，解决了Chan-Vese模型因不能采用快速算法而导致分割过慢的重要不足。山于仍然利用图像的全局信息，从而在保持分割精度不变的同时较大地提高了分割速度，使得该模型在复杂的医学图像分割领域更具实际意义。我们对区域(以及区域中像素)的动态变化进行考察。区域的变化是因轮廓曲线(区域边界)运动而产生的。设AQpAQ2分别表示区域q,A增加或减少的部分，符号lQl表示区域Q的面积，对于离散问题来说即区域上的点的个数。设动态区域QcT为有界连通丌集，AQ为其改变的区域，xeW"，若().xgAQ(1)—1,xe△《}且xgn则称^:5T4{-1,0,1}为AD关于Q的符号隶属函数，简称的符号隶属函数。符号隶属函数^的意义是，在区域Q的变化区域AQ中，对于进入D的像素，符号隶属值取l;对于离开Q的像素，符号隶属值取-1。令i(qi为符号函数，即<formula>formulaseeoriginaldocumentpage5</formula>设动态区域Qcr为有界连通开集，AQ为其改变的区域，若则称A:PW为的符号面积。正的符号面积表示区域增大，负的符号面积表示区域縮小，零则说明区域保持大小(面积)不变的情况下改变了形状或位置。利用式(1)和式(2)，AQ,，AQ2上的灰度变化总量可分别表示为<formula>formulaseeoriginaldocumentpage5</formula>其中^为AQ,关于Q,的符号隶属函数，V/2为AQ2关于Q2的符号隶属函数。注意到，若区域Q经过变化区域后更新为^，则Q,上的平均灰度可按下式计算<formula>formulaseeoriginaldocumentpage6</formula>(5)其中^为Q上的平均灰度。将式(3)和(4)代入上式，得到同理可计算Q经过变化AQ2后的平均灰度<formula>formulaseeoriginaldocumentpage6</formula>(6)<formula>formulaseeoriginaldocumentpage6</formula>(7)综上分析，整理得，，<formula>formulaseeoriginaldocumentpage6</formula>(8)我们将上式称为增量式Chan-Vese模型，并将利用该模型进行医学图像分割的方法称为增量式变分水平集快速医学图像分割方法。事实上，在曲线演化过程中，区域的变化是由于边界曲线(即零水平集)C的单步运动产生的，因而变化区域M^，Af^总是在零水平集邻近的小范围内，可以全部包含在一个窄带之中，因此可以采用窄带法或快速歩进法来快速求解。对于两相分割问题，进入Q,巾的像素就是离开^的像素，因此有=-AQ2，进一步，这两个变化区域上的符号面积及其像素值的和分别满足如下关系<formula>formulaseeoriginaldocumentpage6</formula>(9)增量式Chan-Vese模型依然利用了图像全局信息，保留了全局优化的优点，且减少了计算量。一方面，迭代过程中平均灰度的求取大为简化。由式(8)和式(9)可知，仅仅需要在区域的变化的部分计算，而传统的计算公式必须对整个图像区域计算。另一方面，可以采用一般快速算法來加速水平集曲线演化过程，从而进一歩降低计算量。需要指出的是，从增量式Chan-Vese模型的推导过程可知，该模型与传统的Chan-Vese模型在理论上是完全等价的，因此分割精度不受影响。本发明针对Chan-Vese模型在演化过程中每次迭代都必须在整个图像区域上计算而导致无法采用快速算法的缺陷，从基本理论公式的推导出发，提出一种增量式变分水平集快速医学图像分割方法。该方法根据区域和区域中像素的动态变化，采用增量方式迭代求取区域的平均灰度，将其解析公式改进为递进迭代公式，从而可以采用窄带法等快速算法，较大地提高了分割效率，使得该模型更具实际意义。本发明在医学图像分割领域有着广阔的应用前景。(四)图1为增量式变分水平集方法的算法流程图2为脑MR(T0图像；图3-图5为脑组织分割结果。具体实施方式下面结合附图对本发明作进一歩的说明本发明中提出的增量式变分水平集快速医学图像分割方法主要依据公式(9)进行数值计算，实施简洁，以下结合附图对本发明的具体实施方式作进一歩的详细说明。二维图像的两相分割问题为例，原图像如图2所示。首先根据实际分割问题选择合适的模型参数A,v,;^，;i2，本实施例中分别为//=1.0，v=o，；^=;^=1.0，然后参照图i进行如下分割过程(1)给定初始边界Q)，以此边界可得区域Q，Q，分别计算其平均灰度^,《，并将^,W分别作为^,,2;(2)根据公式(9)对轮廓曲线Co进行单步演化，实现上可以采用窄带法或快速歩进法来快速求解；(3)提取上一步演化结果的零水平集，并以此为新的边界，得到两个更新后的区域Q"Q2，检查演化停止条件，满足则转到(7);(4)分别计算区域QpQ2的变化区域AQ,，AQ2;(5)采用递进方式，分别按照公式(6)和公式(7)计算区域Q,，Q2的平均灰度q，^;(6)转到(2);(7)区域QpQ2即为分割结果。图3-图5给出了分割结果，从左至右分别为脑灰质、脑白质和脑脊液。图l中手工选取的初始边界101，根据初始边界计算的子区域及其平均灰度102，采用窄带法等快速算法求解水平集曲线演化过程103，提取零水平集即新的边界104，判断是否满足停止条件，若是则得到分割结果，若否则利用边界的运动导致区域的改变105，根据增量公式(6)和(7)在窄带范围内计算新区域的平均灰度106，再进行窄带法等快速算法求解水平集曲线演化过程103，最终的零水平集即为分割结果107。图3中分别为脑灰质301、脑白质302、脑脊液303。表1给出了本发明方法与原Chan-Vese模型分割的时间比较，表中数据表明本发明方法减少了约22%的分割时间。虽然对于单幅图像来说分割时间以ms计，但实际医学图像通常都具有几十或几百个切片，因此节省的时间是非常可观的。表1分割时间的比较<table>tableseeoriginaldocumentpage8</column></row><table>权利要求1、一种增量式变分水平集快速医学图像分割方法，其特征在于首先选取初始边界，根据初始边界计算的子区域及其平均灰度，采用窄带法等快速算法求解水平集曲线演化过程，提取零水平集即新的边界，判断是否满足停止条件，若是则得到分割结果，若否则利用边界的运动导致区域的改变，根据增量在窄带范围内计算新区域的平均灰度，再进行窄带法等快速算法求解水平集曲线演化过程，最终的零水平集即为分割结果。2、根据权利要求1所述的一种增量式变分水平集快速医学图像分割方法，其特征在于首先根据实际分割问题选择合适的模型参数P,V，^，^，进行如下分割过程设AQ,，Af^分别表示区域QpQ2增加或减少的部分，符号M表示区域Q的面积，对于离散问题来说即区域上的点的个数；设动态区域Qc9r为有界连通开集，AQ为其改变的区域，xe3T，(1)给定初始边界Q)，以此边界可得子区域Q;1，Q，分别计算其平均灰度^1,《，并将W,《分别作为q,ry2;(2)根据公式<formula>formulaseeoriginaldocumentpage2</formula>对轮廓曲线Co进行单歩演化，实现上可以^用窄带法或快速步进法来快速求解；(3)提取上一步演化结果的零水平集，并以此为新的边界，得到两个更新后的区域^,02,检查演化停止条件，满足则转到歩骤(7);(4)分别计算区域Q,,Q2的变化区域AQ,,AQ2;(5)采用递进方式，分别按照公式<formula>formulaseeoriginaldocumentpage2</formula>计算区域Q"Q2的平均灰度^,^;(6)转到步骤(2);(7)区域Q,，Q2即为分割结果。全文摘要本发明提供了一种增量式变分水平集快速医学图像分割方法，首先选取初始边界，根据初始边界计算的子区域及其平均灰度，采用窄带法等快速算法求解水平集曲线演化过程，提取零水平集即新的边界，判断是否满足停止条件，若是则得到分割结果，若否则利用边界的运动导致区域的改变，根据增量在窄带范围内计算新区域的平均灰度，再进行窄带法等快速算法求解水平集曲线演化过程，最终的零水平集即为分割结果。本发明根据区域和区域中像素的动态变化，采用增量方式迭代求取区域的平均灰度，将其解析公式改进为递进迭代公式，从而可以采用窄带法等快速算法，较大地提高了分割效率，使得该模型更具实际意义。文档编号G06T7/00GK101315699SQ200810064898公开日2008年12月3日申请日期2008年7月10日优先权日2008年7月10日发明者毅沈,艳王,郝家胜申请人:哈尔滨工业大学