一种ECG信号处理方法及装置(700)，用以解决单导联心电采集的单臂测量过程中引入较大的干扰，严重影响了心率的精确计算和心律的精准分析的问题。该方法包括：测量设备采集ECG信号(100)；提取ECG信号的第k个有效QRS波群(110)，计算第k个时间差和第(k+1)个时间差(120)，判断第k个时间差与第(k+1)个时间差是否为目标时间差(130a、130b)。因此，提高了心电特征提取的准确率，极大地消除了单臂测量过程中引入的干扰，能够有效保证心率计算精确度和心律分析的精准度。
1.一种心电图ECG信号处理方法，其特征在于，包括：
提取所述ECG信号的第k个有效QRS波群，其中QRS波群包括第i个R波的峰值点R前的第一极值点Q、第i个R波的峰值点R、和第i个R波的峰值点R后的第二极值点S，其中，i≥2，k≥
所述测量设备计算第k个时间差和第(k+1)个时间差，其中，所述第k个时间差为所述第k个有效QRS波群中R波的峰值点R与第(k+1)个有效QRS波群中R波的峰值点R的时间差，所述第(k+1)个时间差为所述第(k+1)个有效QRS波群中的R波的峰值点R与第(k+2)个有效QRS波群中的R波的峰值点R的时间差；
如果所述第k个时间差与所述第(k+1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值小于预设阈值时，所述测量设备确定所述第k个时间差和所述第(k+1)个时间差为目标时间差；或者，如果所述第k个时间差与所述第(k+1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值大于等于预设阈值时，所述测量设备确定所述第k个时间差和所述第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差，所述时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值；
其中，所述目标时间差用于计算所述ECG信号对应的心率值。
2.如权利要求1所述的方法，其特征在于，所述提取所述ECG信号的第k个有效QRS波群，包括：所述测量设备确定所述第一极值点Q和所述第二极值点S；
所述测量设备计算所述第一极值点Q与所述第i个R波的峰值点R之间的时间差为第一时间差，所述第i个R波的峰值点R与所述第二极值点S之间的时间差为第二时间差，所述第一极值点Q与所述第二极值点S之间的时间差为第三时间差；
如果所述第一时间差、第二时间差、第三时间差分别小于对应的第一阈值、第二阈值、第三阈值时，则所述第一极值点Q对应Q波，所述第二极值点S对应S波，所述测量设备确定所述第i个R波与对应的Q波和S波构成第k个有效QRS波群。
3.如权利要求1或2所述的方法，其特征在于，还包括：
所述测量设备根据至少一个所述信号参数计算信号评估参数；
所述测量设备根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级。
4.如权利要求3所述的方法，其特征在于，所述信号参数包括有效信号功率，基线漂移，和带内噪声；
所述信号评估参数包括信号伪迹比、和信号带内噪声比，其中，所述信号伪迹比为所述有效信号功率和所述基线漂移的函数，所述信号带内噪声比为所述有效信号功率和所述带内噪声的函数；
所述测量设备根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级，包括：所述测量设备确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第一等级；
所述测量设备确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第二等级；
所述测量设备确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第三等级；
所述测量设备确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第四等级；
其中，所述第一等级优于所述第二等级，所述第二等级优于所述第三等级，所述第三等级优于所述第四等级。
5.如权利要求1、2、4任一项所述的方法，其特征在于，在所述提取所述ECG信号的第k个有效QRS波群之前，还包括：所述测量设备对所述ECG信号进行滤波处理，并采用三轴加速计拟合用户的运动轨迹；
所述测量设备将滤波后的ECG信号与所述运动轨迹进行比对，将所述运动轨迹中运动幅度值大于预设幅度阈值的时长对应的所述滤波后的ECG信号中的ECG波形删除。
6.一种可穿戴设备，其特征在于，包括：处理器、存储器、心率采集器，电源；
所述处理器，用于调用所述存储器中的指令，执行上述权利要求 1-5任一项所述的方法。
7.如权利要求6所述的设备，其特征在于，还包括：加速传感器；
所述加速传感器，用于检测各个方向上的加速度的大小，拟合运动轨迹。
8.一种ECG信号处理装置，其特征在于，包括：采集单元，提取单元，计算单元，处理单元，其中，所述采集单元，用于采集ECG信号；
所述提取单元，用于提取所述ECG信号的第k个有效QRS波群，其中QRS波群包括第i个R波的峰值点R前的第一极值点Q、第i个R波的峰值点R、和第i个R波的峰值点R后的第二极值点S，其中，i≥2，k≥2，k≤i；
所述计算单元，用于计算第k个时间差和第(k+1)个时间差，其中，所述第k个时间差为所述第k个有效QRS波群中R波的峰值点R与第(k+1)个有效QRS波群中R波的峰值点R的时间差，所述第(k+1)个时间差为所述第(k+1)个有效QRS波群中的R波的峰值点R与第(k+2)个有效QRS波群中的R波的峰值点R的时间差；
所述处理单元，用于如果所述第k个时间差与所述第(k+1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值小于预设阈值时，确定所述第k个时间差和所述第(k+1)个时间差为目标时间差；或者，如果所述第k个时间差与所述第(k+
1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值大于等于预设阈值时，确定所述第k个时间差和所述第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差，所述时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值；
其中，所述目标时间差用于计算所述ECG信号对应的心率值。
9.如权利要求8所述的装置，其特征在于，所述提取单元，用于：
所述提取所述ECG信号的第k个有效QRS波群时，确定所述第一极值点Q和第二极值点S；
计算所述第一极值点Q与所述第i个R波的峰值点R之间的时间差为第一时间差，所述第i个R波的峰值点R与所述第二极值点S之间的时间差为第二时间差，所述第一极值点Q与所述第二极值点S之间的时间差为第三时间差；
如果所述第一时间差、第二时间差、第三时间差分别小于对应的第一阈值、第二阈值、第三阈值时，则所述第一极值点Q对应Q波，所述第二极值点S对应S波，所述第i个R波与对应的Q波和S波构成第k个有效QRS波群。
10.如权利要求8或9所述的装置，其特征在于，所述装置，还包括：信号质量评估单元，用于计算所述ECG信号对应的信号参数；
根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级。
11.如权利要求10所述的装置，其特征在于，所述信号参数包括有效信号功率，基线漂移，和带内噪声；
所述信号评估参数包括信号伪迹比、和信号带内噪声比，其中，所述信号伪迹比为所述有效信号功率和所述基线漂移的函数，所述信号带内噪声比为所述有效信号功率和所述带内噪声的函数；
所述信号质量评估单元，还用于：根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级时，确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第一等级；
确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第二等级；
确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第三等级；
确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第四等级；
其中，所述第一等级优于所述第二等级，所述第二等级优于所述第三等级，所述第三等级优于所述第四等级。
12.如权利要求8、9、11任一项所述的装置，其特征在于，所述装置还包括：运动轨迹分析单元，用于在提取所述ECG信号的第k个有效QRS波群之前，对所述ECG信号进行滤波处理，并采用三轴加速计拟合用户的运动轨迹；
将滤波后的ECG信号与所述运动轨迹进行比对，将所述运动轨迹中运动幅度值大于预设幅度阈值的时长对应的所述滤波后的ECG信号中的ECG波形删除。
一种ECG信号处理方法及装置
[0001] 本申请要求于2016年11月3日提交中国专利局、申请号为201610974448.2、发明名称为“一种单臂ECG的测量方法和终端”的CN专利申请的优先权，其全部内容通过引用结合在本申请中。
技术领域
[0002] 本发明涉及终端领域，尤其涉及一种ECG信号处理方法及装置。
背景技术
[0003] 针对传统的心电图(Electro Cardio Graph，ECG)测量需要跨心区测量，应用场景受限、造成用户体验较差。
[0004] 现有技术中，提出多种便携式单导联心电采集测量方案。具体的，单导联心电采集设备多采用模拟肢体I导联(双手)，或者模拟胸导联(胸区，跨心脏)的测量方式，这些测量方式接近于标准导联体系获取的心电信号，目前在便携式心电分析服务类的产品中得到广泛应用。但是，在单导联心电采集的单臂测量过程中，可能会引入较大的干扰，其中肌电干扰、运动伪迹等较为显著，严重影响了心率的精确计算和心律的精准分析。
发明内容
[0005] 本发明实施例提供一种ECG信号处理方法及装置，用以解决单导联心电采集的单臂测量过程中，引入较大的干扰，严重影响了心率的精确计算和心律的精准分析的问题。
[0006] 第一方面，一种心电图ECG信号处理方法，包括：
[0007] 测量设备采集ECG信号，提取ECG信号的第k个有效QRS波群，其中QRS波群包括第i个R波的峰值点R前的第一极值点Q、第i个R波的峰值点R、和第i个R波的峰值点R后的第二极值点S，其中，i≥2，k≥2，k≤i，进一步计算第k个时间差和第(k+1)个时间差，其中，第k个时间差为第k个有效QRS波群中R波的峰值点R与第(k+1)个有效QRS波群中R波的峰值点R的时间差，第(k+1)个时间差为第(k+1)个有效QRS波群中的R波的峰值点R与第(k+2)个有效QRS波群中的R波的峰值点R的时间差。如果第k个时间差与第(k+1)个时间差处于预设范围，且第(k+1)个时间差与第k个时间差的差值的绝对值小于预设阈值时，测量设备确定第k个时间差和第(k+1)个时间差为目标时间差；或者，如果第k个时间差与第(k+1)个时间差处于预设范围，且第(k+1)个时间差与第k个时间差的差值的绝对值大于等于预设阈值时，测量设备确定第k个时间差和第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差，时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值。其中，目标时间差用于计算ECG信号对应的心率值。
[0008] 因此，本发明实施例中提出一种ECG信号处理方法，首先从采集的ECG信号中提取有效QRS波群，提高了心电特征提取的准确率，然后，计算相邻有效QRS波群中的R波的峰值点之间的时间差，进一步根据获得的时间差筛选出满足条件的目标时间差，从而确定心率值，极大地消除了单臂测量过程中引入的干扰，能够有效保证心率计算精确度和心律分析的精准度。
[0009] 在一种可能的实现方式中，所述提取所述ECG信号的第k个有效QRS波群，包括：测量设备确定第i个R波对应的第一极值点Q和第二极值点S，计算第一极值点Q与第i个R波的峰值点R之间的时间差为第一时间差，第i个R波的峰值点R与第二极值点S之间的时间差为第二时间差，第一极值点Q与第二极值点S之间的时间差为第三时间差。如果第一时间差、第二时间差、第三时间差分别小于对应的第一阈值、第二阈值、第三阈值时，则第一极值点Q对应Q波，第二极值点S对应S波，测量设备确定第i个R波与对应的Q波和S波构成第k个有效QRS波群。
[0010] 因此，采用本发明实施例提供的方法从采集的ECG信号中提取有效QRS波群，提高了心电特征提取的准确率，进而有效保证心率计算精确度和心律分析的精准度。
[0011] 在一种可能的实现方式中，还包括：
[0012] 测量设备计算ECG信号对应的信号参数，根据至少一个信号参数计算信号评估参数，并根据信号评估参数与信号评估参数对应的预设评估阈值，确定ECG信号的信号质量等级。
[0013] 因此，采用本发明实施例提供的方法，测量设备能够评估不同的信号质量等级，通过终端对用户进行提示，当信号质量较差时，提示用户查找原因，例如佩戴位置错误、接触不良等。
[0014] 在一种可能的实现方式中，信号参数包括有效信号功率，基线漂移，和带内噪声；信号评估参数包括信号伪迹比、和信号带内噪声比，其中，信号伪迹比为有效信号功率和基线漂移的函数，信号带内噪声比为有效信号功率和带内噪声的函数；
[0015] 测量设备根据信号评估参数与信号评估参数对应的预设评估阈值，确定ECG信号的信号质量等级，包括：
[0016] 测量设备确定信号伪迹比大于对应的预设评估阈值，且信号带内噪声比大于对应的预设评估阈值时，确定ECG信号的信号质量等级为第一等级；
[0017] 测量设备确定信号伪迹比小于等于对应的预设评估阈值，且信号带内噪声比大于对应的预设评估阈值时，确定ECG信号的信号质量等级为第二等级；
[0018] 测量设备确定信号伪迹比大于对应的预设评估阈值，且信号带内噪声比小于等于对应的预设评估阈值时，确定ECG信号的信号质量等级为第三等级；
[0019] 测量设备确定信号伪迹比小于等于对应的预设评估阈值，且信号带内噪声比小于等于对应的预设评估阈值时，确定ECG信号的信号质量等级为第四等级；
[0020] 其中，第一等级优于第二等级，第二等级优于第三等级，第三等级优于第四等级。
[0021] 因此，采用本发明实施例提供的方法获得的信号质量等级可用于对采集信号的质量进行预判，通过将信号质量等级反馈给用户，使用户在确定当前信号质量等级较差时纠正佩戴位置，保证采集到的ECG信号的有效性，进而有效保障单臂心电测量设备的信号测量准确性。
[0022] 在一种可能的实现方式中，在提取ECG信号的第k个有效QRS波群之前，还包括：测量设备对ECG信号进行滤波处理，并采用三轴加速计拟合用户的运动轨迹，将滤波后的ECG信号与运动轨迹进行比对，将运动轨迹中运动幅度值大于预设幅度阈值的时长对应的滤波后的ECG信号中的ECG波形删除。
[0023] 因此，本发明实施例的测量设备通过将滤波后的ECG信号与运动轨迹进行比对，将运动轨迹中运动幅度值大于预设幅度阈值的时长内对应的滤波后的ECG信号中的ECG波形删除，从而消除运动伪迹对ECG信号采集的影响，提高采集到的ECG信号的有效性。
[0024] 第二方面，一种可穿戴设备，包括：处理器、存储器、心率采集器，电源；其中，心率采集器，用于采集ECG信号；存储器，用于存储指令；处理器，用于调用存储器中的指令，执行第一方面或第一方面中任意一种可能的实现方式。
[0025] 因此，本发明实施例中提出的可穿戴设备，首先从采集的ECG信号中提取有效QRS波群，提高了心电特征提取的准确率，然后，计算相邻有效QRS波群中的R波的峰值点之间的时间差，进一步根据获得的时间差筛选出满足条件的目标时间差，从而确定心率值，极大地消除了单臂测量过程中引入的干扰，能够有效保证心率计算精确度和心律分析的精准度。
[0026] 在一种可能的实现方式中，还包括：加速传感器；
[0027] 加速传感器，用于检测各个方向上的加速度的大小，拟合运动轨迹。
[0028] 因此，通过加速传感器拟合运动轨迹，将滤波后的ECG信号与运动轨迹进行比对，将运动轨迹中运动幅度值大于预设幅度阈值的时长内对应的滤波后的ECG信号中的ECG波形删除，从而消除运动伪迹对ECG信号采集的影响，提高采集到的ECG信号的有效性。
[0029] 第三方面，本申请提供了一种ECG信号处理装置，用于执行上述第一方面或第一方面的任意可能的实现方式中的方法。具体地，该装置包括用于执行上述第一方面或第一方面的任意可能的实现方式中的方法的单元。
附图说明
[0030] 图1为本发明实施例中的ECG信号处理方法的概述流程图；
[0031] 图2为本发明实施例中标准ECG信号的波形特征示意图；
[0032] 图3为本发明实施例中ECG信号滤波后的波形图；
[0033] 图4(a)为本发明实施例中用户佩戴设备有误时电极位置偏移的示意图；
[0034] 图4(b)为本发明实施例中用户佩戴设备无误时电极位置正常的示意图；
[0035] 图5为本发明实施例中无效QRS波群的示意图；
[0036] 图6为本发明实施例中可穿戴设备的结构示意图；
[0037] 图7为本发明实施例中ECG信号处理装置的结构示意图。
具体实施方式
[0038] 下面结合附图，对本发明的实施例进行描述。
[0039] 参阅图1所示，针对单导联心电采集的单臂测量过程中，引入较大的干扰，严重影响了心率的精确计算和心律的精准分析的问题，本发明实施例提供一种ECG信号处理方法，该方法包括：
[0040] 步骤100：测量设备采集ECG信号。
[0041] 本发明实施例中所指的测量设备为单导联心电采集设备。例如，可穿戴设备等。
[0042] 测量设备采集ECG信号，例如可以为：用户佩戴单导联心电采集设备，打开开关，设备进入自检阶段，设备自检通过表明用户佩戴无误，设备工作正常，否则提示用户纠正佩戴位置，在设备进入正常工作状态后，开始采集用户的ECG信号。
[0043] 步骤110：测量设备提取ECG信号的第k个有效QRS波群，其中QRS波群包括第i个R波峰值点R前的第一极值点Q、第i个R波的峰值点R、和第i个R波峰值点R后的第二极值点S，其中，i≥2，k≥2，k≤i。
[0044] 针对ECG信号中的每个R波，测量设备判断当前R波的峰值点R与左右两个极值点能否构成一个有效的QRS波群，具体的，可以采用但不限于以下方法确定有效的QRS波群：
[0045] 图2为标准ECG信号的波形特征示意图。如图2所示，测量设备可以根据标准ECG信号的波形特征确定第i个R波峰值点R对应的第一极值点Q和第二极值点S，其中，第一极值点Q为第i个R波峰值点R前的第一个极值点Q，第二极值点S为第i个R波峰值点R后的第一个极值点Q，即R波峰值点R左右两侧临近的两个极值点。
[0046] 进一步地，测量设备计算第一极值点Q与第i个R波的峰值点R之间的时间差，将该时间差记为第一时间差QS，计算第i个R波的峰值点R与第二极值点S之间的时间差，将该时间差记为第二时间差RS，计算第一极值点Q与第二极值点S之间的时间差，将该时间差记为第三时间差QS。
[0047] 测量设备确定第一时间差、第二时间差、第三时间差分别小于对应的第一阈值、第二阈值、第三阈值时，确定第一极值点Q对应Q波，第二极值点S对应S波，进而确定第i个R波与对应的Q波和S波构成第k个有效QRS波群。
[0048] 具体的，这里的第一阈值、第二阈值、第三阈值可根据经验值设定，只要符合人体生理特征即可，不作为对本发明实施例的限定。本发明实施例中，第一阈值的范围为90～110ms，例如第一阈值可以为100ms。第二阈值的范围为90～110ms，例如第二阈值可以为
100ms。第三阈值的范围为140～160ms，例如第三阈值可以为150ms，则有效QRS波群的判断条件可以为QR＜100ms，且RS＜100ms，且QS＜150ms。
[0049] 当三个时间差中任一时间差大于等于对应的阈值时，该QRS波群为无效的QRS波群，当三个时间差分别小于对应的阈值时，该QRS波群为有效QRS波群。
[0050] 应理解的是，由于可能存在无效的QRS波群，因此，k≤i。
[0051] 步骤120：测量设备计算第k个时间差和第(k+1)个时间差。
[0052] 其中，第k个时间差为第k个有效QRS波群中R波的峰值点R与第(k+1)个有效QRS波群中R波的峰值点R的时间差，第(k+1)个时间差为第(k+1)个有效QRS波群中的R波的峰值点R与第(k+2)个有效QRS波群中的R波的峰值点R的时间差。
[0053] 可选地，第k个有效QRS波群、第(k+1)个有效QRS波群、第(k+2)个有效QRS波群为三个连续的波群。
[0054] 在确定各个有效QRS波群之后，测量设备计算所有相邻有效QRS波群中两个R波的峰值点R的时间差。
[0055] 步骤130a：如果第k个时间差与第(k+1)个时间差处于预设范围，且第(k+1)个时间差与第k个时间差的差值的绝对值小于预设阈值时，测量设备确定将第k个时间差和第(k+1)个时间差为目标时间差。
[0056] 步骤130b：如果第k个时间差与第(k+1)个时间差处于预设范围，且第(k+1)个时间差与第k个时间差的差值的绝对值大于等于预设阈值时，测量设备确定第k个时间差和第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差，时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值。
[0057] 其中，目标时间差用于计算ECG信号对应的心率值。
[0058] 本发明实施例中在确定第k个时间差和第(k+1)个时间差是否为目标时间差时，测量设备首先判断第k个时间差和第(k+1)个时间差是否均处于预设范围，即测量设备判断在步骤120中计算出的哪些时间差处于预设范围。这里的预设范围为0.3s～1.5s，测量设备将不满足预设范围的第k个时间差或第(k+1)个时间差删除，如果第k个时间差和第(k+1)个时间差均不处于预设范围内，则将第k个时间差和第(k+1)个时间差均删除。应理解的是，删除的时间差不为目标时间差。对于上述预设范围的数值可以根据实际需要进行设定和调整，上述数值仅为举例说明，不作为对本发明实施例的限定。
[0059] 在确定第k个时间差与第(k+1)个时间差处于预设范围时，测量设备进一步判断第k个时间差与第(k+1)个时间差的差值是否小于预设阈值，此处的预设阈值的范围可以为0.1s～0.15s，例如，预设阈值可以为0.12s。对于上述预设阈值的数值可以根据实际需要进行设定和调整，上述数值仅为举例说明，不作为对本发明实施例的限定。
[0060] 如果第(k+1)个时间差与第k个时间差的差值小于预设阈值时，测量设备确定第k个时间差和第(k+1)个时间差为目标时间差。
[0061] 如果第(k+1)个时间差与第k个时间差的差值的绝对值大于等于预设阈值时，测量设备确定第k个时间差和第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差。这里的时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值，即之前未被删除的时间差的平均值。
[0062] 进一步地，测量设备根据获得的目标时间差计算心率。具体的，测量设备统计预设时长内的目标时间差，根据这些目标时间差计算目标时间差的平均值，并根据目标时间差的平均值计算心率，其中，心率值＝60/目标时间差的平均值。
[0063] 综上所述，本发明实施例中提出一种ECG信号处理方法，首先从采集的ECG信号中提取有效QRS波群，提高了心电特征提取的准确率，然后，计算相邻有效QRS波群中的R波的峰值点之间的时间差，进一步根据获得的时间差筛选出满足条件的目标时间差，从而确定心率值，极大地消除了单臂测量过程中引入的干扰，能够有效保证心率计算精确度和心律分析的精准度。
[0064] 进一步的，本发明实施例在执行步骤100之后，执行步骤110之前，测量设备还可以执行如下操作：
[0065] 测量设备对ECG信号进行滤波处理，并采用三轴加速记拟合用户的运动轨迹。
[0066] 例如，测量设备采用梳状滤波器对ECG信号进行滤波处理，具体的，阻带中心点可以为1Hz，50Hz，100Hz，150Hz等，滤波后的波形如图3所示。
[0067] 测量设备采用三轴加速度计拟合运动轨迹时，取x、y、z三轴拟合进行拟合获得运动轨迹对应的曲线。
[0068] 本发明实施例的测量设备通过将滤波后的ECG信号与运动轨迹进行比对，将运动轨迹中运动幅度值大于预设幅度阈值的时长内对应的滤波后的ECG信号中的ECG波形删除，从而消除运动伪迹对ECG信号采集的影响，提高采集到的ECG信号的有效性。
[0069] 进一步的，在步骤130a或步骤130b之后，测量设备还可继续执行如下操作，以判断ECG信号的信号质量。
[0070] 首先，测量设备计算ECG信号对应的信号参数。
[0071] 例如，这里的信号参数可以为以下参数但不限于以下参数，包括：有效信号功率(记为PW_ECG)，工频干扰(记为PW_50)，基线漂移(记为PW_1)，带内噪声(记为PW_NB)等。
[0072] 原始信号功率记为ECG_original，PW_50及PW_1可由原始信号经过频域变换获得，PW_ECG可以通过R波幅度值、QRS宽度、T波幅度值和T波宽度计算获得。
[0073] PW_NB＝ECG_original–PW_50–PW_1–PW_ECG。
[0074] 上述信号参数可通过现有技术中提供的方法计算获得。
[0075] 然后，测量设备根据至少一个信号参数计算信号评估参数。
[0076] 这里的信号评估参数可以为信号伪迹比(Signal Artifical Ratio，SAR)和信号带内噪声比(Signal Noise in Band Ratio，SNBR)。
[0077] 其中，SAR＝lg(PW_ECG/PW_1)，SNBR＝lg(PW_ECG/PW_NB)。
[0078] 最后，测量设备根据信号评估参数与信号评估参数对应的预设评估阈值，确定ECG信号的信号质量等级。
[0079] 进一步的，测量设备可以将ECG信号的信号质量等级通知给用户，例如语音提示，或弹出对话框提示等方式，在此不作限定。
[0080] 具体的，当信号评估参数为SAR和SNBR时，如果SAR大于对应的预设评估阈值(记为TH_SAR)，且SNBR大于对应的预设评估阈值(记为TH_SNBR)时，测量设备确定ECG信号的信号质量等级为第一等级。
[0081] 如果SAR小于等于TH_SAR，且SNBR大于TH_SNBR时，测量设备确定ECG信号的信号质量等级为第二等级；
[0082] 如果SAR大于TH_SAR，且SNBR小于等于TH_SNBR时，测量设备确定ECG信号的信号质量等级为第三等级；
[0083] 如果SAR小于等于TH_SAR，且SNBR小于等于TH_SNBR时，测量设备确定ECG信号的信号质量等级为第四等级；
[0084] 其中，第一等级优于第二等级，第二等级优于第三等级，第三等级优于第四等级。
[0085] 表1为信号质量等级对应表，具体如表1所示，分别对SAR和SNBR按照阈值组合对信号质量进行分级。信号质量等级(SQL)分级为1～4，其中，1为优，依次递减，4为噪声大特征无法提取，对于不同的信号质量等级，可以通过终端对用户进行提示，当信号质量较差时，提示用户查找原因，例如佩戴位置错误、接触不良等。
[0086] 表1
[0087]信号等级 信号质量评估
1 SAR>TH_SAR&&SNBR>TH_SNBR
2 SARTH_SNBR
3 SAR>TH_SAR&&SNBR4 SAR
[0088] 下面以一个具体的心电测量流程为例说明本发明实施例提供的ECG信号处理方法。
[0089] 步骤(Step)1，用户佩戴单导联心电采集设备，打开开关，设备进入自检阶段，进行采集信号质量预判。
[0090] 具体的，设备执行一遍完整的step1～step4后，判断1分钟内连续测量的目标时间差数量是否满足预设条件。若1分钟内目标时间差数量满足预设条件，则表明用户佩戴无误，设备正常工作，否则提示用户纠正佩戴位置。
[0091] 例如，当用户佩戴有误时，电极位置出现偏移，此时的波形如图4(a)所示，当用户佩戴无误时，电极位置正常，此时的波形如图4(b)所示。
[0092] Step 2，正常工作后，设备将获取到的ECG信号进行滤波和伪迹处理。
[0093] 具体的，采用梳状滤波器对采集到的ECG信号进行滤波，并根据三轴加速度计拟合运动轨迹(取x、y、z三轴拟合)，判断运动轨迹对应曲线中大于幅度阈值的时长内对应的滤波后的ECG信号中的ECG波形，删除对应的ECG波形，以增强ECG信号的准确性。
[0094] Step 3，ECG有效特征提取，主要是有效QRS波群以及目标时间差的提取。
[0095] 具体的，采用QRS联合检测的方式，首先检测出R波后，依次按照R波前后的极值点确定Q点和S点，根据QR、RS以及QS间距对对应的阈值进行比较，例如，当QR<100ms,且RS<100ms,且QS<150ms时，该QRS波群被确认为是有效QRS波群，否则为无效QRS波群。对于上述QR、RS、QS的数值可以根据实际需要进行设定和调整，上述数值仅为举例说明，不作为对本发明实施例的限定。
[0096] 如图5所示，椭圆标出来的波形也是局部峰值点，但是不满足QR<100ms,且RS<100ms,且QS<150ms，因此椭圆标出来的波形为无效QRS波群，不会被误判为有效QRS波群。
[0097] 进一步地，计算每两个连续的有效QRS波群中两个R波的峰值点R之间的时间差，即RR间隔，如果RR间隔落入0.3s～1.5s，则满足条件，继续判断是否满足下一个条件，否则排除该RR间隔。进一步地，在满足前述条件后，判断相邻RR间隔差是否小于0.12s，若是则保留这两个RR间隔，即这两个RR间隔均为目标时间差，否则，排除偏离RR平均值较大的一个，保留另一个RR间隔作为目标时间差。特别地，RR平均值为之前保留的目标时间差的平均值。对于上述RR间隔、相邻RR间隔差的数值可以根据实际需要进行设定和调整，上述数值仅为举例说明，不作为对本发明实施例的限定。
[0098] 进一步地，将筛选出的所有目标时间差求RR平均值，例如通过上述过程统计得到一分钟内的目标时间差为30个，然后计算这30个目标时间差的平均值，然后根据30个目标时间差的平均值计算心率，其中，心率＝60/30个目标时间差的平均值。
[0099] 因此，采用上述方法去掉可能影响测量结果的波形以提高心电特征提取的准确率，有效保障单臂心电测量设备信号测量的有效性。
[0100] Step 4，ECG信号的信号质量实时评估。
[0101] 首先计算信号参数，包括有效信号功率(PW_ECG)，工频干扰(PW_50)、基线漂移(PW_1)、带内噪声(PW_NB)等，进一步分别求取信号评估参数，SAR＝lg(PW_ECG/PW_1)，以及SNBR＝lg(PW_ECG/PW_NB)，假设TH_SAR＝0dB,TH_SNBR＝10dB，进一步根据表1确定当前ECG信号的信号质量等级，以提示用户采集到的ECG信号的有效性。对于上述TH_SAR、TH_SNBR的数值可以根据实际需要进行设定和调整，上述数值仅为举例说明，不作为对本发明实施例的限定。
[0102] 例如，Step1可根据Step4得到的结果确定当前ECG信号的信号质量等级，若当前信号质量等级为4，可提示用户查找原因，重新测量。
[0103] Step4获得的信号质量等级可作为对Step1中采集信号质量进行预判的标准，通过将信号质量等级反馈给用户，使用户在确定当前信号质量等级较差时纠正佩戴位置，保证采集到的ECG信号的有效性，进而有效保障单臂心电测量设备的信号测量准确性。
[0104] 本发明实施例提供一种可穿戴设备，包括：处理器、存储器、心率采集器，电源。例如，这里的可穿戴设备可以为运动臂带、心率带或心率贴，在可穿戴设备为心率贴时，心率贴可通过蓝牙模块与其他终端(例如，手机)进行通信。
[0105] 下面结合图6对可穿戴设备600的各个构成部件进行具体的介绍：
[0106] 所述心率采集器610，用于采集ECG信号；
[0107] 所述存储器620，用于存储指令；
[0108] 所述处理器630，用于调用所述存储器中的指令，执行上述ECG信号处理过程；
[0109] 所述电源640，用于为所述可穿戴设备供电。
[0110] 在一种可能的设计中，所述可穿戴设备还包括：运动传感器650；
[0111] 所述运动传感器650，用于检测各个方向上的加速度的大小，拟合运动轨迹。
[0112] 在一种可能的设计中，所述可穿戴设备还包括显示单元660，蓝牙模块670等。
[0113] 采用本发明实施例提供的可穿戴设备用于执行上述图1对应的方法实施例，因此本发明实施例提供的可穿戴设备的实施方式可以参见该方法的实施方式，重复之处不再赘述。
[0114] 综上所述，采用本发明实施例提供的方法，结合运动生理变化规律，提高心电特征提取的准确率，有效保障单臂心电测量设备的信号测量准确性。
[0115] 基于同一发明构思，本发明实施例中还提供了一种ECG信号处理装置，该装置可以用于执行上述图1对应的方法实施例，因此本发明实施例提供的ECG信号处理装置的实施方式可以参见该方法的实施方式，重复之处不再赘述。
[0116] 参阅图7所示，本发明实施例提供一种ECG信号处理装置，包括：采集单元710，提取单元720，计算单元730，处理单元740，其中，
[0117] 采集单元710，用于采集ECG信号；
[0118] 提取单元720，用于提取所述ECG信号的第k个有效QRS波群，其中QRS波群包括第i个R波的峰值点R前的第一极值点Q、第i个R波的峰值点R、和第i个R波的峰值点R后的第二极值点S，其中，i≥2，k≥2，k≤i；
[0119] 计算单元730，用于计算第k个时间差和第(k+1)个时间差，其中，所述第k个时间差为所述第k个有效QRS波群中R波的峰值点R与第(k+1)个有效QRS波群中R波的峰值点R的时间差，所述第(k+1)个时间差为所述第(k+1)个有效QRS波群中的R波的峰值点R与第(k+2)个有效QRS波群中的R波的峰值点R的时间差；
[0120] 处理单元740，用于如果所述第k个时间差与所述第(k+1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值小于预设阈值时，确定所述第k个时间差和所述第(k+1)个时间差为目标时间差；或者，如果所述第k个时间差与所述第(k+1)个时间差处于预设范围，且所述第(k+1)个时间差与所述第k个时间差的差值的绝对值大于等于预设阈值时，确定所述第k个时间差和所述第(k+1)个时间差中偏离时间差平均值较少的一个时间差为目标时间差，所述时间差平均值为第1个时间差至第(k-1)个时间差中所有目标时间差的平均值；
[0121] 其中，所述目标时间差用于计算所述ECG信号对应的心率值。
[0122] 在一种可能的实现方式中，所述提取单元720，用于：所述提取所述ECG信号的第k个有效QRS波群时，确定所述第i个R波对应的第一极值点Q和第二极值点S；
[0123] 计算所述第一极值点Q与所述第i个R波的峰值点R之间的时间差为第一时间差，所述第i个R波的峰值点R与所述第二极值点S之间的时间差为第二时间差，所述第一极值点Q与所述第二极值点S之间的时间差为第三时间差；
[0124] 如果所述第一时间差、第二时间差、第三时间差分别小于对应的第一阈值、第二阈值、第三阈值时，则所述第一极值点Q对应Q波，所述第二极值点S对应S波，所述第i个R波与对应的Q波和S波构成第k个有效QRS波群。
[0125] 在一种可能的实现方式中，所述装置，还包括：
[0126] 信号质量评估单元750，用于计算所述ECG信号对应的信号参数；
[0127] 根据至少一个所述信号参数计算信号评估参数；
[0128] 根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级。
[0129] 在一种可能的实现方式中，所述信号参数包括有效信号功率，基线漂移，和带内噪声；
[0130] 所述信号评估参数包括信号伪迹比、和信号带内噪声比，其中，所述信号伪迹比为所述有效信号功率和所述基线漂移的函数，所述信号带内噪声比为所述有效信号功率和所述带内噪声的函数；
[0131] 所述信号质量评估单元750，用于：根据所述信号评估参数与所述信号评估参数对应的预设评估阈值，确定所述ECG信号的信号质量等级时，确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第一等级；
[0132] 确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比大于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第二等级；
[0133] 确定所述信号伪迹比大于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第三等级；
[0134] 确定所述信号伪迹比小于等于对应的预设评估阈值，且所述信号带内噪声比小于等于对应的预设评估阈值时，确定所述ECG信号的信号质量等级为第四等级；
[0135] 其中，所述第一等级优于所述第二等级，所述第二等级优于所述第三等级，所述第三等级优于所述第四等级。
[0136] 在一种可能的实现方式中，所述装置还包括：运动轨迹分析单元760，用于在提取所述ECG信号的第k个有效QRS波群之前，对所述ECG信号进行滤波处理，并采用三轴加速计拟合用户的运动轨迹；
[0137] 将滤波后的ECG信号与所述运动轨迹进行比对，将所述运动轨迹中运动幅度值大于预设幅度阈值的时长对应的所述滤波后的ECG信号中的ECG波形删除。
[0138] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令处理器完成，所述的程序可以存储于计算机可读存储介质中，所述存储介质是非短暂性(英文：non-transitory)介质，例如随机存取存储器，只读存储器，快闪存储器，硬盘，固态硬盘，磁带(英文：magnetic tape)，软盘(英文：floppy disk)，光盘(英文：optical disc)及其任意组合。
[0139] 本发明是参照本发明实施例的方法和设备各自的流程图和方框图来描述的。应理解可由计算机程序指令实现流程图和方框图中的每一流程和方框、以及流程图和方框图中的流程和方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和方框图一个方框或多个方框中指定的功能的装置。
[0140] 以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
