Nothing Special   »   [go: up one dir, main page]

CN104739399A - 一种检测脉搏波心率计算中运动抵消的方法 - Google Patents

一种检测脉搏波心率计算中运动抵消的方法 Download PDF

Info

Publication number
CN104739399A
CN104739399A CN201410046189.8A CN201410046189A CN104739399A CN 104739399 A CN104739399 A CN 104739399A CN 201410046189 A CN201410046189 A CN 201410046189A CN 104739399 A CN104739399 A CN 104739399A
Authority
CN
China
Prior art keywords
heart rate
fourier
pulse wave
signals
fourier transform
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
CN201410046189.8A
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.)
Beijing Scientific & Trading Co Ltd Of Jin Jitong
Original Assignee
Beijing Scientific & Trading Co Ltd Of Jin Jitong
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 Beijing Scientific & Trading Co Ltd Of Jin Jitong filed Critical Beijing Scientific & Trading Co Ltd Of Jin Jitong
Priority to CN201410046189.8A priority Critical patent/CN104739399A/zh
Publication of CN104739399A publication Critical patent/CN104739399A/zh
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Mathematical Physics (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

一种检测脉搏波心率计算中运动抵消的方法,包括以下步骤:利用光电容积法测量得脉搏波信号x(n),利用另外一种波长或者加速度传感器获得运动信号y(n);截取一定时间长度N的信号,利用快速傅立叶变换分别获得两种信号的傅立叶变换X(k)和Y(k),对运动信号的傅立叶谱进行平滑处理,获得对运动干扰傅立叶谱的估计Pyk;把脉搏信号的傅立叶谱与归一化后的运动傅立叶谱进行如下所示运算,得到第一次加权后的脉搏信号的傅立叶谱,利用之前心率的值来进行第二次功率谱加权,假设之前心率所对应的频率点值为k0,利用高斯函数构建一个加权曲线。

Description

一种检测脉搏波心率计算中运动抵消的方法
技术领域
本发明涉及一种搏波心率计算方法,特别是涉及一种脉搏波心率计算中运动抵消的方法。
背景技术
光电容积法利用血液对光的吸收原理,来探测血管中血液体积的变化情况,获得使用者的脉搏波,从而实现心率和血氧等参数的测量。然而利用光电容积法获得的脉搏波很容易受到运动的干扰,从而导致脉搏波测量的重大误差。许多研究提出采用一个独立的运动信号测量通道,来实现脉搏波中运动干扰的抵消。
利用双波长的方法来实现运动的抵消,该传感器具有绿色LED和红外LED两种波长的光源。其中绿色LED用来测量脉搏波信号,而红外LED用来测量运动信号。对脉搏波信号和运动信号进行快速傅立叶变换后,对傅立叶光谱进行标准化。然后从脉搏波信号的光谱中减去运动信号的光谱从而是实现运动干扰的去除,并且计算人的脉搏。然而,该方法中光谱的强度必须良好适配以允许完全去除在光波的傅立叶光谱的减法运算之后的运动影响,这是一种缺陷,此外该方法进行测量时,观测到高的耗电量。这构成了另一种缺陷。
专利CN102008299B提出采用光谱相干函数来实现运动伪迹的去除。然而该方法中,相干函数的计算要求对多段的傅立叶光谱进行平均,当要求较高的心率精度时,需要较长时间的信号进行运算。不能满足检测快速的心率变化。此外该方法要求在一定时间内,运动信号较为平稳,因此该方法在用于一些无规则运动所造成的非平稳干扰抵消时,其效果较差。
相关文献采用加速度传感器来测量运动信号,同时利用自适应滤波器来实现运动抵消。然而加速度测量的运动信号和脉搏波中的运动干扰差异较大,自适应抵消方法没有获得较好的性能。
发明内容
本发明的目标是通过一种功率谱加权的方法,来克服以上所提到的现有技术的缺陷,更好的降低运动干扰对心率测量的影响。
本发明的技术方案如下:一种脉搏波心率计算中运动抵消的方法,包括以下步骤:
1)利用光电容积法测量得脉搏波信号x(n),利用另外一种波长或者加速度传感器获得运动信号y(n);
2)截取一定时间长度N的信号,利用快速傅立叶变换分别获得两种信号的傅立叶变换X(k)和Y(k),其中X(k)=FFT(X)为脉搏波信号的傅立叶变换,Y(k)=FFT(Y)为运动信号的傅立叶变换,随后对傅立叶变换进行如下运算,获得两个信号的归一化的傅立叶谱Pxk、Py(k),所采用的公式如下:
Pxk=Xk2Xk2,k=0,1…N/2
Pyk=Yk2Yk2,k=0,1…N/2
3)对运动信号的傅立叶谱进行平滑处理,获得对运动干扰傅立叶谱的估计Pyk;
Pyk=12*mn=-mn=m-1Pyk+n,k=0,1…N/2
4)把脉搏信号的傅立叶谱与归一化后的运动傅立叶谱进行如下所示运算,得到第一次加权后的脉搏信号的傅立叶谱Px_weight_1
Pxweight_1(k)=Pyk*Pxk,k=0,1…N/2
5)利用之前心率的值来进行第二次功率谱加权,假设之前心率所对应的频率点值为k0,利用高斯函数构建一个加权曲线Wk:
Wk=e-(k-k0)2δ2,k=0,1…N/2;6
6)最终得到的加权光谱曲线为:
Pxweight_2(k)=Pxweight_1(k)*Wk,k=0,1…N/2
7)最终通过选择Pxweight_2(k)曲线中的峰值点作为最终的心率值,对每一段数据利用上述方法获得心率值后,可以把心率值的变化画成一条曲线,表明运动过程中的心率变化。
附图说明
图1为该方法的实现示意图。
图2为脉搏波信号被运动信号所干扰示意图
图3为脉搏波信号和运动信号的傅立叶谱示意图
图4为加权后的傅立叶谱
图5为利用本方法获得的心率曲线与基于心电方法获得的心率曲线对比
具体实施方式
下面结合附图对本发明作进一步的详细说明。
图1表明了该方法的实现框图,有两个输入信号,一个是利用光电容积法测量得到的脉搏波信号x(n),一个是利用另外一种波长或者加速度传感器获得的运动信号y(n)。图2显示了一段数据,实线为脉搏波信号,虚线为运动信号。可以看到脉搏波信号被运动信号所干扰,导致脉搏跳动的信息完全被掩盖在运动干扰中。
(1)截取一定时间长度N的信号,利用快速傅立叶变换分别获得两种信号的傅立叶变换X(k)和Y(k)。其中X(k)=FFT(X)为脉搏波信号的傅立叶变换,Y(k)=FFT(Y)为运动信号的傅立叶变换;随后我们对傅立叶变换进行如下运算,获得两个信号的归一化的傅立叶谱Pxk、Py(k),所采用的公式如下:
Pxk=Xk2Xk2,k=0,1…N/2
Pyk=Yk2Yk2,k=0,1…N/2
考虑到心率精度的要求,一般选择数据长度为16-32s。分别对应3.75bpm和1.875bpm的心率精度。图2显示了两个信号的归一化傅立叶光谱。从图中可以看出,脉搏波的傅立叶谱在80-120心跳/分钟和110-130心跳每分钟处有两个峰值,而运动信号的傅立叶谱在80-120心跳/分钟范围内也有一个峰值。
(2)随后我们对运动信号的傅立叶谱进行平滑处理,获得对运动干扰傅立叶谱的估计Pyk。
Pyk=12*mn=-mn=m-1Pyk+n,k=0,1…N/2
(3)把脉搏信号的傅立叶谱与归一化后的运动傅立叶谱进行如下所示运算,得到第一次加权后的脉搏信号的傅立叶谱Px_weight_1。图3显示了图1所示信号的第一次加权傅立叶谱。可以看到,运算后的傅立叶谱中,运动成分已经被抵消,而心率所对应的谱能够较好的显示出来。
Pxweight_1(k)=Pyk*Pxk,k=0,1…N/2
(4)基于如下实际情况:心率值的变化一般是比较平稳的,不会出现突然的剧烈的变化。因此我们可以利用之前心率的值来进行第二次功率谱加权。假设之前心率所对应的频率点值为k0。我们利用高斯函数构建一个加权曲线Wk
Wk=e-(k-k0)2δ2,k=0,1…N/2
(5)最终得到的加权光谱曲线为
Pxweight_2(k)=Pxweight_1(k)*Wk,k=0,1…N/2
(6)最终通过选择Pxweight_2(k)曲线中的峰值点作为最终的心率值。对每一段数据利用上述方法获得心率值后,可以把心率值的变化画成一条曲线,表明运动过程中的心率变化。图4显示了某一受试者运动过程中,采用本发明所述方法获得的心率曲线和利用心电图获得的标准心率的对比,可以看出本方法能够很好的追踪运动过程中的心率变化。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此领域技术的人士能够了解本发明的内容并加以实施,并不能以此限制本发明的保护范围,凡根据本发明精神实质所作的等效变化或修饰,效果相似,都应涵盖在本发明的保护范围内。

