CN111096740B - 心电信号分析方法及装置、信号记录仪和三维标测系统 - Google Patents

Abstract

本申请涉及一种心电信号分析方法及装置、信号记录仪、三维标测系统和计算机设备。所述方法包括：获取感兴趣区域或参考点/区域内的N个第一心内电信号以获得N个第一心内波形；根据N个第一心内波形生成第一比较信号频谱图；获取感兴趣点的M个第二心内电信号以获得M个第二心内波形；根据M个第二心内波形生成第二比较信号频谱图；对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生分析结果；根据分析结果输出关于感兴趣点是否为病灶点的提示信息。由此，通过对心内电信号进行频谱分析，并通过对不同区域的频谱进行比较判断，可以使得病灶点的判断更加准确，从而给医生提供更有效的信息，使得心律失常疾病的治愈率得到有效提高。

Description

心电信号分析方法及装置、信号记录仪和三维标测系统

技术领域

本申请涉及医疗技术领域，特别是涉及一种心电信号分析方法及装置、电生理信号记录仪、三维标测系统和计算机设备。

背景技术

目前，对于复杂性心律失常疾病，如房颤，非特异性室上速等，由于通常存在多个局灶或多条异常传播路径，造成心内信号混乱，导致医生难以对病灶点进行判断，现有计算机病灶辅助判断方法和系统能够提供的分析并不足以解决该问题，从而给疾病的治疗带来了很大的挑战。

发明内容

基于此，有必要针对在解决复杂性心律失常疾病时，由于心内信号混乱导致的病灶判断难的技术问题，提供一种心电信号分析方法及装置、电生理信号记录仪、三维标测系统和计算机设备。

本发明提供了一种心电信号分析方法，包括：

获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数；

根据N个第一心内波形生成第一比较信号频谱图；

获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数；

根据M个第二心内波形生成第二比较信号频谱图；

对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果；以及

根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。

本发明还提供了一种心电信号分析装置，包括：

第一获取模块，用于获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数；

第一生成模块，用于根据N个第一心内波形生成第一比较信号频谱图；

第二获取模块，用于获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数；

第二生成模块，用于根据M个第二心内波形生成第二比较信号频谱图；

第一分析模块，用于对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果；以及

输出模块，用于根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。

本发明还提供一种电生理信号记录仪，其包括上述的心电信号分析装置。

本发明还提供一种三维标测系统，其包括上述的心电信号分析装置以及显示模块，显示模块用于显示一心脏三维模型；其中输出模块将提示信息输出于显示模块，并显示于心脏三维模型上。

本发明还提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现上述的心电信号分析方法。

上述心电信号分析方法及装置、电生理信号记录仪、三维标测系统和计算机设备，通过获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，并根据N个第一心内波形生成第一比较信号频谱图，并获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，并根据M个第二心内波形生成第二比较信号频谱图，以及对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果，以及根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。由此，通过利用计算机系统对心内电信号进行频谱分析，并通过对感兴趣区域或参考点/区域和感兴趣点的频谱进行比较判断，可以使得病灶点的判断更加准确，从而给医生提供更有效的信息，使得复杂心律失常疾病的治愈率得到有效提高。

附图说明

图1为一个实施例中心电信号分析方法的应用环境图；

图2为一个实施例中心电信号分析方法的流程示意图；

图3为一个实施例中第一比较信号频谱图获取的流程示意图；

图4为一个实施例中第一单信号频谱图获取的流程示意图；

图5为一个实施例中对第一心内波形进行预处理的示意图；

图6为一个实施例中对信号频谱原图进行平滑处理的示意图；

图7为一个实施例中采用曲线累加生成第一比较信号频谱图的流程示意图；

图8为一个实施例中采用曲线累加生成第一比较信号频谱图的示意图；

图9为一个实施例中采用概率统计法生成第一比较信号频谱图的流程示意图；

图10为另一个实施例中第一比较信号频谱图获取的流程示意图；

图11为一个实施例中第一比较信号频谱图的示意图；

图12为一个实施例中感兴趣点判定的流程示意图；

图13为另一个实施例中感兴趣点判定的流程示意图；

图14为一个实施例中诊疗效果判定的流程示意图；

图15为一个实施例中对房颤进行心电信号分析的示意图；

图16为一个实施例中心电信号分析装置的结构框图；

图17为另一个实施例中心电信号分析装置的结构框图；

图18为又一个实施例中心电信号分析装置的结构框图；

图19为再一个实施例中心电信号分析装置的结构框图；

图20为一个实施例中电生理信号记录仪的结构框图；

图21为一个实施例中三维标测系统的结构框图；

图22为一个实施例中三维标测系统中三维心脏显示示意图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

在传统心律失常疾病诊疗中，医生会将与电生理信号记录仪或三维标测系统连接的标测导管放置在感兴趣区域，此时电生理信号记录仪或三维标测系统会通过标测导管取得感兴趣区域内的心内电信号，并将其呈现给医生，以便医生基于该心内电信号来判断疾病的灶局或异常传播路径。但是对于复杂心律失常疾病来说，通常具有多个灶局或多条异常传播路径，也就是说复杂心律失常疾病是多个异常灶局或异常传播路径共同作用的结果，目前的技术只关注于通过了解某一点或该区域的心内电信号状况来做判断，这常常导致判断结果不准确或不全面，消融区域广却不准，增加了病人的痛苦并降低了手术的成功率。而本申请的发明人通过对心内电信号的深入分析，发现需要有一种技术能够对这些异常进行识别和区分，进而进行分开对待，给医生提供更多量的信息，同时发明人又进一步创造性地发现通过对这些心内电信号进行频域处理，并且对比不同区域或不同点的频域或频率，可以在复杂的心内电信号中获得有规律的数据，进而根据这些数据可提高对心脏电生理行为的认识，进一步提高异常灶局或异常传播路径的判断的准确性和全面性。基于此，本申请提供了一种心电信号分析方法及装置、电生理信号记录仪、三维标测系统和计算机设备。

本申请提供的心电信号分析方法，可以应用于如图1所示的应用环境中。其中，控制器101可与一外部信号采集装置如电生理标测导管相连接，通过该标测导管对心内电信号进行采集，然后控制器101通过标测导管获取心内电信号。具体地，控制器101可通过第一标测导管102获取感兴趣区域内的N个第一心内电信号(图中的*表示采集点)，从而获得N个第一心内波形，并根据N个第一心内波形生成第一比较信号频谱图，或者，通过第一标测导管102获取参考点/区域内的N个第一心内电信号(图中的黑线线条表示采集点)，从而获得N个第一心内波形，并根据N个第一心内波形生成第一比较信号频谱图。控制器101还可通过第二标测导管104获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，并根据M个第二心内波形生成第二比较信号频谱图。最后，控制器101对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果，并根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。其中，N和M为大于等于1的整数，N和M可以相同也可以不同。

