一种基于形状相关性活动轮廓模型的超声图像分割方法技术

本发明专利技术公开了一种基于形状相关性活动轮廓模型的超声图像分割方法，包括以下步骤：步骤1：针对连续超声图像中病灶区域形状变化之间的相关性信息进行挖掘并构建形状相关性的低秩模型；步骤2：构建基于低秩约束的活动轮廓模型；步骤3：构建基于增广拉格朗日的优化算法用于快速计算分割结果。本发明专利技术解决了传统的基于有监督统计学习的分割方法在训练集不足和面对超声图像中病灶区域边缘模糊和病灶区域形状形变的情况下分割结果不准确的结果。

【技术实现步骤摘要】
一种基于形状相关性活动轮廓模型的超声图像分割方法
本专利技术属于计算机辅助治疗中有关超声图像处理领域，具体涉及一种基于形状相关性活动轮廓模型的超声图像分割方法，是针对软组织的超声图像序列分割方法。
技术介绍
一个准确且高效的超声图像分割方法对于提高计算机辅助治疗的效率具有重要意义。目前面临的挑战主要有两个:(1)相较于CT和MRI图像，超声图像普遍具有低信噪比和图像灰度分布不均匀的特点从而造成目标边缘轮廓模糊不清；(2)另一个潜在挑战是软组织受挤压造成形变，反映在超声图像上就是目标区域的形状变化较大；这些原因造成了一些传统的分割方法都不能在以上情况下得到理想的结果。
技术实现思路
为了解决上述的技术问题，本专利技术提供了一种基于形状相关性活动轮廓模型的超声图像分割方法。本专利技术所采用的技术方案是：一种基于形状相关性活动轮廓模型的超声图像分割方法，其特征在于，包括以下步骤：步骤1：针对连续超声图像序列中病灶区域形状变化之间的相关性信息进行挖掘并构建形状相关性的低秩模型；形状的相关性信息主要指在一个超声图像序列中病灶区域的轮廓形状的变化在统一的坐标系统下具有一定的相似性；步骤2：构建基于低秩约束的活动轮廓模型；步骤3：构建基于增广拉格朗日的优化算法用于快速计算分割结果。作为优选，步骤1中所述的低秩模型为：其中(xi,yi)表示轮廓中标志点的坐标，M是线性变换矩阵，该矩阵的元素包含了：缩放、旋转、切错和翻转变换，T表示平移向量；公式(1)转换成：[C1,…,Cn]T＝[Φ1,…,Φn]Ψ(2)其中：和分别表示M和T的元素，和表示向量C1的前一半和后一半；因矩阵Ψ的列空间的维度至多为6，故矩阵[C1,…,Cn]为低秩矩阵。作为优选，步骤2中所述的构建基于低秩约束的活动轮廓模型，其具体实现过程为将图像序列中病灶形状相关性的低秩性，将其低秩约束融入到公式(3)所示的传统的活动轮廓模型中，其中,I表示图像,C表示分割的区域轮廓,Ω1表示分割区域内部,Ω2表示分割区域之外的区域,β表示系数；公式(3)的优化方程如下：其中X＝[f1(C1),…,fn(Cn)]，λ表示形似系数，E表示分割误差的轮廓信息，D表示初始化的分割轮廓，fi(Ci)为每一张超声图像中活动轮廓模型的能量方程用于轮廓的演化，如下式所示：其中Ω1和Ω2分别表示轮廓内外的区域，μ1和μ2分别表示区域内外的强度的均值；利用公式(6)：将公式(4)中的求矩阵X的秩转化为求矩阵X的核范数，||E||1转化为||E||*范数。作为优选，步骤3中所述的构建基于增广拉格朗日的优化算法用于快速计算分割结果，是利用公式(6)构建基于增广拉格朗日的优化算法用于快速计算分割结果，其具体实现过程包括以下子步骤：步骤3.1：手动定义超声图像序列中目标区域的初始化轮廓；步骤3.2：根据公式(5)计算每帧超声图像中目标区域的分割轮廓；步骤3.3：根据公式(6)采用拉格朗日乘子法对步骤3.2中的分割轮廓进行优化；步骤3.4：矩阵X中的图像利用公式(6)计算的分割结果小于设定的阈值时，返回步骤3.3；步骤3.5：收敛以后的轮廓即为目标区域分割的轮廓区域。本专利技术相对于现有技术，本专利技术解决了传统的基于有监督统计学习的分割方法在训练集不足和面对超声图像中病灶区域边缘模糊和病灶区域形状形变的情况下分割结果不准确的结果。具体实施方式为了便于本领域普通技术人员理解和实施本专利技术，下面结合实施例对本专利技术作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本专利技术，并不用于限定本专利技术。本专利技术所采用的技术方案是：一种基于形状相关性活动轮廓模型的超声图像分割方法，包括以下步骤：步骤1：针对连续超声图像序列中病灶区域形状变化之间的相关性信息(形状的相关性信息主要指在一个超声图像序列中病灶区域的轮廓形状的变化在统一的坐标系统下具有一定的相似性)进行挖掘并构建形状相关性的低秩模型；其低秩模型为：其中(xi,yi)表示轮廓中标志点的坐标，M是线性变换矩阵，该矩阵的元素包含了：缩放、旋转、切错和翻转变换，T表示平移向量；公式(1)转换成：[C1,…,Cn]T＝[Φ1,…,Φn]Ψ(2)其中：和分别表示M和T的元素，和表示向量C1的前一半和后一半；因矩阵Ψ的列空间的维度至多为6，故矩阵[C1,…,Cn]为低秩矩阵。步骤2：构建基于低秩约束的活动轮廓模型；其具体实现过程为将图像序列中病灶形状相关性的低秩性，将其低秩约束融入到公式(3)所示的传统的活动轮廓模型中，其中,I表示图像,C表示分割的区域轮廓,Ω1表示分割区域内部,Ω2表示分割区域之外的区域,β表示系数；公式(3)的优化方程如下：其中X＝[f1(C1),…,fn(Cn)]，λ表示形似系数，E表示分割误差的轮廓信息，D表示初始化的分割轮廓，fi(Ci)为每一张超声图像中活动轮廓模型的能量方程用于轮廓的演化，如下式所示：其中Ω1和Ω2分别表示轮廓内外的区域，μ1和μ2分别表示区域内外的强度的均值；利用公式(6)：将公式(4)中的求矩阵X的秩转化为求矩阵X的核范数，||E||1转化为||E||*范数。步骤3：是利用公式(6)构建基于增广拉格朗日的优化算法用于快速计算分割结果，其具体实现过程包括以下子步骤：步骤3.1：手动定义超声图像序列中目标区域的初始化轮廓；步骤3.2：根据公式(5)计算每帧超声图像中目标区域的分割轮廓；步骤3.3：根据公式(6)采用拉格朗日乘子法对步骤3.2中的分割轮廓进行优化；步骤3.4：矩阵X中的图像利用公式(6)计算的分割结果小于设定的阈值时，返回步骤3.3；步骤3.5：收敛以后的轮廓即为目标区域分割的轮廓区域。应当理解的是，本说明书未详细阐述的部分均属于现有技术。应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本专利技术专利保护范围的限制，本领域的普通技术人员在本专利技术的启示下，在不脱离本专利技术权利要求所保护的范围情况下，还可以做出替换或变形，均落入本专利技术的保护范围之内，本专利技术的请求保护范围应以所附权利要求为准。本文档来自技高网...