Claims (1)

1.一种检测脉搏波心率计算中运动抵消的方法,其特征在于,包括以下步骤:
1)利用光电容积法测量得脉搏波信号x(n),利用另外一种波长或者加速度传感器获得运动信号y(n);
2)截取一定时间长度N的信号,利用快速傅立叶变换分别获得两种信号的傅立叶变换X(k)和Y(k),其中X(k)=FFT(X)为脉搏波信号的傅立叶变换,Y(k)=FFT(Y)为运动信号的傅立叶变换,随后对傅立叶变换进行如下运算,获得两个信号的归一化的傅立叶谱Pxk、Py(k),所采用的公式如下:
Pxk=Xk2Xk2,k=0,1…N/2
Pyk=Yk2Yk2,k=0,1…N/2
3)对运动信号的傅立叶谱进行平滑处理,获得对运动干扰傅立叶谱的估计Pyk;
Pyk=12*mn=-mn=m-1Pyk+n,k=0,1…N/2
4)把脉搏信号的傅立叶谱与归一化后的运动傅立叶谱进行如下所示运算,得到第一次加权后的脉搏信号的傅立叶谱Px_weight_1
Pxweight_1(k)=Pyk*Pxk,k=0,1…N/2
5)利用之前心率的值来进行第二次功率谱加权,假设之前心率所对应的频率点值为k0,利用高斯函数构建一个加权曲线Wk:
Wk=e-(k-k0)2δ2, k=0,1…N/2;6
6)最终得到的加权光谱曲线为:
Pxweight_2(k)=Pxweight_1(k)*Wk,k=0,1…N/2
7)最终通过选择Pxweight_2(k)曲线中的峰值点作为最终的心率值,对每一段数据利用上述方法获得心率值后,可以把心率值的变化画成一条曲线,表明运动过程中的心率变化。
CN201410046189.8A 2014-02-10 2014-02-10 一种检测脉搏波心率计算中运动抵消的方法 Pending CN104739399A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410046189.8A CN104739399A (zh) 2014-02-10 2014-02-10 一种检测脉搏波心率计算中运动抵消的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410046189.8A CN104739399A (zh) 2014-02-10 2014-02-10 一种检测脉搏波心率计算中运动抵消的方法