在一个实施例中，如图2所示，提供了一种心电信号分析方法，以该方法应用于图1中的控制器为例进行说明，包括以下步骤：

步骤202，获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数。

针对不同的心律失常疾病，感兴趣区域也不同，例如，房颤的感兴趣区域一般为左心房，室性早搏的感兴趣区域一般为右心室或左心室，室上速的感兴趣区域一般为右心房。参考区域一般为固定放置标测导管且在整个治疗过程中不会发生位置变动的区域，参考点可以为该区域内的某个点，一般房颤的参考区域为冠状窦，其参考点则为冠状窦中的一个点；室上速的参考区域为右心房，其参考点则为右心房中的一个点；房扑的区域为左心房/右心房，其参考点则为左心房/右心房中的一个点，具体医生可以根据病人的具体情况，选择在房颤中放置冠状窦导管，在室上速中放置右房标测导管，或在房扑中放置房扑大环标测导管，以获取参考点/区域内的心内电信号。

在获取感兴趣区域或参考点/区域内的第一心内电信号时，可获取一个或多个第一心内电信号，从而获得一个或多个第一心内波形，具体个数可根据感兴趣区域或参考点/区域的大小进行选择，其中，第一心内电信号的个数越多，越能从整体上了解感兴趣区域或参考点/区域内的心内情况。举例来说，房颤的感兴趣区域一般在左心房，要想从整体上了解整个左心房的特征，可在左心房设置多个采集点(如图1中的多个采集点*)，然后对每个采集点进行第一心内电信号的采集以获得心内波形，记为第一心内波形，假设采集点为50，那么可获得50个第一心内波形。

具体地，可通过标测导管进行心内电信号的采集。在电生理标测技术发展至今，标测导管已经发展到了具有多个电极的配置，可以在同一时间同时采集多个地方的心内电信号，例如，同一时间可以同时采集64路心内电信号，从而在一次采集中可获得64个心内波形。这样，当需要采集N个第一心内波形时，可通过电极数为1的标测导管进行N次采集，也可以通过电极数为a的标测导管进行N/a次采集，具体可根据实时处理的要求进行选择。举例来说，当需要了解整个左心房的特征时，如果使用电极数为1的标测导管，则可能需要采集50-200次，而如果使用电极数为10的标测导管，可能只需要采集5-20次，对于后者，虽然采集的次数较少，但是同时处理的负载高，数据处理的实时性就相对差一些，所以，优选采用电极数为4-10的标测导管，以在满足实时性要求的同时，减少采集的次数。

以电极数为5的第一标测导管为例，在进行一次采集时，医生可将第一标测导管放置在左心房的某一位置，此时控制器通过第一标测导管开始采集该位置处的心内电信号，并在采集时间达到预设时间时停止采集，至此完成一次心内电信号的采集。由于第一标测导管的电极数为5，所以本次采集完，可获得5个第一心内波形，即5个不同采集点的心内电信号的时序波形。接着，医生可将第一标测导管放置在左心房的另一位置，此时控制器通过第一标测导管进行第二次采集，并再次获得5个第一心内波形，依次类推，直至完成所有的采集。

其中，采集时间为分析心内电信号时序的时间，即第一心内波形的时长，范围可以为1s-10s，优选为10s，即以每10s的心内电信号作为一个单元，以形成一个第一心内波形，用于后续频谱分析。在实际应用中，可在控制器中设置一个默认采集时间如10s，同时可设置一个分析窗口，以便医生通过该分析窗口对采集时间进行调整，以满足实际需求。

另外，标测导管采集的心内电信号可以为双极信号，也可以为单极信号。由于双极信号可以有效抑制共模信号干扰，所以一般优选采用双极信号，这样可以保证采集的心内电信号不受环境噪声干扰，使得采集的心内波形更加准确。

步骤204，根据N个第一心内波形生成第一比较信号频谱图。

具体地，在完成感兴趣区域或参考点/区域的第一心内波形的采集后，可对采集的N个第一心内波形进行频域转换、组合处理等，以生成感兴趣区域或参考点/区域的第一比较信号频谱图。

步骤206，获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数。

感兴趣点一般为疾病治疗上通过经验总结得到的容易出现问题的地方，例如，房颤的感兴趣点可优选为肺静脉中的一个点，房扑的感兴趣点可优选为三尖瓣峡部中的一个点。

在获取感兴趣点的第二心内电信号时，可获取一个或多个第二心内电信号，从而获得一个或多个第二心内波形，具体个数可根据感兴趣点的大小进行选择。具体地，仍以房颤为例。在获得感兴趣区域或参考点/区域的第一比较信号频谱图之后，医生可将第二标测导管放置在肺静脉，此时控制器通过第二标测导管开始采集肺静脉处的心内电信号，并在采集时间达到预设时间时停止采集，至此完成一次心内电信号的采集，其中，当第二标测导管的电极数为1时，一次可以获得1个第二心内波形；当第二标测导管的电极数为5时，一次可以获得5个第二心内波形。

步骤208，根据M个第二心内波形生成第二比较信号频谱图。

具体地，在完成感兴趣点的第二心内波形的采集后，可对采集的M个第二心内波形进行频域转换、组合处理等，以生成感兴趣点的第二比较信号频谱图。

步骤210，对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果。

步骤212，根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。

具体地，在获得感兴趣区域或参考点/区域的第一比较信号频谱图和感兴趣点的第二比较信号频谱图之后，可对两个频谱图进行比较分析，例如，对频谱图的能量以及曲线形状等进行比较分析，以输出关于感兴趣点是否为病灶点的提示信息，例如，可以以文字、图像或语音等方式输出关于感兴趣点是否为病灶点的提示信息，具体采用何种方式，可根据实际需求选择。

本实施例中，通过对心内电信号进行频谱分析，并通过对感兴趣区域或参考点/区域和感兴趣点的频谱进行比较判断，即对全局和局部的频谱进行比较，给出感兴趣点是否为病灶点的提示信息，可以有效解决由于心内信号混乱导致的疾病判断难的问题，而且相较于仅通过感兴趣点的频谱给出病灶点的判断结果，其准确度更高，从而给疾病的治疗提供更有利的信息。

在一个实施例中，如图3所示，根据N个第一心内波形生成第一比较信号频谱图的步骤进一步包括：

