基于确定学习理论的心肌缺血辅助检测方法

摘要

本发明公开了一种基于确定学习理论的心肌缺血辅助检测方法，该方法包括如下步骤：数据预处理，获取心电图ST-T环数据；采用确定学习算法对ST-T环的内在系统动态进行局部准确建模，将所建模的系统动态进行三维可视化显示得到CDVG；根据CDVG形态并结合高危因素与临床信息进行分析，得到辅助检测结果。该方法的技术要点是基于确定学习算法的CDVG的获得与分析。本发明方法适用于心电图尚未发生明显改变时的心肌缺血检测，可有效提示心血管疾病的发生，还可用于心血管疾病治疗效果的评估，具有检测过程简单、经济、无创，减轻了患者负担等优点。

说明书

技术领域

本发明属于医学检测技术领域，具体涉及一种基于确定学习理论的心肌 缺血辅助检测方法，该方法是一种针对心电图的ST-T段或者心向量图的 ST-T环的内在系统动态信息提取和检测方法。

背景技术

长期以来，心血管疾病己被公认为是危害人类生命健康最严重的疾病之 一，其中由于心肌缺血导致的心肌梗塞(MI)的发病率和死亡率更是高居各类 疾病之首。由于部分心肌缺血患者在发病早期并没有明显的临床症状或者病 症轻微，使得病情十分容易被忽视。如果在日常保健检查中采用对心肌缺血 更为敏感的检测手段，就可以密切关注其先兆症状，能在心电图尚未发生显 著变化时察觉到心肌缺血病症，使得这类病人在症状早期阶段得到更积极的 干预治疗，从而挽救患者生命。

心肌缺血在理论上有一定的可预测性，体表心电图（ECG）是临床上常用 的心脏特征观察与检测手段，临床上认为心电信号中的ST-T段的偏移反应心脏 冠脉机能的异常，对心肌缺血的诊断有重要的参考价值。对许多患者而言，特 别是ECG的ST-T段无明显改变的心肌缺血患者，现有心电图检测技术的准确率 确实不高，并且难以对ST-T段的微小改变敏感，进而难以把握心血管的器质性 变化。实际上，至今的ST-T段检测方法对微小改变不敏感和准确率不高的原因 有：信号本身微弱，存在各种干扰；以ST-T段对比等电位值的偏移量作为特征 参数的选择依据，灵敏性和准确性并不能令人满意；对ST-T段没有统一的测试 标准和公认的测试数据库，甚至其形态的起始点没有严格统一的定量规定；大 多数病患病情不稳定，使得检测信号呈现出明显的动态变化，不便于检测判定。 对于心肌缺血的诊断，心电向量图（VCG）是另一个工具，它包含了丰富的 动态信息，对了解心脏电激动的顺序及各瞬间电位的改变方面比ECG更为优 越。采用ECG和VCG联合诊断心肌缺血能有效的提高诊断率，但在临床中由 于分析难度相对较大使得VCG的应用并不广泛。其它广泛应用于临床的检查 方法包括平板运动试验、冠脉CTA、核素心肌灌注显像和冠脉造影等。

本发明方法与传统ECG、VCG和平板运动试验相比，对心肌缺血/心肌 梗塞具有更为敏感的检测能力；与冠脉CTA、核素心肌灌注显像和冠脉造影 等检查方法相比，更加简单、经济,无创。本发明方法更适用于心电图尚未 发生明显改变时的心肌缺血检测，可有效提示心血管疾病的发生，还可用于 心血管疾病治疗效果的评估，为心肌缺血/心肌梗塞的准确诊断提供一种新方 法。

发明内容

本发明的目的是为了克服现有技术的缺点与不足，提供一种基于确定学 习理论的心肌缺血辅助检测方法。

基于ECG或者VCG进行心肌缺血检测，本质上可以看作是一个动态模 式的辨识与识别问题，而动态模式问题本身就是模式识别领域的一个难题。 在对径向基函数（RBF）神经网络的持续激励特性研究的基础上，C.Wang 等提出了确定学习理论，其中包括对由非线性动力学系统产生的动态模式的 辨识、表达和快速识别方法，即通过确定学习获得的动态模式内在系统动态 的局部准确神经网络建模，把随时间变化的动态模式以时不变且空间分布的 方式有效地表达，进一步利用动态模式内在的动力学拓扑相似给出动态模式 之间相似性定义，并提出了对动态模式进行快速识别的一套新方法。在实际 工程应用中的系统具有离散时间特性，所以需要对确定学习理论进行扩展以 满足实际工程需求。在中国发明专利申请：一种基于确定学习理论的海量时 态数据建模方法（申请号：201210560790.X）中，王聪等提出了一种针对离 散非线性系统产生的海量时态数据建模的确定学习算法，将确定学习理论在 时态数据建模领域进行了拓展。

