CN107811612A - 一种心率与呼吸率数据处理方法 - Google Patents

一种心率与呼吸率数据处理方法 Download PDF

Info

Publication number
CN107811612A
CN107811612A CN201711232535.1A CN201711232535A CN107811612A CN 107811612 A CN107811612 A CN 107811612A CN 201711232535 A CN201711232535 A CN 201711232535A CN 107811612 A CN107811612 A CN 107811612A
Authority
CN
China
Prior art keywords
waveform
signal
data
frequency domain
time domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201711232535.1A
Other languages
English (en)
Inventor
张法全
罗明会
毛学港
王国富
叶金才
王小红
贾小波
王晶晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN201711232535.1A priority Critical patent/CN107811612A/zh
Publication of CN107811612A publication Critical patent/CN107811612A/zh
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816—Measuring devices for examining respiratory frequency
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physiology (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Cardiology (AREA)
  • Pulmonology (AREA)
  • Power Engineering (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开了一种心率与呼吸率数据处理方法，根据人体体征数据的特点提取采样信号的时域及频域特征重构信号波形并通过不断波形拟合方式进行修正，包括如下步骤：1）让一个采样周期的数据通过滤波器进行降噪处理得到预处理数据；2）对预处理数据进行峰值检测以获得信号的时域幅值，同时对预处理数据作FFT变换以获取其频域特征；3）利用信号频域特征与时域幅值采用内插法重构波形；4）与下一个采样周期的重构波形拟合得到信号的修正波形；5）修正后的波形又将与其之后的波形拟合得到新的重构波形，重复若干次，预测误差逐渐减小，精度逐渐提高。本发明具有预测误差小，精度高的优点，且数据处理流程计算量适中，适用于算力受限的便携式设备。

Description

一种心率与呼吸率数据处理方法
技术领域
本发明涉及心率、呼吸率信号的处理技术领域，尤其涉及一种心率与呼吸率数据处理方法。
背景技术
目前对于心率、呼吸率信号的处理普遍采用傅里叶变换进行滤波处理后，检测数据的峰值，再根据信号峰值点估计采样周期内信号的波形数，最终利用此估计值除以采样周期得到近似的心率与呼吸率。这种算法的优点是直观且运算量小，但是由于采样时间窗口内的波形数不可能是整数值，因而其估计值与实际偏差较大。
发明内容
针对现有技术的不足，本发明所解决的问题是提供一种预测误差小，精度高的心率与呼吸率数据处理方法。
为解决上述技术问题，本发明采用的技术方案是一种心率与呼吸率数据处理方法，根据人体体征数据的特点提取采样信号的时域及频域特征重构信号波形并通过不断波形拟合方式进行修正，包括如下步骤：
1）让一个采样周期时间窗口内的原始采样数据通过一个二阶巴特沃斯带通滤波器进行降噪处理的到预处理数据；
2）对预处理数据一方面进行峰值检测以获得信号的时域幅值，另一方面对预处理数据作FFT变换以获取其频域特征；
3）利用当前信号的频域特征与获得信号的时域幅值采用内插法重构信号波形；
4）与下一个采样周期的重构波形拟合修正得到目标信号重构修正波形；
5）修正后的波形又将与其之后的波形拟合得到新的重构波形，重复拟合次数若干次，预测误差将逐渐减小，精度逐渐提高。
采用本发明的技术方案有益效果：
1）利用信号的时域与频域特征重构信号波形，可以有效的减小短时采样窗口内对信号波形的估计误差；
2）在信号的重构拟合过程中可以通过设定阈值有选择性的舍去如体动干扰信号等偏差明显的数据点，可以最大限度减小无用信号的误差累积，提高计算结果的可靠性；
3）整个数据的处理流程复杂度不高且计算量适中，适用于运算能力受限的便携式处理设备。
附图说明
图1为本发明流程示意图。
具体实施方式
下面结合附图对本发明作进一步的说明，但不是对本发明的限定。
图1示出了一种心率与呼吸率数据处理方法，根据人体体征数据的特点提取采样信号的时域及频域特征重构信号波形并通过不断波形拟合方式进行修正，包括如下步骤：
1）让一个采样周期时间窗口内的原始采样数据通过一个二阶巴特沃斯带通滤波器进行降噪处理的到预处理数据；
2）对预处理数据一方面进行峰值检测以获得信号的时域幅值，另一方面对预处理数据作FFT变换以获取其频域特征；
参考成人静息状态下的心率为每分钟60～100次/分，呼吸率为16～20次/分，从而前述带通滤波器的截止频率可根据当前处理的信号类型可能出现的频率范围自动调整，如呼吸率信号处理时滤波器的截止频率可设定为0.1Hz～0.7Hz，心率信号处理时滤波器的截止频率可设定为0.7Hz～3.4Hz。
在时域分析中对每个数据点进行峰值检测，得到有效的波峰与波谷数据的幅值；在频域分析中作FFT变换以获取其频域特征。
3）利用当前信号的频域特征与获得信号的时域幅值采用内插法重构信号波形；
根据此波形特征对心率或者呼吸率作第一次预估。
4）与下一个采样周期的重构波形拟合修正得到目标信号重构修正波形；此时根据新的波形对心率或者呼吸率作第二次预估。
5）修正后的波形又将与其之后的波形拟合得到新的重构波形，重复拟合次数若干次，预测误差将逐渐减小，精度逐渐提高。
采用本发明的技术方案有益效果：
1）利用信号的时域与频域特征重构信号波形，可以有效的减小短时采样窗口内对信号波形的估计误差；
2）在信号的重构拟合过程中可以通过设定阈值有选择性的舍去如体动干扰信号等偏差明显的数据点，可以最大限度减小无用信号的误差累积，提高计算结果的可靠性；
3）整个数据的处理流程复杂度不高且计算量适中，适用于运算能力受限的便携式处理设备。
以上结合附图对本发明的实施方式做出了详细说明，但本发明不局限于所描述的实施方式。对于本领域技术人员而言，在不脱离本发明的原理和精神的情况下，对这些实施方式进行各种变化、修改、替换和变型仍落入本发明的保护范围内。

Claims (1)

1.一种心率与呼吸率数据处理方法，其特征在于：根据人体体征数据的特点提取采样信号的时域及频域特征重构信号波形并通过不断波形拟合方式进行修正，包括如下步骤：
1）让一个采样周期时间窗口内的原始采样数据通过一个二阶巴特沃斯带通滤波器进 行降噪处理的到预处理数据；
2）对预处理数据一方面进行峰值检测以获得信号的时域幅值，另一方面对预处理数据作FFT变换以获取其频域特征；
3）利用当前信号的频域特征与获得信号的时域幅值采用内插法重构信号波形；
4）与下一个采样周期的重构波形拟合修正得到目标信号重构修正波形；
5）修正后的波形又将与其之后的波形拟合得到新的重构波形，重复拟合次数若干次，预测误差将逐渐减小，精度逐渐提高。
CN201711232535.1A 2017-11-30 2017-11-30 一种心率与呼吸率数据处理方法 Pending CN107811612A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711232535.1A CN107811612A (zh) 2017-11-30 2017-11-30 一种心率与呼吸率数据处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711232535.1A CN107811612A (zh) 2017-11-30 2017-11-30 一种心率与呼吸率数据处理方法

Publications (1)

Publication Number Publication Date
CN107811612A true CN107811612A (zh) 2018-03-20

Family

ID=61605456

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711232535.1A Pending CN107811612A (zh) 2017-11-30 2017-11-30 一种心率与呼吸率数据处理方法

Country Status (1)

Country Link
CN (1) CN107811612A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113208574A (zh) * 2021-04-30 2021-08-06 中科院计算所泛在智能研究院 一种人体特征震动波形提取分析装置及其使用方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101437009A (zh) * 2007-11-15 2009-05-20 华为技术有限公司 丢包隐藏的方法及其系统
US20170049402A1 (en) * 2011-05-02 2017-02-23 The Regents Of The University Of California System and method for reconstructing cardiac activation information

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101437009A (zh) * 2007-11-15 2009-05-20 华为技术有限公司 丢包隐藏的方法及其系统
US20170049402A1 (en) * 2011-05-02 2017-02-23 The Regents Of The University Of California System and method for reconstructing cardiac activation information

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周礼杲，等: "基于匹配滤波的胎儿心电监护仪的研制", 《中国医疗器械杂志》 *
郑南宁: "《数字信号处理简明教程》", 30 September 2015, 西安交通大学出版社 *
龙园，等: "基于太赫兹技术的植物叶片水分检测初步研究", 《光谱学与光谱分析》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113208574A (zh) * 2021-04-30 2021-08-06 中科院计算所泛在智能研究院 一种人体特征震动波形提取分析装置及其使用方法

Similar Documents

Publication Publication Date Title
CN107184203B (zh) 基于自适应集合经验模态分解的心电信号特征点识别方法
CN109907752B (zh) 一种去除运动伪影干扰与心电特征检测的心电诊断与监护系统
CN109522826B (zh) 一种基于fmcw毫米波雷达的生命信号检测方法和系统
CN101138507B (zh) 多普勒血流声音信号的处理方法及装置
CN105508276B (zh) 一种城区变电站主变压器风扇噪声评价方法
CN104644165B (zh) 一种佩戴式的脑电采集装置
CN110680308B (zh) 基于改进emd与阈值法融合的心电信号去噪方法
CN106889984B (zh) 一种心电信号自动降噪方法
Kumar et al. Performance analysis of different filters for power line interface reduction in ECG signal
Soerensen et al. A comparison of IIR and wavelet filtering for noise reduction of the ECG
CN102928517A (zh) 一种基于小波分解阈值去噪的瓷绝缘子振动声学检测数据降噪的方法
CN105078505A (zh) 生理信号处理方法及处理装置
Hu et al. Removal of baseline wander from ECG signal based on a statistical weighted moving average filter
CN104367316A (zh) 基于形态学滤波与提升小波变换的心电信号去噪方法
CN103110422A (zh) 基于生物雷达检测的呼吸和心跳实时分离方法
ATE416670T1 (de) Vorrichtung zur rauschverminderung eines fötalen ekg-signals
CN103417206B (zh) Ecg去除工频干扰的方法及系统
CN105030232A (zh) 一种心电信号的基线漂移校正方法
Taralunga et al. Fetal ECG enhancement: adaptive power line interference cancellation based on Hilbert Huang Transform
CN104702244A (zh) 基于eemd算法滤除肌电信号检测中工频干扰的自适应滤波器
CN102955889B (zh) 一种用于提取时域特征点的脉搏波重构方法
CN108272451A (zh) 一种基于改进小波变换的qrs波识别方法
CN107361764B (zh) 一种心电信号特征波形r波的快速提取方法
Lin et al. Wheeze recognition based on 2D bilateral filtering of spectrogram
TW201225912A (en) Method for measuring physiological parameters

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180320