一种基于房室平面泵模型的心功能诊断方法技术

一种基于房室平面泵模型的心功能诊断方法，包括：通过临床获取的心脏核磁共振图像，在图像上对左心室区域进行分割后计算左心室容积和流量，通过对房室特征点运动轨迹的跟踪计算出房室平面位移和速度，基于心脏动力学实现对左心室的模型构建，通过无迹卡尔曼滤波器进行系统辨识，最后通过非线性优化的方法对模型参数进行优化，利用心脏参数实现心功能诊断。本发明专利技术的方法可以实现心脏磁共振短轴图像的准确分割，有效地解决了需要人工描绘轮廓的费时费力问题；并克服传统归一化互相关匹配法跟踪不准确的问题，所测得的房室平面位移曲线更符合生理学表现；构建了新的心脏模型，可以更为准确地刻画心脏的动力学特征，再通过无迹卡尔曼滤波器确保辨识参数的可靠性，方法简单且辨识准确度高，最后通过非线性优化对左心室模型进行优化处理，所得参数稳定，准确，能够协助医生进行心脏疾病诊断。医生进行心脏疾病诊断。

【技术实现步骤摘要】
一种基于房室平面泵模型的心功能诊断方法


[0001]本专利技术涉及本专利技术涉及医学图像处理技术，心脏动力学模型和系统辨识技术，具体涉及一种基于房室平面泵模型的心功能诊断方法。

技术介绍

[0002]近年来，心血管疾病的死亡率居高不下，提高心脏病的检出率更是迫在眉睫。临床发现，传统的介入式心功能参数的测量方法，不仅测量方法复杂，且对人体伤害较大；此外，射血分数保留型的心脏疾病的发现，也在说明传统的心脏疾病评判指标的局限性。为此，迫切需要一种参数容易测量的、非介入式的、可量化的、具有普遍性的综合性评判指标。
[0003]心脏核磁共振成像具有非侵入式的特点，它可以清晰的显示人体内脏器、肌肉、骨骼、血液、脂肪等的物理结构。基于病患心功能参数的测量结果，医生便可以快速确定病变的位置、炎症的类型、异常物质起因等，从而尽早、准确的实现心脏疾病的诊断，提升心脏病诊断的精度和效率。
[0004]为了实现对心功能参数的量化分析，需要实现对心脏的建模分析。现有的心脏模型中，大量的非线性方程往往会导致其参数辨识困难。在此文中，我们提出了一种基于房室平面位移的建模方法，通过无迹卡尔曼滤波器进行参数辨识，加以非线性优化，实现心功能参数的辨识，协助医生进行心脏疾病诊断。

技术实现思路