本发明将确定学习理论应用于对ECG的ST-T段和VCG的ST-T环的内 在系统动态的局部准确建模和全息特征提取，并将所学习到的系统动态知识 进行三维可视化显示得到心电动力学向量图（Cardiodynamics vectorgram，简 称：CDVG），根据CDVG形态，并结合高危因素与临床信息进行心肌缺血 检测分析。本发明方法适用于心肌缺血早期检测和冠脉健康监测，为日常心 脏保健提供新的依据。

本发明的目的通过如下技术方案实现：一种基于确定学习理论的心肌缺 血辅助检测方法，包含如下步骤：

（1）ST-T环数据获取：对12导联心电图ECG数值数据进行预处理之 后，通过变换法则将ECG数据转换为三维的VCG数据，再截取三维数据中 的ST-T段，从而获取VCG的ST-T环，并在空间中进行归一化处理；

（2）建模并显示CDVG：采用确定学习算法对ST-T环的内在系统动态 进行局部准确RBF神经网络建模，将所建模的系统动态进行三维可视化显示 得到心电动力学向量图CDVG；

（3）检测分析：根据CDVG形态，并结合年龄、病史和血压这些高危 因素与临床症状等信息进行分析，得到辅助检测结果。

上述方法中，步骤（1）中所述12导联心电图ECG数值数据是指传统 12导联Ⅰ、Ⅱ、Ⅲ、aVR、aVL、aVF、V1、V2、V3、V4、V5和V6的数 值形式数据。

上述方法中，步骤（1）中所述的预处理是指常规滤波操作，以使随后的 操作充分稳定和可靠，使用何种滤波方法并不影响本发明方法的有效性。

上述方法中，步骤（1）中所述的通过变换法则将ECG数据转化为三维 的VCG数据是指将12导联心电图（ECG）数据通过既定的公认有效的变换 法则转换为三维的心电向量图（VCG）数据，记为：

V(k)=[v_x(k),v_y(k),v_z(k)]^T∈R³，

式中，k=1,…,m为采样时刻。其中的变换法则根据研究重点不同存在多 种方法，使用何种变换法则并不影响本发明方法的有效性。

上述方法中，步骤（1）中所述的截取三维数据中的ST-T段是为了获得 VCG数据中的ST-T环，所截取的ST-T段的三维数据可表述为：

$V_{\mathrm{ST}}=V\left(k\right){|}_{k=k_S}^{k=k_{\mathrm{TE}}},$

式中，k_S和k_TE分别代表心电图ST-T段的起始点（ST段的起点）和结束 点（T波的终点）。在截取准确的前提下，使用何种截取方法并不影响本发明 方法的有效性。

上述方法中，步骤（1）中所述的归一化是指对VCG的ST-T环三维数据进 行单位球内的归一化，具体操作是：计算VCG截取出的ST-T环的三维数据VST的 范数：

$\left|\right|V\left(k\right)\left|\right|=\sqrt{v_x^2\left(k\right)+v_y^2\left(k\right)+v_z^2\left(k\right),}$

式中，k∈[k_S,k_TE]，找到其中的最大范数||V(k)||_max，将整个ST-T环上所有数 据点坐标除以该最大范数，获得新的数据点位置，即将ST-T环归一化到一个 以原点为球心的单位球内。

上述方法中，步骤（2）中所述的对ST-T环的内在系统动态进行局部准 确建模是指对VCG的ST-T环的内在系统动态利用确定学习算法进行局部准 确RBF神经网络逼近。由于：

（i）心脏系统可以被建模为复杂的高维连续非线性动力学系统，并可进 一步被一个三维非线性动力学系统近似表达，其中：

V(t)=[v_x(t),v_y(t),v_z(t)]^T∈R³是系统状态，代表产生的三维VCG信号，

F(V(t))=[f₁(V(t)),f₂(V(t)),f₃(V(t))]^T是未知的非线性函数向量；

（ii）当采用计算机进行心电图分析时，该系统可以通过采样变换为近 似的离散非线性动力学系统：

V(k+1)=V(k)+T_sF(V(k))，

式中，V(k)=[v_x(k),v_y(k),v_z(k)]^T∈R³是系统状态，代表产生的三维VCG数 据，来源于对心脏系统VCG信号的采样，k=1,…,m为采样时刻，T_s是采样周 期，F(V(k))=[f₁(V(k)),f₂(V(k)),f₃(V(k))]^T表示内在系统动态，当T_s无限小时，此 离散非线性系统无限近似原连续非线性系统；

因此，采用确定学习理论中的动态RBF神经网络辨识器，实现对VCG 的ST-T环的内在系统动态的局部准确神 经网络逼近，即：

