具体实施方式

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。可以理解的是，此处所描述的具体实施例仅用于解释本申请，而非对本申请的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本申请相关的部分而非全部结构。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

本申请中的术语“第一”、“第二”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。

在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。

参阅图1，图1是本申请提供的心电图P波检测方法第一实施例的流程示意图，该方法包括：

步骤11：从心电图中确定满足预设要求的多个心拍。

参阅图2，对心电图片段进行介绍。如图2所示，心电图的横坐标表示时间，心电图的纵坐标表示检测信号值（如电压值）（图2中未示出横纵坐标）。其中，心电图包括多个心电图片段，图2示出了一个心电图片段。具体地，一个心电图片段主要包括PR间期和QT间期，QT间期又包括QRS间期和JT间期。

P波是心房除极波，代表左右二心房的激动。由于窦房结位于右心房内膜下，所以激动首先传到右心房，较晚传到左心房。右心房的除极作用因此也比左心房略早完毕。临床上为了实用起见，P波的前部代表右心房的激动，后部代表左心房的激动。分析P波对心律失常的诊断与鉴别诊断具有重要意义。

其中，P波参数包括：P波时限、P波形态、P波幅值、P波发生时间、是否有对应的QRS波中的至少一个。

P波时限是指P波的起始时间和终止时间之间的时间段，P波幅值包括P波端的最大信号值、最小信号值或最大信号值与最小信号值的差值，P波发生时间表示P波起始时间，是否有对应的QRS波是指P波后面是否有紧跟的QRS波。

进一步，根据上述P波参数还可以确定P波形态，P波形态可以根据不同的分类方式进行分类，例如可以将P波分为直立P波、倒置P波、正负双向P波和负正双向P波。

在一些实施例中，可以按照上述分类方式对心拍进行分类。

在一些实施例中，使用移动心电设备采集心电信号，将采集的心电信号进行滤波处理，对滤波处理后的心电信号执行步骤11。

在一些实施例中，可以从心电图中确定同一类型且RR间期长度满足预设间期阈值的多个心拍。

具体地，在获取到心电信号后，对心电信号进行心拍划分，以得到多个心拍，并为多个心拍确定心拍类型。在确定每个心拍类型时，还为每个心拍获取对应的RR间期长度。在一应用场景中，可以将每一心拍对应的R波与上一心拍对应的R波之间的RR间期作为此心拍的RR间期。在另一应用场景中，可以将每一心拍对应的R波与下一心拍对应的R波之间的RR间期作为此心拍的RR间期。则根据每一心拍对应类型和RR间期长度来确定对应每一类型的多个心拍。

步骤12：确定多个心拍的平均心拍模板。

参阅图3，步骤12可以是如下步骤：

步骤121：将每一心拍对应的第一R波之前的第一数量个采样点，与下一心拍对应的第二R波之后的第二数量个采样点之间采样点区域，作为一个心电片段，以得到多个心电片段。

在确定满足预设要求的多个心拍后，对多个心拍中的每一心拍再次按照步骤121进行处理，得到对应的心电片段。可以理解，此时的心电片段的长度大于其对应的心拍长度。通过这种方式可以获取每一心拍与下一心拍之间的所有波形。

步骤122：对多个心电片段进行平均处理，得到平均心拍模板。

在一些实施例中，因心电片段实际是由多个心电信号的采样点组成，每一采样点对应心电片段中的一个坐标点。可以将每个心电片段中对应的采样点进行平均处理，则可以得到多个平均采样点，将平均采样点形成的心电片段作为平均心拍模板。

步骤13：根据平均心拍模板确定参考P波。

在一些实施例中，参阅图4，步骤13可以具体是如下步骤：

步骤131：在平均心拍模板中的第一R波之后确定T波终点，以及在平均心拍模板中的第二R波之前确定QRS波起点。

在一些实施中，在平均心拍模板中，第一R波之后出现的最大的峰值点则确定为T波峰值点，在T波峰值点之后的第一次出现的最低点为T波终点。如，在第一R波之后，连续出现了峰值点A、峰值点B、峰值点C、峰值点D。其中，峰值点A小于峰值点B，峰值点B小于峰值点C，峰值点D小于峰值点C，则将峰值点C确定为T波峰值点，峰值点C与峰值点D之间的最低点为T波终点。

在第二R波之前确定QRS波起点，Q点为QRS波中的一波谷，则确定在Q点之前最近的峰值点为QRS波起点。如，在Q点之前，连续出现了点E、点F、点G、点H。其中，点E小于点F，点F小于点G，点H小于点G，则点H为离Q点最近的峰值点，将点C确定为QRS波起点。