步骤302，生成每个第一心内波形的第一单信号频谱图，以获得N个第一单信号频谱图。

在一个实施例中，生成每个第一心内波形的第一单信号频谱图，以获得N个第一单信号频谱图的步骤进一步包括：对每个第一心内波形进行频域转换，以获得N个第一单信号频谱图。例如，可以采用快速傅里叶变换或者小波变换等对每个第一心内波形进行频域转换，以获得N个第一单信号频谱图。

具体地，可以在完成对N个第一心内电信号的采集之后，再统一生成N个第一心内波形以及对应N个第一心内波形的N个单信号频谱图，记为N个第一单信号频谱图(即，每个第一心内波形对应的信号频谱图称之为第一单信号频谱图)，也可以在每次采集完成后，即对采集到的一个第一心内电信号进行处理，生成第一心内波形以及对应第一心内波形的第一单信号频谱图。对于后者，可以有效减少数据集中处理导致的整个分析时间长的问题，同时可保证信号频谱图的实时性，所以优选采用后者。

举例来说，当第一标测导管的电极数为5时，在通过该标测导管完成一次采集并获得5个第一心内波形后，开始对获得的5个第一心内波形的每个第一心内波形进行频域转换，以获得5个第一单信号频谱图，并在下一次采集完成后，对下一次采集的5个第一心内波形进行频域转换，以获得5个第一单信号频谱图，依次类推，直至完成所有的采集，最终可获得N个第一单信号频谱图。

在一个实施例中，在对每个第一心内波形进行频域转换之前还包括：对每个第一心内波形进行预处理，以使得每个第一心内波形在预设频率范围内，其中，预设频率范围可根据实际情况进行选择设置，例如可以为0-20Hz；和/或，在对每个第一心内波形进行频域转换之后还包括：对每个第一单信号频谱图进行平滑处理。

也就是说，在对第一心内波形进行频域转换之前，还可以对第一心内波形进行预处理，例如，预处理可包括QRS波减除、带通或高通滤波处理、取模运算和低通滤波处理中的任一种或任意组合，然后再对预处理后的第一心内波形进行频域转换，以获得第一单信号频谱图。另外，在获得第一单信号频谱图之后，还可对获得的第一单信号频谱图进行平滑处理，例如，可通过低通滤波、互补滤波或者卡尔曼滤波等方式进行平滑处理，以获得较为平滑的第一单信号频谱图。对于预处理和平滑处理可以择一使用，也可以同时使用。

作为一个具体示例，如图4所示，生成每个第一心内波形的第一单信号频谱图，以获得N个第一单信号频谱图的步骤可包括：

步骤402，对每个第一心内波形进行预处理，以获得预设频率范围内的N个经预处理的第一心内波形，其中预处理可包括：QRS波减除、带通或高通滤波处理、取模运算和低通滤波处理中的任一种或任意组合。

具体地，在获得第一心内波形之后，可先对第一心内波形进行预处理，例如，参考图5所示，可先对第一心内波形进行频率为40-250Hz的带通滤波处理，然后对带通滤波处理后的第一心内波形进行取模运算，如将负心内电信号取反，以使所有心内电信号为正，最后再对取模运算后的第一心内波形进行0-20Hz的低通滤波处理，以获得预设频率范围0-20Hz内的经预处理的第一心内波形。

步骤404，对每个经预处理的第一心内波形进行频域转换，以获得N个第一单信号频谱原图，其中可以采用快速傅里叶变换对预处理后的第一心内波形进行频域转换，这里将频域转换后的信号频谱图记为第一单信号频谱原图。

步骤406，对每个第一单信号频谱原图进行平滑处理，以获得N个第一单信号频谱图，例如，参考图6所示，可以通过低通滤波方式对第一单信号频谱原图进行平滑处理，以获得较为平滑的信号频谱图，记为第一单信号频谱图。

步骤304，对N个第一单信号频谱图进行组合，生成第一比较信号频谱图，或者，自N个第一单信号频谱图中选择一个第一单信号频谱图，作为第一比较信号频谱图。

也就是说，可以直接从N个第一单信号频谱图中获取一个第一单信号频谱图，以作为第一比较信号频谱图(特别的，当仅有一个第一单信号频谱图时，该第一单信号频谱图即作为第一比较信号频谱图)，也可以将N个第一单信号频谱图进行组合，并将组合后的信号频谱图作为第一比较信号频谱图。对于后者，由于能够从整体上体现感兴趣区域或参考点/区域的心内情况，所以优选采用后者。

具体地，当采用对N个第一单信号频谱图进行组合以生成第一比较信号频谱图时，可以在完成对N个第一单信号频谱图的获取之后，开始对N个第一单信号频谱图进行组合，以生成第一比较信号频谱图，也可以在每次采集完成后，在获得相应的第一单信号频谱图之后，即对已经获得的第一单信号频谱图进行组合。对于后者，可以有效减少数据集中处理导致的整个分析时间长的问题，同时可保证信号频谱图的实时性，所以优选采用后者。

举例来说，当第一标测导管的电极数为5时，在通过该标测导管完成一次采集并获得5个第一心内波形后，开始对获得的5个第一心内波形的每个第一心内波形进行频域转换，以获得5个第一单信号频谱图，并对获得的5个第一单信号频谱图进行组合处理，以获得一个信号频谱图。然后，在下一次采集完成后，先对采集的5个第一心内波形进行频域转换，以获得5个第一单信号频谱图，再将这5个第一单信号频谱图组合到上一次获得的信号频谱图上，以获得一个新的信号频谱图。依次类推，直至完成对N个第一单信号频谱图的组合，并最终获得第一比较信号频谱图。

在一个实施例中，对N个第一单信号频谱图进行组合，生成第一比较信号频谱图的步骤进一步包括：对N个第一单信号频谱图进行累加处理，其中，累加处理可包括曲线累加法或概率统计法。也就是说，可以采用曲线累加法或概率统计法对N个第一单信号频谱图进行累加处理，以生成第一比较信号频谱图。

在一个实施例中，如图7所示，采用曲线累加法对N个第一单信号频谱图进行累加处理，生成第一比较信号频谱图的步骤进一步包括：

步骤702，对每个第一单信号频谱图进行能量归一化处理。

在一个实施例中，在对每个第一单信号频谱图进行能量归一化处理时，可先对每个第一单信号频谱图进行特征提取，如频率以及与频率对应的能量，并根据频率对应的能量从大到小对频率进行排序，并将能量最大的频率定义为第一频率(主频)，对应的能量定义为第一频能量(主频能量)，同时将能量较小的频率去除，然后再进行能量归一化处理。