$\left(\begin{array}{c}{F}_{\mathrm{ST}}=(f_1\left(V\left(k\right)\right),f_2\left(V\left(k\right)\right),f_3\left(V\left(k\right)\right)){|}_{{V\left(k\right)\in V}_{\mathrm{ST}}}\\ \cong ({\bar{w}}_1^{T}S\left(V\left(k\right)\right),{\bar{W}}_2^{T}S\left(V\left(k\right)\right),{\bar{W}}_3^{T}S\left(V\left(k\right)\right)){|}_{V\left(k\right)\in V_{\mathrm{ST}}},\end{array}\right)$

式中，是常值神经网络权值向量，S(V(k))是高斯型径向基函数。

上述方法中，步骤（2）中所述的CDVG是对VCG的ST-T环的内在系统动 态$F_{\mathrm{ST}}=(f_1\left(V\left(k\right)\right),f_2\left(V\left(k\right)\right),f_3\left(V\left(k\right)\right)){|}_{V\left(k\right)\in T_{\mathrm{ST}}}$的近似准确RBF神经网络建模结果

$({\bar{W}}_1^{T}S\left(V\left(k\right)\right),{\bar{W}}_2^{T}S\left(V\left(k\right)\right),{\bar{W}}_3^{T}S\left(V\left(k\right)\right)){|}_{V\left(k\right){\in V}_{\mathrm{ST}}}$的三维可视化显示。

上述方法中，步骤（3）中所述的CDVG形态是指三维空间中显示出的 图像形态，对正常心电数据其CDVG形态为规整环形或略散环形，对心肌缺 血患者的心电数据其CDVG形态为散乱环形或无环形。

上述方法中，步骤（3）中所述的分析是指将经本发明方法获得的CDVG 进行形态判断，再结合患者临床信息进行综合分析诊断。临床信息指患者的 年龄、病史、血压、生活习惯等高危因素与就诊时的症状等信息。

本发明与现有技术相比，具有如下的优点和有益效果：

1、本发明能够克服目前心电图系统对心肌缺血诊断大多只用ST段抬高 或压低的幅度和T波的形态变化为判别的参考所引起的灵敏性和准确性不高 的问题，通过采用确定学习算法对心电信号ST-T环内在的系统动态在局部 区域内的准确建模，产生CDVG，能较为敏感地反映心肌缺血的程度，为心 肌缺血/心肌梗塞的准确诊断提供一种新方法。

2、本发明首次提出CDVG的概念，并将其进行三维可视化动态显示， 与传统ECG、VCG和平板运动试验相比，对心肌缺血/心肌梗塞具有更为敏 感的检测能力，检测结果更为直观，检测准确率高；与冠脉CTA、核素心肌 灌注显像和冠脉造影等检查方法相比，更加简单、经济,无创。

3、本发明所提一种基于确定学习理论的心肌缺血辅助检测方法，在ECG 尚未发生明显改变的情况下，对心肌缺血相关情况的检测具有明显优势，可 作为相关病症诊断的重要参考依据，还可用于心血管疾病治疗效果的评估。

附图说明

图1是实施例中PTB数据库P117的预处理后ECG。

图2是实施例中PTB数据库P117的ST-T环内在的系统动态的CDVG。

图3是实施例中PTB数据库P072的预处理后ECG。

图4是实施例中PTB数据库P072的ST-T环内在的系统动态的CDVG。

具体实施方式

下面结合实施例及附图，对本发明的具体实施方式作进一步详细说明， 但本发明的实施和保护范围不限于此。

实施例

本发明的具体实施例选取标准心电数据库PTB（Physikalisch-Technische  Bundesanstalt）数据库的正常数据P117和心肌梗塞数据P072，说明本发明方 法对心肌缺血相关病症检测的有效性。

采用基于确定学习理论的心肌缺血辅助检测方法的步骤如下：

（1）ST-T环数据获取：

利用心电数据采集设备获取12导联心电图ECG的数值数据：Ⅰ、Ⅱ、 Ⅲ、aVR、aVL、aVF、V1、V2、V3、V4、V5和V6。由于干扰的影响和心 电信号本身微弱的特性，心电信号往往埋藏在噪声中，很容易造成信息丢失 或产生虚假信息，所以要对信号进行预处理来降低噪声、增强信息。本实施 例中进行两项基本滤波：基线矫正和去除50HZ工频干扰，属于常规滤波操 作。滤波预处理选择常用的信号滤波方法皆可，以使随后的操作充分稳定和 可靠。本实施例中以选取的PTB数据库数据P117为例，经过数据预处理之 后的ECG信号如图1所示。