步骤132：在T波终点和QRS波起点之间，确定参考P波。

在一些实施例中，在T波终点和QRS波起点之间确定所有采样点中最大的值为参考P波。

在一些实施例中，因T波终点和QRS波起点之间建立的坐标系存在负值的情况，则对每一采样点求绝对值，并将绝对值中的最大的值作为参考P波的峰值。

步骤14：根据参考P波确定多个心拍中每一心拍的P波。

在一些实施例中，参阅图5，步骤14可以具体是如下步骤：

步骤141：根据参考P波在平均心拍模板中的位置，在每一心拍中确定候选P波。

在本实施例中可将上述步骤121中的心电片段作为此时用于确定候选P波的心拍。

因平均心拍模板基于心电片段得到，则参考P波在平均心拍模板的位置，在每一心拍上也有对应的位置，将该位置对应的采样点作为该心拍的候选P波的峰值。

步骤142：在候选P波对应的时限和幅值满足预设条件时，确定候选P波为对应心拍的P波。

在一些实施例中，根据P波的峰值求出P波对应的幅值、时限，在幅值和时限满足预设条件时，确定候选P波为对应的心拍的P波。

区别于现有技术，本实施例提供的心电图P波检测方法包括：从心电图中确定满足预设要求的多个心拍；确定多个心拍的平均心拍模板；根据平均心拍模板确定参考P波；根据参考P波确定多个心拍中每一心拍的P波。通过上述方式，一方面利用求平均心拍模板的方式，能够均衡每一心拍中的干扰，使对应大干扰的心拍的干扰减少，在确定每一心拍的P波时，能够降低单一心拍中的干扰对于P波检测的影响，提高对P波的检测准确性，另一方面能够避免对每一心拍中的所有采样点进行检测，提高对每一心拍中P波的检测效率。

参阅图6，图6是本申请提供的心电图分析方法第二实施例的流程示意图，该方法包括：

步骤601：对心电图中的所有心拍进行分类。

在一些实施例中，在步骤601之前，对心电图中的所有心拍进行滤波处理。如，在检测P波过程中，为了抑制基线干扰、QRS波、T波等其他信号成分，预处理通过对心电信号进行高通、平滑滤波，减少其他信号成分对于P波的干扰。在其他实施例中，还可使用其他高通、低通、带通滤波器对心电信号进行滤波处理。

在一应用场景中，采集心电信号，对心电信号进行0.5Hz高通滤波处理，抑制基线等低频干扰。高通滤波器采用II阶Butterworth（巴特沃斯）滤波器，其实现方式为：

。

其中，y表示滤波输出；x表示输入的心电信号；x[n-k]表示前k个时刻的输入；y[n-k]表示前k个时刻的输出；b_k和a_k分别为Butterworth滤波器对应的滤波系数。

然后对心电信号进行M点滑动平均滤波，抑制心电信号中的毛刺干扰使波形平滑。其实现方式为：

y[n]=(x[n-M/2]+…+x[n-1]+x[n]+x[n+1]+…x[n+M/2])/M 。

其中，y[n]表示滤波输出；x[n-M/2]表示前M/2个时刻的输入；x[n-1]表示前一个时刻的输入；x[n]表示当前时刻的输入；x[n+1]表示后一个时刻的输入；x[n+M/2]表示后M/2个时刻的输入；M表示滑动窗口对应的采样点数量。

基于滤波后的心电信号生成心电图，对心电图中的心电信号按照心拍进行划分，得到多个心拍，对所有心拍进行分类。具体地，可根据每一心拍的参数进行分类。

步骤602：判断每一类心拍中候选RR间期是否满足：

RR-T≤RR _i ≤RR+T。

其中，RR为目标RR间期时间长度，T为设定时间长度，RR _i 为候选RR间期时间长度。若满足，执行步骤603，若不满足，则舍弃该候选RR间期对应的心拍。

可以理解，根据分类情况，则对每一类中的多个心拍进行判断。

步骤603：确定候选RR间期对应的心拍为满足预设间期阈值的多个心拍中的一个。

当候选RR间期满足时，则确定该候选RR间期与目标RR间期大小相近。

步骤604：将每一心拍对应的第一R波之前的第一数量个采样点，与下一心拍对应的第二R波之后的第二数量个采样点之间采样点区域，作为一个心电片段，以得到多个心电片段。

步骤605：对多个心电片段中，多个相应采样点的值进行求中值运算，以得到平均心拍模板中对应采样点的值。

在一些实施例中，参阅图7，对步骤604-605进行说明：