Publications (1)

Publication Number Publication Date
CN104739399A true CN104739399A (zh) 2015-07-01

Family

ID=53580065

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410046189.8A Pending CN104739399A (zh) 2014-02-10 2014-02-10 一种检测脉搏波心率计算中运动抵消的方法

Country Status (1)

Country Link
CN (1) CN104739399A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105286845A (zh) * 2015-11-29 2016-02-03 浙江师范大学 一种适用于可穿戴式心率测量设备的运动噪声消除方法
CN105796078A (zh) * 2016-05-12 2016-07-27 中世泓利(北京)健康科技有限公司 一种身体生理与运动参数采集与传输系统的方法
CN106383808A (zh) * 2016-09-18 2017-02-08 时瑞科技(深圳)有限公司 心率心电通用的快速计算系统及方法
CN106618542A (zh) * 2015-10-28 2017-05-10 中国科学院上海高等研究院 一种去噪心率信号检测装置及方法
WO2017185846A1 (en) * 2016-04-29 2017-11-02 Boe Technology Group Co., Ltd. Apparatus and method for determining a health parameter of a subject
CN107669253A (zh) * 2017-11-15 2018-02-09 中国科学院光电研究院 基于光谱成像仪的心率和呼吸率测量方法
CN114305371A (zh) * 2021-12-24 2022-04-12 青岛迈金智能科技股份有限公司 一种骑行心率检测稳定算法及心率计

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1195277A (zh) * 1996-04-17 1998-10-07 精工爱普生株式会社 心律失常检测设备
US20050075553A1 (en) * 2003-10-07 2005-04-07 Denso Corporation Portable biological information monitor apparatus and information management apparatus
CN102988036A (zh) * 2012-12-26 2013-03-27 中国科学院自动化研究所 测量脉率的方法
CN103230267A (zh) * 2013-05-14 2013-08-07 北京理工大学 一种抗运动干扰的脉率提取方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1195277A (zh) * 1996-04-17 1998-10-07 精工爱普生株式会社 心律失常检测设备
US20050075553A1 (en) * 2003-10-07 2005-04-07 Denso Corporation Portable biological information monitor apparatus and information management apparatus
CN102988036A (zh) * 2012-12-26 2013-03-27 中国科学院自动化研究所 测量脉率的方法
CN103230267A (zh) * 2013-05-14 2013-08-07 北京理工大学 一种抗运动干扰的脉率提取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
鄂冬,叶树明,周乐川: "耳部光电容积脉搏波的去运动干扰设计", 《传感器与微系统》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106618542A (zh) * 2015-10-28 2017-05-10 中国科学院上海高等研究院 一种去噪心率信号检测装置及方法
CN105286845A (zh) * 2015-11-29 2016-02-03 浙江师范大学 一种适用于可穿戴式心率测量设备的运动噪声消除方法
WO2017185846A1 (en) * 2016-04-29 2017-11-02 Boe Technology Group Co., Ltd. Apparatus and method for determining a health parameter of a subject
US10524676B2 (en) 2016-04-29 2020-01-07 Boe Technology Group Co., Ltd. Apparatus and method for determining a health parameter of a subject
CN105796078A (zh) * 2016-05-12 2016-07-27 中世泓利(北京)健康科技有限公司 一种身体生理与运动参数采集与传输系统的方法
CN106383808A (zh) * 2016-09-18 2017-02-08 时瑞科技(深圳)有限公司 心率心电通用的快速计算系统及方法
CN106383808B (zh) * 2016-09-18 2019-08-02 时瑞科技(深圳)有限公司 心率心电信号的处理系统及方法
CN107669253A (zh) * 2017-11-15 2018-02-09 中国科学院光电研究院 基于光谱成像仪的心率和呼吸率测量方法
CN114305371A (zh) * 2021-12-24 2022-04-12 青岛迈金智能科技股份有限公司 一种骑行心率检测稳定算法及心率计
CN114305371B (zh) * 2021-12-24 2023-11-28 青岛迈金智能科技股份有限公司 一种骑行心率检测稳定算法及心率计