举例来说，在获得N个第一单信号频谱图之后，可对每个第一单信号频谱图的特征进行提取，如第一频率(主频)、第一频能量(主频能量)、第二频率(次频)、第二频能量(次频能量)、第三频率、第三频能量、...、第i频率、第i频能量。其中，i按能量的大小依次进行标记，最大能量处的频率为第一频率即主频，第二大能量处的频率为第二频率即次频，依次类推，并且i根据疾病的复杂程度以及相关能量进行综合考虑，一般疾病越复杂i越大，同时有些频率的能量过低，如低于第一频能量的10％，则可以忽略不计。然后，对每个第一单信号频谱图的第一频能量、第二频能量、第三频能量、...、第i频能量进行能量归一化处理，即将能量限制在1以内。

步骤704，对能量归一化处理后的N个第一单信号频谱图进行同频率能量累加，以获得第一比较信号频谱原图。

其中，同频率能量累加是指不同第一单信号频谱图之间、相同频率对应的能量累加。例如，在N个第一单信号频谱图中，有两个信号频谱图的第一频率为5Hz，那么同频率能量累加后，频率5Hz对应的能量是这两者的能量之和；又如，在N个第一单信号频谱图中，有三个信号频谱图的第一频率为5Hz，那么同频率能量累加后，频率5Hz对应的能量是这三者的能量之和。按照该方式，完成所有频率的能量累加，最终可以形成一个新的信号频谱图，记为第一比较信号频谱原图。

步骤706，对第一比较信号频谱原图进行能量归一化处理，以获得第一比较信号频谱图。

具体地，在获得第一比较信号频谱原图之后，还对该第一比较信号频谱原图进行能量归一化处理，具体过程与对每个第一单信号频谱图的能量归一化处理过程相同。例如，可先对第一比较信号频谱原图的特征进行提取，如第一频率(主频)、第一频能量(主频能量)、第二频率(次频)、第二频能量(次频能量)、第三频率、第三频能量、...、第j频率、第j频能量。其中，j按能量的大小依次进行标记，最大能量处的频率为第一频率即主频，第二大能量处的频率为第二频率即次频，依次类推。然后，对提取的第一频能量、第二频能量、第三频能量、...、第j频能量进行能量归一化处理，以获得第一比较信号频谱图。

作为一个具体示例，当第一标测导管的电极数为5时，在通过该标测导管完成一次采集并获得5个第一心内波形后，开始对获得的5个第一心内波形的每个第一心内波形进行频域转换，以获得5个第一单信号频谱图，并对获得的5个第一单信号频谱图的每个第一单信号频谱图进行能量归一化处理，以及对归一化处理后5个第一单信号频谱图进行同频率能量累加并进行能量归一化处理，以获得一个信号频谱图。接着，在下一次采集完成后，又获得5个第一单信号频谱图，同样进行5个第一单信号频谱图的能量归一化处理，并与上一次获得的信号频谱图进行同频率能量累加并进行能量归一化处理，以获得一个新的信号频谱图。依次类推，直至完成对N个第一单信号频谱图的同频率能量累加并进行归一化处理，最终获得第一比较信号频谱图，如图8所示。至此，完成了感兴趣区域或参考点/区域的第一比较信号频谱图的获取。

在另一个实施例中，如图9所示，采用概率统计法对N个第一单信号频谱图进行累加处理，生成第一比较信号频谱图的步骤进一步包括：

步骤902，获取每个第一单信号频谱图的频率，并获取每个频率在对应第一单信号频谱图中的权重。

其中，频率的获取过程与图7和图8对应实施例的获取过程相同，具体这里不再赘述。通过上述方式在获得每个第一单信号频谱图的频率之后，如第一频率(主频)、第二频率(次频)、第三频率、...、第i频率，获取每个频率在对应第一单信号频谱图中的权重，其中，权重可以根据频率面积、频率的谐频面积以及整个频谱图面积进行确定，例如，第一频率的权重等于第一频率面积与第一频率的谐频面积之和再除以整个频谱面积，如表1所示，在第一个第一单信号频谱图中，经计算后，第一频率在第一单信号频谱图中的权重为0.33。

步骤904，计算N个第一单信号频谱图中相同频率下的权重平均值。

其中，相同频率下的权重平均值是指不同第一单信号频谱图之间、相同频率对应的权重的平均值。参考表1所示，针对第一频率，其权重平均值＝(0.33+0.35+...+0.24)/15＝0.28，其中，在第一频率中，频率1.8对应的权重忽略不计，也就是说，在求取权重平均值时，如果频率相同(如，均为第一频率)但实际频率值(如6.4、1.8)不同，则相应的权重用零取代。

步骤906，根据权重平均值，获得第一比较信号频谱图。

具体地，如表1所示，在进行权重平均值求取后，可获得第一比较信号频谱图的第一频率为6.40、第一频率对应的权重为0.28、第二频率为13.20、第二频率对应的权重为0.10、...，由此可获得第一比较信号频谱图。

表1

作为一个具体示例，参考表1所示，当第一标测导管的电极数为5时，在通过该标测导管完成一次采集并获得5个第一心内波形后，开始对获得的5个第一心内波形的每个第一心内波形进行频域转换，以获得5个第一单信号频谱图，并对获得的5个第一单信号频谱图的每个第一单信号频谱图的特征进行提取，如对第一频率、第二频率、第三频率、第四频率和第五频率进行提取，并计算每个频率在对应第一单信号频谱图中的权重，例如，在采集点1处，第一频率6.40Hz对应的权重为0.33，第二频率9.2Hz对应的权重为0.18。然后，将5个第一频率6.4Hz对应的权重相加并进行平均。接着，在下一次采集完成后，又获得5个第一单信号频谱图，同样处理后进行累加平均，依次类推，最终得到一个能量由大到小分布的第一比较信号频谱图。至此，完成了感兴趣区域或参考点/区域的第一比较信号频谱图的获取。

图10为一个实施例中第一比较信号频谱图的获取流程图，如图10所示，第一比较信号频谱图的获取过程可包括以下步骤：

步骤1002，设置分析窗口，如分析窗口为10s。

步骤1004，获取感兴趣区域或参考点/区域内的单点或多点心内电信号，以获得单个或多个第一心内波形。

步骤1006，对单个或多个第一心内波形进行预处理，如QRS波减除、带通或高通滤波处理、取模运算和低通滤波处理中的任一种或任意组合。

步骤1008，对预处理后的第一心内波形进行频域转换，如进行快速傅里叶变换，以获得第一单信号频谱原图。

步骤1010，对第一单信号频谱原图进行平滑处理以获得单个或多个第一单信号频谱图，如采用低通滤波方式进行平滑处理。

