一种呼吸率提取方法及装置的制造方法
【专利摘要】本发明公开了一种呼吸率提取方法,包括以下步骤：获取原始心电信号，并根据所述原始心电信号得到待处理心电信号；对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一呼吸信号计算得到当前时刻的第一呼吸率；对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据所述第二呼吸信号计算得到当前时刻的第二呼吸率；基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进行分析，计算得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二权重因子，并计算得到当前时刻的呼吸率。本发明还公开了一种提出呼吸率提取装置，本发明能够便捷有效提取呼吸信号，计算得到准确稳定的呼吸率。
【专利说明】
一种呼吸率提取方法及装置
技术领域
[0001 ]本发明涉及呼吸检测领域，尤其涉及一种呼吸率提取方法及装置。
【背景技术】
[0002] 呼吸是人体重要的生理过程，对人体呼吸的监护检测也是现代医学监护技术的一 个重要组成部分。患者不论是呼吸系统本身的病变，或是其他重要脏器的病变，发展到一定 程度都会影响呼吸中枢。多脏器系统功能衰竭，往往累及呼吸功能的衰竭，呼吸功能的衰竭 又导致其他脏器功能的衰竭，互为因果。
[0003] 现有技术对呼吸运动主要使用下列方法检测：阻抗容积法，用高频恒流源测量胸 部阻抗的变化来提取呼吸信息;传感器法，使用温度、压力、湿度和气流传感器作为鼻孔传 感器；电容法，当呼吸时导致电容值产生相应的变化;呼吸音法，通过拾取呼吸音识别呼吸； 超声法，利用超声波产生多谱勒现象，检测出呼吸频率。使用这些方法不但需要增加信号采 集部件，而且受到运动和环境的影晌，不适合用于日常监护。
[0004] 大量临床资料显示，呼吸运动会引起心电图的变化。通过心电图，我们可以观察到 在呼吸周期内由胸部运动和心脏位置变化所引起的心电波形峰峰值的改变。这是由于呼吸 周期内，描述心脏电波主要传播方向的心脏电轴旋转造成QRS波群形态发生了变化。QRS波 是指正常心电图中幅度最大的波群，反映心室除极的全过程。正常心室除极始于室间隔中 部，自左向右方向除极，故QRS波群先呈现一个小向下的q波。正常胸导联QRS波群形态较恒 定。从心电信号中提取呼吸信号(ECG - DerivedRespirationjDR)是一种呼吸信号检测技 术，这种技术不需要专用传感器和硬件模块检测呼吸信号，只需要用心电监护仪获取心电 信号，避免了上述两种检测方法对人体的束缚，使动态呼吸检测成为可能。
[0005] 但现有从心电信号中提取呼吸信号的技术，在计算时主要采用波形法，该方法通 过一段时间内波形的平均值（即基线值），来判定当前呼吸波处于上升或下降趋势，用极值 的方法求得波形的波峰、波谷。根据一定的阈值条件来判定有效的波峰或波谷，再根据有效 波峰或波谷的周期计算波形周期，从而得到呼吸率。这种算法虽然具有比较直观、运算量小 的优点，但在实际过程中获取的呼吸波形或多或少会受到心电活动的影响，当波形出现基 线漂移时，计算的基线值无法很快更新，会导致波形漏检致使呼吸率值偏低，其结果会有较 大偏差。
【发明内容】
[0006] 针对上述问题，本发明的目的在于提供一种呼吸提取方法及装置，可在不受运动 和环境的影响的情况下，实现呼吸率准确稳定的检测。
[0007] 本发明提供了 一种呼吸率提取方法，所述方法包括：
[0008] 获取原始心电信号，并对所述原始心电信号进行工频陷后得到待处理心电信号；
[0009] 对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一呼吸信号计 算得到当前时刻的第一呼吸率；
[0010] 对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据所述第二呼吸信 号计算得到当前时刻的第二呼吸率；
[0011] 基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进行分析，计算得到 与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二权重因子，以根据 所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率。
[0012] 优选地，在对所述待处理心电信号进行小波变换得到第二呼吸信号之前，还包括：
[0013] 对所述待处理心电信号进行降采样。
[0014] 优选地，所述对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第 一呼吸信号计算得到当前时刻的第一呼吸率具体包括：
[0015] 对所述待处理心电信号进行滤波，获取所述待处理心电信号位于预定的通带频率 之间的数据，得到第一呼吸信号；
[0016] 从所述第一呼吸信号中提取最近生成的两个波峰之间的时间间隔；
[0017] 根据所述时间间隔进行采样率换算得到当前时刻的第一呼吸率。
[0018] 优选地，所述对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据所 述第二呼吸信号计算得到当前时刻的第二呼吸率具体包括：
[0019] 根据香农-奈奎斯特采样原理及所述待处理心电信号的采样频率进行频段分层， 计算得到每层的频率范围；
[0020] 依据所述频段分层中每层的频率范围及预置的通带频率确定小波分解和重构所 需的层数；
[0021] 根据与所述小波分解所需的层数及预先选择的母小波进行信号分解，得到按频段 划分的多层波形；
[0022] 根据与所述小波重构所需的层数对应的系数及分解得到的所述多层波形进行信 号重构，得到第二呼吸信号；
[0023] 从所述第二呼吸信号中提取最近生成的两个波峰之间的时间间隔；
[0024]根据所述时间间隔进行采样率换算得到当前时刻的第二呼吸率。
[0025]优选地，所述通带频率的上限截止频率为0.8Hz，所述通带频率的下限截止频率为 O-IHz0
[0026] 优选地，所述基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进行分 析，计算得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二权重 因子，以根据所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率，具体包括：
[0027] 对所述第一呼吸信号和所述第二呼吸信号进行功率谱分析，分析所述第一呼吸信 号和所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应的第一权重因子和与第二呼 吸率对应的第二权重因子；
[0028] 当判断所述第一权重因子大于预设的基准值且所述第二权重因子小于所述基准 值时，将所述第一呼吸率设置为当前时刻的呼吸率；
[0029] 当判断所述第一权重因子小于预设的基准值且所述第二权重因子大于所述基准 值时，将所述第二呼吸率设置为当前时刻的呼吸率；
[0030] 当判断所述第一权重因子计所述第二权重因子均大于预设的基准值时，根据所述 第一权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率进行加权求和，计算得 到当前时刻的呼吸率。
[0031 ]本发明还提供了 一种呼吸率提取装置，包括：
[0032]工频陷波单元，用于获取原始心电信号，并对所述原始心电信号进行工频陷波后 得到待处理心电信号；
[0033] 滤波单元，对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一 呼吸信号计算得到当前时刻的第一呼吸率；
[0034] 小波变换单元，对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据 预设的呼吸算法及所述第二呼吸信号计算得到当前时刻的第二呼吸率；
[0035] 呼吸率计算单元，用于基于信号质量指数对所述第一呼吸信号及所述第二呼吸信 号进行分析，计算得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的 第二权重因子，以根据所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率。
[0036] 优选地，所述呼吸率提取装置还包括：
[0037]降采样单元，用于对所述待处理心电信号进行降采样。
[0038]优选地，所述滤波单元具体包括：
[0039]滤波子模块，用于对所述待处理心电信号进行滤波，获取所述待处理心电信号位 于预定的通带频率之间的数据，得到第一呼吸信号；
[0040] 第一时间间隔提取模块，用于从所述第一呼吸信号中提取最近生成的两个波峰之 间的时间间隔；
[0041] 第一呼吸率计算模块，用于根据所述时间间隔进行采样率换算得到当前时刻的第 一呼吸率。
[0042] 优选地，所述小波变换单元具体包括：
[0043]频段分成模块，用于根据香农-奈奎斯特采样原理及所述待处理心电信号的采样 频率进行频段分层，计算得到每层的频率范围；
[0044] 层数确定模块，用于依据所述频段分层每层的频率范围及通带频率确定小波分解 和重构所需的层数；
[0045] 信号分解模块，用于根据与所述小波分解所需的层数及预先选择的母小波进行信 号分解，得到按频段划分的多层波形；
[0046] 信号重构模块，用于根据与所述小波重构所需的层数对应的系数及分解得到的所 述多层波形进行信号重构，得到第二呼吸信号；
[0047] 第二时间间隔提取模块，用于从所述第二呼吸信号中提取最近生成的两个波峰之 间的时间间隔；
[0048] 第二呼吸率计算模块，用于根据所述时间间隔进行采样率换算得到当前时刻的第 二呼吸率。
[0049] 优选地，所述呼吸率计算单元具体包括：
[0050] 权重因子计算模块，用于对所述第一呼吸信号和所述第二呼吸信号进行功率谱分 析，分析所述第一呼吸信号和所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应的 第一权重因子和与第二呼吸率对应的第二权重因子；
[0051] 第一判断模块，用于当判断所述第一权重因子大于预设的基准值且所述第二权重 因子小于所述基准值时，将所述第一呼吸率设置为当前时刻的呼吸率；
[0052]第二判断模块，用于当判断所述第一权重因子小于预设的基准值且所述第二权重 因子大于所述基准值时，将所述第二呼吸率设置为当前时刻的呼吸率；
[0053]第三判断模块，用于当判断所述第一权重因子计所述第二权重因子均大于预设的 基准值时，根据所述第一权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率进 行加权求和，计算得到当前时刻的呼吸率。
[0054] 本发明提供的呼吸率提取方法及装置，通过利用小波变换技术与滤波技术相结合 的方式处理心电信号得到第一呼吸率及第二呼吸率，并根据与所述第一呼吸率对应的第一 权重因子和与所述第二呼吸率对应的第二权重因子得到当前时刻的呼吸率，相比于现有由 单一技术从心电信号得到呼吸信号的方案，计算结果更准确可靠，并可减轻由于外界或环 境的干扰而引起的测量波动或误差，从而能够得到更为准确稳定的测量结果。
【附图说明】
[0055] 图1是根据本发明一实施例的呼吸率提取方法的流程图；
[0056] 图2是根据本发明一实施例的原始心电信号的波形图；
[0057]图3是根据本发明一实施例的待处理心电信号的波形图；
[0058]图4是根据本发明一实施例的滤波后提取得到的第一呼吸信号的波形图；
[0059]图5是根据本发明一实施例的小波变换提取得到的第二呼吸信号的波形图。
[0060]图6是根据本发明一实施例的呼吸率提取装置的结构示意图。
[0061 ]图7是根据本发明一实施例的滤波单元的结构示意图。
[0062] 图8是根据本发明一实施例的小波变换单元的结构示意图。
[0063] 图9是根据本发明一实施例的呼吸率计算单元的结构示意图。
【具体实施方式】
[0064]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0065] 本发明提供了一种呼吸率提取方法，用于从心电信号中提取出呼吸信息，由于呼 吸作用引起的心电图中的基线漂移，把呼吸信息看作是心电信号的低频成分，通过去除呼 吸频率以外的信号，从而得到所需提取的呼吸信息。
[0066] 参见图1，本发明实施例提供了一种呼吸率提取方法，其包括如下步骤：
[0067] Sl，获取原始心电信号，并对所述原始心电信号进行工频陷波后得到待处理心电 信号。
[0068]参见图2,为根据本发明一实施例的原始心电信号的波形图。所述原始心电信号包 含大量的工频干扰，需要进行50Hz工频陷波，以滤除工频干扰，参见图3，为根据本发明一实 施例的工频陷波后的待处理心电信号的波形图。
[0069]需要说明的是，在本发明实施例中，由于需要对待处理心电信号进行小波变换，若 所述原始心电信号的采用频率较大(通常为500Hz)，则会导致小波分解的层数太高，计算量 增大，不利于后续的小波变换。因此，需对工频陷波后待处理心电信号进行降采样处理。
[0070] 在本发明实施例中，可将所述待处理心电信号降采样至100Hz。当然，也可根据实 际的需要降采样至其他频率，如50Hz，200Hz等，本发明不做具体限定。
[0071] S2,对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一呼吸信 号计算得到当前时刻的第一呼吸率。
[0072]在本发明实施例中，由于呼吸信号的频段范围通常为0.1~0.4Hz，考虑到呼吸急 促的情况，可将频段扩展为0.1~0.8HZ。参见图4,为根据本发明一实施例的滤波器提取得 到的第一呼吸信号的波形图，其中，滤波时使用的通带频率的下限截止频率Fmin = O. 1Hz， 通带频率的上限截止频率Fmax = O. 8Hz。当然，应当理解的是，在本发明的其他实施例中，也 可根据测量用户的年龄、职业或身体状况等对所述通带频率的上限截止频率和下限截止进 行调整，这些方案均在本发明的保护范围之内。
[0073] 在本发明实施例中，在获得所述第一呼吸信号后，即可根据预设的呼吸算法获得 第一呼吸率Rl，具体为：
[0074] 通过求极值法在所述第一呼吸信号的波形图中寻找第一呼吸信号的波峰(或者波 谷），参见图4中的点标记；
[0075] 通过提取最近生成的两个波峰之间的时间间隔，以得到当前时刻的周期T。
[0076] 对所述周期进行采样率换算即可得到当前时刻的第一呼吸率Rl。
[0077] 例如，Rl=6〇/Tl。
[0078] S3,对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据所述第二呼 吸信号计算得到当前时刻的第二呼吸率。
[0079] 在本发明实施例中，步骤S3具体包括：
[0080] S31，根据香农-奈奎斯特采样原理及所述待处理心电信号的采样频率进行频段分 层，计算得到每层的频率范围。
[0081] 根据香农-奈奎斯特采样原理，设待处理心电信号的采样频率为fS，目标频段为 fl-f2(Hz)，运用小波变换分解的层数为N，由奈奎斯特定律可知，则有：
[0082] fl = (fs/2)/ (1)
[0083] f2 = (fs/2)/ (2)
[0084] N>N1(N1>N2) (3)
[0085] 即需要重构的小波的层数为N2~NI层。
[0086] S32,依据所述频段分层中每层的频率范围及预置的通带频率确定小波分解和重 构所需的层数。
[0087] 假设所述待处理心电信号被降采样至IOOHz，则f s为IOOHz，信号最高频率为50Hz， 根据公式(1 )、（2)、（3)可知，每一层对应的频段如下：
[0089] 由于呼吸信号的频段范围通常为0.1~0.4Hz，考虑到呼吸急促的情况，将频段扩 展为0.1~0.8Hz，心电的频率范围为0.9~6Hz，因此，能够很好地分离呼吸信号和心电信 号，所以，选用第9、8、7层的近似系数(D9/D8/D7)来重构信号。
[0090] S33,根据与所述小波分解所需的层数及预先选择的母小波进行信号分解，得到按 频段划分的多层波形。
[0091] 在本发明实施例中，经验证，COifN小波和dmey小波的提取效果较佳，且优选地，以 c〇if3小波基作为母小波时，具有最佳的提取效果。因而本发明实施例采用coif3小波基作 为母小波进行小波分解。当然，可以理解的是，在本发明的其他实施例中，也可选取其他的 母小波，如db小波等，本发明不做具体限定。
[0092]在本发明实施例中，在分解时，可得到N层波形，此时，可提取N2~Nl层对应的波形 来进行重构。
[0093] S34,根据与所述小波重构所需的层数对应的系数及分解得到的所述多层波形进 行信号重构，得到第二呼吸信号。
[0094] 由步骤S33可知，小波重构所需的层数为第9、8、7层，此时，即可根据与所述小波重 构所需的层数对应的系数及分解得到的所述多层波形进行信号重构，得到第二呼吸信号Z ⑴。
[0095] 参见图5,为根据本发明一实施例的小波变换提取得到的第二呼吸信号的波形图。
[0096] 在本发明实施例中，在获得所述第二呼吸信号后，即可根据预设的呼吸算法获得 第二呼吸率R2,具体为
[0097] 通过求极值法在所述第二呼吸信号的波形图中寻找第二呼吸信号的波峰(或者波 谷），参见图5中的点标记；
[0098] 通过提取最近生成的两个波峰之间的时间间隔，以得到周期T2。
[0099] 根据采样率换算即可得到实时的第二呼吸率R2。
[0100] S4,基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进行分析，计算 得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二权重因子，以 根据所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率。
[0101] 具体地，在本发明实施例中，计算得到当前时刻的呼吸率具体包括：
[0102] S41，对所述第一呼吸信号和所述第二呼吸信号进行功率谱分析，分析所述第一呼 吸信号和所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应的第一权重因子μL和与 第二呼吸率对应的第二权重因子μ2。
[0103] S42,当判断所述第一权重因子大于预设的基准值且所述第二权重因子小于所述 基准值时，将所述第一呼吸率设置为当前时刻的呼吸率。
[0104] 当所述第二权重因子小于所述基准值时，可以认为第二呼吸信号的信号质量较 差，此时，直接将所述第一呼吸率Rl设置为当前时刻的呼吸率R。
[0105] S43,当判断所述第一权重因子小于预设的基准值且所述第二权重因子大于所述 基准值时，将所述第二呼吸率设置为当前时刻的呼吸率。
[0106] 当所述第一权重因子小于所述基准值时，可以认为第一呼吸信号的信号质量较 差，此时，直接将所述第一呼吸率Rl设置为当前时刻的呼吸率R。
[0107] S44,当判断所述第一权重因子计所述第二权重因子均大于预设的基准值时，根据 所述第一权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率进行加权求和，计 算得到当前时刻的呼吸率。
[0108] 此时，当前时刻的呼吸率R=yl*Rl+y2*R2。
[0109] 本发明提供的呼吸率提取方法，通过利用小波变换技术与滤波技术相结合的方式 处理心电信号得到第一呼吸率及第二呼吸率，并根据与所述第一呼吸率对应的第一权重因 子和与所述第二呼吸率对应的第二权重因子得到当前时刻的呼吸率，相比于现有由单一技 术从心电信号得到呼吸信号的方案，计算结果更准确可靠，并可减轻由于外界或环境的干 扰而引起的测量波动或误差，从而能够得到更为稳定的测量结果。
[0110] 请参阅图6,本发明实施例还提供一种呼吸率提取装置100,包括：
[0111] 工频陷波单元10,用于获取原始心电信号，并根据所述原始心电信号得到待处理 心电信号。
[0112]滤波单元20,对所述待处理心电信号进行滤波得到第一呼吸信号，并根据预设的 呼吸算法及所述第一呼吸信号计算得到当前时刻的第一呼吸率。
[0113] 小波变换单元30,对所述待处理心电信号进行小波变换得到第二呼吸信号，并根 据预设的呼吸算法及所述第二呼吸信号计算得到当前时刻的第二呼吸率。
[0114] 呼吸率计算单元40,用于基于信号质量指数对所述第一呼吸信号及所述第二呼吸 信号进行分析，计算得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应 的第二权重因子，以根据所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率。
[0115] 其中，优选地，所述呼吸率提取装置100还包括：
[0116] 降采样单元，用于对所述待处理心电信号进行降采样。
[0117]其中，如图7所示，优选地，所述滤波单元20具体包括：
[0118]滤波子模块21，用于对所述待处理心电信号进行滤波，获取所述待处理心电信号 位于预定的通带频率之间的数据，得到第一呼吸信号；
[0119] 第一时间间隔提取模块22,用于从所述第一呼吸信号中提取最近生成的两个波峰 之间的时间间隔；
[0120] 第一呼吸率计算模块23,用于根据所述时间间隔进行采样率换算得到当前时刻的 第一呼吸率。
[0121] 其中，如图8所示，优选地，所述小波变换单元30具体包括：
[0122] 频段分成模块31，用于根据香农-奈奎斯特采样原理及所述待处理心电信号的采 样频率进行频段分层，计算得到每层的频率范围；
[0123] 层数确定模块32,用于依据所述频段分层每层的频率范围及通带频率确定小波分 解和重构所需的层数；
[0124] 信号分解模块33,根据与所述小波分解所需的层数及预先选择的母小波进行信号 分解，得到按频段划分的多层波形。
[0125] 信号重构模块34,用于根据与所述小波重构所需的层数对应的系数及分解得到的 所述多层波形进行信号重构，得到第二呼吸信号。
[0126] 第二时间间隔提取模块35,用于从所述第二呼吸信号中提取最近生成的两个波峰 之间的时间间隔；
[0127] 第二呼吸率计算模块36,用于根据所述时间间隔进行采样率换算得到当前时刻的 第二呼吸率。
[0128] 其中，如图9所示，优选地，所述呼吸率计算单元40具体包括：
[0129] 权重因子计算模块41，用于对所述第一呼吸信号和所述第二呼吸信号进行功率谱 分析，分析所述第一呼吸信号和所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应 的第一权重因子和与第二呼吸率对应的第二权重因子；
[0130]第一判断模块42,用于当判断所述第一权重因子大于预设的基准值且所述第二权 重因子小于所述基准值时，将所述第一呼吸率设置为当前时刻的呼吸率；
[0131]第二判断模块43,用于当判断所述第一权重因子小于预设的基准值且所述第二权 重因子大于所述基准值时，将所述第二呼吸率设置为当前时刻的呼吸率；
[0132] 第三判断模块44,用于当判断所述第一权重因子计所述第二权重因子均大于预设 的基准值时，根据所述第一权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率 进行加权求和，计算得到当前时刻的呼吸率。
[0133] 本发明提供的呼吸率提取装置100,通过利用小波变换技术与滤波技术相结合的 方式处理心电信号得到第一呼吸率及第二呼吸率，并通过加权求和得到当前时刻的呼吸 率，相比于现有由单一技术从心电信号得到呼吸信号的方案，计算结果更准确可靠，并可减 轻由于外界或环境的干扰而引起的测量波动或误差，从而能够得到更为准确稳定的测量结 果。
[0134] 以上所揭露的仅为本发明一种较佳实施例而已，当然不能以此来限定本发明之权 利范围，本领域普通技术人员可以理解实现上述实施例的全部或部分流程，并依本发明权 利要求所作的等同变化，仍属于发明所涵盖的范围。
[0135] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以 通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质 中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁 碟、光盘、只读存储记忆体（Read-Only Memory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。
【主权项】
1. 一种呼吸率提取方法，其特征在于，所述方法包括： 获取原始心电信号，并对所述原始心电信号进行工频陷波后得到待处理心电信号； 对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一呼吸信号计算得 到当前时刻的第一呼吸率； 对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据所述第二呼吸信号计 算得到当前时刻的第二呼吸率； 基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进行分析，计算得到与所 述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二权重因子，以根据所述 第一权重因子及所述第二权重因子得到当前时刻的呼吸率。2. 根据权利要求1所述的呼吸率提取方法，其特征在于，在对所述待处理心电信号进行 小波变换得到第二呼吸信号之前，还包括： 对所述待处理心电信号进行降采样。3. 根据权利要求1所述的呼吸率提取方法，其特征在于，所述对所述待处理心电信号进 行滤波得到第一呼吸信号，并根据所述第一呼吸信号计算得到当前时刻的第一呼吸率具体 包括： 对所述待处理心电信号进行滤波，获取所述待处理心电信号位于预定的通带频率之间 的数据，得到第一呼吸信号； 从所述第一呼吸信号中提取最近生成的两个波峰之间的时间间隔； 根据所述时间间隔进行采样率换算得到当前时刻的第一呼吸率。4. 根据权利要求1所述的呼吸率提取方法，其特征在于，所述对所述待处理心电信号进 行小波变换得到第二呼吸信号，并根据所述第二呼吸信号计算得到当前时刻的第二呼吸率 具体包括： 根据香农-奈奎斯特采样原理及所述待处理心电信号的采样频率进行频段分层，计算 得到每层的频率范围； 依据所述频段分层中每层的频率范围及预置的通带频率确定小波分解和重构所需的 层数； 根据与所述小波分解所需的层数及预先选择的母小波进行信号分解，得到按频段划分 的多层波形； 根据与所述小波重构所需的层数对应的系数及分解得到的所述多层波形进行信号重 构，得到第二呼吸信号； 从所述第二呼吸信号中提取最近生成的两个波峰之间的时间间隔； 根据所述时间间隔进行采样率换算得到当前时刻的第二呼吸率。5. 根据权利要求3或4所述的呼吸率提取方法，其特征在于，所述通带频率的上限截止 频率为0.8Hz，所述通带频率的下限截止频率为0.1Hz。6. 根据权利要求1所述的呼吸率提取方法，其特征在于，所述基于信号质量指数对所述 第一呼吸信号及所述第二呼吸信号进行分析，计算得到与所述第一呼吸率对应的第一权重 因子和与所述第二呼吸率对应的第二权重因子，以根据所述第一权重因子及所述第二权重 因子得到当前时刻的呼吸率，具体包括： 对所述第一呼吸信号和所述第二呼吸信号进行功率谱分析，分析所述第一呼吸信号和 所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应的第一权重因子和与第二呼吸率 对应的第二权重因子； 当判断所述第一权重因子大于预设的基准值且所述第二权重因子小于所述基准值时， 将所述第一呼吸率设置为当前时刻的呼吸率； 当判断所述第一权重因子小于预设的基准值且所述第二权重因子大于所述基准值时， 将所述第二呼吸率设置为当前时刻的呼吸率； 当判断所述第一权重因子计所述第二权重因子均大于预设的基准值时，根据所述第一 权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率进行加权求和，计算得到当 前时刻的呼吸率。7. -种呼吸率提取装置，其特征在于，包括： 工频陷波单元，用于获取原始心电信号，并对所述原始心电信号进行工频陷波后得到 待处理心电信号； 滤波单元，对所述待处理心电信号进行滤波得到第一呼吸信号，并根据所述第一呼吸 信号计算得到当前时刻的第一呼吸率； 小波变换单元，对所述待处理心电信号进行小波变换得到第二呼吸信号，并根据预设 的呼吸算法及所述第二呼吸信号计算得到当前时刻的第二呼吸率； 呼吸率计算单元，用于基于信号质量指数对所述第一呼吸信号及所述第二呼吸信号进 行分析，计算得到与所述第一呼吸率对应的第一权重因子和与所述第二呼吸率对应的第二 权重因子，以根据所述第一权重因子及所述第二权重因子得到当前时刻的呼吸率。8. 根据权利要求7所述的呼吸率提取装置，其特征在于，所述呼吸率提取装置还包括： 降采样单元，用于对所述待处理心电信号进行降采样。9. 根据权利要求7所述的呼吸率提取装置，其特征在于，所述滤波单元具体包括： 滤波子模块，用于对所述待处理心电信号进行滤波，获取所述待处理心电信号位于预 定的通带频率之间的数据，得到第一呼吸信号； 第一时间间隔提取模块，用于从所述第一呼吸信号中提取最近生成的两个波峰之间的 时间间隔； 第一呼吸率计算模块，用于根据所述时间间隔进行采样率换算得到当前时刻的第一呼 吸率。10. 根据权利要求7所述的呼吸率提取装置，其特征在于，所述小波变换单元具体包括： 频段分成模块，用于根据香农-奈奎斯特采样原理及所述待处理心电信号的采样频率 进行频段分层，计算得到每层的频率范围； 层数确定模块，用于依据所述频段分层每层的频率范围及通带频率确定小波分解和重 构所需的层数； 信号分解模块，用于根据与所述小波分解所需的层数及预先选择的母小波进行信号分 解，得到按频段划分的多层波形； 信号重构模块，用于根据与所述小波重构所需的层数对应的系数及分解得到的所述多 层波形进行信号重构，得到第二呼吸信号； 第二时间间隔提取模块，用于从所述第二呼吸信号中提取最近生成的两个波峰之间的 时间间隔； 第二呼吸率计算模块，用于根据所述时间间隔进行采样率换算得到当前时刻的第二呼 吸率。11.根据权利要求7所述的呼吸率提取装置，其特征在于，所述呼吸率计算单元具体包 括： 权重因子计算模块，用于对所述第一呼吸信号和所述第二呼吸信号进行功率谱分析， 分析所述第一呼吸信号和所述第二呼吸信号的谱分布，得到与所述第一呼吸率对应的第一 权重因子和与第二呼吸率对应的第二权重因子； 第一判断模块，用于当判断所述第一权重因子大于预设的基准值且所述第二权重因子 小于所述基准值时，将所述第一呼吸率设置为当前时刻的呼吸率； 第二判断模块，用于当判断所述第一权重因子小于预设的基准值且所述第二权重因子 大于所述基准值时，将所述第二呼吸率设置为当前时刻的呼吸率； 第三判断模块，用于当判断所述第一权重因子计所述第二权重因子均大于预设的基准 值时，根据所述第一权重因子及所述第二权重因子对所述第一呼吸率和第二呼吸率进行加 权求和，计算得到当前时刻的呼吸率。
【文档编号】A61B5/08GK106073784SQ201610679973
【公开日】2016年11月9日
【申请日】2016年8月17日 公开号201610679973.1, CN 106073784 A, CN 106073784A, CN 201610679973, CN-A-106073784, CN106073784 A, CN106073784A, CN201610679973, CN201610679973.1
【发明人】胡静
【申请人】广州视源电子科技股份有限公司