图7中包括A、B、C、D和E五个心电片段，其中，A、B、C和D为按照上述步骤得到的同一类型且满足预设间期阈值的心电片段。其中，A={a1,a2,a3……,an}，B={b1,b2,b3……,bn}，C={c1,c2,c3……,cn}，D={d1,d2,d3……,dn}， E为A，B，C，D对应生成的平均心拍模板，E={e1,e2,e3……,en}，其中e_i的计算方式如下：

e_i=median(a_i，b_i，c_i，d_i)。

其中，median为取中值的计算方法，即对a_i，b_i，c_i，d_i 进行排序，取中间排序的数赋值给e_i；i的范围为1~n。如，a_i=9，b_i=10，c_i=8，d_i=9，通过排序后则为c_i，a_i， d_i，b_i，因此处心电片段为偶数个，则求a_i和 d_i 的平均值，将平均值赋值给e_i。在心电片段为奇数个时，直接取中间的值赋值给e_i。通过上述方式，将具有相似特性的片段聚合在一起生成平均心拍模板，能够有效降低单一心拍中存在的干扰对于P波检测的影响。由于检测只关注平均心拍模板，相比逐一的心拍检测，计算量减少，能够进一步提升算法的运行效率。在生成平均心拍模板时，采用中值的方法替代平均值，当单一片段中存在高尖干扰或显著基线干扰时，不会对平均心拍模板造成影响，有效降低P波误检的概率。

步骤606：在平均心拍模板中的第一R波之后确定T波终点，以及在平均心拍模板中的第二R波之前确定QRS波起点。

参阅图8，对步骤606进行说明：

如图8所示，在平均心拍模板中，第一R波之后出现的最大的峰值点则确定为T波峰值点，如图8所示的“T波”，在T波峰值点之后的第一次出现的最低点为T波终点。在第二R波之前确定QRS波起点，Q点为QRS波中的一波谷，则确定在Q点之前最近的峰值点为QRS波起点，如图8所示的“QRS波起点”。

步骤607：确定T波终点和QRS波起点之间所有采样点对应的绝对值的平均值。

步骤608：获取T波终点和QRS波起点之间，大于平均值的最值，作为参考P波的峰值。

参阅图9，将图8所示的T波终点至QRS波起点之间的采样点作为寻找参考P波的采样点范围。计算该范围内所有采样点绝对值的平均值作为寻找P波峰值点的阈值。遍历图9中的所有采样点，遇到绝对值大于阈值的最值点时则作为参考P波保存（如图9所示，出现两个参考P波）。基于绝对值计算阈值，不仅考虑了正向的P波，同时也可对负向的P波进行检测。相比传统方法中基于最大值计算阈值的方法，平均值能够有效规避高尖干扰对阈值计算的影响，从而降低阈值过高导致P波漏检的概率。通过先确定T波的终点和QRS波的起点，进一步收缩检测范围减少计算量，同时，降低T波和QRS波本身对于P波检测的干扰。由于是对整个搜索范围内所有的最值点进行检测，实现了对该范围内所有P波的检测。

步骤609：根据参考P波在平均心拍模板中的位置，在每一心拍中确定候选P波。

步骤610：在候选P波对应的时限和幅值满足预设条件时，确定候选P波为对应心拍的P波。

参阅图10，对步骤609-610进行说明：

因平均心拍模板基于心电片段得到，则参考P波在平均心拍模板的位置，在每一心拍上也有对应的位置，将该位置对应的采样点作为该心拍的候选P波的峰值。根据P波的峰值求出P波对应的幅值、时限，在幅值和时限满足预设条件时，确定候选P波为对应的心拍的P波。

可根据图10中的平均心拍模板中的参考P波来确定对应的心拍1、心拍2和心拍n中的候选P波。

在一些实施例中，参考P波对应在心拍上的采样点并不是P波的峰值点，则可根据心拍上的采样点周围的采样点来确定峰值点。如，参考P波在一心拍上对应的第一采样点的前一个采样点的绝对值大于第一采样点的绝对值，后一采样点的绝对值小于第一采样点的绝对值，则说明第一采样点的右侧存在一峰值，则在第一采样点的右侧确定该峰值点，将该峰值点确定为候选P波。

又如，参考P波在一心拍上对应的第二采样点的前一个采样点的绝对值小于第二采样点的绝对值，后一采样点的绝对值大于第二采样点的绝对值，则说明第二采样点的左侧存在一峰值，则在第二采样点的左侧确定该峰值点，将该峰值点确定为候选P波。

又如，参考P波在一心拍上对应的第三采样点的前一个采样点的绝对值小于第三采样点的绝对值，后一采样点的绝对值小于第三采样点的绝对值，则说明第三采样点为峰值，将该峰值点确定为候选P波。