步骤1012，提取每个第一单信号频谱图的特征：如主频、次频。

步骤1014，计算感兴趣区域或参考点/区域的整体特征：如主频、次频。

步骤1016，是否继续采样。如果继续采样，则返回执行步骤S1004；否则，采样结束，至此完成对感兴趣区域或参考点/区域的第一比较信号频谱图的获取。

在一个实施例中，第一比较信号频谱图和第二比较信号频谱图的生成方式相同。

具体地，在根据M个第二心内波形生成第二比较信号频谱图时，其步骤可进一步包括：生成每个第二心内波形的第二单信号频谱图，以获得M个第二单信号频谱图；以及对M个第二单信号频谱图进行组合，生成第二比较信号频谱图，或者，自M个第二单信号频谱图中选择一个第二单信号频谱图，作为第二比较信号频谱图。其中，生成每个第二心内波形的第二单信号频谱图的步骤与生成每个第一心内波形的第一单信号频谱图的步骤相同；对M个第二单信号频谱图进行组合，生成第二比较信号频谱图的步骤与对N个第一单信号频谱图进行组合，生成第一比较信号频谱图的步骤相同。特别的，当仅有一个第二单信号频谱图时，该第二单信号频谱图即作为第二比较信号频谱图。

也就是说，在获取感兴趣点的第二比较信号频谱图时，其获取过程可参照感兴趣区域或参考点/区域的第一比较信号频谱图的获取。例如，在有多个第二心内波形的情况下，对于感兴趣点的第二比较信号频谱图的获取，可参照步骤1002-步骤1016，具体这里就不再赘述。

在一个实施例中，对N个第一单信号频谱图进行组合，生成第一比较信号频谱图的步骤进一步包括：对N个第一单信号频谱图进行叠加处理。其中，叠加处理是指将N个第一单信号频谱图进行简单叠加，例如，图11为将5个第一单信号频谱图进行叠加后产生的信号频谱图，记为第一比较信号频谱图。

上述实施例中，既可以直接从获取的单信号频谱图中选择一个单信号频谱图以作为比较信号频谱图，具体方式可以是由医生做出选择，给到计算机系统一个选择指令，从而选择出一个单信号频谱图；也可以通过累加或叠加的方式对单信号频谱图进行组合，以获得比较信号频谱图，而具体选用哪种方式可根据实际需要进行选择。

在一个实施例中，如图12所示，对第一比较信号频谱图和第二比较信号频谱图进行比较分析的步骤进一步包括：

步骤1202，对第二比较信号频谱图和第一比较信号频谱图进行相似度分析，以获得相似度值，例如，在进行相似度分析时，可以采用皮尔森相关系数表示相似度值，也可以采用其他方式进行相似度分析。

步骤1204，判断相似度值与预设相似度阈值的大小关系，其中，预设相似度阈值可以根据实际情况进行选择，例如相似度阈值可以为0.6。

具体地，可以通过比较第一比较信号频谱图和第二比较信号频谱图的曲线相似度来判断感兴趣点是否为病灶点，例如，可以对第一比较信号频谱图和第二比较信号频谱图进行曲线相似度分析，相似度分析可以采用皮尔森相关系数表示。如果分析的相似度值超过预设相似度阈值如0.6，则可以判定感兴趣点为病灶点，此时可以以文字、图像或语音的方式输出关于感兴趣点为病灶点的提示信息，例如，可将相似度值以及初步判断结果以文本方式进行显示，也可以以语音方式进行播放，还可以在第二比较信号频谱图中将相似度较高的部分进行标识，以对医生进行提醒，具体采用何种方式可根据实际需求选择；如果分析的相似度值未超过预设相似度阈值如0.6，则可以判定感兴趣点为非病灶点，此时可以以文字、图像或语音的方式输出关于感兴趣点为非病灶点的提示信息。

在一个实施例中，感兴趣点为多个，从而生成多个第二比较信号频谱图，对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果的步骤还包括将第一比较信号频谱图和多个第二比较信号频谱图分别进行比较分析，从而产生比较分析结果。

可以理解的是，也可以在通过上述方式获得一个感兴趣点的第二比较信号频谱图之后，继续获取其它感兴趣点的第二比较信号频谱图，然后从获得的多个第二比较信号频谱图中获取一个与第一比较信号频谱图相似度最高的第二比较信号频谱图，并将该第二比较信号频谱图对应的感兴趣点作为病灶点。

本实施例中，通过分析感兴趣区域或参考点/区域内的信号频谱图与感兴趣点的信号频谱图的对应关系，即通过分析全局与局部的频谱对应关系，找出相关规律，从而从整体上给出病灶点的判断结果，给医生提供医疗上的建议，帮助医生获得更有价值的心内信息，有利于提高消融的准确性，从而提高心律失常疾病的治愈率，同时也避免了患者的不必要的损伤。

在一个实施例中，如图13所示，对第一比较信号频谱图和第二比较信号频谱图进行比较分析的步骤进一步包括：

步骤1302，计算第二比较信号频谱图中每个频率对应的能量与第一比较信号频谱图中相同频率对应的能量之间的能量差值，以获得多个能量差值。

步骤1304，判断多个能量差值与预设能量阈值的大小关系，其中，预设能量阈值可根据实际情况进行选择，例如，预设能量阈值可以为感兴趣区域归一化后第一频率对应的能量。

具体地，可以通过比较第一比较信号频谱图和第二比较信号频谱图的频率对应的能量来判断感兴趣点是否为病灶点，例如，当第二比较信号频谱图中某一频率(如主频)对应的能量与第一比较信号频谱图中相同频率(如主频)对应的能量之间的差值超过预设能量阈值，则可以判定感兴趣点为病灶点，此时可以以文字、图像或语音的方式输出关于感兴趣点为病灶点的提示信息，例如，可将感兴趣点的频率能量以及初步判断结果以文本方式进行显示，也可以以语音方式进行播放，还可以在第二比较信号频谱图中将能量较高的部分进行标识，以对医生进行提醒，具体采用何种方式可根据实际需求选择；否则，判定感兴趣点为非病灶点，此时可以以文字、图像或语音的方式输出关于感兴趣点为非病灶点的提示信息。

需要说明的是，在判断感兴趣点是否为病灶点时，也可以将上述两种方式进行结合，即，既考虑一定的频率能量，也考虑频谱曲线的相似度，例如，可以通过给予一定的权重对两者进行综合分析，具体这里不再详述。