【技术保护点】
一种基于形状相关性活动轮廓模型的超声图像分割方法，其特征在于，包括以下步骤：步骤1：针对连续超声图像序列中病灶区域形状变化之间的相关性信息进行挖掘并构建形状相关性的低秩模型；步骤2：构建基于低秩约束的活动轮廓模型；步骤3：构建基于增广拉格朗日的优化算法用于快速计算分割结果。

【技术特征摘要】
1.一种基于形状相关性活动轮廓模型的超声图像分割方法，其特征在于，包括以下步骤：步骤1：构建基于低秩约束的活动轮廓模型；具体实现过程为将图像序列中病灶形状相关性的低秩约束融入到公式(1)：其中,I表示图像,C表示分割的区域轮廓,Ω1表示分割区域内部,Ω2表示分割区域之外的区域,β表示系数；公式(1)的优化方程如下：其中X＝[f1(C1),…,fn(Cn)]，λ表示形似系数，E表示分割误差的轮廓信息，D表示初始化的分割轮廓，fi(Ci)为每一张超声图像中活动轮廓模型的能量方程用于轮廓的演化，如下式所示：其中Ω1和Ω2分别表示轮廓内外的区域，μ1和μ2分别表示区域内外的强度的均值；

【专利技术属性】
技术研发人员：何发智，倪波，陈晓，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42