[0005]本专利技术的目的是提供一种简单且准确度高的心功能诊断的方法，通过临床获取的心脏核磁共振图像，在图像上对左心室区域进行分割后计算左心室容积和流量，通过对房室特征点运动轨迹的跟踪计算出房室平面位移，基于心脏动力学实现对左心室的模型构建，通过无迹卡尔曼滤波器进行系统辨识，最后通过非线性优化的方法对模型进行优化，实现心功能诊断。
[0006]为了实现上述专利技术目的，本专利技术提供了一种基于房室平面泵模型的心功能诊断方法，它包含以下步骤：
[0007]S1：通过临床获取心脏核磁共振图像，完成左心室的分割，实现心室容积的计算；
[0008]S2：在左心室的核磁共振图像中，选定左心室的特征点，利用跟踪算法跟踪特征点得出房室平面位移；
[0009]S3：基于心脏动力学原理，建立左心室模型：
[0010]S4：通过无迹卡尔曼滤波器对心脏的生理参数进行辨识；
[0011]S5：使用非线性规划的方法对模型进行优化。
[0012]进一步的，步骤S1中，所述心室容积计算步骤如下：
[0013]基于临床采集的左心室核磁共振图像以及医学影像专家手动标注的左心室轮廓构建左心室分割训练集；对训练集进行数据增强，获取足量的训练集和测试集用以U
‑
NET网络模型的训练；采用训练好的U
‑
NET分割左心室，得到每张切片中左心室像素意义上的面
积；结合每张切片中左心室的面积及相邻切片间的物理间距等数据，实现对心脏个体左心室容积及其流量的计算。
[0017]进一步的，步骤S2中，所述房室平面位移计算步骤如下：
[0018]通过临床测得的心室长轴图像，制作房室特征点数据集并进行标注，划分训练集和验证集；构建模型网络结构，选择跟踪目标图像并确定搜索区域，提取特征后将不同尺度的特征图分别送入SiamRPN网络，直接对其输出采用加权和，由合并的SiamRPN输出包含分类分支和回归分支；使用L1损失来计算锚框的偏移量，和一个交叉熵损失来计算置信度的损失，实现对心室特征点的自动跟踪并计算出房室平面位移。相比于传统归一化互相关匹配法跟踪房室平面关键点位置有许多异常值的情况，所测得的房室平面位移曲线更符合生理学表现。
[0021]进一步的，步骤S3中，心脏模型构建原理如下：
[0023]将左心室外壁近似为一个圆柱和一个半球，利用微元法进行分析，建立轴向函数模型，然而这种模型建立时只考虑了心脏纵向运动的因素，根据心脏造影显示，实际上在心脏纵向位移时也伴随着少量周向位移。
[0024]为了能够更加准确的描述心脏运动特性，在房室平面模型中，将横向位移等效成一个系数k
RAD
。建立的心脏模型为：
[0025][0026][0027]其中，V
AVP
表示房室平面位移速度，表示房室平面位移变化率，A
LV
是心室侧等效活塞面积，A
LA
是心房侧等效活塞面积，F
C
是心肌力，Q
LV
是左心流量，P
LV
表示左心室的压力，表示左心室压力变化率，P
LA
表示左心房的压力，C
LV
表示左心室顺应性，k
RAD
表示径向函数系数。
[0028]进一步的，步骤S4中，基于无迹卡尔曼滤波器的参数辨识方法如下：
[0029]基于心肌的张力与收缩关系，将心室力等效成梯形，选取的待辨识的参数如下：
[0030][0031]其中，A1和A2分别表示心室收缩力F
VC
和心房收缩力F
AC
的幅值，R
AVP
表示AVP的阻尼，L
AVP
表示AVP的惯性，以左心室流量作为输入量，房室平面位移速度作为观测量，利用无迹卡尔曼滤波器实现对房室平面泵模型参数的辨识；基于所测的阻尼、惯性、响应性、心肌力的参数信息，以及心室的几何性及周向扩张等信息，可实现心脏状态的评估。
[0032]进一步的，步骤S5中，基于非线性优化的方法对模型进行优化。以梯形心肌力作为初始值，基于离散化的思想，将心肌力曲线中的每个点都作为待优化的变量，以动力学方程作为约束，房室平面位移速度估计误差最小化作为评判指标，对其进行非线性优化，使用内点法求解，实现对心脏参数和心肌力的形状进行进一步优化。临床发现，心脏疾病的产生，往往伴随着相关心脏参数的大幅改变，因此，通过对心脏参数的测量，可以协助医生进行心功能诊断。
附图说明
[0033]图1是本专利技术整体方案流程图；
[0034]图2是基于U
‑
NET网络的分割效果；
[0035]图3是房室特征点示意图；
[0036]图4是SiamRPN网络框架；
[0037]图5是房室平面的心脏模型；
[0038]图6是最大面积/纵向位移与射血分数关系示意图；
[0039]图7是心肌力F
C
示意图。
具体实施方式：
[0040]下面将结合附图和具体实施例，对本专利技术的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0041]通过临床获取的心脏核磁共振图像，在图像上对左心室区域进行分割后计算左心室容积和流量，通过对房室特征点运动轨迹的跟踪计算出房室平面位移，基于心脏动力学实现对左心室的模型构建，通过无迹卡尔曼滤波器进行系统辨识，最后通过非线性优化的方法对模型进行优化，实现心本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于房室平面泵模型的心功能诊断方法，其特征在于，包括：S1：通过临床获取心脏核磁共振图像，完成左心室的分割，实现左心室容积的计算；S2：在左心室的核磁共振图像中，选定左心室的特征点，利用跟踪算法跟踪特征点得出房室平面位移；S3：基于心脏动力学原理建立左心室模型；S4：通过无迹卡尔曼滤波器对心脏的生理参数进行辨识；S5：使用非线性优化的方法对模型进行优化。2.根据权利要求1所述的一种基于房室平面泵模型的心功能诊断方法，其特征在于，步骤S1中，左心室容积计算步骤包括：基于临床采集的左心室核磁共振图像以及医学影像专家手动标注的左心室轮廓构建左心室分割训练集；对训练集进行数据增强，获取足量的训练集和测试集用以U
‑
NET网络模型的训练；采用训练好的U
‑
NET分割左心室，得到每张切片中左心室像素意义上的面积；结合每张切片中左心室的面积及相邻切片间的物理间距等数据，实现对心脏个体左心室容积及其流量的计算。3.根据权利要求1所述的一种基于房室平面泵模型的心功能诊断方法，其特征在于，步骤S2中，所述房室平面位移计算步骤包括：通过临床测得的左心室长轴图像，制作房室特征点数据集并进行标注，划分训练集和验证集；构建模型网络结构，选择跟踪目标图像并确定搜索区域，提取特征后将不同尺度的特征图分别送入SiamRPN网络，直接对其输出采用加权和，由合并的SiamRPN输出包含分类分支和回归分支；使用L1损失来计算锚框的偏移量，和一个交叉熵损失来计算置信度的损失，实现对心室特征点的自动跟踪并计算出房室平面位移。相比于传统归一化互相关匹配法算法跟踪房室平面关键点位置有许多异常值的情况，所测得的房室平面位移曲线更符合生理学表现。4.根据权利要求1所述的一种基于房室平面泵模型的心功能诊断方法，其特征在于，步骤S3中，所述房室平面位移计算步骤如下：1)将左心室外壁近似为一个圆柱和一个半球，利用微元法进行分析，建立轴向函数模型，然而这种模型建立时只考虑了心脏纵向运动的因素，根据心脏造影显示，实际上在心脏纵向位...

【专利技术属性】
技术研发人员：黄欢，杨跃，
申请(专利权)人：江苏师范大学，
类型：发明
国别省市：