本实施例中，通过分析感兴趣区域或参考点/区域内的信号频谱图与感兴趣点的信号频谱图的对应关系，即通过分析全局与局部的频谱对应关系，找出相关规律，从而从整体上给出病灶点的判断结果，给医生提供医疗上的建议，帮助医生获得更有价值的心内信息，有利于提高心律失常疾病的治愈率。

进一步地，在对感兴趣点判断完成后，可将判断结果呈现给医生，这样医生可以基于判断结果，结合第一比较信号频谱图和第二比较信号频谱图做进一步的比较判断，以确认是否将其作为病灶进行治疗。通常，心律失常疾病的治疗手段是使心肌发生物理性破坏，例如，可以通过射频消融能量使得心肌失活，或者通过冷冻、激光以及超声等方式使心肌失活，从而达到心律失常疾病治疗的目的。

在一个实施例中，如图14所示，心电信号分析方法还包括：

步骤1402，获取病灶点被治疗后感兴趣区域或参考点/区域内的K个第三心内电信号，从而获得K个第三心内波形，其中K为大于等于1的整数。

步骤1404，根据K个第三心内波形生成第三比较信号频谱图(将治疗后感兴趣区域或参考点/区域对应的信号频谱图称之为第三比较信号频谱图)。

在一个实施例中，第一比较信号频谱图和第三比较信号频谱图的生成方式相同。

具体地，在根据K个第三心内波形生成第三比较信号频谱图时，其步骤可进一步包括：生成每个第三心内波形的第三单信号频谱图，以获得K个第三单信号频谱图；以及对K个第三单信号频谱图进行组合，生成第三比较信号频谱图，或者，自K个第三单信号频谱图中选择一个第三单信号频谱图，作为第三比较信号频谱图。其中，生成每个第三心内波形的第三单信号频谱图的步骤与生成每个第一心内波形的第一单信号频谱图的步骤相同；对K个第三单信号频谱图进行组合，生成第三比较信号频谱图的步骤与对N个第一单信号频谱图进行组合，生成第一比较信号频谱图的步骤相同。特别的，当仅有一个第三单信号频谱图时，该第三单信号频谱图即作为第三比较信号频谱图。

具体生成过程可参考第一比较信号频谱图的生成过程，这里就不再赘述

步骤1406，将第一比较信号频谱图和第二比较信号频谱图中的至少一个与第三比较信号频谱图进行比较分析，以提供病灶点的治疗效果的信息，具体可参考对感兴趣点的分析过程。

具体而言，在病灶点治疗完成后，可通过比较治疗前后感兴趣区域或参考点/区域的比较信号频谱图来确定治疗效果，由于参考点/区域在整个治疗过程中位置不会变化，所以在一次治疗完成后，优选地，对治疗前后参考点/区域的比较信号频谱图进行比较判断，以判断治疗是否有效，并根据治疗是否有效做进一步处理。一般情况下，在进行治疗后，会出现以下几种情况：1)在参考点/区域的第三比较信号频谱图上，某一频率能量(如主频能量)消失，心律失常恢复窦律；2)在参考点/区域的第三比较信号频谱图上，某一频率能量(如主频能量)减弱，没有恢复窦律，可能还有其他感兴趣点具有频率能量，继续使用标测导管进行查找；3)在参考点/区域的第三比较信号频谱图上，某一频率能量(如主频能量)没有变化，没有恢复窦律，继续使用标测导管进行查找，或者治疗没有其他隔离作用，继续优化消融。不断循环，直至完成对心律失常疾病的治疗。

也就是说，复杂心律失常疾病通常可能存在多个驱动频率，所以可以先解决主要矛盾，即主频，如果疾病仍然复发，次频就成为主要矛盾，此时次频上升为主频，进入下一循环，依次类推，直至心律失常疾病治疗完成，恢复窦律。

以房颤治疗为例。通常房颤的感兴趣区域在左心房，如图15所示，可对左心房进行标测，并获取感兴趣区域的比较信号频谱图。在房颤手术中，冠状窦放置的标测导管在整个手术过程中保持不动，所以可以采用某一心内电信号作为参考点，并获取参考点的比较信号频谱图。房颤的发病局灶有很大概率来自于4条肺静脉，所以可以优先对这些区域进行标测，即获得感兴趣点的比较信号频谱图。

在本实施例中，可以采用电极数为20的环状标测导管，该标测导管通过20个电极可以同时采集10组心内电信号，在肺静脉区域可分成10个点，分别进行频谱分析，得到10个感兴趣点的比较信号频谱图。然后，分别将10个感兴趣点的比较信号频谱图与感兴趣区域的比较信号频谱图进行比较判断，如进行曲线相似度比较或频率能量比较，以确定感兴趣点是否为病灶点，如图15所示，虚线内的感兴趣点的比较信号频谱图与感兴趣区域的比较信号频谱图相似度最高，所以该点可判定为病灶点。此时，医生可根据判断结果对该点进行治疗，并在治疗完成后，再次获取参考点的比较信号频谱图，并与治疗前的参考点的比较信号频谱图进行比较判断，如进行曲线相似度比较或频率能量比较，以确定治疗效果。

在一次治疗完成后，如果治疗效果欠佳，则更换感兴趣点，并继续治疗，直至治疗后感兴趣区域或参考点/区域的比较信号频谱图上的频率能量相对于治疗前感兴趣区域或参考点/区域的比较信号频谱图的能量消失，心律失常恢复窦律。

本实施例中，通过对全局和局部的频谱关系进行对比分析，可以从整体上了解心内频谱情况，从而可以更加精确的找到病灶点的位置，在提高心律失常疾病治愈率的同时，可以有效减少对健康心肌组织造成的损害。

应该理解的是，虽然图1-14的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1-14中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

在一个实施例中，如图16所示，提供了一种心电信号分析装置100，包括：第一获取模块110、第一生成模块120、第二获取模块130、第二生成模块140、第一分析模块150和输出模块160，其中：

第一获取模块110可与一外部信号采集装置(如电生理标测导管)相连接，用于获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数；

第一生成模块120用于根据N个第一心内波形生成第一比较信号频谱图；

第二获取模块130可与一外部信号采集装置(如电生理标测导管)相连接，用于获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数；

第二生成模块140用于根据M个第二心内波形生成第二比较信号频谱图；

第一分析模块150用于对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果；

输出模块160用于根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。

在一个实施例中，如图17所示，第一生成模块120进一步包括：

第一处理单元121，用于生成每个第一心内波形的第一单信号频谱图，以获得N个第一单信号频谱图；

第二处理单元122，用于对N个第一单信号频谱图进行组合，生成第一比较信号频谱图，或者，自N个第一单信号频谱图中选择一个第一单信号频谱图，作为第一比较信号频谱图。