通过既定的公认有效的变换法则将ECG数据转换为三维的VCG数据。 变换法则根据研究重点不同存在多种方法，使用何种变换法则并不影响本实 施例的最终检测结果。除非另外特别指明，否则随后的操作都在经过变换后 的VCG三维数据上进行。由于VCG的ST-T环更能反应心肌缺血病症，包 含了心肌梗塞早期症状的主要动态，所以需要截取ST-T段，获得ST-T环。 在截取准确的前提下，使用何种截取方法并不影响本发明方法的有效性。本 实施例中的截取操作是：以ECG转换得到的VCG的X导联为截取标准，选 取在整个心电周期中用最高点R波定位，使用最大函数方法搜索到R波后， 沿后寻找J波，也就是ST段的起点，再根据面积法推断ST-T段。VCG三维 数据的ST-T段的空间表达即为ST-T环，再对ST-T环进行单位球内的归一 化，具体归一化操作如步骤说明所述。

（2）建模并显示CDVG：

如步骤说明中所述，当采用计算机进行心电图分析时，产生三维VCG 数据的心脏系统可近似为离散非线性动力学系统：

V(k+1)=V(k)+T_sF(V(k))，

式中，V(k)=[v_x(k),v_y(k),v_z(k)]^T∈R³是系统状态，代表产生的三维VCG数 据，来源于对心脏系统VCG信号的采样，k=1,…,m为采样时刻，T_s是采样周 期，F(V(k))=[f₁(V(k)),f₂(V(k)),f₃(V(k))]^T表示系统动态。VCG的ST-T环的内 在系统动态记为：

$F_{\mathrm{ST}}=(f_1\left(V\left(k\right)\right),f_2\left(V\left(k\right)\right),f_3\left(V\left(k\right)\right)){|}_{V\left(k\right)\in T_{\mathrm{ST}}}$

式中，k_S和k_TE分别代表心电图ST-T段的起始点（ST段 的起点）和结束点（T波的终点）。采用确定学习理论中的动态RBF神经网 络辨识器，可以实现对VCG的ST-T环的内在系统动态 $F_{\mathrm{ST}}=(f_1\left(V\left(k\right)\right),f_2\left(V\left(k\right)\right),f_3\left(V\left(k\right)\right)){|}_{V\left(k\right)\in T_{\mathrm{ST}}}$的局部准确神经网络逼近，即：

$\left(\begin{array}{c}{F}_{\mathrm{ST}}=(f_1\left(V\left(k\right)\right),f_2\left(V\left(k\right)\right),f_3\left(V\left(k\right)\right)){|}_{{V\left(k\right)\in V}_{\mathrm{ST}}}\\ \cong ({\bar{w}}_1^{T}S\left(V\left(k\right)\right),{\bar{W}}_2^{T}S\left(V\left(k\right)\right),{\bar{W}}_3^{T}S\left(V\left(k\right)\right)){|}_{V\left(k\right)\in V_{\mathrm{ST}}},\end{array}\right)$

式中，是常值神经网络权值向量，S(V(k))是高斯型径向基函数。

上述常值神经网络权值在本实施例中是指在一段时间内权值收敛至最优 值，选取收敛的部分权值取数值平均，获取相应的神经网络常值权值便于存储。将进行三维立体空间显示即 得到CDVG。CDVG是对学习到的系统动态知识的再利用，可反映心脏冠脉 血液动力学状况，对冠脉狭窄所导致的心肌缺血/心肌梗塞有较为敏感的检测 能力。本实施例中的PTB数据库数据P117的可视化建模结果CDVG如图2 所示。

（3）检测分析：

根据CDVG在三维空间中显示出的图像形态，结合年龄、病史和血压、 这些高危因素与就诊时的临床症状等信息进行综合分析，得到辅助检测结果。

本实施例中以选取的PTB数据库数据P117为例，PTB数据库内头文件 记载数据信息为：男性，37岁，无急性梗塞，无前期梗塞；正常。其CDVG 如图2所示，三维空间中的图像形态为规整环形。CDVG形态分析结果与头 文件内记载的正常结果一致。

本实施例中再以PTB数据库数据P072为例，PTB数据库内头文件记载 数据信息为：男性，40岁，数据采集时间为1992年1月17日，下侧壁心肌 梗塞，无前期梗塞；血脂蛋白II级，有抽烟史；血压120/70mmHg。采用基 于确定学习理论的心肌缺血辅助检测方法的步骤（1）（2），其滤波预处理之 后的ECG信号如图3所示，可视化建模结果CDVG如图4所示。利用步骤 （3），根据本发明方法观察其CDVG形态为无环散乱。CDVG形态分析结果 与头文件内记载的心肌梗塞诊断结果一致。

实施例验证了本发明方法的有效性和检测结果的可靠性，本发明的辅助 检查结果可作为心血管疾病诊断的重要参考，为此患者下一步病症诊断提供 依据。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述 实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、 修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护 范围之内。