Similar Documents

Publication Publication Date Title
CN104739399A (zh) 一种检测脉搏波心率计算中运动抵消的方法
CN107949321B (zh) 时域干扰去除和改进的心率测量跟踪机理
US11744520B2 (en) Accuracy of heart rate estimation from photoplethysmographic (PPG) signals
EP3476282A1 (en) Heart rate measuring method and apparatus, and electronic terminal
CN108903929B (zh) 心率检测修正的方法、装置、存储介质和系统
CN103989462B (zh) 一种脉搏波形第一特征点和第二特征点的提取方法
KR101849263B1 (ko) 견고한 심박수 추정
CN104706336B (zh) 一种光电式脉搏信号测量方法、装置及测量设备
WO2016057781A1 (en) Method and apparatus for non-contact fast vital sign acquisition based on radar signal
CN109414203A (zh) 基于光学测量的在线心率估计
CN103315728B (zh) 心率检测和显示方法及其装置
CN102460356A (zh) 用于确定包括两个或更多个换能器的物体的表面上的一个或多个触击或触摸的位置的方法
JP6056389B2 (ja) 心拍推定装置、心拍推定方法及びプログラム
CN106443178A (zh) 一种基于IQuinn‑Rife综合的正弦信号频率估计方法
CN105249925B (zh) 一种中医脉象采集装置及降噪系统和降噪方法
RU2015146506A (ru) Устройство и способ для устранения артефактов движения при проведении экг
CN105588577A (zh) 一种用于运动监测装置的异常计步的检测方法和装置
US20230091163A1 (en) Heart activity monitoring during physical exercise
TW201438669A (zh) 血氧飽和度檢測方法及裝置
JP6519344B2 (ja) 心拍間隔特定プログラム、心拍間隔特定装置、及び心拍間隔特定方法
CN104352234A (zh) 一种生理电信号尖峰奇异点检测方法
CN107526064A (zh) 基于二维特征的自适应lfm信号参数估计方法
Ma et al. A multichannel nonlinear adaptive noise canceller based on generalized FLANN for fetal ECG extraction
US10980485B2 (en) Measuring apparatus, measuring method and non-transitory computer readable medium
CN103027692B (zh) 一种基于不确定度的动态光谱数据处理方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150701

WD01 Invention patent application deemed withdrawn after publication