在一个实施例中，如图17所示，第一处理单元121进一步包括：

转换单元1211，用于对每个第一心内波形进行频域转换，以获得N个第一单信号频谱图。

在一个实施例中，如图17所示，第一处理单元121还进一步包括：

预处理单元1212，用于在转换单元1211对每个第一心内波形进行频域转换之前，对每个第一心内波形进行预处理，以使得每个第一心内波形在预设频率范围内；

平滑处理单元1213，用于在转换单元1211对每个第一心内波形进行频域转换之后，对每个第一单信号频谱图进行平滑处理。

在一个实施例中，如图17所示，第二处理单元122进一步包括：

累加处理单元1221，用于对N个第一单信号频谱图进行累加处理，生成第一比较信号频谱图，其中，累加处理包括曲线累加法或概率统计法。

在一个实施例中，如图17所示，第二处理单元122进一步包括：

叠加处理单元1222，用于对N个第一单信号频谱图进行叠加处理，生成第一比较信号频谱图。

在一个实施例中，如图17所示，第一分析模块150进一步包括：

相似度分析单元151，用于对第二比较信号频谱图和第一比较信号频谱图进行相似度分析，以获得相似度值；

第一判断单元152，用于判断相似度值与预设相似度阈值的大小关系。

在一个实施例中，如图18所示，第一分析模块150进一步包括：

计算单元153，用于计算第二比较信号频谱图中每个频率对应的能量与第一比较信号频谱图中相同频率对应的能量之间的能量差值，以获得多个能量差值；

第二判断单元154，用于判断多个能量差值与预设能量阈值的大小关系。

在一个实施例中，输出模块160以文字、图像或语音的方式输出关于感兴趣点是否为病灶点的提示信息。

在一个实施例中，如图19所示，心电信号分析装置100还包括：

第三获取模块170，用于获取病灶点被治疗后感兴趣区域或参考点/区域内的K个第三心内电信号，从而获得K个第三心内波形，其中K为大于等于1的整数；

第三生成模块180，用于根据K个第三心内波形生成第三比较信号频谱图；

第二分析模块190，用于将第一比较信号频谱图和第二比较信号频谱图中的至少一个与第三比较信号频谱图进行比较分析，以提供病灶点的治疗效果的信息。

在一个实施例中，第二生成模块140和/或第三生成模块180与第一生成模块120为同一功能模块；第二获取模块130和/或第三获取模块170与第一获取模块110为同一功能模块；第一分析模块150与第二分析模块190为同一功能模块。

关于心电信号分析装置100的具体限定可以参见上文中对于心电信号分析方法的限定，在此不再赘述。上述心电信号分析装置100中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中，提供了一种电生理信号记录仪1000，如图20所示，电生理信号记录仪1000包括上述的心电信号分析装置100。

在一个实施例中，提供了一种三维标测系统2000，如图21所示，三维标测系统2000包括上述的心电信号分析装置100和显示模块200，显示模块200用于显示一心脏三维模型，其中，输出模块160将提示信息输出于显示模块200，并显示于心脏三维模型上。

在一个实施例中，可将心电信号分析装置100嵌入至三维标测系统2000中，以心脏三维图形区域化显示的方式给医生提供更有利的信息。参考图22所示，心电信号分析装置100获取并保存感兴趣区域内所有采集点的第一单信号频谱图，并进行特征提取，例如，特征优先为第一频率、第一频能量、第二频率、第二频能量、第三频率、第三频能量、第四频率、第四频能量、第五频率、第五频能量，并对提取的能量进行归一化处理。然后，医生选择所需频率，此时心电信号分析装置100将所有采集点中，只要出现过该频率的采集点对应的心内波形，按照能量的高低，显示在左心房的三维构建图上。例如，当医生选择的频率为6.4Hz时，所有采集过的点，只要出现6.4Hz频率的心内波形，就按照能量的高低，共同显示在三维标测系统的显示界面上，以便医生方便地了解不同频率在心脏内的分布情况，有利于做出准确的判断，进一步地，心电信号分析装置的输出模块将提示信息输出于所述显示界面，并以颜色不同或深浅等方式将可能被病灶区显示于所述心脏三维模型上。

本实施例中，通过将心电信号分析装置嵌入到三维标测系统中，可以增加三维标测系统中的频谱分析功能，从而增加产品的核心竞争力。

在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：获取感兴趣区域或参考点/区域内的第一心内电信号，从而获得第一心内波形；根据第一心内波形生成第一比较信号频谱图；获取感兴趣点的第二心内电信号，以获得第二心内波形；根据第二心内波形生成第二比较信号频谱图；对第一比较信号频谱图和第二比较信号频谱图进行比较分析并产生比较分析结果；以及根据比较分析的结果，输出关于感兴趣点是否为病灶点的提示信息。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (21)

1.一种心电信号分析方法，其特征在于，所述的心电信号分析方法包括：

获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数；

根据N个所述第一心内波形生成第一比较信号频谱图；

获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数；

根据M个所述第二心内波形生成第二比较信号频谱图；

对所述第一比较信号频谱图和所述第二比较信号频谱图进行比较分析并产生比较分析结果；以及

根据比较分析的结果，输出关于所述感兴趣点是否为病灶点的提示信息；

所述根据N个所述第一心内波形生成第一比较信号频谱图，包括：

生成每个所述第一心内波形的第一单信号频谱图，以获得N个所述第一单信号频谱图；

对N个所述第一单信号频谱图进行累加处理或叠加处理，生成所述第一比较信号频谱图，或者，自N个所述第一单信号频谱图中选择一个所述第一单信号频谱图，作为所述第一比较信号频谱图。

2.根据权利要求1所述的心电信号分析方法，其特征在于，所述生成每个所述第一心内波形的第一单信号频谱图，以获得N个所述第一单信号频谱图的步骤进一步包括：对每个所述第一心内波形进行频域转换，以获得N个所述第一单信号频谱图。

3.根据权利要求2所述的心电信号分析方法，其特征在于，其中，

在对每个所述第一心内波形进行频域转换之前还包括：对每个所述第一心内波形进行预处理，以使得每个所述第一心内波形在预设频率范围内；和/或，

在对每个所述第一心内波形进行频域转换之后还包括：对每个所述第一单信号频谱图进行平滑处理。

4.根据权利要求3所述的心电信号分析方法，其特征在于，所述预处理包括：QRS波减除、带通或高通滤波处理、取模运算和低通滤波处理中的任一种或任意组合。