在一应用场景中，在确定候选P波峰值点后，需要对该候选P波进行幅值和时限的计算，并将计算出来的幅值和时限与阈值进行比较，满足阈值时，确定该候选P波为对应心拍的P波。

具体地，根据候选P波波峰位置，寻找对应的左侧和右侧的两个转折点，利用两个转折点和峰值点进行时限和幅值的计算。

计算方法如下：

时限= P _end - P _start 。

幅值=(|ECG[P _peak ]- ECG[P _start ]|+ |ECG[P _peak ]- ECG[P _end ]|)/2。

其中，P _end 表示候选P波峰值右侧转折点，P _start 表示候选P波峰值左侧转折点，ECG [P _peak ]表示在ECG信号中候选P波峰值处对应的值；ECG[P _start ]表示在ECG信号中候选P波峰值左侧转折点处对应的值；ECG[P _end ]表示在ECG信号中候选P波峰值右侧转折点处对应的值。

在计算出时限后，与参考P波对应的时限和幅值进行比较，当候选P波的时限大于参考P波的时限，且候选P波的幅值大于参考P波的幅值时，确定该候选P波为对应心拍的P波。

通过上述方式，能够避免人工设定固定阈值不当导致P波漏检和误检。通过计算每个心拍中所有候选P波对应的幅值和时限，与参考P波的幅值和时限进行比较，保留满足条件的P波作为最终的P波检测结果。通过对每个心拍中的候选P波的幅值、时限进行判定，避免平均心拍模板中存在P波而部分心拍无对应P波而误检的情况。

在本实施例中，在平均心拍模板中进行检测时，检测范围不再局限于R波之前的一段范围，而是对T波终点到R波起点内的所有采样点进行判断，从而实现所有P波的同步检测。能够方便的应用于信号干扰较大的长时间心电信号P波定位，对心电信号的完整分析起到良好的改善效果。

参阅图11，图11为利用本申请实施例的方法对一个心电信号的P波检测结果，在图中可以看出，在RR间期内存在一个P波和多个P波的情况，本申请实施例的方法能够针对连续出现的单个和多个P波，均能有效准确的定位到P波的峰值点。

参阅图12，图12是本申请提供的心电图分析装置一实施例的结构示意图，该心电图分析装置200包括处理器201和存储器202，存储器202用于存储程序数据，处理器201用于执行程序数据，以实现如下的方法：

从心电图中确定满足预设要求的多个心拍；确定多个心拍的平均心拍模板；根据平均心拍模板确定参考P波；根据参考P波确定多个心拍中每一心拍的P波。

可以理解地，本实施例中的处理器201还可以实现上述实施例中任一方法步骤，这里不再赘述。

另外，可以理解地，该心电图分析装置200还包括采集器接口（图未示），采集器接口连接外部心电采集器，心电采集器进一步通过多个导联连接电极，以在将电极贴敷于人体时，进行心电数据采集。

另外，该心电图分析装置200还包括显示屏（图未示），显示屏用于显示心电图以及对应的分析数据。进一步，该显示屏为触摸显示屏，用于获取用户的触摸指令对心电图及其分析过程进行操作。

在一实施例中，心电图分析装置200为移动心电图机，移动心电图机用于动态心电图检查，作为一种在患者日常生活中连续记录24小时或更长时间的心电活动的检查，并辅助计算机进行分析处理的新技术，越来越受到欢迎。较之普通心电图，动态心电图在24小时内可连续记录多达10万次左右的心动周期，这样，可以发现常规体表心电检查不易发现的心律失常和心肌缺血等，能够提高对非持续性心律失常，尤其是对一过性心律失常及短暂的心肌缺血发作的检出率，很好地扩大了心电图临床运用的范围。

参阅图13，图13是本申请提供的计算机可读存储介质一实施例的结构示意图，该计算机可读存储介质300中存储有程序数据301，该程序数据301在被处理器执行时，用以实现如下的方法：

从心电图中确定满足预设要求的多个心拍；确定多个心拍的平均心拍模板；根据平均心拍模板确定参考P波；根据参考P波确定多个心拍中每一心拍的P波。

可以理解地，本实施例中的计算机可读存储介质300还可以实现上述实施例的任一方法，这里不再赘述。

在本申请所提供的几个实施方式中，应该理解到，所揭露的方法以及设备，可以通过其它的方式实现。例如，以上所描述的设备实施方式仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施方式方案的目的。

另外，在本申请各个实施方式中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

以上所述仅为本申请的实施方式，并非因此限制本申请的专利范围，凡是根据本申请说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本申请的专利保护范围内。