5.根据权利要求1所述的心电信号分析方法，其特征在于，所述对N个所述第一单信号频谱图进行累加处理，其中，所述累加处理包括曲线累加法或概率统计法。

6.根据权利要求5所述的心电信号分析方法，其特征在于，采用所述曲线累加法对N个所述第一单信号频谱图进行累加处理，生成所述第一比较信号频谱图的步骤进一步包括：

对每个所述第一单信号频谱图进行能量归一化处理；

对能量归一化处理后的N个所述第一单信号频谱图进行同频率能量累加，以获得第一比较信号频谱原图；

对所述第一比较信号频谱原图进行能量归一化处理，以获得所述第一比较信号频谱图。

7.根据权利要求5所述的心电信号分析方法，其特征在于，采用所述概率统计法对N个所述第一单信号频谱图进行累加处理，生成所述第一比较信号频谱图的步骤进一步包括：

获取每个所述第一单信号频谱图的频率，并获取每个所述频率在对应第一单信号频谱图中的权重；

计算N个所述第一单信号频谱图中相同频率下的权重平均值；

根据所述权重平均值，获得所述第一比较信号频谱图。

8.根据权利要求1-7中任一项所述的心电信号分析方法，其特征在于，所述对所述第一比较信号频谱图和所述第二比较信号频谱图进行比较分析的步骤进一步包括：

对所述第二比较信号频谱图和所述第一比较信号频谱图进行相似度分析，以获得相似度值；以及判断所述相似度值与预设相似度阈值的大小关系；

和/或，

计算所述第二比较信号频谱图中每个频率对应的能量与所述第一比较信号频谱图中相同频率对应的能量之间的能量差值，以获得多个能量差值；以及判断所述多个能量差值与预设能量阈值的大小关系。

9.根据权利要求1所述的心电信号分析方法，其特征在于，根据比较分析的结果，输出关于所述感兴趣点是否为病灶点的提示信息的步骤进一步包括：以文字、图像或语音的方式输出关于所述感兴趣点是否为所述病灶点的所述提示信息。

10.根据权利要求1所述的心电信号分析方法，其特征在于，其中所述感兴趣点为多个，从而生成多个第二比较信号频谱图，对所述第一比较信号频谱图和所述第二比较信号频谱图进行比较分析并产生比较分析结果的步骤还包括将所述第一比较信号频谱图和多个所述第二比较信号频谱图分别进行比较分析，从而产生所述比较分析结果。

11.一种心电信号分析装置，其特征在于，所述的心电信号分析装置包括：

第一获取模块，用于获取感兴趣区域或参考点/区域内的N个第一心内电信号，从而获得N个第一心内波形，其中N为大于等于1的整数；

第一生成模块，用于根据N个所述第一心内波形生成第一比较信号频谱图；

第二获取模块，用于获取感兴趣点的M个第二心内电信号，从而获得M个第二心内波形，其中M为大于等于1的整数；

第二生成模块，用于根据M个所述第二心内波形生成第二比较信号频谱图；

第一分析模块，用于对所述第一比较信号频谱图和所述第二比较信号频谱图进行比较分析并产生比较分析结果；以及

输出模块，用于根据比较分析的结果，输出关于所述感兴趣点是否为病灶点的提示信息；

所述第一生成模块包括：

第一处理单元，用于生成每个所述第一心内波形的第一单信号频谱图，以获得N个所述第一单信号频谱图；

第二处理单元，用于对N个所述第一单信号频谱图进行累加处理或叠加处理，生成所述第一比较信号频谱图，或者，自N个所述第一单信号频谱图中选择一个所述第一单信号频谱图，作为所述第一比较信号频谱图。

12.根据权利要求11所述的心电信号分析装置，其特征在于，所述第一处理单元进一步包括：

转换单元，用于对每个所述第一心内波形进行频域转换，以获得N个所述第一单信号频谱图。

13.根据权利要求12所述的心电信号分析装置，其特征在于，所述第一处理单元还进一步包括：

预处理单元，用于在所述转换单元对每个所述第一心内波形进行频域转换之前，对每个所述第一心内波形进行预处理，以使得每个所述第一心内波形在预设频率范围内；和/或，

平滑处理单元，用于在所述转换单元对每个所述第一心内波形进行频域转换之后，对每个所述第一单信号频谱图进行平滑处理。

14.根据权利要求11所述的心电信号分析装置，其特征在于，所述对N个所述第一单信号频谱图进行累加处理，其中，所述累加处理包括曲线累加法或概率统计法。

15.根据权利要求11-14中任一项所述的心电信号分析装置，其特征在于，所述第一分析模块进一步包括：

相似度分析单元，用于对所述第二比较信号频谱图和所述第一比较信号频谱图进行相似度分析，以获得相似度值；以及第一判断单元，用于判断所述相似度值与预设相似度阈值的大小关系；

和/或，

计算单元，用于计算所述第二比较信号频谱图中每个频率对应的能量与所述第一比较信号频谱图中相同频率对应的能量之间的能量差值，以获得多个能量差值；以及第二判断单元，用于判断所述多个能量差值与预设能量阈值的大小关系。

16.根据权利要求11所述的心电信号分析装置，其特征在于，所述输出模块以文字、图像或语音的方式输出关于所述感兴趣点是否为所述病灶点的所述提示信息。

17.根据权利要求11所述的心电信号分析装置，其特征在于，所述的心电信号分析装置还包括：

第三获取模块，用于获取所述病灶点被治疗后所述感兴趣区域或所述参考点/区域内的K个第三心内电信号，从而获得K个第三心内波形，其中K为大于等于1的整数；

第三生成模块，用于根据K个所述第三心内波形生成第三比较信号频谱图；

第二分析模块，用于将所述第一比较信号频谱图和所述第二比较信号频谱图中的至少一个与所述第三比较信号频谱图进行比较分析，以提供所述病灶点的治疗效果的信息。

18.根据权利要求17所述的心电信号分析装置，其特征在于，所述第二生成模块和/或所述第三生成模块与所述第一生成模块为同一功能模块；或者，所述第二获取模块和/或所述第三获取模块与所述第一获取模块为同一功能模块；或者，所述第一分析模块与所述第二分析模块为同一功能模块。

19.一种电生理信号记录仪，其特征在于，包括如权利要求11-18中任一项所述的心电信号分析装置。

20.一种三维标测系统，其特征在于，包括如权利要求11-18中任一项所述的心电信号分析装置，以及显示模块，所述显示模块用于显示一心脏三维模型；其中所述输出模块将所述提示信息输出于所述显示模块，并显示于所述心脏三维模型上。

21.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1-10中任一项所述的心电信号分析方法。

Images